In scalapack-2.0.1 when compiling pslaqr4 and pdlaqr4 you get the following warning:
"pslaqr4.f", line 457.27: 1513-029 (W) The number of arguments to "sgemm" differs from the number of arguments in a previous reference. You should use the OPTIONAL attribute and an explicit interface to define a procedure with optional arguments.
** pslaqr4 === End of Compilation 1 ===
Going to line 457 shows, that SGEMM (DGEMM for pdlasqr4) is called with
CALL SGEMM( 'N', 'N', KLN, NH, NH, ONE,
$ A( KKROW+(ICOL-1)*LDA ), LDA, V, LDV, ZERO,
$ ZERO, WORK, KLN )
The parameter beta of SGEMM is called twice, zero at the end of line 2 and at the beginning of line 3
of the call. This is obviously a typo.
Institute for Advanced Simulation
Juelich Supercomputing Centre
Forschungszentrum Juelich GmbH
phone: +49 2461 61 3135