aboutsummaryrefslogtreecommitdiff
path: root/gcc
diff options
context:
space:
mode:
authorPaul Thomas <pault@gcc.gnu.org>2006-09-10 17:21:44 +0000
committerPaul Thomas <pault@gcc.gnu.org>2006-09-10 17:21:44 +0000
commitf0e871d647e8ccdf048be36ba9196f410adc677a (patch)
tree8042aa969d343b56b5ee82785eb2c28922abe06c /gcc
parent380bfbbd61795428e53826d379ce1f4bfe1768f0 (diff)
downloadgcc-f0e871d647e8ccdf048be36ba9196f410adc677a.zip
gcc-f0e871d647e8ccdf048be36ba9196f410adc677a.tar.gz
gcc-f0e871d647e8ccdf048be36ba9196f410adc677a.tar.bz2
re PR fortran/28947 (Double MATMUL() uses wrong array elements)
2006-09-10 Paul Thomas <pault@gcc.gnu.org> 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 <pault@gcc.gnu.org> PR libfortran/28947 gfortran.dg/matmul_4.f90: New test. From-SVN: r116817
Diffstat (limited to 'gcc')
-rw-r--r--gcc/testsuite/ChangeLog5
-rw-r--r--gcc/testsuite/gfortran.dg/matmul_4.f9022
2 files changed, 27 insertions, 0 deletions
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 <pault@gcc.gnu.org>
+ PR libfortran/28947
+ gfortran.dg/matmul_4.f90: New test.
+
+2006-09-10 Paul Thomas <pault@gcc.gnu.org>
+
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 <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