[petsc-dev] reading vector from HDF5 with DMPlex

[Matthew Knepley]

On Mon, Feb 15, 2016 at 4:22 PM, Adrian Croucher <a.croucher at auckland.ac.nz>
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: a.croucher at auckland.ac.nz
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20160216/a4eed479/attachment.html>