g95 vs. gfortran

[Dominique Dhumieres]

> Yes, the algorithm could be unstable, not sure why though.

One thing you can do with g95 in order to estimate how unstable
is your algorithm is to change the rounding mode (through
compilation options and/or environment setting as documented
in the g95 manual). Looking at how many digits remain unchanged
give you a rough estimate of the accuracy of the results.

Cheers

Dominique

PS it would be nice to have that in gfortran!-)