anchor.diff:     integer, parameter :: n = 10
   
      integer :: i, j
  +   logical :: test_ok
   
      real(8), allocatable :: a(:, :)
      real(8), allocatable :: b(:, :)
  
positive.code: 
     call dgemm('n', 'n', n, n, n, 1.0d0, a, n, b, n, 0.0d0, c, n)
  
     do i = 1, n
        do j = 1, n
           if (dabs(c(i, j) - 20.0d0) > tiny(0.0d0)) return
        end do
     end do
  
  
negative.code:    deallocate(b)
     deallocate(c)
  
     if (test_ok) then
        print *, 'dgemm test ok'
     else
        print *, 'dgemm test failed'
     end if
  
  end program
  