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 PetscErrorCode ierr; 12c4762a1bSJed Brown DM da; 13c4762a1bSJed Brown MatFDColoring matfdcoloring = 0; 14c4762a1bSJed Brown ISColoring iscoloring; 15c4762a1bSJed Brown 16c4762a1bSJed Brown ierr = PetscInitialize(&argc,&argv,(char*)0,help);if (ierr) return ierr; 17c4762a1bSJed Brown /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 18c4762a1bSJed Brown Create distributed array (DMDA) to manage parallel grid and vectors 19c4762a1bSJed Brown - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ 20*5f80ce2aSJacob 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)); 21*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetFromOptions(da)); 22*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetUp(da)); 23*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetMatType(da,MATAIJ)); 24*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateMatrix(da,&J)); 25*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateColoring(da,IS_COLORING_LOCAL,&iscoloring)); 26*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatFDColoringCreate(J,iscoloring,&matfdcoloring)); 27*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatFDColoringSetUp(J,iscoloring,matfdcoloring)); 28*5f80ce2aSJacob Faibussowitsch CHKERRQ(ISColoringDestroy(&iscoloring)); 29c4762a1bSJed Brown 30c4762a1bSJed Brown /* free spaces */ 31*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatDestroy(&J)); 32*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatFDColoringDestroy(&matfdcoloring)); 33*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDestroy(&da)); 34c4762a1bSJed Brown ierr = PetscFinalize(); 35c4762a1bSJed Brown return ierr; 36c4762a1bSJed Brown } 37