At 09:35 09-02-98 +0930, you wrote:
>Phillip rightly points out that double complex is not part of the F77
>standard. The LAPACK authors, in their thoroughness, also point this
>out in their documentation, but indicate that in all other respects
>the double complex routines are standard conforming.
>
>However, from a practical point of view I have the following
>question (which is not rhetorical); I would be interested if anybody
>could point out any "serious" compilers that do not support the
>double complex extension?
>
>In Fortran90, is the following a standard conforming way of
>implementing the double complex type?
>
>DP = KIND(0.0D0)
^^^^^
>COMPLEX (KIND = DP) :: X
In my opinion you should write:
Integer, Parameter :: DP = KIND(0.0D0)
COMPLEX (KIND = DP) :: X
---
Meilleures Salutations,
Best Greetings,
/---
Jan van Oosterwijk
Computing Centre
Delft University of Technology
Postbus 354
2600 AJ Delft
Netherlands / Pays Bas
Phone: +31 15 278 50 17
Fax: +31 15 278 37 87
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|