------- Comment #2 from kargl at gcc dot gnu dot org 2006-05-18 20:41 -------
Actually, transpose by itself does the right thing. It is the
combination of matmul and transpose that is screwing up indexing.
A slightly modified version of your program gives
troutmask:kargl[290] ./z
1.0 0.0 0.0 1.0 <-- x
1.0 0.0 0.0 1.0 <-- a = test()
1.0 0.0 0.0 1.0 <-- a = transpose(test())
1.0 0.0 0.0 1.0 <-- z = matmul(x,a)
1.0 0.0 1.0 0.0 <-- z = matmul(x,transpose(test()))
Abort (core dumped)
Here the program.
program foo
implicit none
real(kind=kind(1.0d0)), dimension(2, 2):: x, z, a
integer i, j
x(1,1) = 1.d0
x(2,1) = 0.d0
x(1,2) = 0.d0
x(2,2) = 1.d0
!
! Uncomment to see the output.
!
print '(4F4.1,1X,A)', x(1,1), x(2,1), x(1,2), x(2,2), '<-- x'
a = test()
print '(4F4.1,1X,A)', a(1,1), a(2,1), a(1,2), a(2,2), '<-- a = test()'
a = transpose(test())
print '(4F4.1,1X,A)', a(1,1), a(2,1), a(1,2), a(2,2), &
& '<-- a = transpose(test())'
z = matmul(x, a)
print '(4F4.1,1X,A)', z(1,1), z(2,1), z(1,2), z(2,2), '<-- z = matmul(x,a)'
z = matmul(x, transpose(test()))
print '(4F4.1,1X,A)', z(1,1), z(2,1), z(1,2), z(2,2), &
& '<-- z = matmul(x,transpose(test()))'
z = matmul(x, transpose(test()))
do i = 1, 2
do j = 1, 2
if (x(i,j) .ne. z(i,j)) call abort()
end do
end do
close (10)
contains
function test() result (res)
real(kind=kind(1.0d0)), dimension(2,2) :: res
res(1,1) = 1.d0
res(2,1) = 0.d0
res(1,2) = 0.d0
res(2,2) = 1.d0
end function
end
--
kargl at gcc dot gnu dot org changed:
What |Removed |Added
----------------------------------------------------------------------------
Status|UNCONFIRMED |NEW
Ever Confirmed|0 |1
Last reconfirmed|0000-00-00 00:00:00 |2006-05-18 20:41:45
date| |
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=27662
