Gideon Simpson <[email protected]> writes: > Suppose I have a 1D dmda structure distributed across N processors, > where there are some number of degrees of freedom (dof) associated > with each element. Now, suppose I want to have a vec also distributed > across the N processors, but distributed "in the same way" as the > dmda. That is to say, if I call > > DMDAGetCorners(da, &idx_lo, NULL, NULL, &idx_width, NULL, NULL); > > idx_lo to idx_lo+ idx_width > > are the same indices I'd get if I called > > VecGetOwnershipRange(vec, &idx_low, &idx_hi); > > Can I do this? How?
What's wrong with DMCreateGlobalVector()?
pgpMrvxoxD2oT.pgp
Description: PGP signature
