On Mon, Feb 15, 2016 at 4:22 PM, Adrian Croucher <[email protected]> wrote:

On 15/02/16 20:36, Matthew Knepley wrote:

Sorry I took too long. It has been hectic this semester. The logic is somewhat complicated with regard to global and local here, but it is to ensure that boundary values are put in on output. I am not sure how to make it more transparent, but I am willing to take suggestions.

Thanks Matt.

I see now the documentation for VecLoad() does say "The input file must contain the full global vector...".

Maybe it'd be useful to add a note to the VecView() documentation saying that the vectors to view should be global rather than local?

The problem is that VecView knows nothing about global/local. Its at the algebraic level. The DM is using extra information to organize the HDF5 file. It should be documented in the manual at least. Thanks, Matt

- Adrian

-- Dr Adrian Croucher Senior Research Fellow Department of Engineering Science University of Auckland, New Zealand email: [email protected] tel: +64 (0)9 923 84611

-- 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