Plex: Unified interpolated submesh creation
- Now both versions create an SF
- Moved DMPlexFilterPoint_Internal()
- DMPlexMarkCohesiveSubmesh_Interpolated() now takes a DMLabel instead of a string

Plex: Allow NULL Label in DMPlexMarkSubmesh_Interpolated()

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c
src/sys/threadcomm/interface/threadcomm.c

show more ...

Merge branch 'master' into pr252/dmay/PetscViewerBinarySkipHeaderVecViewLoad
Manual merge of makefile

Conflicts:
src/sys/classes/viewer/examples/tutorials/makefile

Merge branch 'tisaac/feature-gamg-extensible' into barry/experiment

Conflicts:
src/ksp/pc/impls/gamg/gamg.h

Merge branch 'knepley/feature-combined-tabulation'

* knepley/feature-combined-tabulation: (38 commits)
Plex: Fix compiler warning
Plex+SNES: Fixes for nonzero cStart
Plex: Preserve periodicity

Merge branch 'knepley/feature-combined-tabulation'

* knepley/feature-combined-tabulation: (38 commits)
Plex: Fix compiler warning
Plex+SNES: Fixes for nonzero cStart
Plex: Preserve periodicity when adding ghost cells
Plex: Fixed DMPlexGetGlobalFieldOffset_Private() in optimized mode
Plex+SNES: Offset should be into field components, not basis functions
Plex+SNES: Add guards for optional arrays
DS: Added PetscDSGetComponentOffset(), Added to petsc.h
Plex+TS: FVM RHS computation uses a global F
Plex+FV: Corrected accumulation of fluxes to cells and added time derivative term
Plex: Added DMPlexPointGlobalFieldRead/Ref()
Plex+SNES: Fixed DMPlexGetFaceFields() - It returns field values, not coefficients - FE evaluation is now done
FV: In PetscFVIntegrateRHSFunction(), the arguments are field values, not coefficients
FE: Added EvaluateFaceFields() to FE fields at face centroids
Plex: Added DMPlexPointLocalFieldRead()
FE: Have a standard tabulation at face centroids
Plex: Copy BC for hierarchical refinement
FE: Allow P0 over 1D cells
Plex+FV: Fixed BC with gradient calculation
Plex ex1: Move overlap computation after SetFromOptions
Plex: DMPlexDistribute() no longer calls DMSetFromOptions() - Fixed examples
...

show more ...

Merge branch 'knepley/feature-combined-tabulation' into knepley/feature-ts-fem+fvm-example

* knepley/feature-combined-tabulation: (41 commits)
Plex: Fix distribution of uninterpolated meshes - Doe

Merge branch 'knepley/feature-combined-tabulation' into knepley/feature-ts-fem+fvm-example

* knepley/feature-combined-tabulation: (41 commits)
Plex: Fix distribution of uninterpolated meshes - Does SF stratification by depth, not by dimension - The old code numbered cells after vertices
Plex+SNES: Fixes for nonzero cStart
cannot call PetscObjectGetName() on a subset of the MPI processes
MatView() with some matrix formats and ASCII_MATLAB did not print correct matrix name
Plex: Preserve periodicity when adding ghost cells
Plex: Fixed DMPlexGetGlobalFieldOffset_Private() in optimized mode
Plex+SNES: Offset should be into field components, not basis functions
Plex+SNES: Add guards for optional arrays
DS: Added PetscDSGetComponentOffset(), Added to petsc.h
Plex+TS: FVM RHS computation uses a global F
Plex+FV: Corrected accumulation of fluxes to cells and added time derivative term
Plex: Added DMPlexPointGlobalFieldRead/Ref()
Plex+SNES: Fixed DMPlexGetFaceFields() - It returns field values, not coefficients - FE evaluation is now done
FV: In PetscFVIntegrateRHSFunction(), the arguments are field values, not coefficients
FE: Added EvaluateFaceFields() to FE fields at face centroids
Plex: Added DMPlexPointLocalFieldRead()
FE: Have a standard tabulation at face centroids
Plex: Copy BC for hierarchical refinement
FE: Allow P0 over 1D cells
Plex+FV: Fixed BC with gradient calculation
...

