Merged in knepley/feature-plex-wedge-geometry (pull request #1537)

Knepley/feature plex wedge geometry

Approved-by: BarryFSmith <bsmith@mcs.anl.gov>

Plex: Added support for hybrid (prism) cell geometry
- There are 2 differences from normal cells:
1) The bottom face is inverted, since the normal points toward the top face, instead of outward
2

Plex: Added support for hybrid (prism) cell geometry
- There are 2 differences from normal cells:
1) The bottom face is inverted, since the normal points toward the top face, instead of outward
2) The side faces are ordered as tensor products of intervals, not as quadrilaterals

show more ...

Merge branch 'master' of bitbucket.org:petsc/petsc

Merge branch 'knepley/fix-ds-init' [PR #1415]

* knepley/fix-ds-init:
DM: Fix DS initialization, so that we can setup the discretization multiple times

Revert "Merge branch 'knepley/fix-ds-init' [PR #1415]"

This reverts commit 37bb5b94b6914f9a726e34e5bf9c9c8b40dda9f0, reversing
changes made to 4a21182031185fd5d7353432f9213cf374069ae2.

The reason f

Revert "Merge branch 'knepley/fix-ds-init' [PR #1415]"

This reverts commit 37bb5b94b6914f9a726e34e5bf9c9c8b40dda9f0, reversing
changes made to 4a21182031185fd5d7353432f9213cf374069ae2.

The reason for this revert is an incorrect merge of a local copy of knepley/fix-ds-init.

show more ...

Merge branch 'knepley/fix-ds-init' [PR #1415]

* knepley/fix-ds-init:
Fix DS initialization, so that we can setup the discretization multiple times

Plex: Change slightly the check for higher order coordinates

DM: Fix DS initialization, so that we can setup the discretization multiple times
- Clear out DS when changing DM fields
- Do not error if we cannot find a DS, just return NULL
- If there is no defau

DM: Fix DS initialization, so that we can setup the discretization multiple times
- Clear out DS when changing DM fields
- Do not error if we cannot find a DS, just return NULL
- If there is no default DS, make one in DMCreateDS()
- When checking for a DS in the coordinateDM, allow NULL

show more ...

Plex: Change slightly the check for higher order coordinates

DM: Fix DS initialization, so that we can setup the discretization multiple times
- Clear out DS when changing DM fields
- Do not error if we cannot find a DS, just return NULL
- If there is no defau

DM: Fix DS initialization, so that we can setup the discretization multiple times
- Clear out DS when changing DM fields
- Do not error if we cannot find a DS, just return NULL
- If there is no default DS, make one in DMCreateDS()
- When checking for a DS in the coordinateDM, allow NULL

show more ...

Merge branch 'master' into wence/feature-patch-all-at-once

Merge branch 'knepley/doc-added-seealso-for-pair-of-functions-1550777000537' [PR #1391]

Doc: Added seealso for pair of functions

Merge branch 'master' into jpusztay/feature-swarm-symplectic-example

Merge branch 'master' into jpusztay/feature-swarm-symplectic-example

Merged petsc/petsc into master

Merged petsc/petsc into master

Merge branch 'knepley/feature-multiple-ds' [PR #1309]

Merge branch 'master' into knepley/feature-snes-patch

DM+DS: Reorganized DS handling completely do that multiple regions can be supported
- A region is labeled portion of the mesh which has the same approximation space for the solution. Note that this i

DM+DS: Reorganized DS handling completely do that multiple regions can be supported
- A region is labeled portion of the mesh which has the same approximation space for the solution. Note that this is NOT the same as the label indicating the support of a field. Regions are the refinement of the field supports into unique approximation space pieces. Right now, we only support a single strategy, namely one space on the regular mesh and one on the hybrid cells, but we only have to change DMCreateDS() to extend the capability.
- Now instead of using PetscDSSetDiscretization() and DMSetDS(), you use DMSetField() and DMCreateDS(). Also, you now must use DMCopyDisc(), instead of DMSetDS(), which copies both the field and DS structure to another DM. These changes have been made in all examples.
- DMGetCellDS() returns the approximation space for that cell
- DMGetDS() return the default space since most people will not use regions. There is still a lot of work left t generalize all methods to multiple DSes, but all old code still functions.
- DMGetRegion/RegionNumDS() gets spaces by label and region number
- DMCreateDS() is the setup method which makes the regions and corresponding DSes. With some work, this could be converted to a lazy initialization.
- Cleaned up access to DS since now we can have more than one
- Projection was improved to extract the DS for the particular cells it is iterating over
- The DS are set from options during Section creation. The idea is that people using DS would also use automatic Section creation, or they would manage things themselves and call SetFromOptions manually.
- I updated the subDM and superDM constructors, but they only work for a single region

show more ...

DM: Reorganized DM fields
- A field is defined by a discretization (PetscFE or PetscFV) and a region (DMLabel or NULL)
- This will feed into the creation of PetscDS for every different cell system
-

DM: Reorganized DM fields
- A field is defined by a discretization (PetscFE or PetscFV) and a region (DMLabel or NULL)
- This will feed into the creation of PetscDS for every different cell system
- Changed interface for Get/SetField to add the region

show more ...

Merge branch 'knepley/feature-label-object' [PR #1284]

* knepley/feature-label-object:
DMLabel: Make it a PetscObject

Merge remote-tracking branch 'origin/knepley/feature-label-object' into knepley/feature-section-label

* origin/knepley/feature-label-object:
DMLabel: When distributing, do not need a valid Label f

Merge remote-tracking branch 'origin/knepley/feature-label-object' into knepley/feature-section-label

* origin/knepley/feature-label-object:
DMLabel: When distributing, do not need a valid Label from processes which we are sending to
Plex: Cannot assume label "markers" exists
Forest ex2: Fixed DMLabelCreate()
DMLabel: Remove Fortran binding for DMLabelGetName()
TS ex11: Fixed interface for DMLabelCreate()
DMLabel: Make it a PetscObject - Add a CLASSID - Constructor takes an MPI_Comm - Remove interface that is in Object - Update interface in library - Fixed location of destruct in DMPlexCreateSubmeshGeneric_Interpolated() - Added DMLabelReset()

show more ...

Merge branch 'master' into knepley/fix-plex-partition-parallel

* master: (414 commits)
very small updates to users manual
syntax is define(XXX) in examples tests not defined(XXX)
Increase patc

Merge branch 'master' into knepley/fix-plex-partition-parallel

* master: (414 commits)
very small updates to users manual
syntax is define(XXX) in examples tests not defined(XXX)
Increase patchlevel to 3.10.3
petsc4py: updat to latest maint snapshot with SNESTEST updates
Fix link for CUDA example in PETSc from GPUs.html page Minor fixes to installation instructions
Mat: Fixed types in Pardiso interface
Converted example comment to latex
This example crashes (for unknown reasons, probably related to the handling of extended derived types passed as arguments) with the PGI Fortran compiler
Remove outdated requires: from top of several examples
remove outdated test requires: command from top of source file
Fix the documentation for -ksp_gmres_cgs_refinement_type
Commit f20c2d659b168217cb455989f28afe052d71a64c broke lex.py so the users manual was no longer being built. There reverted the changes to lex.py
Plex: Add PetscPartitioner to fort list Suggested-by: Jiaoyan Li <Jiaoyan.Li@inl.gov>
Doc: Fix spelling Suggested-by: Manuel Colera Rico <m.colera@upm.es>
Mat ex88: check in double or quad precision only
Fix error with PetscComplex when PETSc is compiled with C and used in C++
configure: fix typo endwith -> endswith
Sys ex23: Migrated to ts_max_time in yaml file as well.
AO ex5: Fixed missing CHKERRQ();
Mat test ex88: add tests for MatMultAdd and MatMultTranposeAdd
...

Conflicts:
src/dm/impls/plex/examples/tests/ex12.c

show more ...

DMLabel: Make it a PetscObject
- Add a CLASSID
- Constructor takes an MPI_Comm
- Remove interface that is in Object
- Update interface in library
- Fixed location of destruct in DMPlexCreateSubmeshGe

DMLabel: Make it a PetscObject
- Add a CLASSID
- Constructor takes an MPI_Comm
- Remove interface that is in Object
- Update interface in library
- Fixed location of destruct in DMPlexCreateSubmeshGeneric_Interpolated()
- Added DMLabelReset()

show more ...