[petsc-dev] Weird interaction between VecDuplicate and VecSetBlockSize

[Jed Brown]

The PetscLayout is shared. It should probably not be allowed to change
block size after duplicating.
On Mar 30, 2012 9:02 AM, "Blaise Bourdin" <bourdin at lsu.edu> wrote:

> Hi,
>
> I am trying to make use of the more advanced features of the hdf5 viewer.
>
> I get the following error when running
> src/vec/vec/examples/tutorials/ex19.c:
> [0]PETSC ERROR: VecLoad_HDF5() line 251 in
> /opt/HPC/petsc-dev/src/vec/vec/utils/vecio.c Block size 1 specified for
> vector does not match blocksize in file 2
>
> Upon further investigating, I noticed that if  Vec x2 is obtained by
> duplicating x1, then setting the block size for x2 also sets it for x1. Is
> this really the desired behavior?
>
> I am attaching an trivial example.
>
> Blaise
>
>
>
>
> --
> Department of Mathematics and Center for Computation & Technology
> Louisiana State University, Baton Rouge, LA 70803, USA
> Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276
> http://www.math.lsu.edu/~bourdin
>
>
>
>
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120330/890dd56a/attachment.html>