pslaqr4/pdlaqr4

Post here if you want to report a bug to the LAPACK team

pslaqr4/pdlaqr4

Postby gutheil » Thu Feb 23, 2012 6:42 am

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.
--
Inge Gutheil
Institute for Advanced Simulation
Juelich Supercomputing Centre
Forschungszentrum Juelich GmbH
52425 Juelich
Germany

phone: +49 2461 61 3135
e-mail: i.gutheil@fz-juelich.de
gutheil
 
Posts: 3
Joined: Thu Feb 23, 2012 5:51 am

Re: pslaqr4/pdlaqr4

Postby admin » Thu Feb 23, 2012 7:55 pm

Indeed! nice catch.
This just has been corrected in the ScaLAPACK repository, rev 168.
Thank you for reporting the problem
admin
Site Admin
 
Posts: 501
Joined: Wed Dec 08, 2004 7:07 pm


Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest