[petsc-users] (no subject)

[Barry Smith]

  ierr = DMDACreateNaturalVector(da,&natural);CHKERRQ(ierr);
    ierr = DMDAGlobalToNaturalBegin(da,xin,INSERT_VALUES,natural);CHKERRQ(ierr);
    ierr = DMDAGlobalToNaturalEnd(da,xin,INSERT_VALUES,natural);CHKERRQ(ierr);

Now do VecScatterCreateToZero() from natural and the vector will be in the natural ordering on process zero with the dof interlaced.


   Barry

On Dec 27, 2012, at 12:22 AM, amlan barua <abarua at iit.edu> wrote:

> Hi,
> Is there an analogue of VecScatterCreateToZero for DA vectors? The DMDA object has more than one degrees of freedom.
> If there isn't any, should I use an IS object to do the scattering?
> Amlan