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?

- Adrian

-- 
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.croucher@auckland.ac.nz
tel: +64 (0)9 923 84611