refactor bsi and add the stormer verlet method

wikipedia is wrong about the order of updating p and p

Add the first symplectic integrator SIEuler

The interface is designed particularly for separable Hamiltonian systems

Merge branch 'master' into knepley/feature-snes-ex13-spectral

* master: (2350 commits)
migrate warning fix from next-tmp
Revert "Remove dead code."
clarifying comment in the code
Remove dead

Merge branch 'master' into knepley/feature-snes-ex13-spectral

* master: (2350 commits)
migrate warning fix from next-tmp
Revert "Remove dead code."
clarifying comment in the code
Remove dead code.
Add JuMP citation
Copyedits and updates to developers manual
Add more protection to intrinsics
MPI: protect arguments in MPI macros
Add the missing conditionals for 64-bit build
Plex+VTK: Fix field type determination
Fix bug in name cleanup
Bib: Added ref
Plex+VTK: Fix field type determination
Plex: Updated test output - I am now hating this shit
forgot to save file after moving prototypes to petscsys.h
update Apple example project settings per Xcode
Two functions should actually be public, used, for example from Apple MacOS and iOS programs
includes should always be included relative to include directory
Apple removed old locations of docsets
SNES ex12: Updated test output
...

Conflicts:
src/snes/examples/tutorials/ex13.c

show more ...

Merge branch 'barry/correct-radau5'

Start implementation of Radau5 interface

Reproduce Radau5 result for simple test case

Commit-type: feature

Merge commit 'd47bf9aaf1e5266cc3f0ff499f934c85788965a9' into jed/fix-matcreatempibaij/maint

Obtain changed (v3.5) handling of MatMPIBAIJSetPreallocationCSR to match
documentation.

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

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 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 'stefano_zampini/feature-convert-matnest-matis' into stefano_zampini/feature-pcbddc-saddlepoint

I have to merge it here, since there's an API addition for MatIS in stefano_zampini/featu

Merge branch 'stefano_zampini/feature-convert-matnest-matis' into stefano_zampini/feature-pcbddc-saddlepoint

I have to merge it here, since there's an API addition for MatIS in stefano_zampini/feature-pcbddc-saddlepoint
that prevents me to smoothly merge stefano_zampini/feature-convert-matnest-matis into next.

This will also facilitate merging stefano_zampini/feature-pcbddc-saddlepoint into master

Conflicts:
src/docs/website/documentation/changes/dev.html

show more ...

Merge branch 'emil/ts-glee'

TSGLEE: Implemented the General Linear methods with Error Estimation

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 mark/snes-ex56

* master: (93 commits)
Web: Added a tutorial
typo
rename X -> U
sowing: update to version that creates manpage URLs for 'enum' types/values
Increa

Merge branch 'master' into mark/snes-ex56

* master: (93 commits)
Web: Added a tutorial
typo
rename X -> U
sowing: update to version that creates manpage URLs for 'enum' types/values
Increase patchlevel to 3.7.3
move runex3_3 to NOTSINGLE [ partial backport 9ed0bd76d5c329eb91401ebabfbd851446214e20 ]
DMPlex: local coordinate vector should be a sequential vector
DM: Name coordinate vectors in subdomain and restrict hooks
tests: disable a couple of TS examples from single precision tests
reduce TS ex12_2 output
fix PetscFormatRealArray() to cast PetscReal in print statement to double so works for __float128 precision
Vec: enable Fortran interfaces for VecGetLocalVectorRead and VecRestoreLocalVectorRead
Bib: Added refs
PCShell: man page fixes
update handling of latex in example code to use https
fix typo in isltog comment
test: Avoid conversion of object <-> function pointer type [-Wpedantic]
fortran: enable interface for VecGetLocalVector(), VecRestoreLocalVector()
update citation for cmz2016
PCBDDC: silent compilation warnings from nightlybuilds
...

# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.

show more ...

Merge branch 'emil/rename-gl-to-glle'

TS: Changed TSGL to TSGLLE because TSGL is too general and we want to emphasize that TSGLLE is a type of GL with local error estimation.

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (180 commits)
Bib: Added ref
-random_view caused error because PetscRandomView() was not provided in PetscRandomCreate()

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (180 commits)
Bib: Added ref
-random_view caused error because PetscRandomView() was not provided in PetscRandomCreate()
Update the xSDKTrilinos test message to use interactive session and ctest
gcc6: fix comple warnings with -Wmisleading-indentation
fix compile warning
bugfix for MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable()
minor formatting cleanup of gamg.c
add ml example to tutorials page for xsdk
bugfix of MatDuplicate for mpiaij matrix product Reported-by: Miorelli, Federico <Federico.Miorelli@cgg.com>
The deprecated -log_summary_exclude was still being used when it should have been -log_exclude
cpardiso: requires mkl_blacs so check for that. Also common code to detect MKL for both pardiso and cpardiso
cpardiso: fix segv and memory leaks
resume MatGetSubMatrix_SeqSBAIJ() because it cannot share routine with BAIJ due to cast
fix bad relative links
replace MatSeqBAIJSetPreallocation_SeqBAIJ with MatSeqBAIJSetPreallocation to skip this call when input mat is sbaij type
fix a variety of set but not used variables in examples detected by clang static analyzer
removed declared but not used variable
ANL Report formatting requires listing authors by affiliation
DMP4est: uncrustify
DMP4est: remove memset()
...

show more ...

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

Conflicts:
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcgraph.c
src/ksp/pc/impls/bddc/bddcprivate.c

Merge branch 'maint'

Merge branch 'dalcinl/tsbdf' into maint

Merge branch 'dalcinl/ts-2nd-order-systems' into maint

TSALPHA2: New implementation for second order systems

* Support for event detection.
* Support for time adaptivity based on approximating the local truncation
error using backward difference formu

TSALPHA2: New implementation for second order systems

* Support for event detection.
* Support for time adaptivity based on approximating the local truncation
error using backward difference formulas.
* Use an step-halving approach in the initial step to compute a consistent
initial state as well as approximating the local truncation error.

show more ...