On Mar 1, 2012, at 11:11 AM, Matthew Knepley wrote: > On Thu, Mar 1, 2012 at 11:03 AM, Barry Smith <bsmith at 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.
I have pushed this. > > 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
