Gianluca Meneghello <[email protected]> writes: > Hi again, > > I am now trying to write a vector to an HDF5, this time using a DM > structure to impose the vector layout. > > Creating the vector with (as in src/dm/examples/tutorials/ex10.c) > > DMCreateGlobalVector(da2D,&gauss); > > 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: > > VecCreate(PETSC_COMM_WORLD,&gauss); > VecSetSizes(gauss,PETSC_DECIDE,Nx*Ny); > VecSetDM(gauss,da2D);
VecSetDM only associates a DM with the Vec. It does not interpose itself into VecView. Should we change this so that creating a Vec of the correct size/block size/layout and associating a DM is equivalent? Perhaps, but it's more nuanced than it might seem. Why don't you want to just call DMCreateGlobalVector()?
signature.asc
Description: PGP signature
