Field propagator

MaxwellWeakAmpere

Advances e and b.

What it solves

FEEC discretization of the following equations:

find $E \in H (curl)$ and $B \in H (div)$ such that

\int_{Ω} \frac{\partial E}{\partial t} \cdot F d x - \int_{Ω} B \cdot \nabla \times F d x \frac{\partial B}{\partial t} + \nabla \times E = 0, \forall F \in H (curl) = 0 .

Time discretization:

System size reduction via struphy.linear_algebra.schur_solver.SchurSolver.

The variables this step updates over a time step dt, and the discrete spaces each one accepts.

Name	Kind	Accepted spaces
`e`	FEECVariable	`Hcurl`
`b`	FEECVariable	`Hdiv`

Set these on the propagator in a parameter file, for example MaxwellWeakAmpere.Options(algo=...).

Option	Default	Choices	Description
`algo`	`implicit`	`implicitexplicit`	Time stepping scheme to use for the Maxwell equations.
`solver`	`pcg`	`pcgcg`	Name of the symmetric iterative solver passed to `inverse`.
`precond`	`MassMatrixPreconditioner`	`MassMatrixPreconditionerMassMatrixDiagonalPreconditionerNone`	Name of the preconditioner configuration. Currently used only for implicit time stepping.
`solver_params`	`None`	SolverParameters	Iterative-solver controls (for example `tol`, `maxiter`, `verbose`, `info`, `recycle`). If `None`, defaults to `SolverParameters()`.
`butcher`	`None`	ButcherTableau	Butcher tableau for explicit Runge-Kutta methods. Only used when `algo="explicit"`. If `None`, defaults to `ButcherTableau()`.

Each model splits its time step into a sequence of propagators. This one appears in 5 models, at the position shown.