pyacs.lib.glinalg.make_normal_system_tarantola

Form the Tarantola-style normal system (with prior).

pyacs.lib.glinalg.make_normal_system_tarantola.make_normal_system_tarantola(G, d, m0, inv_Cd, inv_Cm)

Form the Tarantola-style normal system (with prior).

Returns N = G.T inv(Cd) G + inv(Cm), b = G.T inv(Cd) d + inv(Cm) m0.

Parameters:

• G (numpy.ndarray) – Design matrix.

• d (array_like) – Observation vector.

• m0 (array_like or float) – Prior model.

• inv_Cd (numpy.ndarray) – Inverse data covariance.

• inv_Cm (numpy.ndarray or float) – Inverse prior covariance.

Returns:

(N, b) normal matrix and right-hand side.

Return type:

tuple