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" <[email protected]> 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