1c4762a1bSJed Brown /* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. 2c4762a1bSJed Brown 3c4762a1bSJed Brown Detected bug in DMCreateMatrix() for skinny domains with periodic boundary conditions in overestimating nonzero preallocation 4c4762a1bSJed Brown 50e3d61c9SBarry Smith Creation Date : 08-12-2016 6c4762a1bSJed Brown 70e3d61c9SBarry Smith Last Modified : Thu 08 Dec 2016 10:46:02 AM CET 8c4762a1bSJed Brown 90e3d61c9SBarry Smith Created By : Davide Monsorno 10c4762a1bSJed Brown 11c4762a1bSJed Brown _._._._._._._._._._._._._._._._._._._._._.*/ 12c4762a1bSJed Brown 13c4762a1bSJed Brown #include <petscdmda.h> 14c4762a1bSJed Brown 15c4762a1bSJed Brown int main(int argc, char *argv[]) 16c4762a1bSJed Brown { 17c4762a1bSJed Brown PetscInt nx = 2; 18c4762a1bSJed Brown PetscInt ny = 2; 19c4762a1bSJed Brown PetscInt nz = 128; 20c4762a1bSJed Brown DM da; 21c4762a1bSJed Brown Mat A; 22c4762a1bSJed Brown 23*327415f7SBarry Smith PetscFunctionBeginUser; 249566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&argv,NULL,NULL)); 25d0609cedSBarry Smith PetscCall(DMDACreate3d(PETSC_COMM_WORLD,DM_BOUNDARY_PERIODIC,DM_BOUNDARY_PERIODIC,DM_BOUNDARY_GHOSTED,DMDA_STENCIL_BOX,nx,ny,nz, 26d0609cedSBarry Smith PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,1,2,NULL,NULL,NULL,&da)); 279566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(da)); 289566063dSJacob Faibussowitsch PetscCall(DMSetUp(da)); 299566063dSJacob Faibussowitsch PetscCall(DMView(da,PETSC_VIEWER_STDOUT_WORLD)); 309566063dSJacob Faibussowitsch PetscCall(DMCreateMatrix(da,&A)); 31c4762a1bSJed Brown 329566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A)); 339566063dSJacob Faibussowitsch PetscCall(DMDestroy(&da)); 349566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 35b122ec5aSJacob Faibussowitsch return 0; 36c4762a1bSJed Brown } 37c4762a1bSJed Brown 38c4762a1bSJed Brown /*TEST 39c4762a1bSJed Brown 40c4762a1bSJed Brown test: 41c4762a1bSJed Brown nsize: 5 42c4762a1bSJed Brown 43c4762a1bSJed Brown TEST*/ 44