SimPEG.electromagnetics.frequency_domain.Simulation3DCurrentDensity#

class SimPEG.electromagnetics.frequency_domain.Simulation3DCurrentDensity(mesh, survey=None, forward_only=False, permittivity=None, **kwargs)[source]#

Bases: BaseFDEMSimulation

We eliminate \(mathbf{h}\) using

\[\mathbf{h} = \frac{1}{i \omega} \mathbf{M_{\mu}^e}^{-1} \left(-\mathbf{C}^{\top} \mathbf{M_{\rho}^f} \mathbf{j} + \mathbf{M^e} \mathbf{s_m} \right)\]

and solve for \(mathbf{j}\) using

\[\left(\mathbf{C} \mathbf{M_{\mu}^e}^{-1} \mathbf{C}^{\top} \mathbf{M_{\rho}^f} + i \omega\right)\mathbf{j} = \mathbf{C} \mathbf{M_{\mu}^e}^{-1} \mathbf{M^e} \mathbf{s_m} - i\omega\mathbf{s_e}\]

Note

This implementation does not yet work with full anisotropy!!

Parameters:

mesh (discretize.base.BaseMesh) – mesh

Attributes

permittivity

Dielectric permittivity (F/m)

Methods

fieldsPair

alias of Fields3DCurrentDensity

getA(freq)

System matrix

getADeriv_rho(freq, u, v[, adjoint])

Product of the derivative of our system matrix with respect to the model and a vector

getRHS(freq)

Right hand side for the system

getRHSDeriv(freq, src, v[, adjoint])

Derivative of the right hand side with respect to the model

getADeriv

getADeriv_mu