Plex: Renamed the functions which insert boundary values and made them public

Merge remote-tracking branch 'origin/knepley/feature-fe-dual-apply' into tisaac/feature-fe-dual-apply

* origin/knepley/feature-fe-dual-apply: (75 commits)
Plex: Screwed up last fix
Plex: Fix ove

Merge remote-tracking branch 'origin/knepley/feature-fe-dual-apply' into tisaac/feature-fe-dual-apply

* origin/knepley/feature-fe-dual-apply: (75 commits)
Plex: Screwed up last fix
Plex: Fix overeager protection of function
PetscFE: Fixed declaration
Travis-CI: Attempt to solve issues with gfortran
Travis-CI: Use newer macOS image
clarify damping and other issues related to SNESLINESEARCHBASIC
change example in new harness to use simple partitioner so gets consistent output on all systems
remove unused function from test file
cleanup formatting in example, plus myrank not rank
all requires need to go on the same line in test definition otherwise they are skipped
missing - makes the manual page formatting wrong
The developer manual clearly states that MPI_Comm_size() returns size it doesn't return numProcs or other nonsense. Pretty soon you'll start using mype instead of rank?
fix typo and cut and paste error in docs
bugfix for !HYPRE
fix warning for ';;'
may switch to scalable only when B->cmap->N > 100000
MatLRC: store in the context the parallel matrices rather than the sequential ones
Bib: Add ref
set default algorithm based on B->cmap->N for MatMatMult_MPIAIJ_MPIAIJ()
Error in VecScatterView_MPI() had the meaning of from and to backwards
...

show more ...

Merge branch 'knepley/feature-plex-tests-new' into knepley/feature-fe-dual-apply

* knepley/feature-plex-tests-new: (73 commits)
Plex: Fix overeager protection of function
Travis-CI: Attempt to s

Merge branch 'knepley/feature-plex-tests-new' into knepley/feature-fe-dual-apply

* knepley/feature-plex-tests-new: (73 commits)
Plex: Fix overeager protection of function
Travis-CI: Attempt to solve issues with gfortran
Travis-CI: Use newer macOS image
clarify damping and other issues related to SNESLINESEARCHBASIC
change example in new harness to use simple partitioner so gets consistent output on all systems
remove unused function from test file
cleanup formatting in example, plus myrank not rank
all requires need to go on the same line in test definition otherwise they are skipped
missing - makes the manual page formatting wrong
The developer manual clearly states that MPI_Comm_size() returns size it doesn't return numProcs or other nonsense. Pretty soon you'll start using mype instead of rank?
fix typo and cut and paste error in docs
bugfix for !HYPRE
fix warning for ';;'
may switch to scalable only when B->cmap->N > 100000
MatLRC: store in the context the parallel matrices rather than the sequential ones
Bib: Add ref
set default algorithm based on B->cmap->N for MatMatMult_MPIAIJ_MPIAIJ()
Error in VecScatterView_MPI() had the meaning of from and to backwards
fix build warning as reported by the nightly tests
Update alquimia.py to use Glenn's branch that updates alquimia for pflotran xsdk
...

show more ...

Merge branch 'master' into knepley/feature-plex-tests-new

* master: (94 commits)
Plex: Fix overeager protection of function
Travis-CI: Attempt to solve issues with gfortran
Travis-CI: Use newe

Merge branch 'master' into knepley/feature-plex-tests-new

* master: (94 commits)
Plex: Fix overeager protection of function
Travis-CI: Attempt to solve issues with gfortran
Travis-CI: Use newer macOS image
clarify damping and other issues related to SNESLINESEARCHBASIC
change example in new harness to use simple partitioner so gets consistent output on all systems
remove unused function from test file
cleanup formatting in example, plus myrank not rank
missing - makes the manual page formatting wrong
The developer manual clearly states that MPI_Comm_size() returns size it doesn't return numProcs or other nonsense. Pretty soon you'll start using mype instead of rank?
fix typo and cut and paste error in docs
bugfix for !HYPRE
fix warning for ';;'
may switch to scalable only when B->cmap->N > 100000
MatLRC: store in the context the parallel matrices rather than the sequential ones
Bib: Add ref
set default algorithm based on B->cmap->N for MatMatMult_MPIAIJ_MPIAIJ()
Error in VecScatterView_MPI() had the meaning of from and to backwards
fix build warning as reported by the nightly tests
Update alquimia.py to use Glenn's branch that updates alquimia for pflotran xsdk
fix divide by zero exposed by 8534686b5096cd4c91a8f3895408b386c0628431 - and produced by petsc4py testsuite
...

