| 6f33a894 | 14-Aug-2015 |
Barry Smith <bsmith@mcs.anl.gov> |
MatShift_MPI/SeqXAIJ() could hang if some processes had no entries on a process while others had entries because some processes would attempt a parallel preallocation and the others would not.
Fixed
MatShift_MPI/SeqXAIJ() could hang if some processes had no entries on a process while others had entries because some processes would attempt a parallel preallocation and the others would not.
Fixed by first checking if no preallocation was done, and if not doing. Otherwise preallocation is only done if approprate by each process on the diagonal block portion of the matrix, thus not requiring all processes that share the matrix to call the parallel preallocation routine
Reported-by: Patrick Lacasse <patrick.m.lacasse@gmail.com>
show more ...
|