1c4762a1bSJed Brown static char help[] = "Tests for periodic mesh output\n\n"; 2c4762a1bSJed Brown 3c4762a1bSJed Brown #include <petscdmplex.h> 4c4762a1bSJed Brown 59371c9d4SSatish 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 259371c9d4SSatish 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 349371c9d4SSatish 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 71*1fe9daf4SMatthew G. Knepley # This checks that the SF with extra root for periodic cut still checks 72*1fe9daf4SMatthew G. Knepley test: 73*1fe9daf4SMatthew G. Knepley suffix: 5_hdf5 74*1fe9daf4SMatthew G. Knepley requires: hdf5 75*1fe9daf4SMatthew G. Knepley nsize: 2 76*1fe9daf4SMatthew 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 hdf5:mesh.h5 77c4762a1bSJed Brown test: 78c4762a1bSJed Brown suffix: 6 79c4762a1bSJed Brown nsize: 4 80e600fa54SMatthew 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 81c4762a1bSJed Brown 82c4762a1bSJed Brown TEST*/ 83