Field propagator

CurrentCoupling6DDensity

Advances u.

What it solves

FEEC discretization of the following equations:

find $U ~ \in {H (curl), H (div), (H^{1})^{3}}$ such that

\int_{Ω} ρ_{0} \frac{\partial U ~}{\partial t} \cdot V d x = \frac{A _{h}}{A _{b}} \frac{1}{ε} \int_{Ω} n_{h} U ~ \times (B_{0} + B ~) \cdot V d x \forall V \in {H (curl), H (div), (H^{1})^{3}}, n_{h} = \int_{R^{3}} f_{h} d^{3} v .

time_discret: Crank-Nicolson (implicit mid-point).

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

Name	Kind	Accepted spaces
`u`	FEECVariable	`HcurlHdivH1vec`

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

Option	Default	Choices	Description
`u_space`	`Hdiv`	`HcurlHdivH1vec`	FEEC space used for the unknown `u` variable.
`solver`	`pcg`	`pcgcg`	Symmetric iterative solver used for the linear system.
`precond`	`MassMatrixPreconditioner`	`MassMatrixPreconditionerMassMatrixDiagonalPreconditionerNone`	Preconditioner applied to the FEEC mass matrix.
`solver_params`	`None`	SolverParameters	Iterative-solver controls. If `None`, defaults to `SolverParameters()`.
`filter_params`	`None`	FilterParameters	Particle-to-grid filtering parameters used by the accumulator.
`boundary_cut`	`(0.0, 0.0, 0.0)`	tuple	Boundary clipping parameters forwarded to the accumulation kernel.

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