Actual source code: ex50.c

  1: static char help[] = "Test global numbering\n\n";

  3: #include <petscdmplex.h>
  4: #include <petscsf.h>

  6: int main(int argc, char **argv)
  7: {
  8:   DM      dm;
  9:   IS      point_numbering;
 10:   PetscSF point_sf;

 13:   PetscInitialize(&argc, &argv, NULL, help);
 14:   DMCreate(PETSC_COMM_WORLD, &dm);
 15:   DMSetType(dm, DMPLEX);
 16:   DMSetFromOptions(dm);
 17:   DMViewFromOptions(dm, NULL, "-dm_view");

 19:   DMPlexCreatePointNumbering(dm, &point_numbering);
 20:   ISViewFromOptions(point_numbering, NULL, "-point_numbering_view");
 21:   ISDestroy(&point_numbering);

 23:   DMGetPointSF(dm, &point_sf);
 24:   PetscSFViewFromOptions(point_sf, NULL, "-point_sf_view");

 26:   DMDestroy(&dm);
 27:   PetscFinalize();
 28:   return 0;
 29: }

 31: /*TEST

 33:   test:
 34:     nsize: 2
 35:     args: -dm_plex_simplex 0 -dm_plex_box_faces 2,2 -dm_view -point_numbering_view -petscpartitioner_type simple
 36: TEST*/