plexdistribute.c (488f980df2802144bd4c89a2c07257da386f1916) plexdistribute.c (fef8a7d2fedf2f44b3883e5ff4b3df3e69e8cfa7)
1#include <petsc/private/dmpleximpl.h> /*I "petscdmplex.h" I*/
2#include <petsc/private/dmlabelimpl.h> /*I "petscdmlabel.h" I*/
3
4/*@C
5 DMPlexSetAdjacencyUser - Define adjacency in the mesh using a user-provided callback
6
7 Input Parameters:
8+ dm - The DM object

--- 1697 unchanged lines hidden (view full) ---

1706 }
1707
1708 /* Create non-overlapping parallel DM and migrate internal data */
1709 ierr = DMPlexCreate(comm, dmParallel);CHKERRQ(ierr);
1710 ierr = PetscObjectSetName((PetscObject) *dmParallel, "Parallel Mesh");CHKERRQ(ierr);
1711 ierr = DMPlexMigrate(dm, sfMigration, *dmParallel);CHKERRQ(ierr);
1712
1713 /* Build the point SF without overlap */
1#include <petsc/private/dmpleximpl.h> /*I "petscdmplex.h" I*/
2#include <petsc/private/dmlabelimpl.h> /*I "petscdmlabel.h" I*/
3
4/*@C
5 DMPlexSetAdjacencyUser - Define adjacency in the mesh using a user-provided callback
6
7 Input Parameters:
8+ dm - The DM object

--- 1697 unchanged lines hidden (view full) ---

1706 }
1707
1708 /* Create non-overlapping parallel DM and migrate internal data */
1709 ierr = DMPlexCreate(comm, dmParallel);CHKERRQ(ierr);
1710 ierr = PetscObjectSetName((PetscObject) *dmParallel, "Parallel Mesh");CHKERRQ(ierr);
1711 ierr = DMPlexMigrate(dm, sfMigration, *dmParallel);CHKERRQ(ierr);
1712
1713 /* Build the point SF without overlap */
1714 ((DM_Plex*) (*dmParallel)->data)->partitionBalance = ((DM_Plex*) dm->data)->partitionBalance;
1714 ierr = DMPlexCreatePointSF(*dmParallel, sfMigration, PETSC_TRUE, &sfPoint);CHKERRQ(ierr);
1715 ierr = DMSetPointSF(*dmParallel, sfPoint);CHKERRQ(ierr);
1716 ierr = DMGetCoordinateDM(*dmParallel, &dmCoord);CHKERRQ(ierr);
1717 if (dmCoord) {ierr = DMSetPointSF(dmCoord, sfPoint);CHKERRQ(ierr);}
1718 if (flg) {ierr = PetscSFView(sfPoint, PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);}
1719
1720 if (overlap > 0) {
1721 DM dmOverlap;

--- 237 unchanged lines hidden --- 		1715 ierr = DMPlexCreatePointSF(*dmParallel, sfMigration, PETSC_TRUE, &sfPoint);CHKERRQ(ierr);
1716 ierr = DMSetPointSF(*dmParallel, sfPoint);CHKERRQ(ierr);
1717 ierr = DMGetCoordinateDM(*dmParallel, &dmCoord);CHKERRQ(ierr);
1718 if (dmCoord) {ierr = DMSetPointSF(dmCoord, sfPoint);CHKERRQ(ierr);}
1719 if (flg) {ierr = PetscSFView(sfPoint, PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);}
1720
1721 if (overlap > 0) {
1722 DM dmOverlap;

--- 237 unchanged lines hidden ---