PetscDAEnsembleComputeAnomalies

Forms the state-space anomalies matrix for a PetscDA.

Synopsis

PetscErrorCode PetscDAEnsembleComputeAnomalies(PetscDA da, Vec mean_in, Mat *anomalies_out)

Collective

Input Parameters

• da - the PetscDA context

• mean_in - optional mean state vector (pass NULL to compute internally)

Output Parameter

• anomalies_out - location to store the newly created anomalies matrix

Notes

If mean is NULL, the function will create a temporary vector and compute the ensemble mean using PetscDAEnsembleComputeMean(). If mean is provided, it will be used directly, which can improve performance when the mean has already been computed.

See Also

PetscDA: Data Assimilation, PetscDA, PETSCDAETKF, PETSCDALETKF, PetscDAEnsembleComputeMean()

Level

intermediate

Location

src/ml/da/impls/ensemble/daensemble.c

Examples

src/ml/da/tutorials/ex2.c

---

Index of all PetscDA routines
Table of Contents for all manual pages
Index of all manual pages