TS: shrink time step if t0 + dt is greater the final time with matchstep

TS: small bug fix

Fix typos in man pages

TS: Add TSRestartStep() to flag restart the next step

doc: fixes

TS: Move call to TSPostEvaluate() call from TSRollBack() to TSSolve()

Merge remote-tracking branch 'origin/maint'

TS: need to flag the restart of the step if PreStep or PostStep change the solution

TS: Add command line option -ts_monitor_lg_timestep_log

TS: Deprecate TSGetTotalSteps()

Fix extra semicolon after function

TS: Deprecate TSSetInitialTimeStep()

TS: Change default maximum time and steps to unlimited

* Error if users do not specify either final time or max steps
* Remove unnecessary "zero" initializations in TSCreate()
* Add command line opt

TS: Change default maximum time and steps to unlimited

* Error if users do not specify either final time or max steps
* Remove unnecessary "zero" initializations in TSCreate()
* Add command line option -ts_max_time (but keep -ts_final_time)

show more ...

TS: Deprecate TSGetTimeStepNumber()

TS: Deprecate TS{Set|Get}Duration()

TS: Add TSGetExactFinalTime()

TS: Fix handling of step counter

* Proper support for successive calls to TSSolve()
* Reset various counters only if ts->steps is zero (first step)
* Merge ts->steps with ts->total_steps, remove the

TS: Fix handling of step counter

* Proper support for successive calls to TSSolve()
* Reset various counters only if ts->steps is zero (first step)
* Merge ts->steps with ts->total_steps, remove the second
* Add ts->adjoint_steps to track steps in adjoint solve

show more ...

TS: Add TSGetStepNumber() and TSSetStepNumber()

* Remove TSGetTotalSteps()
* Remove TSGetTimeStepNumber()

TS: Add TS{Set|Get}MaxSteps() and TS{Set|Get}MaxTime()

Fix many dead assignments, mostly due to missing CHKERRQ(ierr);

Commit-type: style-fix
Reported-by: clang static analyzer

Fix memory leak for duplicate setup calls

TSSetCostIntegrand() creates new vecs. If called multiple times, there would be memory leak.

Refactor forward SA

Renamed a TS public function (TSAdjointComputeCostIntegrand->TSComputeCostIntegrand); implementations have been simplifed.

Roll back forward sensitivities in TSRollBack

Without this fix, the forward sensitivity (FS) analysis may generate wrong results when a step is rolled back due to TSEvent. Note that we use one rollb

Roll back forward sensitivities in TSRollBack

Without this fix, the forward sensitivity (FS) analysis may generate wrong results when a step is rolled back due to TSEvent. Note that we use one rollback function for both TS and TSForward.

show more ...

add the forward SA capability for Theta methods

-Support the calculation of the gradients of integrals.
-The integrals are treated as additional differential equations.
-Trajectory sensitivities

add the forward SA capability for Theta methods

-Support the calculation of the gradients of integrals.
-The integrals are treated as additional differential equations.
-Trajectory sensitivities to parameters and initial values are handled separatedly.
-Two examples are included: van der Pol and a power grid problem.
-The forward SA modle is implemented in a separate file because ts.c is already a monster and it is still growing.
-TSAdjointComputeDRDYFunction is always required.

show more ...

Jacobians for adjoints for explicit methods and constant Jacobians never got initialized

Since this is the copy of the Jacobian used in the explicit methods it must be initialized
with the Jacobian

Jacobians for adjoints for explicit methods and constant Jacobians never got initialized

Since this is the copy of the Jacobian used in the explicit methods it must be initialized
with the Jacobian values since the values are never reset for the constant Jacobian case.

We really need to refactor the management of the Jacobians, it is overly convoluted with
many corner cases and likely still buggy

Commit-type: bug-fix
Funded-by: IMEX
Project: adjoints
Time: 8 hours
Reported-by: Oana Marin <oanam@mcs.anl.gov>
Development Tools: Emacs

show more ...