adding PetscGaussLobattoLegendre function to dt, removing GLL, updating examples

Plex+FE: Better error message for nonexistent label

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

Merge branch 'knepley/feature-snes-patch' [PR #1332]

* knepley/feature-snes-patch:
Implement SNESPATCH, a SNES for the topological construction of nonlinear relaxation methods.

SNES PATCH: Removed unused FVM-related variables.

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

Merged petsc/petsc into master

Merged petsc/petsc into master

Some cleanups of Matt's FVM code

Remove some of Matt's dead code

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

Fix some Jenkins errors with careful compilers

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 ...

PCPatch: Adding in support for forming residuals

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

Merge branch 'master' into jpustay/feature-swarm-example

Merged petsc/petsc into master

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (393 commits)
Bib: Update reference
Mat: Doc fix
Bib: Updated ref
PetscDS: Doc fixes
PC+LU: Do not try to refactor an

Merge branch 'master' into arcowie-rem/feature-error-logging

* master: (393 commits)
Bib: Update reference
Mat: Doc fix
Bib: Updated ref
PetscDS: Doc fixes
PC+LU: Do not try to refactor an already factored matrix
Mat: Small fix for checking and docs
Mat: Added MatSetFactorType() - Needed it when making a shell matrix look factored
PetscDS: Added PetscDSUpdateBoundary() - Lets the user change the boundary condition
single precision produces different convergence history
p4est: has a dependency on zlib - so handle it correctly
Add -mat_mffd_complex to use the Lyness complex number trick to compute J_u * v instead of differencing.
Replace VecWAXPY by VecAXPY if needed
Revert "Replace VecWAXPY by VecAXPY if needed"
Replace VecWAXPY by VecAXPY if needed
MatHermitianTransposeGetMat and MatCreateVecs
ex19: Updated test filter to avoid false positives.
DM: Improved Global-To-Natural docs Suggested-by: Josh L <ysjosh.lo@gmail.com>
Do not build shared openblas when doing static build
Upgrade OpenBLAS to 0.3.3
test: add alt output file for changes in OSX-10.14 and Xcode-10.0 [perhaps related to ML]
...

Conflicts:
src/vec/is/utils/vsectionis.c

show more ...

Merged petsc/petsc into master

Merge branch 'pr1149/m.diehl/more-modern-fortran-interfaces/master' [PR #1149]

* pr1149/m.diehl/more-modern-fortran-interfaces/master:
This is my first attempt to provide interfaces for modern Fortr

Merge branch 'pr1149/m.diehl/more-modern-fortran-interfaces/master' [PR #1149]

* pr1149/m.diehl/more-modern-fortran-interfaces/master:
This is my first attempt to provide interfaces for modern Fortran.
copied interface from PetscOptionsClearValue
Fortran module files enable type checking at compile time
Fixed wrong intent(in/out) statement
Enables error detection in fortran

show more ...

Merge branch 'knepley/feature-dm-nearnullspace' [PR #1128]

* knepley/feature-dm-nearnullspace:
DM: Added interface for (near) nullspace constructors
DM: Now DMCreateMatrix() adds a (near) nullspace

Merge branch 'knepley/feature-dm-nearnullspace' [PR #1128]

* knepley/feature-dm-nearnullspace:
DM: Added interface for (near) nullspace constructors
DM: Now DMCreateMatrix() adds a (near) nullspace if we have only 1 field and the constructor set
Plex+FEM: Fix corner case where space is smaller than the proposed number of null modes
Plex+FEM: Fix another corner case in CreateRigidBody
SNES ex17: Now use rigid body modes as the near null space

show more ...

Plex+FEM: Must have all prods agree how many rigid body modes to make

Merged petsc/petsc into master

autogeneration should work

the only problematic argument could be MatNullSpace which is a struct