Thanks for the quick response. I have altered my code so that it
does not depend on the way extended precision is handled. It
gives the right answer no mater how it is compiled.
I guess the moral is: don't try to check for exact equality between two
double precision numbers that have been calcul
On Oct 26, 2004, at 3:01 PM, Steven E. Williamson wrote:
Two other symptoms that I have noticed: if the -ffloat-store option is
used, then the program always works correctly (does only 2 iterations).
And, if the write statement immediately preceding the check for
equal averages is un-commented, the
I have run into a strange problem with g77 when run with the -O or -O1
option but NOT with -O2, -O3, or no optimization. Below is piece
of Fortran code (removed from a larger program) which illustrates the
problem. The basic function of the code is to go through a list
of data (the RAW array) and
I have run into a strange problem with g77 when run with the -O or -O1
option but NOT with -O2, -O3, or no optimization. Below is piece
of Fortran code (removed from a larger program) which illustrates the
problem. The basic function of the code is to go through a list
of data (the RAW array) and
