1c4762a1bSJed Brown static char help[] = "Tests for coarsening\n\n"; 2c4762a1bSJed Brown 3c4762a1bSJed Brown #include <petscdmplex.h> 4c4762a1bSJed Brown 5c4762a1bSJed Brown typedef struct { 6c4762a1bSJed Brown PetscBool uninterpolate; /* Uninterpolate the mesh at the end */ 7c4762a1bSJed Brown } AppCtx; 8c4762a1bSJed Brown 9d71ae5a4SJacob Faibussowitsch PetscErrorCode ProcessOptions(MPI_Comm comm, AppCtx *options) 10d71ae5a4SJacob Faibussowitsch { 11c4762a1bSJed Brown PetscFunctionBegin; 12c4762a1bSJed Brown options->uninterpolate = PETSC_FALSE; 13c4762a1bSJed Brown 14d0609cedSBarry Smith PetscOptionsBegin(comm, "", "Meshing Problem Options", "DMPLEX"); 159566063dSJacob Faibussowitsch PetscCall(PetscOptionsBool("-uninterpolate", "Uninterpolate the mesh at the end", "ex14.c", options->uninterpolate, &options->uninterpolate, NULL)); 16d0609cedSBarry Smith PetscOptionsEnd(); 17*3ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 18c4762a1bSJed Brown } 19c4762a1bSJed Brown 20d71ae5a4SJacob Faibussowitsch PetscErrorCode CreateMesh(MPI_Comm comm, AppCtx *user, DM *dm) 21d71ae5a4SJacob Faibussowitsch { 22c4762a1bSJed Brown PetscFunctionBegin; 239566063dSJacob Faibussowitsch PetscCall(DMCreate(comm, dm)); 249566063dSJacob Faibussowitsch PetscCall(DMSetType(*dm, DMPLEX)); 259566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(*dm)); 269566063dSJacob Faibussowitsch PetscCall(DMViewFromOptions(*dm, NULL, "-orig_dm_view")); 27c4762a1bSJed Brown if (user->uninterpolate) { 28c4762a1bSJed Brown DM udm = NULL; 29c4762a1bSJed Brown 309566063dSJacob Faibussowitsch PetscCall(DMPlexUninterpolate(*dm, &udm)); 319566063dSJacob Faibussowitsch PetscCall(DMDestroy(dm)); 32c4762a1bSJed Brown *dm = udm; 339566063dSJacob Faibussowitsch PetscCall(DMViewFromOptions(*dm, NULL, "-un_dm_view")); 34c4762a1bSJed Brown } 359566063dSJacob Faibussowitsch PetscCall(DMViewFromOptions(*dm, NULL, "-dm_view")); 36*3ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS); 37c4762a1bSJed Brown } 38c4762a1bSJed Brown 39d71ae5a4SJacob Faibussowitsch int main(int argc, char **argv) 40d71ae5a4SJacob Faibussowitsch { 41c4762a1bSJed Brown DM dm; 42c4762a1bSJed Brown AppCtx user; 43c4762a1bSJed Brown 44327415f7SBarry Smith PetscFunctionBeginUser; 459566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &argv, NULL, help)); 469566063dSJacob Faibussowitsch PetscCall(ProcessOptions(PETSC_COMM_WORLD, &user)); 479566063dSJacob Faibussowitsch PetscCall(CreateMesh(PETSC_COMM_WORLD, &user, &dm)); 489566063dSJacob Faibussowitsch PetscCall(DMDestroy(&dm)); 499566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 50b122ec5aSJacob Faibussowitsch return 0; 51c4762a1bSJed Brown } 52c4762a1bSJed Brown 53c4762a1bSJed Brown /*TEST 54c4762a1bSJed Brown test: 55c4762a1bSJed Brown suffix: 0 56c4762a1bSJed Brown requires: triangle 57c4762a1bSJed Brown args: -dm_view -dm_refine 1 -dm_coarsen -dm_plex_check_symmetry -dm_plex_check_skeleton -dm_plex_check_faces 58c4762a1bSJed Brown 59c4762a1bSJed Brown TEST*/ 60