DMPlexBuildFromCellSectionParallel#
Build distributed DMPLEX topology from a list of vertices for each cell (common mesh generator output) allowing multiple celltypes
Synopsis#
#include "petscdmplex.h"   
#include "petscdmplextransform.h"   
PetscErrorCode DMPlexBuildFromCellSectionParallel(DM dm, PetscInt numCells, PetscInt numVertices, PetscInt NVertices, PetscSection cellSection, const PetscInt cells[], PetscSF *vertexSF, PetscInt **verticesAdjSaved)
Collective; No Fortran Support
Input Parameters#
- dm - The - DM
- numCells - The number of cells owned by this process 
- numVertices - The number of vertices to be owned by this process, or - PETSC_DECIDE
- NVertices - The global number of vertices, or - PETSC_DETERMINE
- cellSection - The - PetscSectiongiving the number of vertices for each cell (layout of cells)
- cells - An array of the global vertex numbers for each cell 
Output Parameters#
- vertexSF - (Optional) - PetscSFdescribing complete vertex ownership
- verticesAdjSaved - (Optional) vertex adjacency array 
Notes#
A triangle and quadrilateral sharing a face
        2----------3
      / |          |
     /  |          |
    /   |          |
   0  0 |     1    |
    \   |          |
     \  |          |
      \ |          |
        1----------4
would have input
  numCells = 2, numVertices = 5
  cells = [0 1 2  1 4 3 2]
which would result in the DMPLEX
        4----------5
      / |          |
     /  |          |
    /   |          |
   2  0 |     1    |
    \   |          |
     \  |          |
      \ |          |
        3----------6
Vertices are implicitly numbered consecutively 0,…,NVertices.
Each rank owns a chunk of numVertices consecutive vertices.
If numVertices is PETSC_DECIDE, PETSc will distribute them as evenly as possible using PetscLayout.
If NVertices is PETSC_DETERMINE and numVertices is PETSC_DECIDE, NVertices is computed by PETSc as the maximum vertex index in cells + 1.
If only NVertices is PETSC_DETERMINE, it is computed as the sum of numVertices over all ranks.
The cell distribution is arbitrary non-overlapping, independent of the vertex distribution.
See Also#
DMPlex: Unstructured Grids, DM, DMPLEX, DMPlexBuildFromCellListParallel(), DMPlexCreateFromCellSectionParallel(), DMPlexBuildCoordinatesFromCellListParallel(),
PetscSF
Level#
advanced
Location#
src/dm/impls/plex/plexcreate.c
Index of all DMPlex routines
Table of Contents for all manual pages
Index of all manual pages