From f0e871d647e8ccdf048be36ba9196f410adc677a Mon Sep 17 00:00:00 2001 From: Paul Thomas Date: Sun, 10 Sep 2006 17:21:44 +0000 Subject: re PR fortran/28947 (Double MATMUL() uses wrong array elements) 2006-09-10 Paul Thomas PR libfortran/28947 * m4/matmul.m4: For the case where the second input argument is transposed, ensure that the case with rank (a) == 1 is correctly calculated. * generated/matmul_r4.c: Regenerate. * generated/matmul_r8.c: Regenerate. * generated/matmul_r10.c: Regenerate. * generated/matmul_r16.c: Regenerate. * generated/matmul_c4.c: Regenerate. * generated/matmul_c8.c: Regenerate. * generated/matmul_c10.c: Regenerate. * generated/matmul_c16.c: Regenerate. * generated/matmul_i4.c: Regenerate. * generated/matmul_i8.c: Regenerate. * generated/matmul_i16.c: Regenerate. 2006-09-10 Paul Thomas PR libfortran/28947 gfortran.dg/matmul_4.f90: New test. From-SVN: r116817 --- gcc/testsuite/ChangeLog | 5 +++++ gcc/testsuite/gfortran.dg/matmul_4.f90 | 22 ++++++++++++++++++++++ 2 files changed, 27 insertions(+) create mode 100644 gcc/testsuite/gfortran.dg/matmul_4.f90 (limited to 'gcc') diff --git a/gcc/testsuite/ChangeLog b/gcc/testsuite/ChangeLog index c482122..327a048 100644 --- a/gcc/testsuite/ChangeLog +++ b/gcc/testsuite/ChangeLog @@ -1,5 +1,10 @@ 2006-09-10 Paul Thomas + PR libfortran/28947 + gfortran.dg/matmul_4.f90: New test. + +2006-09-10 Paul Thomas + PR fortran/28959 gfortran.dg/used_types_10: New test. diff --git a/gcc/testsuite/gfortran.dg/matmul_4.f90 b/gcc/testsuite/gfortran.dg/matmul_4.f90 new file mode 100644 index 0000000..8bbaef9 --- /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 +! +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 -- cgit v1.1