------- Comment #2 from jv244 at cam dot ac dot uk 2007-06-01 11:06 ------- reduced:
SUBROUTINE polint(n) INTEGER, PARAMETER :: dp=KIND(0.0D0) REAL(dp) :: ya(n), xa(n), x DO m = 1,n-1 DO i = 1,n-m ho=xa(i)-x hp=xa(i+m)-x den=ho-hp IF(den.eq.0.) STOP 'POLINT' END DO END DO END SUBROUTINE polint -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=32176