conjg(transpose(a)) produces the wrong answers. Both absoft and ifort
produce correct answers.
program main
implicit none
complex (kind=4),dimension(2,2)::a,b,c
a(1,1) = (1.,1.)
a(2,1) = (2.,2.)
a(1,2) = (3.,3.)
a(2,2) = (4.,4.)
print *,"original",a
b=conjg(transpose(a))
print *,"H(a)-wrong",b
c=transpose(a)
c=conjg(c)
print *,"H(a)-right",c
END program main
--
Summary: conjg(transpose(a)) produces wrong answer.
Product: gcc
Version: 4.3.0
Status: UNCONFIRMED
Severity: critical
Priority: P3
Component: fortran
AssignedTo: unassigned at gcc dot gnu dot org
ReportedBy: elizabeth dot l dot yip at boeing dot com
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31994
