Hi Mikael, I had only a quick glance at this patch, but am a little concerned about the tests involving nans.
E.g.:
+ subroutine check_all_nans() + real, allocatable :: a(:,:,:) + real :: nan + integer, allocatable :: m(:) + nan = 0 + nan = nan / nan + allocate(a(3,3,3), source = nan) + m = maxloc(a) + if (size(m, dim=1) /= 3) stop 161 + if (any(m /= (/ 1, 1, 1 /))) stop 162 + end subroutineT
Is there a reason you do not use the ieee intrinsic module way to get a quiet nan? Otherwise, how do you prevent exceptions to happen, possibly leading to a failing test? (The test cases need a workaround to run with NAG). Thanks, Harald
