Possible Buglet in gfortran Release 143193 (and 4.3.1 and 4.3.2)

[Tom Browder]

On Sat, Jan 24, 2009 at 11:08 AM, Thomas Koenig <tkoenig@netcologne.de> wrote:
...
> I'll give you some stock answers before I even know the question :-)
>
> - If you use -ffast-math, leave it off and see if the problem persists.
>
> - For i386, try "-ffloat-store" or "-march=native -mfpmath=sse2" (if you
>  have a system which supports SSE2) to make sure you don't run into the
>  extended precision problem.
>
> - If your code is standard-conforming f77, running it through ftncheck
>  tends to spot a lot of mistakes.

Thanks for the hints, Thomas--I'm saving those for when I get to more
extensive tests for the full-up program.

Right now I'm testing some functions I had to write to replace the
pgf77 trig functions accepting arguments in degrees or returning
degrees instead of radians.

I'm getting the same answers to 6 decimal places for real*4 and 15
decimal places for real*8.  Does that sound about right?

-Tom