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 24c4762a1bSJed Brown ierr = PetscInitialize(&argc,&argv,NULL,NULL);if (ierr) return ierr; 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 29*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetFromOptions(da)); 30*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetUp(da)); 31*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMView(da,PETSC_VIEWER_STDOUT_WORLD)); 32*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateMatrix(da,&A)); 33c4762a1bSJed Brown 34*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatDestroy(&A)); 35*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDestroy(&da)); 36c4762a1bSJed Brown ierr = PetscFinalize(); 37c4762a1bSJed Brown return ierr; 38c4762a1bSJed Brown } 39c4762a1bSJed Brown 40c4762a1bSJed Brown /*TEST 41c4762a1bSJed Brown 42c4762a1bSJed Brown test: 43c4762a1bSJed Brown nsize: 5 44c4762a1bSJed Brown 45c4762a1bSJed Brown TEST*/ 46