| ! { dg-do run } |
| ! { dg-options "-finline-matmul-limit=0" } |
| ! Stress-test the matmul blocking code with sizes close to or |
| ! equal to powers ot two. |
| |
| program main |
| implicit none |
| integer, dimension(*), parameter :: nn = & |
| & [2,3,4,5, 7,8,9, 15,16,17, 31,32,33, 63,64,65, & |
| 127 ,228,129, 255,256,257]; |
| integer, parameter :: s = size(nn) |
| real, dimension(:,:),allocatable :: a, b, c |
| integer :: i1, i2, i3 |
| integer :: nx, ny, count |
| real :: sm |
| |
| sm = 0.0 |
| do i1=1, s |
| nx = nn(i1) |
| do i2=1,s |
| ny = nn(i2) |
| do i3=1,s |
| count = nn(i3) |
| allocate (a(nx,ny), b(ny,count), c(nx,count)) |
| call random_number(a) |
| call random_number(b) |
| c = matmul(a,b) |
| sm = sm + sum(c) |
| deallocate(a,b,c) |
| end do |
| end do |
| end do |
| |
| end program main |
