1 /* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. 2 3 Detected bug in DMCreateMatrix() for skinny domains with periodic boundary conditions in overestimating nonzero preallocation 4 5 Creation Date : 08-12-2016 6 7 Last Modified : Thu 08 Dec 2016 10:46:02 AM CET 8 9 Created By : Davide Monsorno 10 11 _._._._._._._._._._._._._._._._._._._._._.*/ 12 13 #include <petscdmda.h> 14 15 int main(int argc, char *argv[]) 16 { 17 PetscInt nx = 2; 18 PetscInt ny = 2; 19 PetscInt nz = 128; 20 DM da; 21 Mat A; 22 23 PetscFunctionBeginUser; 24 PetscCall(PetscInitialize(&argc,&argv,NULL,NULL)); 25 PetscCall(DMDACreate3d(PETSC_COMM_WORLD,DM_BOUNDARY_PERIODIC,DM_BOUNDARY_PERIODIC,DM_BOUNDARY_GHOSTED,DMDA_STENCIL_BOX,nx,ny,nz, 26 PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,1,2,NULL,NULL,NULL,&da)); 27 PetscCall(DMSetFromOptions(da)); 28 PetscCall(DMSetUp(da)); 29 PetscCall(DMView(da,PETSC_VIEWER_STDOUT_WORLD)); 30 PetscCall(DMCreateMatrix(da,&A)); 31 32 PetscCall(MatDestroy(&A)); 33 PetscCall(DMDestroy(&da)); 34 PetscCall(PetscFinalize()); 35 return 0; 36 } 37 38 /*TEST 39 40 test: 41 nsize: 5 42 43 TEST*/ 44