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
main(int argc,char * argv[])15*d71ae5a4SJacob Faibussowitsch int main(int argc, char *argv[])
16*d71ae5a4SJacob Faibussowitsch {
17c4762a1bSJed Brown PetscInt nx = 2;
18c4762a1bSJed Brown PetscInt ny = 2;
19c4762a1bSJed Brown PetscInt nz = 128;
20c4762a1bSJed Brown DM da;
21c4762a1bSJed Brown Mat A;
22c4762a1bSJed Brown
23327415f7SBarry Smith PetscFunctionBeginUser;
249566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &argv, NULL, NULL));
259371c9d4SSatish 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));
269566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(da));
279566063dSJacob Faibussowitsch PetscCall(DMSetUp(da));
289566063dSJacob Faibussowitsch PetscCall(DMView(da, PETSC_VIEWER_STDOUT_WORLD));
299566063dSJacob Faibussowitsch PetscCall(DMCreateMatrix(da, &A));
30c4762a1bSJed Brown
319566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A));
329566063dSJacob Faibussowitsch PetscCall(DMDestroy(&da));
339566063dSJacob Faibussowitsch PetscCall(PetscFinalize());
34b122ec5aSJacob Faibussowitsch return 0;
35c4762a1bSJed Brown }
36c4762a1bSJed Brown
37c4762a1bSJed Brown /*TEST
38c4762a1bSJed Brown
39c4762a1bSJed Brown test:
40c4762a1bSJed Brown nsize: 5
41c4762a1bSJed Brown
42c4762a1bSJed Brown TEST*/
43