1c4762a1bSJed Brown static char help[] = "Tests for periodic mesh output\n\n"; 2c4762a1bSJed Brown 3c4762a1bSJed Brown #include <petscdmplex.h> 4c4762a1bSJed Brown 5*9371c9d4SSatish Balay PetscErrorCode CheckMesh(DM dm) { 65f80ce2aSJacob Faibussowitsch PetscReal detJ, J[9]; 7c4762a1bSJed Brown PetscReal vol; 85f80ce2aSJacob Faibussowitsch PetscInt dim, depth, cStart, cEnd, c; 9c4762a1bSJed Brown 10c4762a1bSJed Brown PetscFunctionBegin; 119566063dSJacob Faibussowitsch PetscCall(DMGetDimension(dm, &dim)); 129566063dSJacob Faibussowitsch PetscCall(DMPlexGetDepth(dm, &depth)); 139566063dSJacob Faibussowitsch PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); 14c4762a1bSJed Brown for (c = cStart; c < cEnd; ++c) { 159566063dSJacob Faibussowitsch PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, NULL, J, NULL, &detJ)); 1663a3b9bcSJacob Faibussowitsch PetscCheck(detJ > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Mesh cell %" PetscInt_FMT " is inverted, |J| = %g", c, (double)detJ); 17c4762a1bSJed Brown if (depth > 1) { 189566063dSJacob Faibussowitsch PetscCall(DMPlexComputeCellGeometryFVM(dm, c, &vol, NULL, NULL)); 1963a3b9bcSJacob Faibussowitsch PetscCheck(vol > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Mesh cell %" PetscInt_FMT " is inverted, vol = %g", c, (double)vol); 20c4762a1bSJed Brown } 21c4762a1bSJed Brown } 22c4762a1bSJed Brown PetscFunctionReturn(0); 23c4762a1bSJed Brown } 24c4762a1bSJed Brown 25*9371c9d4SSatish Balay PetscErrorCode CreateMesh(MPI_Comm comm, DM *dm) { 26c4762a1bSJed Brown PetscFunctionBegin; 279566063dSJacob Faibussowitsch PetscCall(DMCreate(comm, dm)); 289566063dSJacob Faibussowitsch PetscCall(DMSetType(*dm, DMPLEX)); 299566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(*dm)); 309566063dSJacob Faibussowitsch PetscCall(DMViewFromOptions(*dm, NULL, "-dm_view")); 31c4762a1bSJed Brown PetscFunctionReturn(0); 32c4762a1bSJed Brown } 33c4762a1bSJed Brown 34*9371c9d4SSatish Balay int main(int argc, char **argv) { 3530602db0SMatthew G. Knepley DM dm; 36c4762a1bSJed Brown 37327415f7SBarry Smith PetscFunctionBeginUser; 389566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &argv, NULL, help)); 399566063dSJacob Faibussowitsch PetscCall(CreateMesh(PETSC_COMM_WORLD, &dm)); 409566063dSJacob Faibussowitsch PetscCall(CheckMesh(dm)); 419566063dSJacob Faibussowitsch PetscCall(DMDestroy(&dm)); 429566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 43b122ec5aSJacob Faibussowitsch return 0; 44c4762a1bSJed Brown } 45c4762a1bSJed Brown 46c4762a1bSJed Brown /*TEST 47c4762a1bSJed Brown 48c4762a1bSJed Brown test: 49c4762a1bSJed Brown suffix: 0 5030602db0SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 3,1,0 -dm_plex_box_bd periodic,none -dm_view ::ascii_info_detail 51c4762a1bSJed Brown test: 52c4762a1bSJed Brown suffix: 1 53c4762a1bSJed Brown nsize: 2 54e600fa54SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 3,1,0 -dm_plex_box_bd periodic,none -petscpartitioner_type simple -dm_view ::ascii_info_detail 55c4762a1bSJed Brown test: 56c4762a1bSJed Brown suffix: 2 57c4762a1bSJed Brown nsize: 2 58e600fa54SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 6,2,0 -dm_plex_box_bd periodic,none -petscpartitioner_type simple -dm_view ::ascii_info_detail 59c4762a1bSJed Brown test: 60c4762a1bSJed Brown suffix: 3 61c4762a1bSJed Brown nsize: 4 62e600fa54SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 6,2,0 -dm_plex_box_bd periodic,none -petscpartitioner_type simple -dm_view ::ascii_info_detail 63c4762a1bSJed Brown test: 64c4762a1bSJed Brown suffix: 4 65c4762a1bSJed Brown nsize: 2 66e600fa54SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 3,1,0 -dm_plex_box_bd periodic,none -dm_plex_periodic_cut -petscpartitioner_type simple -dm_view ::ascii_info_detail 67c4762a1bSJed Brown test: 68c4762a1bSJed Brown suffix: 5 69c4762a1bSJed Brown nsize: 2 70e600fa54SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 6,2,0 -dm_plex_box_bd periodic,none -dm_plex_periodic_cut -petscpartitioner_type simple -dm_view ::ascii_info_detail 71c4762a1bSJed Brown test: 72c4762a1bSJed Brown suffix: 6 73c4762a1bSJed Brown nsize: 4 74e600fa54SMatthew G. Knepley args: -dm_plex_simplex 0 -dm_plex_box_faces 6,2,0 -dm_plex_box_bd periodic,none -dm_plex_periodic_cut -petscpartitioner_type simple -dm_view ::ascii_info_detail 75c4762a1bSJed Brown 76c4762a1bSJed Brown TEST*/ 77