Conflicts:
src/dm/impls/plex/examples/tests/ex13.c
src/dm/impls/plex/plexgmsh.c
src/snes/examples/tests/ex2.c
src/snes/examples/tests/output/ex2_3.out
src/snes/examples/tests/output/ex2_4.out
src/snes/examples/tests/output/ex2_5.out

show more ...

Merge branch 'knepley/fix-plex-docs'

* knepley/fix-plex-docs:
Plex: Fix the fix
fix -Wmaybe-uninitialized
Doc+Plex: Add PetscUnit to Fortran
Doc+Plex: Fix function name
Plex: Fix declarati

Merge branch 'knepley/fix-plex-docs'

* knepley/fix-plex-docs:
Plex: Fix the fix
fix -Wmaybe-uninitialized
Doc+Plex: Add PetscUnit to Fortran
Doc+Plex: Fix function name
Plex: Fix declaration
Doc+Plex: Added function docs
Doc+Plex: More function docs - Make functions static - Move internal declarations
Doc+Plex: More function docs
Doc+Plex: Document DMPlexGetOrientedFace() and DMPlexSetSubpointMap() - Make functions static
Doc+Plex: Document DMPlexCreateOverlapMigrationSF() - Make functions static
Plex: Make functions static - Move internal declarations to private header
Plex: Make functions internal
Plex+HDF5: Make functions internal
Plex: Make functions static
Plex: Make functions static
Doc+Plex: Document DMPlexCreateNeighborCSR() - Make functions static
Doc+Plex: Document DMPlexPreallocateOperator() - Make functions static
Plex: Make function static
Plex+VTK: Make functions static
Plex+VTK: Make DMPlexVTKGetCellType() internal

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

show more ...

FE+DT: Now scalar quadrature can be used for vector systems
- You do the obvious thing
- Fixed minor field indexing bugs

Plex: Fixed general interpolation

PetscFE+FV: Now discretization can be vector spaces
- Space and DualSpace have a number of components
- Now the dimension of a space incorporates the components, so all size determinations changed
-

PetscFE+FV: Now discretization can be vector spaces
- Space and DualSpace have a number of components
- Now the dimension of a space incorporates the components, so all size determinations changed
- The default is to form direct sums of scalar spaces, e.g. for constants (1, 0) and (0, 1) so every
basis vector is a vector with only one component nonzero
- DualSpaceApply() is now handled by the implementation
- FE now uses numDof directly from the DualSpace
- Changed SpaceDG --> SpacePoint
- Updated ex3 tests (DA and p4est is still broken)

show more ...

Doc+Plex: More function docs
- Make functions static
- Move internal declarations

Merge branch 'tisaac/dmplex-feature-quadrature-in-geom'

* tisaac/dmplex-feature-quadrature-in-geom:
DMPlex: use new geometry in L2 calculations
DMPlex: fixes to recent geometry calculation chang

Merge branch 'tisaac/dmplex-feature-quadrature-in-geom'

* tisaac/dmplex-feature-quadrature-in-geom:
DMPlex: use new geometry in L2 calculations
DMPlex: fixes to recent geometry calculation changes
DMPlexComputeCellGeometryFEM(): change interface to take quadrature points

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

show more ...

Merge branch 'barry/remove-funct'

Remove the use and definition of __FUNCT__ throughout the code

Since all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETSc
Time: 1.5 hours
Thanks-t

Remove the use and definition of __FUNCT__ throughout the code

Since all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETSc
Time: 1.5 hours
Thanks-to: Andreas Mang <andreas@ices.utexas.edu>

show more ...

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

This is needed to deal with some changes in MatNullSpace from commit 9d2471e

DMPlex: use new geometry in L2 calculations

Merge branch 'master' into tisaac/dmplex-feature-real-to-ref

I need updated PetscFE interfaces in what I'm going to do

* master: (603 commits)
parmetis: fix typo in testsuite
Added PETSC_HAVE_A

Merge branch 'master' into tisaac/dmplex-feature-real-to-ref

I need updated PetscFE interfaces in what I'm going to do

