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 PetscErrorCode ierr; 18c4762a1bSJed Brown PetscInt nx = 2; 19c4762a1bSJed Brown PetscInt ny = 2; 20c4762a1bSJed Brown PetscInt nz = 128; 21c4762a1bSJed Brown DM da; 22c4762a1bSJed Brown Mat A; 23c4762a1bSJed Brown 24*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscInitialize(&argc,&argv,NULL,NULL)); 25c4762a1bSJed Brown 26c4762a1bSJed Brown ierr = DMDACreate3d(PETSC_COMM_WORLD,DM_BOUNDARY_PERIODIC,DM_BOUNDARY_PERIODIC,DM_BOUNDARY_GHOSTED,DMDA_STENCIL_BOX,nx,ny,nz, 27c4762a1bSJed Brown PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,1,2,NULL,NULL,NULL,&da);CHKERRQ(ierr); 28c4762a1bSJed Brown 295f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetFromOptions(da)); 305f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetUp(da)); 315f80ce2aSJacob Faibussowitsch CHKERRQ(DMView(da,PETSC_VIEWER_STDOUT_WORLD)); 325f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateMatrix(da,&A)); 33c4762a1bSJed Brown 345f80ce2aSJacob Faibussowitsch CHKERRQ(MatDestroy(&A)); 355f80ce2aSJacob Faibussowitsch CHKERRQ(DMDestroy(&da)); 36*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscFinalize()); 37*b122ec5aSJacob Faibussowitsch return 0; 38c4762a1bSJed Brown } 39c4762a1bSJed Brown 40c4762a1bSJed Brown /*TEST 41c4762a1bSJed Brown 42c4762a1bSJed Brown test: 43c4762a1bSJed Brown nsize: 5 44c4762a1bSJed Brown 45c4762a1bSJed Brown TEST*/ 46