Do some housekeeping and add missing CHKERRQ calls

Remove memory leaks of XXXRegister() is called but XXX object is never used

This means adding YYYInitialPackage() at the beginning of XXXRegister() routines

Minor fixes of missing prototypes

Commi

Remove memory leaks of XXXRegister() is called but XXX object is never used

This means adding YYYInitialPackage() at the beginning of XXXRegister() routines

Minor fixes of missing prototypes

Commit-type: style-fix
Reported-by: Jed Brown <jed@jedbrown.org>
Thanks-to: Fande Kong <fande.kong@inl.gov>

show more ...

Merge branch 'jed/doc-xsetdm' into maint [PR #1121]

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:

should not pass user contexts through CHKFORTRANNULLOBJECT in Fortran stubs

Commit-type: bug-fix
Reported-by: jenkins

Remove CHKFORTRANNULLOBJECTDEREFERENCE since it is no longer needed with Fortran null equal to C null

Commit-type: feature

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

Use linked list for TSRHSSplit

Funded-by:
Project: PETSc time integration
Time: 4.0 hours
Reported-by:
Thanks-to:

Allow RHSSplit to return sub TS

Funded-by:
Project: PETSc time integration
Time: 6.0 hours
Reported-by:
Thanks-to:

Redesign TS API for split methods

Some explanation of the description on the wikipedia:
- It is correct, but different from the implementations in common codes
- The general framework for the repres

Redesign TS API for split methods

Some explanation of the description on the wikipedia:
- It is correct, but different from the implementations in common codes
- The general framework for the representation orginates from Ronald Ruth's paper on 3rd- and 4th-order symplectic methods
- The order in which velocity and position get updated does not matter; it only leads to variants of the same algorithm

Funded-by:
Project: PETSc time integration
Time: 4.0 hours
Reported-by:
Thanks-to:

show more ...

missing newline in -ts_monitor_extreme

Commit-type: bug-fix
Thanks-to: Satish Balay <balay@mcs.anl.gov>

Add -ts_monitor_extreme

This is useful for monitoring ODEs where the solutions should not grow (for example advection)

Commit-type: feature

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 remote-tracking branch 'origin/maint'

Restore Push/Pop model for ASCII viewers

Log: various fixes and cleanups

- Reorganize declarations in petsclog.h
- Better -log_exclude, now applies to all registered stages
- Minor fixes to implementation of -info_exclude
- PCASM: Log appl

Log: various fixes and cleanups

- Reorganize declarations in petsclog.h
- Better -log_exclude, now applies to all registered stages
- Minor fixes to implementation of -info_exclude
- PCASM: Log apply-on-blocks KSP solves
- Rename PetscLogDestroy() -> PetscLogFinalize() and hide it from the API

show more ...

Move tssen.c into a subsection of TS

Sensitivity related functions should appear in a separate page than TS

Funded-by:
Project: FASTMATH adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:

Merge branch 'maint'

Fix the doc for call-backs drdyf drdpf

Reported-by: Mirko Hahn <hahnm@anl.gov>

Merge branch 'maint'