[Bug fortran/60392] Problem with TRANSPOSE and CONTIGUOUS dummy arguments
mikael at gcc dot gnu.org
gcc-bugzilla@gcc.gnu.org
Fri Mar 7 21:01:00 GMT 2014
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=60392
--- Comment #5 from Mikael Morin <mikael at gcc dot gnu.org> ---
Created attachment 32307
--> http://gcc.gnu.org/bugzilla/attachment.cgi?id=32307&action=edit
preliminary patch
This removes the difference between my_mul/my_mul_cont.
However this is not yet correct, with the patch the program output is:
Normal: 0
7 10 15 22
7 10 15 22
Transposed: 0
5 11 11 25
5 11 11 25
and according to maxima (so rather accurate):
matmul(transpose(a), a) == matrix([10, 14], [14, 20])
and
matrix([5, 11], [11, 25]) == matmul(a, transpose(a))
So there remains a wrong transposition hiding somewhere.
More information about the Gcc-bugs
mailing list