PetscDAEnsembleApplyTInverse#
Apply T^{-1} to a vector [Alg 6.4 line 8]
Synopsis#
PetscErrorCode PetscDAEnsembleApplyTInverse(PetscDA da, Vec sdel, Vec w)
Collective
Input Parameters#
da - the
PetscDAcontextsdel - input vector S^T-delta
Output Parameter#
w - output vector w = T^{-1} * sdel
Notes#
This function applies the inverse of \(T = (1/\rho) I + S^T S\) (with \(\rho\) the inflation factor) using the stored symmetric eigendecomposition: \(T^{-1} = V D^{-1} V^T\).
See Also#
PetscDA: Data Assimilation, PetscDA, PETSCDALETKF, PetscDAEnsembleTFactor(), PetscDAEnsembleApplySqrtTInverse()
Level#
advanced
Location#
src/ml/da/impls/ensemble/daensemble.c
Index of all PetscDA routines
Table of Contents for all manual pages
Index of all manual pages