SimPEG.electromagnetics.frequency_domain.Simulation3DCurrentDensity

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

Bases: BaseFDEMSimulation

We eliminate $mathbfh$ using

$$\mathbf{h}=\frac{1}{i\omega }{{\mathbf{M}}_{\mu }^{\mathbf{e}}}^{-1}(-{\mathbf{C}}^{\mathrm{\top }}{\mathbf{M}}_{\rho }^{\mathbf{f}}\mathbf{j}+{\mathbf{M}}^{\mathbf{e}}{\mathbf{s}}_{\mathbf{m}})$$

and solve for $mathbfj$ using

$$(\mathbf{C}{{\mathbf{M}}_{\mu }^{\mathbf{e}}}^{-1}{\mathbf{C}}^{\mathrm{\top }}{\mathbf{M}}_{\rho }^{\mathbf{f}}+i\omega )\mathbf{j}=\mathbf{C}{{\mathbf{M}}_{\mu }^{\mathbf{e}}}^{-1}{\mathbf{M}}^{\mathbf{e}}{\mathbf{s}}_{\mathbf{m}}-i\omega {\mathbf{s}}_{\mathbf{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