plex.c (c5e4d11f6a865b7e67dc0aae98421b672bf01406) plex.c (df0420ec0cabd052d76b6e285652119a1edfc9b2)
1#include <petsc/private/dmpleximpl.h> /*I "petscdmplex.h" I*/
2#include <petsc/private/isimpl.h>
3#include <petscsf.h>
4#include <petscds.h>
5
6/* Logging support */
7PetscLogEvent DMPLEX_Interpolate, PETSCPARTITIONER_Partition, DMPLEX_Distribute, DMPLEX_DistributeCones, DMPLEX_DistributeLabels, DMPLEX_DistributeSF, DMPLEX_DistributeOverlap, DMPLEX_DistributeField, DMPLEX_DistributeData, DMPLEX_Migrate, DMPLEX_GlobalToNaturalBegin, DMPLEX_GlobalToNaturalEnd, DMPLEX_NaturalToGlobalBegin, DMPLEX_NaturalToGlobalEnd, DMPLEX_Stratify, DMPLEX_Preallocate, DMPLEX_ResidualFEM, DMPLEX_JacobianFEM, DMPLEX_InterpolatorFEM, DMPLEX_InjectorFEM, DMPLEX_IntegralFEM, DMPLEX_CreateGmsh;
8

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

1896 DMPlexStratify() should be called after all calls to DMPlexSymmetrize()
1897
1898 Level: beginner
1899
1900.seealso: DMPlexCreate(), DMPlexSymmetrize()
1901@*/
1902PetscErrorCode DMPlexStratify(DM dm)
1903{
1#include <petsc/private/dmpleximpl.h> /*I "petscdmplex.h" I*/
2#include <petsc/private/isimpl.h>
3#include <petscsf.h>
4#include <petscds.h>
5
6/* Logging support */
7PetscLogEvent DMPLEX_Interpolate, PETSCPARTITIONER_Partition, DMPLEX_Distribute, DMPLEX_DistributeCones, DMPLEX_DistributeLabels, DMPLEX_DistributeSF, DMPLEX_DistributeOverlap, DMPLEX_DistributeField, DMPLEX_DistributeData, DMPLEX_Migrate, DMPLEX_GlobalToNaturalBegin, DMPLEX_GlobalToNaturalEnd, DMPLEX_NaturalToGlobalBegin, DMPLEX_NaturalToGlobalEnd, DMPLEX_Stratify, DMPLEX_Preallocate, DMPLEX_ResidualFEM, DMPLEX_JacobianFEM, DMPLEX_InterpolatorFEM, DMPLEX_InjectorFEM, DMPLEX_IntegralFEM, DMPLEX_CreateGmsh;
8

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

1896 DMPlexStratify() should be called after all calls to DMPlexSymmetrize()
1897
1898 Level: beginner
1899
1900.seealso: DMPlexCreate(), DMPlexSymmetrize()
1901@*/
1902PetscErrorCode DMPlexStratify(DM dm)
1903{
1904 DM_Plex *mesh = (DM_Plex*) dm->data;
1904 DMLabel label;
1905 PetscInt pStart, pEnd, p;
1906 PetscInt numRoots = 0, numLeaves = 0;
1907 PetscErrorCode ierr;
1908
1909 PetscFunctionBegin;
1910 PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1911 ierr = PetscLogEventBegin(DMPLEX_Stratify,dm,0,0,0);CHKERRQ(ierr);

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

1965 ++level;
1966 ierr = ISDestroy(&pointIS);CHKERRQ(ierr);
1967 ierr = DMLabelGetStratumIS(label, level, &pointIS);CHKERRQ(ierr);
1968 if (pointIS) {ierr = ISGetLocalSize(pointIS, &numPoints);CHKERRQ(ierr);}
1969 else {numPoints = 0;}
1970 }
1971 ierr = ISDestroy(&pointIS);CHKERRQ(ierr);
1972 }
1905 DMLabel label;
1906 PetscInt pStart, pEnd, p;
1907 PetscInt numRoots = 0, numLeaves = 0;
1908 PetscErrorCode ierr;
1909
1910 PetscFunctionBegin;
1911 PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
1912 ierr = PetscLogEventBegin(DMPLEX_Stratify,dm,0,0,0);CHKERRQ(ierr);

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

1966 ++level;
1967 ierr = ISDestroy(&pointIS);CHKERRQ(ierr);
1968 ierr = DMLabelGetStratumIS(label, level, &pointIS);CHKERRQ(ierr);
1969 if (pointIS) {ierr = ISGetLocalSize(pointIS, &numPoints);CHKERRQ(ierr);}
1970 else {numPoints = 0;}
1971 }
1972 ierr = ISDestroy(&pointIS);CHKERRQ(ierr);
1973 }
1974 ierr = DMLabelGetState(label, &mesh->depthState);CHKERRQ(ierr);
1973 ierr = PetscLogEventEnd(DMPLEX_Stratify,dm,0,0,0);CHKERRQ(ierr);
1974 PetscFunctionReturn(0);
1975}
1976
1977#undef __FUNCT__
1978#define __FUNCT__ "DMPlexGetJoin"
1979/*@C
1980 DMPlexGetJoin - Get an array for the join of the set of points

--- 4284 unchanged lines hidden --- 		1975 ierr = PetscLogEventEnd(DMPLEX_Stratify,dm,0,0,0);CHKERRQ(ierr);
1976 PetscFunctionReturn(0);
1977}
1978
1979#undef __FUNCT__
1980#define __FUNCT__ "DMPlexGetJoin"
1981/*@C
1982 DMPlexGetJoin - Get an array for the join of the set of points

--- 4284 unchanged lines hidden ---