[Patch, Fortran] Another floating point speedup

Dominique Dhumieres dominiq@lps.ens.fr
Mon Jan 5 15:02:00 GMT 2009


Daniel,

> I don't think this has something to do with it, ...

You are right for the test I have sent. Could you try:

character(15) :: str="+ .339567+2"
real :: x
read(str,'(F15.6)') x
print *, x
read(str,'(G15.7)') x
print *, x
read(str,*) x
print *, x
end

With an unpatched gfortran I get:

   33.956699    
   33.956699    
At line 7 of file read_tst_1.f90
Fortran runtime error: Bad real number in item 1 of list input

I don't have access to the machine on which I have applied your patch,
so I cannot do the test before this evening.

I am not familiar with the format subtilities, so I don't know what
the default format * is supposed to do with "+ .3...".

>  as this error message is not generated by my code ...

The message is generated because there is an error in reading the data.

Dominique



More information about the Fortran mailing list