On Thu, Mar 1, 2012 at 11:03 AM, Barry Smith <bsmith@mcs.anl.gov> 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