<div dir="ltr">Hi again,<div><br></div><div>I am now trying to write a vector to an HDF5, this time using a DM structure to impose the vector layout.</div><div><br></div><div>Creating the vector with (as inĀ src/dm/examples/tutorials/ex10.c)</div><div><br></div><div>DMCreateGlobalVector(da2D,&gauss);<br></div><div><br></div><div>every works fine, i.e. the hdf5 file has the correct layout given by the dimension of the DM. On the other side, if I use:</div><div><br></div><div><div>VecCreate(PETSC_COMM_WORLD,&gauss);</div><div>VecSetSizes(gauss,PETSC_DECIDE,Nx*Ny);</div><div>VecSetDM(gauss,da2D);</div><div>VecSetFromOptions(gauss);</div><div>VecSetUp(gauss);</div></div><div><br></div><div>the hdf5 file contains a single dimension vector, apparently negecting the VecSetDM instruction (I would need this for vectors created by MatCreateVecs, in order to have both right and left vectors available).</div><div><br></div><div>I did some testing but I did not manage to solve this. Is there anything I am doing wrong, or is this the way it is supposer to work?</div><div><br></div><div>Thanks as always for your help,</div><div><br></div><div>Gianluca</div><div><br></div></div>