Coupling propagator

CurrentCoupling6DCurrent

Advances ions and u.

What it solves

FEEC discretization of the following equations:

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

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

time_discret: Crank-Nicolson (implicit mid-point). 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
`ions`	PICVariable	`Particles6D`
`u`	FEECVariable	`HcurlHdivH1vec`

Set these on the propagator in a parameter file, for example CurrentCoupling6DCurrent.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 by `SchurSolver`.
`precond`	`MassMatrixPreconditioner`	`MassMatrixPreconditionerMassMatrixDiagonalPreconditionerNone`	Preconditioner applied to the FEEC mass matrix block.
`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 accumulation/pushing kernels.

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