Struct CylinderResidual 

pub struct CylinderResidual {
    blocks: Vec<Block>,
    solves: Vec<MieCylResidual>,
    kz_count: usize,
    boundaries: usize,
}

What cylinder saves for its pullback: one boundary solve per (kz, m) pair, shared with the mirror pair (-kz, -m) when both are requested.

Fields

blocks: Vec<Block>

solves: Vec<MieCylResidual>

kz_count: usize

boundaries: usize

Implementations

impl CylinderResidual

pub fn state_size( kz_count: usize, mmax: usize, boundaries: usize, ) -> Result<usize>

Byte count of saved derivative state for these input shape parameters.

impl CylinderResidual

pub const fn input_counts(&self) -> (usize, usize)

The counts of recorded axial wavenumbers and layer boundaries.

pub fn shape(&self) -> (usize, usize)

The shape of the T-matrix: two helicities for each (kz, m) block.

pub fn pushforward( &self, kzs: &[f64], k0: f64, radii: &[f64], materials: &[Material], ) -> Result<DMatrix<Complex>>

The T-matrix tangent along changes of the axial wavenumbers, k0, radii and materials. Material entries store changes of epsilon, mu and kappa. Mirror blocks reuse their cached solve, with the sign of their own kz direction reversed; their axial directions need not preserve the symmetry.

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

Gradients of kzs, k0, the radii and the materials from cotangent, the gradient of a real loss with respect to the T-matrix.

The blocks run in parallel and their gradients add in block order, so the Rayon pool does not change them.

Trait Implementations

impl Debug for CylinderResidual

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

Formats the value using the given formatter.

impl SavedState for CylinderResidual

fn save_state(&self) -> Result<Vec<u8>>

Copy the saved numerical state into its fixed byte layout.

fn from_state(bytes: &[u8]) -> Result<Self>

Restore numerical state from bytes produced by this extension.