simpeg.electromagnetics.frequency_domain.Simulation3DMagneticField.getADeriv_mu#

Simulation3DMagneticField.getADeriv_mu(freq, u, v, adjoint=False)[source]#

Permeability derivative operation for the system matrix times a vector.

The system matrix at each frequency is given by:

\[\mathbf{A} = \mathbf{C^T M_{f\rho} C} + i\omega \mathbf{M_{e\mu}}\]

where

  • \(\mathbf{M_{f\rho}}\) is the inner-product matrix for resistivities projected to faces

  • \(\mathbf{M_{e\mu}}\) is the inner-product matrix for permeabilities projected to edges

See the Notes section of the doc strings for Simulation3DMagneticField for a full description of the formulation.

Where \(\mathbf{m}_\boldsymbol{\mu}\) are the set of model parameters defining the permeability, \(\mathbf{v}\) is a vector and \(\mathbf{h}\) is the discrete magnetic field solution, this method assumes the discrete solution is fixed and returns

\[\frac{\partial (\mathbf{A \, h})}{\partial \mathbf{m}_\boldsymbol{\mu}} \, \mathbf{v}\]

Or the adjoint operation

\[\frac{\partial (\mathbf{A \, h})}{\partial \mathbf{m}_\boldsymbol{\mu}}^T \, \mathbf{v}\]
Parameters:
freqfloat

The frequency in Hz.

u(n_edges,) numpy.ndarray

The solution for the fields for the current model at the specified frequency.

vnumpy.ndarray

The vector. (n_param,) for the standard operation. (n_edges,) for the adjoint operation.

adjointbool

Whether to perform the adjoint operation.

Returns:
numpy.ndarray

Derivative of system matrix times a vector. (n_edges,) for the standard operation. (n_param,) for the adjoint operation.