https://gcc.gnu.org/bugzilla/show_bug.cgi?id=94246
--- Comment #8 from anlauf at gcc dot gnu.org ---
Reduced testcase:
implicit none
integer :: i
if (any (BESSEL_YN(0, 10, 0.0) /= [ (BESSEL_YN(i, 0.0), i = 0, 10) ])) &
then
STOP 6
end if
if (any (abs (BESSEL_YN(0, 10, 1.0) &
- [ (BESSEL_YN(i, 1.0), i = 0, 10) ]) &
> epsilon(0.0)*32)) then
STOP 8
end if
end
