WIP: Reorganized PtAP for MPIAIJ and MAIJ
It is the second attempt to reorganize PtAP
for MPIAIJ
and MAIJ
. The movitaiton is to be able to select an PtAP algorithm in PCHMG. I want to use allatonce
by default because we concern about the memory usage of other algorithms. BTW, allatonce
is negligible slower than scalable
for realistic problems.
@BarrySmith I am still not sure which way it the right path to reorganize PtAP. If you do not like the current design, please add suggestions.
Edited by Barry Smith