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
47
|
! { dg-do run }
! PR97063 - Wrong result for vector (step size is negative) * matrix
program p
implicit none
integer, parameter :: m = 3, k = 2*m, l = k-1, n = 4
integer :: i, j, m1, m2, ms
integer :: ai(k), bi(k,n), ci(n), ci_ref(n), c1, c2
real :: ar(k), br(k,n), cr(n), cr_ref(n)
ai(:) = [(i,i=0,k-1)]
bi(:,:) = reshape ([(((5*i+j),i=0,k-1),j=0,n-1)],[k,n])
! Parameters of subscript triplet
m1 = 1; m2 = l; ms = 2
! Reference values for cross-checks: integer variant
c1 = dot_product (ai(m1:m2: ms), bi(m1:m2: ms,1))
c2 = dot_product (ai(m1:m2: ms), bi(m1:m2: ms,2))
ci_ref = matmul (ai(m1:m2: ms), bi(m1:m2: ms,:))
ci = matmul (ai(m2:m1:-ms), bi(m2:m1:-ms,:))
if (ci_ref(1) /= c1 .or. ci_ref(2) /= c2) stop 1
if (any (ci /= ci_ref)) stop 2
! Real variant
ar = real (ai)
br = real (bi)
cr_ref = matmul (ar(m1:m2: ms), br(m1:m2: ms,:))
cr = matmul (ar(m2:m1:-ms), br(m2:m1:-ms,:))
if (any (cr_ref /= real (ci_ref))) stop 3
if (any (cr /= cr_ref )) stop 4
! Mixed variants
cr_ref = matmul (ar(m1:m2: ms), bi(m1:m2: ms,:))
cr = matmul (ar(m2:m1:-ms), bi(m2:m1:-ms,:))
if (any (cr_ref /= real (ci_ref))) stop 5
if (any (cr /= cr_ref )) stop 6
cr_ref = matmul (ai(m1:m2: ms), br(m1:m2: ms,:))
cr = matmul (ai(m2:m1:-ms), br(m2:m1:-ms,:))
if (any (cr_ref /= real (ci_ref))) stop 7
if (any (cr /= cr_ref )) stop 8
end program
|
