All I did was to replace the
"real*8" statements by
"double precision" in the *.f and *.h
and recompile with Intel Fortran.
It is necessary to reduce the
predefined machine epsilon for the
NAG library routines. It is in
x02ajft.f. I simply changed one line
in this file to
X02AJF = 1.e-22
Normally, one would expect that a mantissa
of more than 30 decimal places is supported,
but this turned out not to be the case
on a Pentium M and Intel Fortran V.8
(I experimented a bit to determine
the effective length of the mantissa).
However, on an Athlon64 with Intel Fortran 9.1
the mantissa seems to have the 30+ decimal
places one would expect for 128bit.