| #
7e33834a
|
| 10-May-2019 |
Hong Zhang <hongzhang@anl.gov> |
Merged in hongzh/ts-adjoint2 (pull request #1562)
Second-order adjoint ODE solvers
Approved-by: BarryFSmith <bsmith@mcs.anl.gov>
|
| #
9ffb3502
|
| 30-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
TSAdjoint: set TSTrajectory not to use TSHistory
Funded-by:
Project: FASTMath adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:
|
| #
298bade4
|
| 25-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
Add object reference for JacobianP
- Move the setup from TSAdjont and TSForward to TS.
- TSSetRHSJacobianP has to be called before TSSetUp.
Funded-by:
Project: FASTMath adjoint
Time: 1.5 hours
Repo
Add object reference for JacobianP
- Move the setup from TSAdjont and TSForward to TS.
- TSSetRHSJacobianP has to be called before TSSetUp.
Funded-by:
Project: FASTMath adjoint
Time: 1.5 hours
Reported-by: Stefano Zampini
Thanks-to:
show more ...
|
| #
b1e111eb
|
| 24-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
Change names for Hessian functions
Funded-by:
Project: FASTMath adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:
|
| #
369cf35f
|
| 24-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
Add notes for Hessian callbacks
Funded-by:
Project: FASTMath adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:
|
| #
8e224c67
|
| 24-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
Replace PetscStrcmp with PetscObjectTypeCompare
Funded-by:
Project: FASTMath adjoint
Time:
Reported-by: Stefano Zampini
Thanks-to:
|
| #
881c1a9b
|
| 22-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
TSAdjoint: add check for checkpointing settings
- Error out when -ts_trajectory_type basic and -ts_trajectory_solution_only are used together
- The basic type simply saves everything to disk and doe
TSAdjoint: add check for checkpointing settings
- Error out when -ts_trajectory_type basic and -ts_trajectory_solution_only are used together
- The basic type simply saves everything to disk and does not recompute stage values needed by TSAdjoint
- Ugly workaround. Erroring out from TSTrajectory could be better
Funded-by:
Project: FASTMath adjoint
Time: 2.0 hours
Reported-by:
Thanks-to:
show more ...
|
| #
f63bf25f
|
| 21-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
Use TS communicator instead of PETSC_COMM_WORLD
Funded-by:
Project: FASTMath adjoint
Time: 0.5 hours
Reported-by: Barry Smith
Thanks-to:
|
| #
ecf68647
|
| 09-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
Add TSAdjointSetForward and TSAdjointResetForward
|
| #
7207e0fd
|
| 05-Mar-2019 |
Hong Zhang <hongzhang@anl.gov> |
Use quadrature ts to hold tlm sensitivities
- Need guards for integral sensitivity
- Remove unneeded check since TLM of quadrature TS use a different TS
Funded-by:
Project: FASTMath adjoint
Time: 6
Use quadrature ts to hold tlm sensitivities
- Need guards for integral sensitivity
- Remove unneeded check since TLM of quadrature TS use a different TS
Funded-by:
Project: FASTMath adjoint
Time: 6.0 hours
Reported-by:
Thanks-to:
show more ...
|
| #
cd4cee2d
|
| 11-Dec-2018 |
Hong Zhang <hongzhang@anl.gov> |
Refactor quadrature
- Use a separate TS for quadrature
- Initially focus on Theta and Beuler
- Tested on ex3opt and ex9opt
- RK support is not complete
- Fixed memory leak. log_view for ex3sa shows
Refactor quadrature
- Use a separate TS for quadrature
- Initially focus on Theta and Beuler
- Tested on ex3opt and ex9opt
- RK support is not complete
- Fixed memory leak. log_view for ex3sa shows that there is one matrix that is not destroyed.
Funded-by:
Project: FASTMath adjoint
Time: 9.0 hours
Reported-by:
Thanks-to:
show more ...
|
| #
13af1a74
|
| 24-Jul-2018 |
Hong Zhang <hongzhang@anl.gov> |
Add second-order adjoint for ERK
Funded-by:
Project: FASTMath adjoint
Time: 4.0 hours
Reported-by:
Thanks-to:
|
| #
67633408
|
| 12-Apr-2019 |
Hong Zhang <hongzhang@anl.gov> |
TSTheta: hack Hessian computation for explicit ODE
The adjoint of Theta methods does not work if users provide callbacks for an explicit ODE (RHSFunction,RHSJacobian,RHSHessian,etc.).
TSComputeXXX h
TSTheta: hack Hessian computation for explicit ODE
The adjoint of Theta methods does not work if users provide callbacks for an explicit ODE (RHSFunction,RHSJacobian,RHSHessian,etc.).
TSComputeXXX handles the conversion implicitly, which is convoluted.
Funded-by:
Project: FASTMATH adjoint
Time: 1.0 hours
Reported-by:
Thanks-to:
show more ...
|
| #
33f72c5d
|
| 07-May-2018 |
Hong Zhang <hongzhang@anl.gov> |
Refactor TSAdjoint
- Add TSComputeIJacobianP() and TSSetIJacobianP() to pave the way for IMEX adjoint.
- Working on the second-order adjoint to compute parametric sensitivities using backward Euler.
Refactor TSAdjoint
- Add TSComputeIJacobianP() and TSSetIJacobianP() to pave the way for IMEX adjoint.
- Working on the second-order adjoint to compute parametric sensitivities using backward Euler.
- Remove redundant code in TSReset_Theta. TSReset_Theta should not do what TSForwardReset_Theta or TSAdjointReset_Theta do
Funded-by:
Project: FASTMath adjoint
Time: 7.0 hours
Reported-by:
Thanks-to:
Fix TSAdjointInitializeForward for second-order adjoint
To calculate second-order parametric sensitivity, the TLM variable should be initialized as 2*dIdP*dir.
Funded-by:
Project: FASTMath adjoint
Time: 2.0 hours
Reported-by:
Thanks-to:
show more ...
|
| #
d8979ba4
|
| 24-Jul-2018 |
Hong Zhang <hongzhang@anl.gov> |
TSView should not be called in TSAdjointStep
Funded-by:
Project: FASTMATH adjoint
Time: 0.2 hours
Reported-by:
Thanks-to:
|
| #
05c9950e
|
| 25-Apr-2018 |
Hong Zhang <hongzhang@anl.gov> |
Add check for convergence in adjoint run
Funded-by:
Project: FASTMATH adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:
|
| #
3fca9621
|
| 20-Apr-2018 |
Hong Zhang <hongzhang@anl.gov> |
Add a switch for second-order adjoint
One can pass NULL to TSSetCostHessianProducts() for the 3rd argument lambda2 to disable second-order adjoint while not resetting the TSAdjoint context
Funded-b
Add a switch for second-order adjoint
One can pass NULL to TSSetCostHessianProducts() for the 3rd argument lambda2 to disable second-order adjoint while not resetting the TSAdjoint context
Funded-by:
Project: FASTMATH adjoint
Time: 0.5 hours
Reported-by:
Thanks-to:
show more ...
|
| #
cda2db4b
|
| 20-Apr-2018 |
Hong Zhang <hongzhang@anl.gov> |
Separate TSReset(),TSForwardReset() and TSAdjointReset()
Allows the TLM and adjoint context to be removed from the TS context. Use cases are when one wants to disable sensitivity analysis while reus
Separate TSReset(),TSForwardReset() and TSAdjointReset()
Allows the TLM and adjoint context to be removed from the TS context. Use cases are when one wants to disable sensitivity analysis while reusing TS for the simulation.
Funded-by:
Project: FASTMath adjoint
Time: 0.2 hours
Reported-by:
Thanks-to:
show more ...
|
| #
7adebdde
|
| 20-Apr-2018 |
Hong Zhang <hongzhang@anl.gov> |
Add reset function for TSForward
Funded-by:
Project: FASTMath adjoint
Time: 2.0 hours
Reported-by:
Thanks-to:
|
| #
ece46509
|
| 19-Apr-2018 |
Hong Zhang <hongzhang@anl.gov> |
Add reset function for TSAdjoint
This allows TSAdjoint to be turned on/off between different runs using the same TS
Funded-by:
Project: FASTMATH adjoint
Time: 2.0 hours
Reported-by:
Thanks-to:
|
| #
6affb6f8
|
| 05-Apr-2018 |
Hong Zhang <hongzhang@anl.gov> |
Checkpoint the TLM sensitivity for second-order adjoint
Add TSForwardGetStages to acces the stage values for TLM sensitivities
Funded-by:
Project: FASTMath adjoint
Time: 8 hours
Reported-by:
Thanks
Checkpoint the TLM sensitivity for second-order adjoint
Add TSForwardGetStages to acces the stage values for TLM sensitivities
Funded-by:
Project: FASTMath adjoint
Time: 8 hours
Reported-by:
Thanks-to:
show more ...
|
| #
b5b4017a
|
| 08-Mar-2018 |
Hong Zhang <hongzhang@anl.gov> |
Add second-order adjoint backbone
Compiles, but not tested
Funded-by:
Project: FASTMath adjoint
Time: 10.0 hours
Reported-by:
Thanks-to:
Fix the forward setup for endpoint=0 case
Funded-by:
Proje
Add second-order adjoint backbone
Compiles, but not tested
Funded-by:
Project: FASTMath adjoint
Time: 10.0 hours
Reported-by:
Thanks-to:
Fix the forward setup for endpoint=0 case
Funded-by:
Project: FASTMath adjoint
Time:
Reported-by:
Thanks-to:
show more ...
|
| #
c9ad9de0
|
| 08-Mar-2018 |
Hong Zhang <hongzhang@anl.gov> |
Replace deprecated adjoint functions with new ones
Change y to U for consistency
Funded-by:
Project: FASTMATH adjoint
Time: 1.0 hours
Reported-by:
Thanks-to:
|
| #
8000f006
|
| 15-Mar-2019 |
Barry Smith <bsmith@mcs.anl.gov> |
Merge branch 'master' of bitbucket.org:petsc/petsc
|
| #
553225e2
|
| 14-Mar-2019 |
Satish Balay <balay@mcs.anl.gov> |
Merge remote-tracking branch 'origin/jczhang/vecscat-in-sf' [PR #1396]
|