Struct IterativeSystem 

pub struct IterativeSystem { /* private fields */ }

The assembled, frequency-independent, matrix-free mesh system of a geometry and its ports. See the module documentation.

Implementations

impl IterativeSystem

pub fn assemble( geometry: &Geometry, ports: &[Port], discretization: &Discretization, params: &IterativeParams, ) -> Result<Self, SolveError>

Discretizes geometry, builds the loop basis for ports and the pFFT operator of the filaments, and precomputes the port columns and the preconditioner.

Errors

As MeshSystem::assemble, plus SolveError::Pfft if the operator cannot be built from params.pfft.

pub fn ports(&self) -> &[Port]

The ports, in the order of the rows and columns of Z.

pub fn filaments(&self) -> &[Filament]

Every filament (branch), segment by segment.

pub fn mesh(&self) -> &MeshMatrix

The loop basis.

pub fn counts(&self) -> Counts

Problem size: nodes, segments, filaments, meshes, ports.

pub fn operator_stats(&self) -> PfftStats

Sizes of the pFFT operator.

pub fn gmres_params(&self) -> GmresParams

The GMRES parameters every solve uses.

pub fn impedance( &self, frequency: f64, ) -> Result<DMatrix<Complex<f64>>, SolveError>

The port impedance matrix Z(ω), in ohms, at frequency hertz; the iterative counterpart of MeshSystem::impedance.

Errors

SolveError::InvalidFrequency for a negative or non-finite frequency (reported with index 0), SolveError::NotConverged if a GMRES solve misses its tolerance, SolveError::Singular if the result is not finite.

pub fn solve(&self, frequency: f64) -> Result<IterativeSolution, SolveError>

As impedance, also returning how every port’s GMRES solve went.

Errors

As impedance.

pub fn sweep(&self, frequencies: &[f64]) -> Result<SweepResult, SolveError>

Z(ω) at every one of frequencies (hertz), with provenance; the iterative counterpart of MeshSystem::sweep.

Frequencies are solved concurrently, as many at a time as fit a fixed memory budget; each is an independent deterministic computation, so the result does not depend on the number of threads. Provenance::timing.inductance_s reports the time spent building the pFFT operator.

Errors

As impedance; the first failing frequency in list order is reported.

Trait Implementations

impl Debug for IterativeSystem

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.