On Thu, Mar 1, 2012 at 11:03 AM, Barry Smith <[email protected]> wrote:
MatCreateMPIAIJ(), MatCreateMPIBAIJ(), MatCreateMPISBAIJ() actually create Seq version of the matrices when run on one process.
Meanwhile MatMPIXXXSetPreallocation() always do preallocation for MPI case.
Shouldn't MatCreateMPIXXX() be MatCreateXXX()?
I agree, and we already have MATAIJ so it matches a type.
Also should we add MatXXXSetPreallocation() that calls both MatSeqXXXSetPreallocation() and MatMPIXXXSetPreallocation()?
In fact, that exists: MatXAIJSetPreallocation(). I use it for DMComplex. Matt
Barry
-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener