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 15*9371c9d4SSatish Balay int main(int argc, char *argv[]) { 16c4762a1bSJed Brown PetscInt nx = 2; 17c4762a1bSJed Brown PetscInt ny = 2; 18c4762a1bSJed Brown PetscInt nz = 128; 19c4762a1bSJed Brown DM da; 20c4762a1bSJed Brown Mat A; 21c4762a1bSJed Brown 22327415f7SBarry Smith PetscFunctionBeginUser; 239566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &argv, NULL, NULL)); 24*9371c9d4SSatish Balay PetscCall(DMDACreate3d(PETSC_COMM_WORLD, DM_BOUNDARY_PERIODIC, DM_BOUNDARY_PERIODIC, DM_BOUNDARY_GHOSTED, DMDA_STENCIL_BOX, nx, ny, nz, PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, 1, 2, NULL, NULL, NULL, &da)); 259566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(da)); 269566063dSJacob Faibussowitsch PetscCall(DMSetUp(da)); 279566063dSJacob Faibussowitsch PetscCall(DMView(da, PETSC_VIEWER_STDOUT_WORLD)); 289566063dSJacob Faibussowitsch PetscCall(DMCreateMatrix(da, &A)); 29c4762a1bSJed Brown 309566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A)); 319566063dSJacob Faibussowitsch PetscCall(DMDestroy(&da)); 329566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 33b122ec5aSJacob Faibussowitsch return 0; 34c4762a1bSJed Brown } 35c4762a1bSJed Brown 36c4762a1bSJed Brown /*TEST 37c4762a1bSJed Brown 38c4762a1bSJed Brown test: 39c4762a1bSJed Brown nsize: 5 40c4762a1bSJed Brown 41c4762a1bSJed Brown TEST*/ 42