1c4762a1bSJed Brown /* -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-. 2c4762a1bSJed Brown 3c4762a1bSJed Brown Detected bug in DMCreateMatrix() for skinny domains with periodic boundary conditions in overestimating nonzero preallocation 4c4762a1bSJed Brown 5*0e3d61c9SBarry Smith Creation Date : 08-12-2016 6c4762a1bSJed Brown 7*0e3d61c9SBarry Smith Last Modified : Thu 08 Dec 2016 10:46:02 AM CET 8c4762a1bSJed Brown 9*0e3d61c9SBarry 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 29c4762a1bSJed Brown ierr = DMSetFromOptions(da);CHKERRQ(ierr); 30c4762a1bSJed Brown ierr = DMSetUp(da);CHKERRQ(ierr); 31c4762a1bSJed Brown ierr = DMView(da,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr); 32c4762a1bSJed Brown ierr = DMCreateMatrix(da,&A);CHKERRQ(ierr); 33c4762a1bSJed Brown 34c4762a1bSJed Brown ierr = MatDestroy(&A);CHKERRQ(ierr); 35c4762a1bSJed Brown ierr = DMDestroy(&da);CHKERRQ(ierr); 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