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 8*d71ae5a4SJacob Faibussowitsch int main(int argc, char **argv) 9*d71ae5a4SJacob Faibussowitsch { 10c4762a1bSJed Brown Mat J; 11c4762a1bSJed Brown DM da; 12c4762a1bSJed Brown MatFDColoring matfdcoloring = 0; 13c4762a1bSJed Brown ISColoring iscoloring; 14c4762a1bSJed Brown 15327415f7SBarry 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