[petsc-dev] Weird interaction between VecDuplicate and VecSetBlockSize

[Blaise Bourdin]

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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: TestBlockSize.c
Type: application/octet-stream
Size: 2177 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120330/fdecc7c8/attachment.obj>
-------------- next part --------------


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