DMProject: If dmAux is a sub mesh, get correct mapped point corresponding to p

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

There has unfortunately been a lot of interplay between our work to
change the petscspace/dualspace to accept vector co

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

There has unfortunately been a lot of interplay between our work to
change the petscspace/dualspace to accept vector components, and the
work to transfer the tests in config/builder.py to the new test harness.
I'm merging them together to get them to play nice with each other, and
so that I only have one feature branch to handler when merging into
master/next.

* tisaac/feature-fe-dual-apply: (64 commits)
snes tutorials ex12: fix restart test arguments
nightly build fixes for arch-viennacl
fix run command to match the command in master/makefile. [the diff will be correct in master/next]
PetscFE: fix typo in for loop declaration
PetscPartitioner: use simple partitioning for regression tests
DT tests ex1: fix test arguments to match regression output
DT ex1: fix options
DMPlexCreateHexBoxMesh: no cells in z direction for 2D
PCBDDC: use %D for PetscInt printfs
Plex: helper function used in all cases now
Plex tests ex7: make partitions portable for testing
Plex: use MPIU_2INT instead of MPI_2INT
MatAIJ KLU: remove redefinition of SuiteSparse_long
DMPlexRefine_CTetgen() mimics DMPlexRefine_Tetgen()
DMPlexGenerate_Tetgen(): mimic ctetgen so it can be used equivalently
Plex: Cylinder fixes for complex
Plex ex1: Test for periodic cylinder
Plex: I think periodic cylinder is working
Plex ex1: Updated torus test
Plex: Torus can now be arbitrarily refined
...

show more ...

Merge branch 'master' into knepley/feature-plex-pragmatic-parallel

* master: (1379 commits)
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes

Merge branch 'master' into knepley/feature-plex-pragmatic-parallel

* master: (1379 commits)
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes
add help menu for MatCreateSubMatricesMPI()
rename struct Mat_SubMat -> Mat_SubSppt
add a little information to the manual about using tchem
rename MatCreateSubMatrixComposite -> MatCreateSubMatrixVirtual; Mat_SubMatrix -> Mat_SubVirtual
need to protect the \ passed on the makefile on my Mac
Bib: Updated ref
update src/mat/examples/tests/output/ex23_#.out
rename MatGetSubMat -> MatCreateSubMat
Mark MatCreateMPIAIJWithSplitArrays() and MatMPIAIJSetPreallocationCSR() with MAT_NO_OFF_PROC_ENTRIES
superlu: libray is installed in PREFIX/lib64 - fix this to use PREFIX/lib
backport fixes cc105cc3 2b35c90a - required for 3bf428d - that was backported from master - 80c96bb
'echo -n' fails on opensolaris/osx - so avoid using it.
fix diffs in nightlybuilds due to options in ~/.petscrc
add Kelley03 to src/docs/tex/petsc.bib
build: determine moudles, mpi.h used during buildtime and print in 'info'
fix analyzer warnings
fix build failure on opensolaris - where -lm from flibs is removed due to -lm in cxxlibs. But then - only flib gets used for some configure tests [blaslapack, parms]
config: prefer __restrict over C99 restrict keyword
...

show more ...

Plex: Choose correct dual space for sub meshes

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

DualSpace: Made Simple implementation multicomponent

DMProjectPoint_Func_Private(): remove bad assertion

Sometimes (non-conforming, wie immer) things get evaluated on faces, we can't assume they are cells

Plex: Need to initialize value array in projection

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

Plex: Allow projection to work for submeshes of lower dimension

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

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

Remove unused variables

Plex: In projection, get Nc from DS instead of calculating it again

Plex: Need invJ for field projection

Plex: Fix bug with field projection
- We must tabulate the FE on the dual basis quadrature points

Plex: Correct handling of u_t in projection
- Should allow X_t to be passed in

DM+KSP: Changed projection functions to accept multiple BC function types
- DMProjectField() now takes the time