1! { dg-do run }
2program main
3  integer, parameter :: sz=5, su=3
4  integer, parameter :: l=2
5  integer, parameter :: u=l-1+su
6  integer(kind=4), dimension(sz,sz) :: r,a,b
7  integer :: i,j
8  do i=1,4
9     do j=1,4
10        a(i,j) = i*10+j
11        b(i,j) = 100+i*10+j
12     end do
13  end do
14  r = -1
15  b(l:u,l:u) = reshape([(i,i=1,su*su)],[su,su]);
16  a(l:u,l:u) = reshape([(i,i=1,su*su)],[su,su]);
17
18  r(1:su,1:su) = matmul(a(l:u,l:u),b(l:u,l:u))
19  if (any(reshape(r,[sz*sz]) /= [30, 36, 42, -1, -1, 66, 81, 96, -1, -1,&
20       & 102, 126, 150, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1])) &
21       STOP 1
22end program main
23