simpeg.electromagnetics.static.resistivity.Simulation3DHierarchicalNodal.getJtJdiag#

Simulation3DHierarchicalNodal.getJtJdiag(m, W=None, f=None)[source]#

Return the diagonal of JtJ.

Where \(\mathbf{d}\) are the data, \(\mathbf{m}\) are the model parameters, and the sensitivity matrix is defined as:

\[\mathbf{J} = \dfrac{\partial \mathbf{d}}{\partial \mathbf{m}}\]

this method computes and returns the diagonal elements of \(\mathbf{J^T J}\).

Parameters:
m(n_param,) numpy.ndarray

The model parameters.

W(n_param,) numpy.ndarray, optional

Cell weights.

fresistivity.fields.FieldsDC, optional

Fields solved for all sources.

Returns:
(n_param,) numpy.ndarray

The diagonal of \(\mathbf{J^T J}\).