simpeg.electromagnetics.time_domain.Simulation3DElectricField.getRHS#

Simulation3DElectricField.getRHS(tInd)[source]#

Right-hand sides for the given time index.

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

\[\mathbf{q}_k = -\frac{1}{\Delta t_k} \big [ \mathbf{s}_{\mathbf{e}, k} - \mathbf{s}_{\mathbf{e}, k-1} \big ] - \frac{1}{\Delta t_k} \mathbf{C^T M_{f\frac{1}{\mu}}} \big [ \mathbf{s}_{\mathbf{m}, k} - \mathbf{s}_{\mathbf{m}, k-1} \big ]\]

where

  • \(\Delta t_k\) is the step length

  • \(\mathbf{C}\) is the discrete curl operator

  • \(\mathbf{s_m}\) and \(\mathbf{s_e}\) are the integrated magnetic and electric source terms, respectively

  • \(\mathbf{M_{f\frac{1}{\mu}}}\) is the inner-product matrices for inverse permeabilities projected to faces

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

Parameters:
tIndint

The time index; between [0, n_steps].

Returns:
(n_edges, n_sources) numpy.ndarray

The right-hand sides.