simpeg.electromagnetics.time_domain.Simulation3DMagneticField.getAdiag#
- Simulation3DMagneticField.getAdiag(tInd)[source]#
Diagonal system matrix for the given time-step index.
This method returns the diagonal system matrix for the time-step index provided:
\[\mathbf{A}_k = \mathbf{C^T M_{f\rho} C} + \frac{1}{\Delta t_k} \mathbf{M_{e\mu}}\]where
\(\Delta t_k\) is the step length
\(\mathbf{C}\) is the discrete curl operator
\(\mathbf{M_{f \rho}}\) is the resistivity inner-product matrix on faces
\(\mathbf{M_{e\mu}}\) is the permeability inner-product matrix on edges
See the Notes section of the doc strings for
Simulation3DMagneticField
for a full description of the formulation.- Parameters:
- tInd
int
The time-step index; between
[0, n_steps-1]
.
- tInd
- Returns:
- (
n_edges
,n_edges
)sp.sparse.csr_matrix
The diagonal system matrix.
- (