1! { dg-do run }
2! PR 34566 - logical matmul used to give the wrong result.
3! We check this by running through every permutation in
4! multiplying two 3*3 matrices, and all permutations of multiplying
5! a 3-vector and a 3*3 matrices  and checking against equivalence
6! with integer matrix multiply.
7program main
8  implicit none
9  integer, parameter :: ki=4
10  integer, parameter :: dimen=3
11  integer :: i, j, k
12  real, dimension(dimen,dimen) :: r1, r2
13  integer, dimension(dimen,dimen) :: m1, m2
14  logical(kind=ki), dimension(dimen,dimen) :: l1, l2
15  logical(kind=ki), dimension(dimen*dimen) :: laux
16  logical(kind=ki), dimension(dimen) :: lv
17  integer, dimension(dimen) :: iv
18
19  do i=0,2**(dimen*dimen)-1
20     forall (k=1:dimen*dimen)
21        laux(k) = btest(i, k-1)
22     end forall
23     l1 = reshape(laux,shape(l1))
24     m1 = ltoi(l1)
25
26     ! Check matrix*matrix multiply
27     do j=0,2**(dimen*dimen)-1
28        forall (k=1:dimen*dimen)
29           laux(k) = btest(i, k-1)
30        end forall
31        l2 = reshape(laux,shape(l2))
32        m2 = ltoi(l2)
33        if (any(matmul(l1,l2) .neqv. (matmul(m1,m2) /= 0))) then
34          call abort
35        end if
36     end do
37
38     ! Check vector*matrix and matrix*vector multiply.
39     do j=0,2**dimen-1
40        forall (k=1:dimen)
41           lv(k) = btest(j, k-1)
42        end forall
43        iv = ltoi(lv)
44        if (any(matmul(lv,l1) .neqv. (matmul(iv,m1) /=0))) then
45          call abort
46        end if
47        if (any(matmul(l1,lv) .neqv. (matmul(m1,iv) /= 0))) then
48          call abort
49        end if
50     end do
51  end do
52
53contains
54  elemental function ltoi(v)
55    implicit none
56    integer :: ltoi
57    real :: rtoi
58    logical(kind=4), intent(in) :: v
59    if (v) then
60       ltoi = 1
61    else
62       ltoi = 0
63    end if
64  end function ltoi
65
66end program main
67