gies_base⚓︎
EnRML type schemes
GIESMixIn
⚓︎
Bases: Ensemble
This is a base template for implementating the generalized iterative ensemble smoother (GIES) in the following papers: Luo, Xiaodong. "Novel iterative ensemble smoothers derived from a class of generalized cost functions." Computational Geosciences 25.3 (2021): 1159-1189. Luo, Xiaodong, and William C. Cruz. "Data assimilation with soft constraints (DASC) through a generalized iterative ensemble smoother." Computational Geosciences 26.3 (2022): 571-594.
__init__(keys_da, keys_fwd, sim)
⚓︎
The class is initialized by passing the PIPT init. file upwards in the hierarchy to be read and parsed in
pipt.input_output.pipt_init.ReadInitFile.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
init_file
|
PIPT init. file containing info. to run the inversion algorithm |
required |
calc_analysis()
⚓︎
Calculate the update step in LM-EnRML, which is just the Levenberg-Marquardt update algorithm with the sensitivity matrix approximated by the ensemble.
check_convergence()
⚓︎
Check if LM-EnRML have converged based on evaluation of change sizes of objective function, state and damping parameter.
Returns:
| Name | Type | Description |
|---|---|---|
conv |
bool
|
Logic variable telling if algorithm has converged |
why_stop |
dict
|
Dict. with keys corresponding to conv. criteria, with logical variable telling which of them that has been met |