! { dg-do run }
! { dg-options "-funsigned" }
! Test matrix multiplication
program memain
  implicit none
  call test1
  call test2
contains
  subroutine test1
    integer, parameter :: n = 10, m = 28
    unsigned, dimension(n,n) :: u, v, w
    integer(kind=8), dimension(n,n) :: i, j, k
    real(8), dimension(n,n) :: a, b

    call random_number(a)
    call random_number(b)
    u = uint(a*2.0**m)
    v = uint(b*2.0**m)
    i = int(a*2.0**m,8)
    j = int(b*2.0**m,8)
    w = matmul(u,v)
    k = mod(matmul(i,j),2_8**32)
    if (any(uint(k) /= w)) error stop 1
  end subroutine test1
  subroutine test2
    unsigned, parameter :: u(3,3) = reshape ([1u, uint(-2), 3u, uint(-4), &
         5u, uint(-6), 7u, uint(-8), 9u],[3,3])
    unsigned, parameter :: v(3,3) = 1u - u
    unsigned, parameter :: w(3,3) = matmul(u,v)
    integer(kind=8), dimension(3,3), parameter :: &
         i = int(u,8), j = int(v,8)
    integer(kind=8), dimension(3,3) :: k = matmul(i,j)
    if (any(uint(k) /= w)) error stop 2
  end subroutine test2
end program memain