Conflicts:
config/builder.py

show more ...

Plex: Preserve periodicity when adding ghost cells

Plex: Forgot to set block size for coordinates

Doc: Better comment

Merge branch 'knepley/fix-cohesive-intersection'

* knepley/fix-cohesive-intersection:
Plex: Fix leak
Plex: Fix for cohesive labeling - We now allow two cohesive patches to touch directly - This

Merge branch 'knepley/fix-cohesive-intersection'

* knepley/fix-cohesive-intersection:
Plex: Fix leak
Plex: Fix for cohesive labeling - We now allow two cohesive patches to touch directly - This seems to allow hybrid cells to be in the star of fault vertices for other faults, but we have not reproduced this on a small enough example to check - I just filter out hybrid points when labeling

show more ...

Plex: Fix leak

Merge branch 'barry/remove-unneeded-paren-petscmalloc'

Conflicts:
src/mat/impls/baij/seq/baij.c

convert some PetscMalloc() to PetscMalloc1()
remove unneeded () in first argument to PetscMalloc/Calloc1()
fixes for MATLAB examples

Plex: Fix for cohesive labeling
- We now allow two cohesive patches to touch directly
- This seems to allow hybrid cells to be in the star of fault vertices for other faults, but we have not reproduc

Plex: Fix for cohesive labeling
- We now allow two cohesive patches to touch directly
- This seems to allow hybrid cells to be in the star of fault vertices for other faults, but we have not reproduced this on a small enough example to check
- I just filter out hybrid points when labeling

show more ...

Merge branch 'master' into knepley/fix-plex-examples-refinement

