simpeg.electromagnetics.static.resistivity.simulation.BaseDCSimulation.getJtJdiag#
- BaseDCSimulation.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.
- f
resistivity.fields.FieldsDC,optional Fields solved for all sources.
- m(n_param,)
- Returns:
- (n_param,)
numpy.ndarray The diagonal of \(\mathbf{J^T J}\).
- (n_param,)