On Mon, Feb 6, 2012 at 20:35, Matthew Knepley <[email protected]> wrote:

VecCreate()

...

VecSetType() VecSetGhosted()

Is there a compelling reason for memory to be allocated before VecSetUp()? Deferring it would make it easier to specify what memory to use or specify ghosting in any order. My preference is that ghosting is orthogonal to type, to the extent that it would make sense for the storage format (non-contiguous AMR grids might not be).