Merge branch 'master' into knepley/feature-ex11-simplex

* master: (63 commits)
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Docs+PetscFE: Evaluate() operates in

Merge branch 'master' into knepley/feature-ex11-simplex

* master: (63 commits)
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Docs+PetscFE: Evaluate() operates in reference space
Plex ex7: Added 600 cell test
Plex: Added 600-cell for S^3
Math: Add the golden ratio
Plex ex7: Meshing a sphere
Plex: Changed to DMPlexCreateSphereMesh() - Now can do quads or triangles
Share: Added large 3D mesh file - Its the only one I have gotten to work
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
initialize ap as complained by the nightly tests
Mat: Added example for MatIS
Plex ex3: Update to new point function API
initialize some variables as complained by nightly tests
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
update tests
...

show more ...

Merge branch 'knepley/feature-plex-boundary-values'

* knepley/feature-plex-boundary-values:
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix upd

Merge branch 'knepley/feature-plex-boundary-values'

* knepley/feature-plex-boundary-values:
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix update fund for new API
PetscFE: Derivatives for Point class are all zero
Small fixes
TS ex46: Fix sin and cos
TS ex18: Fix sin and cos
PetscFE: Fixing Point subclass
Plex: Fixed bug with refinement, fixed composed functions for Plex - When we refine, we were calling DMDestroy_Plex(), which wiped out the composed funds, but not calling DMInitialize_Plex() to reinsert them in to the refined mesh
PetscSpacePoint: Fix handling of points
PetscDS: Added update() function for fields - This is intended to directly update the values at points
Plex: DMPlexInsertBoundaryValues() can now be overridden by the user

Conflicts:
include/petsc/private/petscdsimpl.h
include/petscds.h

show more ...

Merge branch 'scott/test-datafiles-fix' of https://bitbucket.org/petsc/petsc into scott/test-datafiles-fix

Merge branch 'master' into knepley/feature-plex-sphere-mesh

