int main(int argc,char **args)
{
PetscErrorCode ierr;
ierr = PetscInitialize(&argc,&args,NULL,help);
ierr = PetscSequentialPhaseBegin(PETSC_COMM_WORLD,1);CHKERRQ(ierr);
ierr = PetscSequentialPhaseEnd(PETSC_COMM_WORLD,1);CHKERRQ(ierr);
ierr = PetscFinalize();
return ierr;
}
int main(int argc,char **argv)
{
PetscErrorCode ierr;
PetscInt M = -2, N = -3, P = 4,stencil_width = 1, dof = 1,m,n,p,xstart,ystart,zstart,i,j,k,c;
DM da;
Vec global,local;
PetscScalar ****vglobal;
ierr = PetscInitialize(&argc,&argv,(char*)0,help);CHKERRQ(ierr);
PetscFunctionBeginUser;
ierr = PetscOptionsGetInt(0,"-stencil_width",&stencil_width,0);CHKERRQ(ierr);
ierr = PetscOptionsGetInt(0,"-dof",&dof,0);CHKERRQ(ierr);
ierr = DMDACreate3d(PETSC_COMM_WORLD,DM_BOUNDARY_MIRROR,DM_BOUNDARY_MIRROR,DM_BOUNDARY_MIRROR,DMDA_STENCIL_STAR,M,N,P,PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,dof,stencil_width,NULL,NULL,NULL,&da);CHKERRQ(ierr);
ierr = DMDAGetCorners(da,&xstart,&ystart,&zstart,&m,&n,&p);CHKERRQ(ierr);
ierr = DMCreateGlobalVector(da,&global);CHKERRQ(ierr);
ierr = DMDAVecGetArrayDOF(da,global,&vglobal);CHKERRQ(ierr);
for (k=zstart; k<zstart+p; k++) {
for (j=ystart; j<ystart+n; j++) {
for (i=xstart; i<xstart+m; i++) {
for (c=0; c<dof; c++) {
vglobal[k][j][i][c] = 1000*k + 100*j + 10*(i+1) + c;
}
}
}
}
ierr = DMDAVecRestoreArrayDOF(da,global,&vglobal);CHKERRQ(ierr);
ierr = DMCreateLocalVector(da,&local);CHKERRQ(ierr);
ierr = DMGlobalToLocalBegin(da,global,INSERT_VALUES,local);CHKERRQ(ierr);
ierr = DMGlobalToLocalEnd(da,global,INSERT_VALUES,local);CHKERRQ(ierr);
ierr = PetscSequentialPhaseBegin(PETSC_COMM_WORLD,1);CHKERRQ(ierr);
ierr = VecView(local,PETSC_VIEWER_STDOUT_SELF);CHKERRQ(ierr);
ierr = PetscSequentialPhaseEnd(PETSC_COMM_WORLD,1);CHKERRQ(ierr);
ierr = VecView(global,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);
ierr = DMDestroy(&da);CHKERRQ(ierr);
ierr = VecDestroy(&local);CHKERRQ(ierr);
ierr = VecDestroy(&global);CHKERRQ(ierr);
ierr = PetscFinalize();
return 0;
}
int main(int argc,char **argv)
{
PetscErrorCode ierr;
PetscInt M = 6,stencil_width = 1, dof = 1,m,xstart,i,j;
DM da;
Vec global,local;
PetscScalar **vglobal;
ierr = PetscInitialize(&argc,&argv,(char*)0,help);CHKERRQ(ierr);
PetscFunctionBeginUser;
ierr = PetscOptionsGetInt(0,"-stencil_width",&stencil_width,0);CHKERRQ(ierr);
ierr = PetscOptionsGetInt(0,"-dof",&dof,0);CHKERRQ(ierr);
ierr = DMDACreate1d(PETSC_COMM_WORLD,DM_BOUNDARY_MIRROR,M,dof,stencil_width,NULL,&da);CHKERRQ(ierr);
ierr = DMDAGetCorners(da,&xstart,0,0,&m,0,0);CHKERRQ(ierr);
ierr = DMCreateGlobalVector(da,&global);CHKERRQ(ierr);
ierr = DMDAVecGetArrayDOF(da,global,&vglobal);CHKERRQ(ierr);
for (i=xstart; i<xstart+m; i++) {
for (j=0; j<dof; j++) {
vglobal[i][j] = 100*(i+1) + j;
}
}
ierr = DMDAVecRestoreArrayDOF(da,global,&vglobal);CHKERRQ(ierr);
ierr = DMCreateLocalVector(da,&local);CHKERRQ(ierr);
ierr = DMGlobalToLocalBegin(da,global,INSERT_VALUES,local);CHKERRQ(ierr);
ierr = DMGlobalToLocalEnd(da,global,INSERT_VALUES,local);CHKERRQ(ierr);
ierr = PetscSequentialPhaseBegin(PETSC_COMM_WORLD,1);CHKERRQ(ierr);
ierr = VecView(local,PETSC_VIEWER_STDOUT_SELF);CHKERRQ(ierr);
ierr = PetscSequentialPhaseEnd(PETSC_COMM_WORLD,1);CHKERRQ(ierr);
ierr = VecView(global,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);
ierr = DMDestroy(&da);CHKERRQ(ierr);
ierr = VecDestroy(&local);CHKERRQ(ierr);
ierr = VecDestroy(&global);CHKERRQ(ierr);
ierr = PetscFinalize();
return 0;
}
void PETSC_STDCALL petscsequentialphasebegin_(MPI_Fint * comm,int *ng, int *__ierr ){
*__ierr = PetscSequentialPhaseBegin(
MPI_Comm_f2c( *(comm) ),*ng);
}
