1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
! { dg-do run }
! { dg-options "-fdump-tree-original" }
!
! PR fortran/56318
!
! Contributed by Alberto Luaces
!
SUBROUTINE mass_matrix
DOUBLE PRECISION,PARAMETER::m1=1.d0
DOUBLE PRECISION,DIMENSION(3,2),PARAMETER::A1=reshape([1.d0,0.d0, 0.d0, &
0.d0,1.d0, 0.d0],[3,2])
DOUBLE PRECISION,DIMENSION(2,2),PARAMETER::Mel=reshape([1.d0/3.d0, 0.d0, &
0.d0, 1.d0/3.d0],[2,2])
DOUBLE PRECISION,DIMENSION(3,3)::MM1
MM1=m1*matmul(A1,matmul(Mel,transpose(A1)))
!print '(3f8.3)', MM1
if (any (abs (MM1 &
- reshape ([1.d0/3.d0, 0.d0, 0.d0, &
0.d0, 1.d0/3.d0, 0.d0, &
0.d0, 0.d0, 0.d0], &
[3,3])) > epsilon(1.0d0))) &
STOP 1
END SUBROUTINE mass_matrix
program name
implicit none
integer, parameter :: A(3,2) = reshape([1,2,3,4,5,6],[3,2])
integer, parameter :: B(2,3) = reshape([3,17,23,31,43,71],[2,3])
integer, parameter :: C(3) = [-5,-7,-21]
integer, parameter :: m1 = 1
! print *, matmul(B,C)
if (any (matmul(B,C) /= [-1079, -1793])) STOP 2
! print *, matmul(C,A)
if (any (matmul(C,A) /= [-82, -181])) STOP 3
! print '(3i5)', m1*matmul(A,B)
if (any (m1*matmul(A,B) /= reshape([71,91,111, 147,201,255, 327,441,555],&
[3,3]))) &
STOP 4
call mass_matrix
end program name
! { dg-final { scan-tree-dump-times "matmul" 0 "original" } }
|
