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
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
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
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.
