Bill and Mike,
Thanks for your help. I think Mike was right about the
80-bit/64-bit compare. As Mike thought, the error occurs at a test for
equality, i.e.
if (x .ge. y) then
I know better than to test reals for equality, but this is a closed
interval test, and it still fails.
on so others don't struggle with this? As far as I
know the same code worked under g77.
Thanks, Dave Roberts
--
David W. Roberts office 406-994-4548
Professo
you all, Dave
Prof Brian Ripley wrote:
> On Mon, 17 Dec 2007, Dave Roberts wrote:
>
>> I have been revising some FORTRAN 77 routines in R packages I have
>> previously submitted. Since R is now using gfortan I experimented with
>> some Fortran 90 code (array intrinsics prim
I have been revising some FORTRAN 77 routines in R packages I have
previously submitted. Since R is now using gfortan I experimented with
some Fortran 90 code (array intrinsics primarily). So far the code is
still in F77 fixed format, in files suffixed .f (not .f90), but
incorporates some F90
