>From Dominique, a potential wrong-code:
$ cat mat-2.f90
program mat
implicit none
complex, allocatable :: a(:,:),b(:,:)
complex :: d(1,1)
allocate(a(4,1),b(4,1))
a = cmplx(0.5,2.0)
b = cmplx(2.0,0.5)
d = matmul(a,transpose(b))
print *, d
print *, matmul(a,transpose(b))
deallocate(a,b)
end program mat
$ gfortran -g mat-2.f90
$ gdb ./a.out
GNU gdb 6.8-debian
Copyright (C) 2008 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law. Type "show copying"
and "show warranty" for details.
This GDB was configured as "i486-linux-gnu"...
(gdb) r
Starting program: /home/ig25/Krempel/Matmul/a.out
( 0.0000000 , 4.2500000 )
a.out: ../../../../gcc/trunk/libgfortran/generated/matmul_c4.c:96: matmul_c4:
Assertion `((a)->dtype & 0x07) == 2 || ((b)->dtype & 0x07) == 2' failed.
Program received signal SIGABRT, Aborted.
0xffffe410 in __kernel_vsyscall ()
(gdb) up
#1 0xb7d41640 in raise () from /lib/i686/cmov/libc.so.6
(gdb) up
#2 0xb7d43018 in abort () from /lib/i686/cmov/libc.so.6
(gdb) up
#3 0xb7d3a5be in __assert_fail () from /lib/i686/cmov/libc.so.6
(gdb) up
#4 0xb7f45d7d in *_gfortran_matmul_c4 (retarray=0xbfd6661c, a=0xbfd6680c,
b=0xbfd667b4,
try_blas=0, blas_limit=0, gemm=0)
at ../../../../gcc/trunk/libgfortran/generated/matmul_c4.c:95
95 assert (GFC_DESCRIPTOR_RANK (a) == 2
(gdb) p/x a->dtype
$1 = 0x40880000
(gdb) p/x b->dtype
$2 = 0x0
(gdb)
--
Summary: [4.4 Regression] array descriptors for matmul messed up?
Product: gcc
Version: 4.4.0
Status: UNCONFIRMED
Severity: normal
Priority: P3
Component: fortran
AssignedTo: unassigned at gcc dot gnu dot org
ReportedBy: tkoenig at gcc dot gnu dot org
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=37802
