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
ProcessOptions(MPI_Comm comm,AppCtx * options)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
CreateMesh(MPI_Comm comm,AppCtx * user,DM * dm)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
main(int argc,char ** argv)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