simpeg.electromagnetics.time_domain.Simulation3DMagneticField.getRHS#

Simulation3DMagneticField.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 = \mathbf{C^T M_{f\rho} s}_{\mathbf{e},k} \; - \frac{1}{\Delta t_k} \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\rho}}\) is the inner-product matrix for resisitivites projected to faces

See the Notes section of the doc strings for Simulation3DMagneticField 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.