Merge branch 'psanan/2021-03-09/ts-sundials2-fix' into 'main'

TSSUNDIALS: update PETSC_HAVE_SUNDIALS to PETSC_HAVE_SUNDIALS2

See merge request petsc/petsc!3692

TSSUNDIALS: update PETSC_HAVE_SUNDIALS to PETSC_HAVE_SUNDIALS2

Merge branch 'danfinn/feature-ts-discrete-gradients' into 'master'

TS: Adding discrete gradients method.

See merge request petsc/petsc!2835

TS: Added new discrete gradients TS
- Added accessors for Discrete Gradient formulation
- Registered TSDISCGRAD

Merged in psanan/man-pages-remove-keywords (pull request #1717)

Man pages: Remove .keywords fields

Approved-by: BarryFSmith <bsmith@mcs.anl.gov>

Man pages: remove .keywords: fields

This field is recognized by doctext from Sowing (http://wgropp.cs.illinois.edu/projects/software/sowing/doctext/node20.htm#Node20),
but it doesn't seem to be unif

Man pages: remove .keywords: fields

This field is recognized by doctext from Sowing (http://wgropp.cs.illinois.edu/projects/software/sowing/doctext/node20.htm#Node20),
but it doesn't seem to be uniformly used or maintained.

Thus, remove all .keywords: fields, and a following blank line, if present.

This is accomplished with GNU sed (gsed on OS X), with the following commands.
*Warning* that this type of command can corrupt a .git directory,
so be cautious in reusing or modifying these commands. They first look
for and delete matching lines with a following line consisting of only whitespace,
and then delete any remaining matching lines.

find src/ -type f -exec gsed -i '/keywords:/ {N; /\n\s*$/d}' {} +
find src/ -type f -exec gsed -i '/keywords:/d' {} +
find include/ -type f -exec gsed -i '/keywords:/ {N; /\n\s*$/d}' {} +
find include/ -type f -exec gsed -i '/keywords:/d' {} +

Hints on the sed command obtained from: https://unix.stackexchange.com/questions/100754/how-to-delete-a-specific-line-and-the-following-blank-line-using-gnu-sed

show more ...

Merged in luzhanghpp/petsc-multirate/luzhang/multirate (pull request #1043)

Multirate time integration methods

Approved-by: BarryFSmith <bsmith@mcs.anl.gov>

Overhaul the MPRK implementation using RHSSplit

add partitioned rk method

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

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

Merged petsc/petsc into master

Merged petsc/petsc into master

Merge from master

Merge branch 'master' into knepley/feature-dm-nearnullspace

* master: (219 commits)
petsc4py: update with VecScatter API changes
remove duplicate code
petsc4py: update to version 3.10.0
fixe

Merge branch 'master' into knepley/feature-dm-nearnullspace

* master: (219 commits)
petsc4py: update with VecScatter API changes
remove duplicate code
petsc4py: update to version 3.10.0
fixed typos in documentation
fixed typos in documentation
Fix some remaining Register()s to call InitializePackage()
Remove memory leaks of XXXRegister() is called but XXX object is never used
Add to PCTFS manual page that it does not work for complex numbers
Stray , in manual page that should be .
Only register MPI3 VecScatter if they are built
clarification of KSPSetComputeInitialGuess()
Add missing makefiles for VecScatter implementations
Changed VecScatterCreate() to take a communicator and return scatter as other first class PETSc objects do. Changed old VecScatterCreate() to VecScatterCreateWithData()
Added VecScatterSetUp() and converted implementations to use it to compute scatter information
Pull "MPI1" and "MPI3" specific code out of common VecScatter impls code This allows future implementations to use the common code also without changing the common code
Seperate VecScatter implementation code into three subdirectories seq, mpi1, and mpi3
doc: fix author list in manual wrt 6402aee
doc: fix LOCDIR wrt broken URL from MatHYPRESetPreallocation.html
ifort: ifx compile error with --with-mpi=0
changes/310.html: merged two TAO sections
...

Conflicts:
src/dm/interface/dm.c

show more ...

Merge branch 'master' into knepley/feature-ksp-ex10-transpose

* master: (219 commits)
petsc4py: update with VecScatter API changes
remove duplicate code
petsc4py: update to version 3.10.0
fi

Merge branch 'master' into knepley/feature-ksp-ex10-transpose

* master: (219 commits)
petsc4py: update with VecScatter API changes
remove duplicate code
petsc4py: update to version 3.10.0
fixed typos in documentation
fixed typos in documentation
Fix some remaining Register()s to call InitializePackage()
Remove memory leaks of XXXRegister() is called but XXX object is never used
Add to PCTFS manual page that it does not work for complex numbers
Stray , in manual page that should be .
Only register MPI3 VecScatter if they are built
clarification of KSPSetComputeInitialGuess()
Add missing makefiles for VecScatter implementations
Changed VecScatterCreate() to take a communicator and return scatter as other first class PETSc objects do. Changed old VecScatterCreate() to VecScatterCreateWithData()
Added VecScatterSetUp() and converted implementations to use it to compute scatter information
Pull "MPI1" and "MPI3" specific code out of common VecScatter impls code This allows future implementations to use the common code also without changing the common code
Seperate VecScatter implementation code into three subdirectories seq, mpi1, and mpi3
doc: fix author list in manual wrt 6402aee
doc: fix LOCDIR wrt broken URL from MatHYPRESetPreallocation.html
ifort: ifx compile error with --with-mpi=0
changes/310.html: merged two TAO sections
...

Conflicts:
src/ksp/ksp/interface/itfunc.c

show more ...

Merged petsc/petsc into master

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

* master: (345 commits)
Update dev notes for the coming release
Increase patchlevel to 3.9.4
doc fixes
Added some AIJSELL tests t

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

* master: (345 commits)
Update dev notes for the coming release
Increase patchlevel to 3.9.4
doc fixes
Added some AIJSELL tests that mirror the SELL tests.
Bib: More refs
Removed use of hardwired -1 to represent PETSC_NULL_XXX object, since the null value is now 0.
fix memleaks flagged by valgrind
ex62: disable long running test and workaround differences in iteration counts on various builds.
Add FAQ we can point to when users get error messages such as: "PETSc was configured with one OpenMPI mpi.h version but now appears to be compiling using a different OpenMPI mpi.h version"
swith test from mpich to openmpi as mpich ch3:sock is unuserable per PetscShmCommGet()
DM: Forgot to move include statement
fix nightlybuild issues - disable hypre [for openmp] and fix uninitialized variables
tests: enable openmp and hwloc with one of the mumps tests
Default BoomerAMG smoothers
Configure: Revert changes from debb90c, so that we are back to the original Libmesh
Configure: Reverted change for Libmesh from 83cd5d7
Split vecimpl.h into vecimpl.h and vecscatterimpl.h
SNES ex62: Turned off some tests for other precisions and shortened output of sensitive tests
Use mmap to allocate the shared pthread barrier if available
Add configure test for functional mmap
...

Conflicts:
include/petscconvest.h
src/snes/examples/tutorials/ex13.c
src/snes/interface/snes.c
src/snes/utils/convest.c

show more ...

Merge branch 'master' into knepley/fix-superdm-is

* master: (123 commits)
Fix a few more Fortran objects to have their NULL value match the C NULL value
remove dead code
Add variable point blo

Merge branch 'master' into knepley/fix-superdm-is

* master: (123 commits)
Fix a few more Fortran objects to have their NULL value match the C NULL value
remove dead code
Add variable point block Jacobi for MPIAIJ matrix.
should not pass user contexts through CHKFORTRANNULLOBJECT in Fortran stubs
fix mpiuni test
Remove CHKFORTRANNULLOBJECTDEREFERENCE since it is no longer needed with Fortran null equal to C null
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Basic as external
Change Fortran null pointer to match C null pointer
one more fix
jenkins: run 'cleantest' after examples are run - to recover diskspace
logical order bug in diagbrdn update fixed, tests updated again
replace new -petscspace_order instances with -petscspace_degree
PETSCSPACETENSOR: clear up naming confusion of subspaces
regression.py: remove configure code related to old test suite - now that SLEPc no longer needs it
PetscSpaceSetup_Poly(): make idempotent
Plex ex3: update tests
Lift changes to defunct dtfe.c from [PR #1074]
PetscSpace_Tensor: clarification of vector-valued
PetscSpace_Tensor: allow height subspaces, fixes tabulation initialization
...

show more ...

Merge remote-tracking branch 'origin/hongzh/rename-bsi'

Rename BSI to BASICSYMPLECTIC

Funded-by:
Project: PETSc time integration
Time: 2.3 hours
Reported-by: Matthew Knepley <knepley@buffalo.edu> Barry Smith <bsmith@mcs.anl.gov>
Thanks-to:

Merge branch 'master' into knepley/feature-pc-patch

* master: (131 commits)
fix mpiuni test
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Ba

Merge branch 'master' into knepley/feature-pc-patch

* master: (131 commits)
fix mpiuni test
PETSCFEOPENCL: declare PetscFEGetTabulation_Basic as external
PETSCFEOPENCL: declare PetscFESetUp_Basic as external
one more fix
jenkins: run 'cleantest' after examples are run - to recover diskspace
logical order bug in diagbrdn update fixed, tests updated again
replace new -petscspace_order instances with -petscspace_degree
PETSCSPACETENSOR: clear up naming confusion of subspaces
regression.py: remove configure code related to old test suite - now that SLEPc no longer needs it
PetscSpaceSetup_Poly(): make idempotent
Plex ex3: update tests
Lift changes to defunct dtfe.c from [PR #1074]
PetscSpace_Tensor: clarification of vector-valued
PetscSpace_Tensor: allow height subspaces, fixes tabulation initialization
PetscSpace: add a dimension member
PetscSpaceSetDegree(): now take maxDegree argument
-petscspacepoly_tensor -> -petscspace_type tensor
PetscSpaceEvaluate_Polynomial(): fix temp array dependencies
PETSCSPACETENSOR: A generic tensor product of spaces.
PetscSpace: -petscspace_order -> -petscspace_degree
...

Conflicts:
src/dm/dt/interface/dtfe.c
src/snes/examples/tutorials/ex13.c
src/snes/examples/tutorials/ex17.c
src/snes/examples/tutorials/ex62.c
src/snes/examples/tutorials/output/ex62_bddc_quad.out
src/snes/examples/tutorials/output/ex62_fetidp_2d_quad.out
src/snes/examples/tutorials/output/ex62_fetidp_2d_tri.out
src/snes/examples/tutorials/output/ex62_fetidp_3d_hex.out
src/snes/examples/tutorials/output/ex62_fetidp_3d_tet.out

show more ...

Merged petsc/petsc into master

Merge branch 'hongzh/symplectic-integrator' [PR #1049]

* hongzh/symplectic-integrator:
This branch implements symplectic time integrators.

We should come up with a unifying framework that works for

Merge branch 'hongzh/symplectic-integrator' [PR #1049]

* hongzh/symplectic-integrator:
This branch implements symplectic time integrators.

We should come up with a unifying framework that works for all multi- methods such as multirate methods (see PR #1043). The idea of this branch is to ask users to provide a two-way splitting of the separable Hamiltonian system.

TSGetRHSFunctionSplit2w(TS ts,Vec *r,TSRHSFunctionSplit2w *fun1,TSRHSFunctionSplit2w *fun2,void **ctx)
TSSetRHSFunctionSplit2w(TS ts,Vec r,TSRHSFunctionSplit2w fun1,TSRHSFunctionSplit2w fun2,void *ctx)

The split/partitioned functions typically have an output vector that has lower DOF. For example, RHSFunction() represents a mapping Rn → Rn , but RHSFunctionSplit2w() provides a mapping Rn → Rn/2 since the subvector is always half of the combined vector. This is also why I use nested vectors for simplicity.

Using IS is a more general way to support multi-way splitting. I am not sure if we should use IS instead of nested vectors for symplectic integrators. Perhaps we can keep both approaches.

show more ...

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

* master: (2931 commits)
remove N in docs for SNESQN
replaced test tolerance with PETSC_SQRT_MACHINE_EPSILON
fortran: PetscComp

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

* master: (2931 commits)
remove N in docs for SNESQN
replaced test tolerance with PETSC_SQRT_MACHINE_EPSILON
fortran: PetscComplex type is missing from fortran
initialize variables so that static analyzer doesn't mistakenly think they are used uninitialized.
Move access of prefix until it is actually used
Remove DMDASplitComm2d(), used in one example that was not tested
Functions that do not return PetscErrorCode should not use PetscFunctionBegin/End paradgm
remove unneeded reassignment of same value
Missing CHKERRQ(); Mostly in Tao code
Remove clang static analyser warning about passing NULL variable to strncat()
missing newline in -ts_monitor_extreme
minor fixes for integer counters
Configure: Deal with '__func__' in C89 builds with newer Clang and GCC.
Jenkinsfile: Increased timeout for tests to 5 minutes.
3.9 changes: further clarify removal of -pc_factor_mat_solver_package
Options deprecations: add "PETSc" to inline warning
removing unused variable
MatShift support added to MATSUBMAT. Mimics MATSHELL functionality. Needed for matrix-free Hessians in TAO Newton-Line-Search methods.
Add -ts_monitor_extreme
Dev manual: add note on deprecating options database keys
...

Conflicts:
src/dm/impls/plex/examples/tests/ex4.c

show more ...