PetscErrorCode DMCreateGlobalVector_DA(DM da,Vec *g)
{
PetscErrorCode ierr;
DM_DA *dd = (DM_DA*)da->data;
PetscFunctionBegin;
PetscValidHeaderSpecific(da,DM_CLASSID,1);
PetscValidPointer(g,2);
if (da->defaultSection) {
ierr = DMCreateGlobalVector_Section_Private(da,g);CHKERRQ(ierr);
/* The view and load functions break for general layouts */
PetscFunctionReturn(0);
} else {
ierr = VecCreate(PetscObjectComm((PetscObject)da),g);CHKERRQ(ierr);
ierr = VecSetSizes(*g,dd->Nlocal,PETSC_DETERMINE);CHKERRQ(ierr);
ierr = VecSetBlockSize(*g,dd->w);CHKERRQ(ierr);
ierr = VecSetType(*g,da->vectype);CHKERRQ(ierr);
ierr = VecSetDM(*g, da);CHKERRQ(ierr);
ierr = VecSetLocalToGlobalMapping(*g,da->ltogmap);CHKERRQ(ierr);
}
ierr = VecSetOperation(*g,VECOP_VIEW,(void (*)(void))VecView_MPI_DA);CHKERRQ(ierr);
ierr = VecSetOperation(*g,VECOP_LOAD,(void (*)(void))VecLoad_Default_DA);CHKERRQ(ierr);
ierr = VecSetOperation(*g,VECOP_DUPLICATE,(void (*)(void))VecDuplicate_MPI_DA);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
static PetscErrorCode DMCreateLocalVector_Plex(DM dm,Vec *vec)
{
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = DMCreateLocalVector_Section_Private(dm,vec);CHKERRQ(ierr);
ierr = VecSetOperation(*vec, VECOP_VIEW, (void(*)(void)) VecView_Plex_Local);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
static PetscErrorCode DMCreateGlobalVector_Plex(DM dm,Vec *vec)
{
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = DMCreateGlobalVector_Section_Private(dm,vec);CHKERRQ(ierr);
/* ierr = VecSetOperation(*vec, VECOP_DUPLICATE, (void(*)(void)) VecDuplicate_MPI_DM);CHKERRQ(ierr); */
ierr = VecSetOperation(*vec, VECOP_VIEW, (void (*)(void))VecView_Plex);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
EXTERN_C_END
#undef __FUNCT__
#define __FUNCT__ "DMCreateGlobalVector_Composite"
PetscErrorCode DMCreateGlobalVector_Composite(DM dm,Vec *gvec)
{
PetscErrorCode ierr;
DM_Composite *com = (DM_Composite*)dm->data;
PetscFunctionBegin;
PetscValidHeaderSpecific(dm,DM_CLASSID,1);
ierr = DMSetUp(dm);CHKERRQ(ierr);
ierr = VecCreateMPI(((PetscObject)dm)->comm,com->n,com->N,gvec);CHKERRQ(ierr);
ierr = VecSetDM(*gvec, dm);CHKERRQ(ierr);
ierr = VecSetOperation(*gvec,VECOP_VIEW,(void(*)(void))VecView_DMComposite);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
