simpeg.regularization.SmoothnessSecondOrder.f_m

SmoothnessSecondOrder.f_m(m)

Evaluate the regularization kernel function.

For second-order smoothness regularization in the x-direction, the regularization kernel function is given by:

\[\mathbf{f_m}(\mathbf{m}) = \mathbf{L_x} \left[ \mu(\mathbf{m}) - \mu(\mathbf{m}^\text{ref}) \right]\]

where \(\mathbf{L_x}\) is the discrete second order x-derivative operator, \(\mathbf{m}\) are the discrete model parameters defined on the mesh, \(\mathbf{m}^{(ref)}\) is the reference model (optional), and \(\mu\) is the mapping function.

Parameters:

mnumpy.ndarray

The model.

Returns:

numpy.ndarray

The regularization kernel function.

Notes

The objective function for second-order smoothness regularization along the x-direction is given by:

\[\phi(\mathbf{m}) = \lVert \mathbf{W} \mathbf{L_x} \left[ \mu(\mathbf{m}) - \mu(\mathbf{m}^\text{ref}) \right] \rVert^2.\]

where \(\mathbf{m}\) are the discrete model parameters (model), \(\mathbf{m}^\text{ref}\) is the reference model, \(\mathbf{L_x}\) is the second-order x-derivative operator, \(\mu\) is the mapping function, and \(\mathbf{W}\) is the weighting matrix. Similar for smoothness along y and z. See the SmoothnessSecondOrder class documentation for more detail.

We define the regularization kernel function \(\mathbf{f_m}\) as:

\[\mathbf{f_m}(\mathbf{m}) = \mathbf{L_x} \left[ \mu(\mathbf{m}) - \mu(\mathbf{m}^\text{ref}) \right]\]

such that

\[\phi_m(\mathbf{m}) = \lVert \mathbf{W \, f_m} \rVert^2.\]