scalapack icon indicating copy to clipboard operation
scalapack copied to clipboard

P?GETRF's documentation is missing contrains on IA and JA

Open grisuthedragon opened this issue 3 years ago • 0 comments

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.

grisuthedragon avatar May 28 '21 07:05 grisuthedragon