[Patch,Fortran,committed] PR 45367/36158 - Increase numeric tolerance in bessel_{6,7}.f90
Tobias Burnus
burnus@net-b.de
Mon Aug 23 03:41:00 GMT 2010
H.J. Lu wrote:
> gfortran.dg/bessel_7.f90 still fails with -m32 -mfpmath=sse:
>
> http://gcc.gnu.org/ml/gcc-testresults/2010-08/msg02282.html
> http://gcc.gnu.org/ml/gcc-testresults/2010-08/msg02284.html
I think that's because we have different libm versions. And because the
recurrence algorithm is somehow not very stable for certain combinations
of N and X.
$ gfortran -m32 -mfpmath=sse bessel_7.f90 && ./a.out
$
Can you uncomment the "print" lines in bessel_7 and comment the if and
abort line - and send me the output?
Tobias
More information about the Gcc-patches
mailing list