Question about Tpetra::RowMatrixTransposer (was "Tpetra: RowMatrixTransposer::createTranspose not communicating rectangular matrices")
Created by: hcedwar
Tpetra: RowMatrixTransposer::createTranspose not communicating rectangular matrices.
Current Behavior
A unit test in a CSCS application @Patrick81 has rectangular sparse matrix with numRow > numCol, first column all 1 and last column all 2. In serial the transpose is correct. With two processors the upper right and lower left submatrices are not communicated to their new owning process.
Let:
- numRow = 22,
- numCol = 10,
- numProc = 2,
- P0 owns row [0..11),
- P1 owns row [11,22),
- A(0..10,0) = 1 on P0
- A(11..21,0) = 1 on P1
- A(0..10,9) = 2 on P0
- A(11..21,9) = 2 on P1
then RowMatrixTransposer result, AT = A^T
- AT(0,0..10) = 1 on P0
- AT(0,11..21) does not exist
- AT(9,0..10) does not exist
- AT(9,11..21) = 2 on P1
@trilinos/tpetra
Definition of Done
New 2+ processor unit test with rectangular row matrix transpose yielding correct result.