diff options
Diffstat (limited to 'gcc/testsuite/gfortran.dg/matmul_4.f90')
| -rw-r--r-- | gcc/testsuite/gfortran.dg/matmul_4.f90 | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/gcc/testsuite/gfortran.dg/matmul_4.f90 b/gcc/testsuite/gfortran.dg/matmul_4.f90 new file mode 100644 index 000000000..8bbaef934 --- /dev/null +++ b/gcc/testsuite/gfortran.dg/matmul_4.f90 @@ -0,0 +1,22 @@ +! { dg-do run } +! Check the fix for PR28947, in which the mechanism for dealing +! with matmul (a, transpose (b)) would cause wrong results for +! a having a rank == 1. +! +! Contributed by Harald Anlauf <anlauf@gmx.de> +! +program gfcbug40 + implicit none + + real :: h(3,3), mat(2,3) + + h(:,:) = - HUGE (1.0)/4 ! Preset unused elements suitably... + + h(3,:) = 0 + h(3,3) = 1 + mat(:,:) = 1 + h(3,:) = h(3,:) + matmul (matmul (h(3,:), transpose (mat)), mat) + + if (any (h(3,:) .ne. (/2.0, 2.0, 3.0/))) call abort () + +end program gfcbug40 |