ZCGESV

[kdallen]

ZCGESV (as well as other routines) call ZLANGE( NORM, M, N, A, LDA, WORK ) where WORK is a DOUBLE PRECISION work array.

In the case of ZCGESV though, the work array used in the call to ZLANGE is COMPLEX, causing a compiler warning.

Is there an error in ZCGESV?

[Julien Langou]

Hello

yes we have observed the same warning recently. This was with the Intel Fortran
compiler + all the flags on (Julie might be able give the exact flags). It turns out that the
compiler managed to check the interface of ZLANGE while compiling ZCGESV, and so
complained since ZCGESV is calling ZLANGE with a COMPLEX*16 WORK array instead
of a DOUBLE PRECISION array.

Great.

ZLANGE is used in ZCGESV once, with the argument 'I' so we really need to provide an
array of DOUBLE PRECISION of size N and, well we provide the DOUBLE COMPLEX
WORK array that is at this point in the routine unused and large enough.

Since we already have in the interface two WORK arrays (mainly for strong typing of the Fortran language), we were not willing to have a third one which would:
- overload the interface
- overload the memory requirement of the routine.

If someone sees a way to avoid this warning, let me know. But technically it is safe and
we know what we are doing.

Note as well that the interface might change as written in the header, so, well, if there is
comments or suggestion, send them out.

Thanks to report the warning,
Julien.