[Bug fortran/31994] conjg(transpose(a)) produces wrong answer.

fxcoudert at gcc dot gnu dot org gcc-bugzilla@gcc.gnu.org
Sun May 20 14:43:00 GMT 2007



------- Comment #6 from fxcoudert at gcc dot gnu dot org  2007-05-20 15:43 -------
(In reply to comment #5)
> What is more
> confusing is that I can't trigger this scalarization bug with other intrinsics

10 seconds after I hit the Commit button, I thought about another testcase, and
it does also trigger the bug:

  integer :: a(1,1)
  real :: b(1,1)
  a = 0
  b = real(transpose(a))
  print *, b(1,1)
  end

This time, the interesting thing is that changing "b = real(transpose(a))" into
"b = transpose(a)" gets rid of the bug. I think it's somewhere hidden in the
intrinsics scalarization, but can't find it. Hopefully someone can go further
with this information!


-- 

fxcoudert at gcc dot gnu dot org changed:

           What    |Removed                     |Added
----------------------------------------------------------------------------
   Last reconfirmed|2007-05-19 03:29:11         |2007-05-20 15:43:41
               date|                            |


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31994



More information about the Gcc-bugs mailing list