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