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
> > you need: ZGETRI.F computes the inverse of a double precision complex
> ^^^^^^^^^^^^^^^^^^^^^^^^
> > general matrix. It is written in standard Fortran77 and so by
> ^^^^^^^^^^^^^^^^^^
> At most one of these can be true, strictly speaking.
>
>
> --
> Phillip Helbig Email .......... [log in to unmask]
> Nuffield Radio Astronomy Laboratories Tel. ..... +44 1477 571 321 (ext. 297)
> Jodrell Bank Fax ................. +44 1477 571 618
> Macclesfield Telex ................. 36149 JODREL G
> UK-Cheshire SK11 9DL Web .... http://www.jb.man.ac.uk/~pjh/
>
> My opinions are not necessarily those of NRAL or the University of Manchester.
>
>
Regards,
David.
----------------------------------------------------------
David Vowles
Research Officer
Department of Electrical and Electronic Engineering
The University of Adelaide
Australia 5005
Voice: +61 8 8303 5416
Fax: +61 8 8303 4360
Email: [log in to unmask]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|