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*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscInitialize(&argc,&argv,(char*)0,help)); 16c4762a1bSJed Brown /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 17c4762a1bSJed Brown Create distributed array (DMDA) to manage parallel grid and vectors 18c4762a1bSJed Brown - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ 195f80ce2aSJacob Faibussowitsch CHKERRQ(DMDACreate2d(PETSC_COMM_WORLD,DM_BOUNDARY_PERIODIC, DM_BOUNDARY_NONE,DMDA_STENCIL_BOX,-5,-5,PETSC_DECIDE,PETSC_DECIDE,1,2,0,0,&da)); 205f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetFromOptions(da)); 215f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetUp(da)); 225f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetMatType(da,MATAIJ)); 235f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateMatrix(da,&J)); 245f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateColoring(da,IS_COLORING_LOCAL,&iscoloring)); 255f80ce2aSJacob Faibussowitsch CHKERRQ(MatFDColoringCreate(J,iscoloring,&matfdcoloring)); 265f80ce2aSJacob Faibussowitsch CHKERRQ(MatFDColoringSetUp(J,iscoloring,matfdcoloring)); 275f80ce2aSJacob Faibussowitsch CHKERRQ(ISColoringDestroy(&iscoloring)); 28c4762a1bSJed Brown 29c4762a1bSJed Brown /* free spaces */ 305f80ce2aSJacob Faibussowitsch CHKERRQ(MatDestroy(&J)); 315f80ce2aSJacob Faibussowitsch CHKERRQ(MatFDColoringDestroy(&matfdcoloring)); 325f80ce2aSJacob Faibussowitsch CHKERRQ(DMDestroy(&da)); 33*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscFinalize()); 34*b122ec5aSJacob Faibussowitsch return 0; 35c4762a1bSJed Brown } 36