* master: (62 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would mean a ton of

Merge branch 'master' into knepley/fix-plex-examples-refinement

* master: (62 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would mean a ton of infrastructure I think
Plex ex12: Fix output after label distribution fix
Plex: Only call SetFromOptions() after Distribute, since we cannot avoid it
TS ex11: Since DMPlexDistribute() calls SetFromOptions(), do not call it again, and cleaned up tests
Plex+FV: If using FVM adjacency, do not create ghost cells for faces in the overlap, or set BC on them
Plex: If using FVM adjacency, we must include the closure in the overlap
Plex+Label: All processes with a valid label should participate in label distribution - Formerly it was just proc 0
Plex: Add check that parallel renumbering fully covers the cone points
Plex: For one level adjacency, include the point itself since this makes sense when looking in parallel from the overlap
Plex: Cannot Malloc2() the arrays passed to PetscSFSetGraph()
Plex: Preallocation fixes for indegree > 1 - Again this should be rewritten
SNES ex62: Changes to output for fixes with quads
Plex+Label: Set default output flag to PETSC_TRUE
Plex ex1: Fixed LaTeX output
Plex ex12: Fixed overlap test
Plex: Fixed incorrect usage of PetscSFGatherBegin/End() in preallocation - This should be rewritten using new style
PetscSF: Incorrect check for computed degrees
Builder: Fix bad search and replace
Configure: Fix for Barry's change to organization - I do not know what happened here
Builder: Fix for Barry's change
...

Conflicts:
src/snes/examples/tutorials/ex12.c

show more ...

Merge branch 'mlange/feature-parallel-distribute'

* mlange/feature-parallel-distribute: (33 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would

Merge branch 'mlange/feature-parallel-distribute'

* mlange/feature-parallel-distribute: (33 commits)
Plex+Label: Fix logic for label distribution - It is restrictive now, but doing it right would mean a ton of infrastructure I think
Plex ex12: Fix output after label distribution fix
Plex: Only call SetFromOptions() after Distribute, since we cannot avoid it
TS ex11: Since DMPlexDistribute() calls SetFromOptions(), do not call it again, and cleaned up tests
Plex+FV: If using FVM adjacency, do not create ghost cells for faces in the overlap, or set BC on them
Plex: If using FVM adjacency, we must include the closure in the overlap
Plex+Label: All processes with a valid label should participate in label distribution - Formerly it was just proc 0
Plex: Add check that parallel renumbering fully covers the cone points
Plex: For one level adjacency, include the point itself since this makes sense when looking in parallel from the overlap
Plex: Cannot Malloc2() the arrays passed to PetscSFSetGraph()
Plex: Preallocation fixes for indegree > 1 - Again this should be rewritten
SNES ex62: Changes to output for fixes with quads
Plex ex1: Fixed LaTeX output
Plex ex12: Fixed overlap test
Plex: Fixed incorrect usage of PetscSFGatherBegin/End() in preallocation - This should be rewritten using new style
PetscSF: Incorrect check for computed degrees
Builder: Fix for Barry's change
PetscSF: Fixes for sparse leaves - When sf->mine exists, the leaves can be numbered non-sequentially. This causes allocation based upon sf->nleaves to underallocate. Replacing sf->nleaves with the maximum leaf number is sufficient in this case.
Plex: Fix leaks
Configure: How did this break?
...

Conflicts:
config/builder.py

show more ...

Plex+FV: If using FVM adjacency, do not create ghost cells for faces in the overlap, or set BC on them

Plex: Cannot Malloc2() the arrays passed to PetscSFSetGraph()

Merge branch 'tisaac/feature-fe-geom-dim' into tisaac/plex-flexible-projection

Use the new PetscFECellGeom

Conflicts:
include/petscdmplex.h
include/petscfe.h
src/dm/impls/plex/plexfem.c

Merge branch 'knepley/fix-submesh-crossedge'

* knepley/fix-submesh-crossedge:
Plex: Do not check hybrid cells when classifying sides of a mesh split - This is not strictly necessary I think, but i

Merge branch 'knepley/fix-submesh-crossedge'

* knepley/fix-submesh-crossedge:
Plex: Do not check hybrid cells when classifying sides of a mesh split - This is not strictly necessary I think, but is true
Plex: Loop edges cannot be cross edges - Loop edges have the same endpoint twice and exist on the endpoints of clamped edges

show more ...

Plex: Do not check hybrid cells when classifying sides of a mesh split
- This is not strictly necessary I think, but is true

Plex: Loop edges cannot be cross edges
- Loop edges have the same endpoint twice and exist on the endpoints of clamped edges

Merge remote-tracking branch 'origin/knepley/fix-plex-submesh-intersection'

* origin/knepley/fix-plex-submesh-intersection:
DMPlex: Fix bug with intersecting submeshes - The label for clamped vert

Merge remote-tracking branch 'origin/knepley/fix-plex-submesh-intersection'

* origin/knepley/fix-plex-submesh-intersection:
DMPlex: Fix bug with intersecting submeshes - The label for clamped vertices for a given submesh can get extraneous points from an adjacent submesh - We filter these out before checking for cross-edges

show more ...

Merge branch 'knepley/feature-dm-dimension-handling'

* knepley/feature-dm-dimension-handling:
DM: Fixed DMDASetDim() in examples
Plex: Fix for isoparametric geometry
Plex: Complex fix
TS ex1

Merge branch 'knepley/feature-dm-dimension-handling'

* knepley/feature-dm-dimension-handling:
DM: Fixed DMDASetDim() in examples
Plex: Fix for isoparametric geometry
Plex: Complex fix
TS ex11: Fixed calls to DMSetCoordinateSection()
DM: FIxed DMClone() to copy dimension
Plex+DA: Now cell geometry routine can take a PetscFE in order to evaluate at quadrature points - We evaluate isoparametrically at the quadrature points - Changed DMDAComputeCellGeometry() --> DMDAComputeCellGeometryFEM() - Changed DMPlexComputeCellGeometry() --> DMPlexComputeCellGeometryFEM() - Added DMPlexComputeCellGeometryAffineFEM()
DM: Added a embedding dimension, different from the topological dimension - This can now be set in DMSetCoordinateSection(), or it can be automatically determined
Plex: Fix output for ex4
DM: Added DMGetDimPoints() - This replicates DMPlex/DAGetDepthStratum(), but uses the dimension which people have said is more intuitive
DM: Move dimension up into DM - Both DA and Plex had this and I think it makes sense

Conflicts:
src/dm/impls/plex/plexcreate.c

show more ...