slamch.f and dlamch.f should be replaced by the end user to fit their
hardware. These functions are meant to return constants and are only
provided as an example. Any modern Fortran implementation has these
constants available in the IEEE module.
I'm glad you found a solution and I hope it matches your floating-point hardware.
Apparently, only a FORTRAN 77 compiler is required to compile PLASMA,
as we can see if we inspect the source file ! Nothing is written in the
documentation... e.g. small examples included in the User Guide use
the F90-style comment ("!") but they use the old F77 syntax:
PARAMETER (N = 10)
should writes now:
integer, parameter :: N = 10
So, I deduced that a FORTRAN 77 is sufficient to compile PLASMA.
Therefore, the old 'dlamch.f' and 'slamch.f' files are required to
compute the machine FP parameters...
Moreover, some FP parameters can be obtained with the new intrinsic
functions of Fortran 90 (tiny, huge, ...), but not all.
The IEEE intrinsic module will be available only with full compliant
F2003 standard and today, IMHO, only one compiler verify that,
i.e. the XLF 12 Fortran compiler of IBM.