Use the new API instead of a private function for adjoint Theta

Funded-by:
Project: FASTMath adjoint
Time: 0.5 hours
Reported-by:
Thanks-to: Barry Smith <bsmith@mcs.anl.gov>

Fix SNESTSFormFunction_Theta in order to use FDColoring

VecAXPBYPCZ complains about aliasing vectors in the input argument. This can be seen by
./ex5adj -ts_max_steps 10 -ts_monitor -ts_adjoint_mon

Fix SNESTSFormFunction_Theta in order to use FDColoring

VecAXPBYPCZ complains about aliasing vectors in the input argument. This can be seen by
./ex5adj -ts_max_steps 10 -ts_monitor -ts_adjoint_monitor -ksp_monitor_short -da_grid_x 15 -da_grid_y 15 -snes_fd -implicitform 0 -ts_type cn

Funded-by:
Project: FASTMath adjoint
Time: 1.0 hours
Reported-by:
Thanks-to:

show more ...

Use FD for the adjoint RK

- Works for both -snes_fd and -snes_fd_color.
- Querying for "SNESMatFDColoring" does not work for explicit methods since SNES is not used in the forward run; so check the

Use FD for the adjoint RK

- Works for both -snes_fd and -snes_fd_color.
- Querying for "SNESMatFDColoring" does not work for explicit methods since SNES is not used in the forward run; so check the Jacobian compute function instead.

Funded-by:
Project: FASTMath adjoint
Time: 2.0 hours
Reported-by:
Thanks-to: Barry Smith <bsmith@mcs.anl.gov>

show more ...

Use FD Jacobian in adjoint Theta methods

Funded-by:
Project: FASTMath adjoint
Time: 1.0 hours
Reported-by:
Thanks-to: Barry Smith <bsmith@mcs.anl.gov>

Add the shift variable to the TSTheta context

So the shift can be changed dynamically for SNES function and Jacobian evaluations. Previously it is fixed.

Funded-by:
Project: FASTMath adjoint
Time:

Add the shift variable to the TSTheta context

So the shift can be changed dynamically for SNES function and Jacobian evaluations. Previously it is fixed.

Funded-by:
Project: FASTMath adjoint
Time: 1.0 hours
Reported-by:
Thanks-to: Barry Smith <bsmith@mcs.anl.gov>

show more ...

Initial attempt to use FD Jacobian in adjoint BEuler

-snes_fd_color has more aggressive performance optimizations than -snes_fd; it takes advantage of the fact that a base vector may be ready before

Initial attempt to use FD Jacobian in adjoint BEuler

-snes_fd_color has more aggressive performance optimizations than -snes_fd; it takes advantage of the fact that a base vector may be ready before the Jacobian evaluation and the function evaluation can be skipped for the base vector. This commit makes -snes_fd_colr work for BEuler.

Funded-by:
Project: FASTMath adjoint
Time: 4.0 hours
Reported-by:
Thanks-to: Barry Smith <bsmith@mcs.anl.gov>

show more ...

Merge branch 'hongzh/improve-adjoint-logging' into 'master'

Hongzh/improve adjoint logging

See merge request petsc/petsc!1984

Performance improvement for adjoint BEuler

Adjoint BEuler can use one Jacobian evaluatation per adjoint time step instead of two for explicit ODEs. But the equation type needs to be set by the user.

Performance improvement for adjoint BEuler

Adjoint BEuler can use one Jacobian evaluatation per adjoint time step instead of two for explicit ODEs. But the equation type needs to be set by the user.

Funded-by:
Project: FASTMath adjoint
Time: 1.5 hours
Reported-by:
Thanks-to:

show more ...

DMGetSection: replace all internal usage with DMGetLocalSection

Same for DMSetSection -> DMSetLocalSection. We retain the interface
for DMGetSection/DMSetSection without deprecation warnings (for n

DMGetSection: replace all internal usage with DMGetLocalSection

Same for DMSetSection -> DMSetLocalSection. We retain the interface
for DMGetSection/DMSetSection without deprecation warnings (for now).

Suggested-by: Barry Smith <bsmith@mcs.anl.gov>

show more ...

TSBDF: Add support for FAS

Merge branch 'gropp/doc-fixes-1' [PR #1838]

Merge branch 'maint'

typo in fixing __PETSCSUNDIALS guard

Commit-type: portability-fix, style-fix
Reported-by: nightly tests

Remove use of _ and __ in front of PETSc include guards. Reason: C99 Reserved Identifiers

Commit-type: portability-fix

Merge from doc-fixes branch

Fixes for doctext update