SimPEG.regularization.SmoothnessSecondOrder.f_m#

SmoothnessSecondOrder.f_m(m)[source]#

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} \big [ \mathbf{m} - \mathbf{m}^{(ref)} \big ]$

where where $$\mathbf{m}$$ are the discrete model parameters (model), $$\mathbf{m}^{(ref)}$$ is the reference model (optional), $$\mathbf{L_x}$$ is the discrete second order x-derivative operator.

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_m (\mathbf{m}) = \frac{1}{2} \Big \| \mathbf{W L_x} \big [ \mathbf{m} - \mathbf{m}^{(ref)} \big ] \Big \|^2$

where $$\mathbf{m}$$ are the discrete model parameters (model), $$\mathbf{m}^{(ref)}$$ is the reference model, $$\mathbf{L_x}$$ is the second-order x-derivative operator, 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} \big [ \mathbf{m} - \mathbf{m}^{(ref)} \big ]$

such that

$\phi_m (\mathbf{m}) = \frac{1}{2} \Big \| \mathbf{W \, f_m} \Big \|^2$