1c4762a1bSJed Brown 2c4762a1bSJed Brown static char help[] = "Tests error message in DMCreateColoring() with periodic boundary conditions. \n\n"; 3c4762a1bSJed Brown 4c4762a1bSJed Brown #include <petscdm.h> 5c4762a1bSJed Brown #include <petscdmda.h> 6c4762a1bSJed Brown #include <petscmat.h> 7c4762a1bSJed Brown 8c4762a1bSJed Brown int main(int argc,char **argv) 9c4762a1bSJed Brown { 10c4762a1bSJed Brown Mat J; 11c4762a1bSJed Brown DM da; 12c4762a1bSJed Brown MatFDColoring matfdcoloring = 0; 13c4762a1bSJed Brown ISColoring iscoloring; 14c4762a1bSJed Brown 15*327415f7SBarry Smith PetscFunctionBeginUser; 169566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&argv,(char*)0,help)); 17c4762a1bSJed Brown /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 18c4762a1bSJed Brown Create distributed array (DMDA) to manage parallel grid and vectors 19c4762a1bSJed Brown - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ 209566063dSJacob Faibussowitsch PetscCall(DMDACreate2d(PETSC_COMM_WORLD,DM_BOUNDARY_PERIODIC, DM_BOUNDARY_NONE,DMDA_STENCIL_BOX,-5,-5,PETSC_DECIDE,PETSC_DECIDE,1,2,0,0,&da)); 219566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(da)); 229566063dSJacob Faibussowitsch PetscCall(DMSetUp(da)); 239566063dSJacob Faibussowitsch PetscCall(DMSetMatType(da,MATAIJ)); 249566063dSJacob Faibussowitsch PetscCall(DMCreateMatrix(da,&J)); 259566063dSJacob Faibussowitsch PetscCall(DMCreateColoring(da,IS_COLORING_LOCAL,&iscoloring)); 269566063dSJacob Faibussowitsch PetscCall(MatFDColoringCreate(J,iscoloring,&matfdcoloring)); 279566063dSJacob Faibussowitsch PetscCall(MatFDColoringSetUp(J,iscoloring,matfdcoloring)); 289566063dSJacob Faibussowitsch PetscCall(ISColoringDestroy(&iscoloring)); 29c4762a1bSJed Brown 30c4762a1bSJed Brown /* free spaces */ 319566063dSJacob Faibussowitsch PetscCall(MatDestroy(&J)); 329566063dSJacob Faibussowitsch PetscCall(MatFDColoringDestroy(&matfdcoloring)); 339566063dSJacob Faibussowitsch PetscCall(DMDestroy(&da)); 349566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 35b122ec5aSJacob Faibussowitsch return 0; 36c4762a1bSJed Brown } 37