* master: (2633 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
ini

Merge branch 'master' into knepley/feature-plex-sphere-mesh

* master: (2633 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
initialize ap as complained by the nightly tests
Plex ex3: Update to new point function API
initialize some variables as complained by nightly tests
SNES ex12: Mix misnaming Nf --> Nc
update tests
add tests
add MatView_SeqXAIJ_ASCII_structonly(); enable A->structure_only in MatSetValuesBlocked_SeqBAIJ
typo
fill-in some references
Bib: Fixed ref
Plex: Added PLY to file support
Plex: Adding PLY support
PetscDS: Turn off Fortran bindings for Constants
Plex tests: Updated to new point function API
SNES ex56: Moved tests to new format
DM: Allow field projection functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields
...

Conflicts:
config/builder.py
include/petscdmplex.h
src/dm/impls/plex/examples/tests/ex1.c
src/dm/impls/plex/plexcreate.c
src/dm/impls/plex/plexgmsh.c
src/snes/examples/tutorials/ex63.c

show more ...

Merge branch 'master' into knepley/feature-convergence-estimator

* master: (27 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables

Merge branch 'master' into knepley/feature-convergence-estimator

* master: (27 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
initialize ap as complained by the nightly tests
Plex ex3: Update to new point function API
initialize some variables as complained by nightly tests
update tests
add tests
add MatView_SeqXAIJ_ASCII_structonly(); enable A->structure_only in MatSetValuesBlocked_SeqBAIJ
Plex: Added PLY to file support
Plex: Adding PLY support
PetscDS: Turn off Fortran bindings for Constants
Plex tests: Updated to new point function API
SNES ex56: Moved tests to new format
DM: Allow field projection functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields
PetscDS: Allow point functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields
dmda matstructure_only
add 'structure_only' to seqbaij matrices -- MatView() does not work yet
PetscFE: Initialize temp array before calling f1
...

Conflicts:
include/petsc/private/petscdsimpl.h
include/petscds.h
src/dm/dt/interface/dtds.c

show more ...

Merge branch 'knepley/feature-fe-kernel-scalars'

* knepley/feature-fe-kernel-scalars:
SNES ex56: Update test output
Plex ex3: Update to new point function API
PetscDS: Turn off Fortran binding

Merge branch 'knepley/feature-fe-kernel-scalars'

* knepley/feature-fe-kernel-scalars:
SNES ex56: Update test output
Plex ex3: Update to new point function API
PetscDS: Turn off Fortran bindings for Constants
Plex tests: Updated to new point function API
SNES ex56: Moved tests to new format
DM: Allow field projection functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields
PetscDS: Allow point functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields

Conflicts:
src/dm/dt/interface/dtfe.c

show more ...

PetscDS: Added exact solution functions
- We will use these to evaluate convergence using MMS

PetscDS: Fix update fund for new API

PetscDS: Allow point functions to take an array of scalars
- This will allow us to pass in user data to kernels without void *, more efficiency than aux fields

Merge branch 'master' into knepley/feature-plex-torus

* master: (174 commits)
Bib: Another ref
Increase patchlevel to 3.7.6
Explicitly listing the LTO library causes difficulties on MacOS so a

Merge branch 'master' into knepley/feature-plex-torus

* master: (174 commits)
Bib: Another ref
Increase patchlevel to 3.7.6
Explicitly listing the LTO library causes difficulties on MacOS so avoid it
TSBDF: Update changes/dev.html
Added Fortran TaoType and TaoLineSearchType definitions
Add Fortran TaoType and TaoLineSearchType definitions
setblocksize for iscol_sub
TSBDF: Use standard TS API and options to control adaptivity
Bib: More refs and a correction
SNESTest: Bug-fix in SNESTest when threshold printing option is turned on.
fix warning and valgrind complains reported by the nightly tests
Section+Doc: Small doc fix
pflotran: switch to using the new upstream git repo
DMPlex: Error if attempt to partition a distributed graph with Chaco
fix complex build error
fix complex build
PetscError() did not properly copy over the specific error message when error type was PETSC_ERROR_INITIAL
TSAdapt: Refactor handling of the always_accept flag
Additional clarification for MatMeshToCellGraph() manual page
optimize MatCreateSubMatrix_MPIAIJ() by merging check of sameXXXDist() into one MPI_AllReduce
...

Conflicts:
config/BuildSystem/config/packages/med.py

show more ...

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

* master: (108 commits)
Bib: Added ref
Bib: Updated reference
Bib: Updated reference
Work around problems with Apple Xcode 8.3 and h

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

* master: (108 commits)
Bib: Added ref
Bib: Updated reference
Bib: Updated reference
Work around problems with Apple Xcode 8.3 and handling of rpath argument with gfortran
help menu update for MatCreateSubMatrices() Reported-by: Brad Aagaard <baagaard@usgs.gov>
fix flops for PCApply_PBJacobi_1 and PCApply_PBJacobi_7 Reported-by: Pierre Jolivet <pierre.jolivet@enseeiht.fr>
TSRK: Fix FSAL optimization to not skip update and check or stages.
Travis-CI: Workaround too long TMPDIR in macOS with Open MPI
test: Fix type mismatch in KSP test.
test: Fix warnings from gfortran -Wconversion (OS X, Homebrew build)
MPIUNI: Silence GCC -Wconversion warnings
Fortran: Add error checking for incorrect PETSC_NULL_XXX argument when expecting function
fix test output per changes in 517c9301b2e9239931bfd644803fd6bb859387ea
MUMPS: update to version 5.1.1 with opensolaris compile fixes
TSRK: Mention 5BS in the changelog.
TSRK: Mention the new 5BS scheme in the manual.
fix 'unused variable' warning
Remove unused code
added hbw_realloc
fixed undefined ptr
...

Conflicts:
config/BuildSystem/config/packages/med.py

show more ...

PetscDS: Added update() function for fields
- This is intended to directly update the values at points

Merge remote-tracking branch 'origin/balay/fix-redundant-decls-warnings-master'

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

fix -Wstrict-prototypes warnings

Reported-by: Till Martens <till.martens@dlr.de>

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

Merge branch 'pr566/knepley/plex_bd_aux_field/master' into knepley/feature-plex-bdfunc

* pr566/knepley/plex_bd_aux_field/master: (22 commits)
PetscDS: reset max nr values for PetscOptionsIntArray(

Merge branch 'pr566/knepley/plex_bd_aux_field/master' into knepley/feature-plex-bdfunc

* pr566/knepley/plex_bd_aux_field/master: (22 commits)
PetscDS: reset max nr values for PetscOptionsIntArray()
SNES+Plex: dm was destroyed too early
Plex: use tabulation based on quadrature points of dual space, not the default tabulation in DMProjectFieldLocal()
SNES ex77: remove zero pointwise functions
PetscFE: only transform pointwise function if it is present
PetscDS: get rid of boundary discretizations
SNES ex77: update test
SNES ex77: simplify residuals/jacobians and delete unnecessary variables and functions
SNES+Plex: use all cell values to interpolate to face
PetscFE: add face quadrature and tabulation
PetscDS: cache number of components and dimensions of discretization space in PetscDS object
SNES ex77: simplify pointwise residual and jacobian for deformation
SNES ex77: move specifications of boundary conditions to command line options
SNES ex77: use auxiliary field for wall pressure
SNES+Plex: allow auxiliary fields in boundary residual/jacobian calculations
PetscDS: get rid of initialization with PetscContainer
SNES+Plex: move boundary jacobian calculation into seperate function
SNES ex77: remove BdResidual for pressure field
SNES+Plex: remove unnecessary (and wrong) loop in boundary residual/jacobian calculations
SNES+Plex: change DMGetField to PetscDSGetBdDiscretization in residual and jacobian calculation
...

Conflicts:
include/petsc/private/petscfeimpl.h
src/dm/dt/interface/dtfe.c
src/dm/impls/plex/plexfem.c
src/snes/examples/tutorials/ex56.c
src/snes/examples/tutorials/ex77.c
src/snes/examples/tutorials/output/ex77_0.out
src/snes/utils/dmplexsnes.c

show more ...

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

PetscDS: get rid of boundary discretizations

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

PetscDS: cache number of components and dimensions of discretization space in PetscDS object

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