ScaLAPACK Archives

[Scalapack] Incorrect global matrix indexing in P{DSY/ZHE}GST for IBTYPE

Summary (and in case formatting fails):  in the case IBTYPE==1 and 
UPPER==.TRUE.  for both PDSYGST and PZHEGST, the call to P{DSY/ZHE}GS2 uses 
"IB+K-1" as the 'JB' argument, and I believe it should be "JB+K-1", especially 
since "JB+K-1" is used in the LOWER case.

In the HTML exploration and available tarbell of ScaLAPACK 2.0.2,

Currently lines 292-3 of SRC/pdsygst.f the code read as:


00292             CALL 
PDSYGS2<http://www.netlib.org/scalapack/explore-html/dc/db9/pdsygs2_8f.html#a5826e0a546b188e7fd378d4b1db3b20d>(
 IBTYPE, UPLO, KB, A, IA+K-1, JA+K-1, DESCA, B,

00293      $                    IB+K-1, IB+K-1, DESCB, INFO )

but I'm fairly certain this should be:


00292             CALL 
PDSYGS2<http://www.netlib.org/scalapack/explore-html/dc/db9/pdsygs2_8f.html#a5826e0a546b188e7fd378d4b1db3b20d>(
 IBTYPE, UPLO, KB, A, IA+K-1, JA+K-1, DESCA, B,

00293      $                    IB+K-1, JB+K-1, DESCB, INFO )

Simlarly, lines 295-6 of SRC/pzhegst.f are:


00295             CALL 
PZHEGS2<http://www.netlib.org/scalapack/explore-html/df/d6f/pzhegs2_8f.html#ae96030ef5ceb2e76b1983518be17c83b>(
 IBTYPE, UPLO, KB, A, IA+K-1, JA+K-1, DESCA, B,

00296      $                    IB+K-1, IB+K-1, DESCB, INFO )

and probably should be:


00295             CALL 
PZHEGS2<http://www.netlib.org/scalapack/explore-html/df/d6f/pzhegs2_8f.html#ae96030ef5ceb2e76b1983518be17c83b>(
 IBTYPE, UPLO, KB, A, IA+K-1, JA+K-1, DESCA, B,

00296      $                    IB+K-1, JB+K-1, DESCB, INFO )

-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
<http://lists.eecs.utk.edu/mailman/private/scalapack/attachments/20131120/9bf7845f/attachment.html>

<Prev in Thread] Current Thread [Next in Thread>
  • [Scalapack] Incorrect global matrix indexing in P{DSY/ZHE}GST for IBTYPE==1, UPPER, Tim Wentz <=


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or