scalapack
scalapack copied to clipboard
P?GETRF's documentation is missing contrains on IA and JA
I ran into an issue while executing pdgetrf
on a submatrix A(2:M, 2:M)
. Therefore I used the following call:
CALL PDGETRF(M-1, N-1, MATA, 2, 2, DESCA, IPIV, INFO)
which caused the following error:
{ 0, 0}: On entry to PDGETRF parameter number 4 had an illegal value
{ 1, 1}: On entry to PDGETRF parameter number 4 had an illegal value
....
but in the documentation there is no word about the contraint, that IA and JA must be a multiple of the block size.
More general, why is this restriction? I currently searching in P?GETRF
, P?GETF2
, .... why this is required? If there is no reason, then the checks should be removed.