simpeg.regularization.SmoothnessSecondOrder.f_m_deriv#

SmoothnessSecondOrder.f_m_deriv(m)[source]#

Derivative of the regularization kernel function.

For second-order smoothness regularization, the derivative of the regularization kernel function with respect to the model is given by:

$$\frac{\partial {\mathbf{f}}_{\mathbf{m}}}{\partial \mathbf{m}}={\mathbf{L}}_{\mathbf{x}}$$

where ${\mathbf{L}}_{\mathbf{x}}$ is the second-order derivative operator with respect to x.

Parameters:

mnumpy.ndarray

The model.

Returns:

scipy.sparse.csr_matrix

The derivative of the regularization kernel function.

Notes

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

$${\varphi }_m(\mathbf{m})=\Vert \mathbf{W}{\mathbf{L}}_{\mathbf{x}}[\mathbf{m}-{\mathbf{m}}^{(ref)}]{\Vert }^2$$

where $\mathbf{m}$ are the discrete model parameters (model), ${\mathbf{m}}^{(ref)}$ is the reference model, ${\mathbf{L}}_{\mathbf{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}}_{\mathbf{m}}$ as:

$${\mathbf{f}}_{\mathbf{m}}(\mathbf{m})={\mathbf{L}}_{\mathbf{x}}[\mathbf{m}-{\mathbf{m}}^{(ref)}]$$

such that

$${\varphi }_m(\mathbf{m})=\Vert \mathbf{W}{\mathbf{f}}_{\mathbf{m}}{\Vert }^2$$

The derivative of the regularization kernel function with respect to the model is:

$$\frac{\partial {\mathbf{f}}_{\mathbf{m}}}{\partial \mathbf{m}}={\mathbf{L}}_{\mathbf{x}}$$