* master: (603 commits)
parmetis: fix typo in testsuite
Added PETSC_HAVE_ATTRIBUTEALIGNED flag to makefile
PCSetUp_SVD() used incorrect matrix for matrix size in parallel, also need to skip LAPACK call if matrix size is zero
strumpack: update to latest version 1.0.4
Remove valgrind warning about MPI_Allreduce() with uninitialized data
sed for test output that has atol convergence with quad precision
updated ex1_nest tests
fixed incorrect key value
fix: bad merge a69d34a
test: move runex78 to TESTEXAMPLES_C_NOTSINGLE
build: fix warnings
sys: update test to produce same result - irrespective of debug/opt build - and presence of ~/.petscrc
sys: src/sys/examples/tests/ex30.c -> src/sys/examples/tests/ex32.c to avoid conflict in master
Sed output, since quad precision reduces below atol on these tests
remove automatically generated matrix names from test output files since can be different
Fixes for MatMPIXXXSetPreallocation(), they were not reentrant, since they kept the information used for matrix-vector products
GAMG: fix PCView_GAMG() format - and update output/ex56_0.out
superlu_dist: 64-bit-indices fix https://github.com/xiaoyeli/superlu_dist/issues/4
Add SNESSetDivergenceTolerance / -snes_divergence_tolerance.
fixed unit read in Plex
...

show more ...

Merge branch 'knepley/feature-plex-bdfunc'

* knepley/feature-plex-bdfunc: (56 commits)
updated snes/ex56
updated snes/ex56
Plex: Goddamn it! Shadowing a variable name
PetscDS: Must be setup

Merge branch 'knepley/feature-plex-bdfunc'

* knepley/feature-plex-bdfunc: (56 commits)
updated snes/ex56
updated snes/ex56
Plex: Goddamn it! Shadowing a variable name
PetscDS: Must be setup in order to access field sizes and components
fix bad merge 44ac993612e1ca111ab14eba1c22039e0ba32a12 wrt src/snes/examples/tutorials/ex56.c [by attempting to redo the merge with kdiff3 - and then comparing the final result with the result in the above merge]
fixed snes/ex56 for PetscDSAddBoundary change
dmplexsnes: Remove more unused variables
Remove unused variables
SNES ex77: Must also set face quadrature
SNES ex12: No longer need boundary discretization orders
SNES ex62: Numerical flutter
SNES ex12: Numerical flutter
Plex: In projection, get Nc from DS instead of calculating it again
PetscFE: Added missing destroy
SNEX ex77: Update boundary args and FE creation
Plex: Update boundary integrals to new API
SNES ex12: Tests of bc from fields
SNES ex12: Added the ability to use finite element fields as boundary conditions
Plex: Need invJ for field projection
Plex: Fix bug with field projection - We must tabulate the FE on the dual basis quadrature points
...

show more ...

Plex: Added field BC values

DM+DS: Allow different kinds of BC functions
- Added DMBoundaryConditionType
- Changed DMGetBoudnary() and DMAddBoundary() and the DS versions

Plex: Fixed name for internal function

Plex: Moved projection code to a new file

Plex: use tabulation based on quadrature points of dual space, not the default tabulation in DMProjectFieldLocal()

SNES+Plex: use all cell values to interpolate to face

Previously the boundary residual/jacobian routines interpolated to face
quadrature points by taking the closure of those faces and using a
separ

SNES+Plex: use all cell values to interpolate to face

Previously the boundary residual/jacobian routines interpolated to face
quadrature points by taking the closure of those faces and using a
separate boundary discretization. This only worked if there were dofs on
faces and the discretization order of the cell and boundary were the
same.
Inverse element jacobians are now calculated from adjacent cells
to also have information on normal derivatives.

show more ...

Merge branch 'master' into barry/fix-successful-numerical-factorization-after-failure/master

The maint patch for the given problem cannot merge cleanly into next due to changes
in the same code in m

Merge branch 'master' into barry/fix-successful-numerical-factorization-after-failure/master

The maint patch for the given problem cannot merge cleanly into next due to changes
in the same code in master (the introduction of MatFactorGetError()) therefor need
to merge in master and use the new construct before merging to next

show more ...

Merge branch 'master' into psanan/pc-telescope-nearnullspace

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddc.c

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcfetidp.c
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcschurs.c
src/mat/impls/is/matis.c
src/mat/impls/is/matis.h

show more ...