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:
- tInd
int
The time index; between
[0, n_steps]
.
- tInd
- Returns:
- (
n_edges
,n_sources
)numpy.ndarray
The right-hand sides.
- (