simpeg.electromagnetics.frequency_domain.Simulation3DCurrentDensity.getRHS#

Simulation3DCurrentDensity.getRHS(freq)[source]#

Right-hand sides for the given frequency.

This method returns the right-hand sides for the frequency provided. The right-hand side for each source is constructed according to:

$$\mathbf{q}=-i\omega {\mathbf{s}}_{\mathbf{e}}-i\omega \mathbf{C}{\mathbf{M}}_{\mathbf{e}\mu }^{\mathbf{-}\mathbf{1}}{\mathbf{s}}_{\mathbf{m}}$$

where

• $\mathbf{C}$ is the discrete curl operator

• ${\mathbf{s}}_{\mathbf{m}}$ and ${\mathbf{s}}_{\mathbf{e}}$ are the integrated magnetic and electric source terms, respectively

• ${\mathbf{M}}_{\mathbf{e}\mu }$ is the inner-product matrices for permeabilities projected to edges

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

Parameters:

freqfloat

The frequency in Hz.

Returns:

(n_faces, n_sources) numpy.ndarray

The right-hand sides.