[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