Opened 11 years ago
Closed 11 years ago
#575 closed upgrade (Done)
Added Hessian Drivers
| Reported by: | arango | Owned by: | arango |
|---|---|---|---|
| Priority: | major | Milestone: | Adjoint Based Algorithms |
| Component: | Adjoint | Version: | 3.6 |
| Keywords: | Cc: |
Description
Added new drivers to compute singular vectors and stochastic optimals:
- hessian_op_ocean.h: Hessian singular vectors driver. Use HESSIAN_SV to activate.
- hessian_so_ocean.h: Hessian stochastic optimals driver. Use STOCHASTIC_OPT and HESSIAN_SO to activate.
- so_ocean.h: Full stochastic optimals driver. Use STOCHASTIC_OPT to activate.
The associated propagators are:
- propagator_hop.h: Hessian singular vectors propagator.
- propagator_hso.h: Hessian stochastic optimals propagator
- propagator_so.h: Full stochastic optimals propagator.
The Hessian algorithms require the Lanczos vectors generated by any of the 4D-Var algorithms, which is stored in the adjoint output file. That is, we need to run any of the 4D-Var algorithms first. All these algorithms will be described in a publication in the near future. They are complex and requires a lot of expertise.
Many thanks to Andy Moore for helping us to design and code these algorithms. They are released now because of substantial changes and corrections to the GST drivers. Currently, we are testing these algorithms and tracking couple of parallel issues with the Hessian Stochastic Optimals algorithm.
