| ca31476a | 28-Oct-2011 |
Jed Brown <jed@59A2.org> |
Fix MatIncreaseOverlap (-pc_type asm) with MPIAIJ and bs>1 by reverting PetscLayoutFindOwner to original API.
In 46250962cc92: eliminate non-scalable "rtable" from MatIncreaseOverlap_MPIBAIJ_Once, u
Fix MatIncreaseOverlap (-pc_type asm) with MPIAIJ and bs>1 by reverting PetscLayoutFindOwner to original API.
In 46250962cc92: eliminate non-scalable "rtable" from MatIncreaseOverlap_MPIBAIJ_Once, using updated PetscLayoutFindOwner(),
the API for PetscLayoutFindOwner() was changed to index by blocks, but the calls from mpiov.c were not updated, so this function was returning erroneous data which broke -pc_type asm for MPIAIJ if bs>1. This semantic change was not accompanied by a documentation update. Since PetscLayoutFindOwner() does not include the word "block", I have reverted it to operate on scalar indices instead of block indices, so now the original calls from mpiov.c are correct once again. This patch updates the calls from baijov.c to use scalar indices.
Hg-commit: 7782710e3f96d2562f32a7b21cb4bef00b1dbc9d
show more ...
|