simpeg.electromagnetics.time_domain.Simulation3DCurrentDensity.getAdiag#
- Simulation3DCurrentDensity.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 M_{e\mu}^{-1} C^T M_{f\rho}} + \frac{1}{\Delta t_k} \mathbf{I}\]where
\(\Delta t_k\) is the step length
\(\mathbf{I}\) is the identity matrix
\(\mathbf{C}\) is the discrete curl operator
\(\mathbf{M_{e \mu}}\) is the permeability inner-product matrix on faces
\(\mathbf{M_{f \rho}}\) is the permeability inner-product matrix on edges
See the Notes section of the doc strings for
Simulation3DCurrentDensity
for a full description of the formulation.- Parameters:
- tInd
int
The time-step index; between
[0, n_steps-1]
.
- tInd
- Returns:
- (
n_faces
,n_faces
)sp.sparse.csr_matrix
The diagonal system matrix.
- (