Struct InteractionFactor 

pub struct InteractionFactor {
    local: LocalMatrix,
    coupling: DMatrix<Complex>,
    lu: Lu,
}

The LU factorization of I - T C, kept with the local T-matrix T and the coupling C, for any number of incident-field solves.

A block-diagonal T stays a list of particle blocks; it is never expanded to a dense matrix. record shares the factor through an Arc, so every recorded solve keeps only its own incident and scattered fields.

treams-rs extension. treams solves the full system once with TMatrix.interaction.solve(); interaction computes that result.

Fields

local: LocalMatrix

coupling: DMatrix<Complex>

lu: Lu

Implementations

impl InteractionFactor

fn write_state(&self, writer: &mut Writer)

fn read_body( reader: &mut Reader<'_>, sizes: &[usize], dimension: usize, ) -> Result<Self>

impl InteractionFactor

pub fn new(local: DMatrix<Complex>, coupling: DMatrix<Complex>) -> Result<Self>

Factor I - T C for a dense local matrix T.

pub fn from_blocks( local: Vec<DMatrix<Complex>>, coupling: DMatrix<Complex>, ) -> Result<Self>

Factor I - T C for a block-diagonal T given by its particle blocks.

fn factorize(local: LocalMatrix, coupling: DMatrix<Complex>) -> Result<Self>

pub fn dimension(&self) -> usize

Number of multipole channels.

pub fn solve(&self, incident: &DMatrix<Complex>) -> Result<DMatrix<Complex>>

Solve (I - T C) scattered = T incident for only the supplied columns.

pub fn pullback_incident( &self, cotangent: &DMatrix<Complex>, ) -> Result<DMatrix<Complex>>

Gradient of incident columns for a fixed factor: Tᴴ (I - T C)⁻ᴴ G. This borrows the saved LU and needs no incident or scattered fields.

pub fn record( self: &Arc<Self>, incident: DMatrix<Complex>, ) -> Result<IlluminateResidual>

Solve like solve and keep B + C X for both derivative directions. This response replaces the incident fields, and the residual shares this factor. Use solve when only the scattered fields are needed.

fn solve_system(&self, rhs: DMatrix<Complex>) -> Result<DMatrix<Complex>>

(I - T C)⁻¹ rhs; a non-finite solution means a numerically singular operator.

fn validate(&self, values: &DMatrix<Complex>) -> Result<()>

fn validate_tangents( &self, local: &[DMatrix<Complex>], coupling: &DMatrix<Complex>, ) -> Result<()>

Trait Implementations

impl Clone for InteractionFactor

fn clone(&self) -> InteractionFactor

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for InteractionFactor

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

Formats the value using the given formatter.