------- Comment #10 from elizabeth dot l dot yip at boeing dot com 2007-05-15 07:48 ------- The following code exposes the problem relating to conjg(transpose(.)), with which gfortran returns the wrong answer. If we break conjg(transpose(.)) into two lines, then we can get the right answer.
Is this related to bugs 31196, 31197 and 31258? Should I resubmit this as a new bug? I have tried versions 4.1, 4.2 and the latest 4.3 (05-14-07). I could'nt find a version that works. Elizabeth Yip --------------------------------------------------------------------- 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 -- elizabeth dot l dot yip at boeing dot com changed: What |Removed |Added ---------------------------------------------------------------------------- CC| |elizabeth dot l dot yip at | |boeing dot com http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31197