Struct ForwardModel 

pub struct ForwardModel<'a> { /* private fields */ }

Allocation-friendly CPU implementation of the image-plane FPM forward model.

Implementations

impl<'a> ForwardModel<'a>

pub fn new(model: &'a ImagePlaneModel) -> Result<Self>

pub fn with_backend( model: &'a ImagePlaneModel, backend: Arc<dyn Backend>, ) -> Result<Self>

pub fn model(&self) -> &ImagePlaneModel

pub fn workspace(&self) -> ForwardWorkspace

pub fn extract_patch( &self, object_spectrum: &Array2<Complex64>, source: usize, ) -> Result<Array2<Complex64>>

pub fn insert_patch_update( &self, object_spectrum: &mut Array2<Complex64>, source: usize, update: &[Complex64], scale: f64, ) -> Result<()>

pub fn apply_pupil( &self, patch: &[Complex64], pupil: &Pupil, destination: &mut [Complex64], ) -> Result<()>

pub fn forward_source_field( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, source: usize, ) -> Result<Array2<Complex64>>

Coherent low-resolution field for one illumination source.

pub fn forward_source_field_into<'b>( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, source: usize, workspace: &'b mut ForwardWorkspace, ) -> Result<&'b [Complex64]>

pub fn forward_field( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, frame: usize, ) -> Result<Array2<Complex64>>

Coherent field for a non-multiplexed frame. A coded frame with exactly one source is represented by a field scaled by the square-root weight.

pub fn forward_intensity( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, frame: usize, ) -> Result<Array2<f64>>

pub fn forward_intensity_into( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, frame: usize, workspace: &mut ForwardWorkspace, destination: &mut [f64], ) -> Result<()>

pub fn forward_intensity_stack_into( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, destination: &mut [f64], worker_count: usize, ) -> Result<()>

Predicts every measured frame into a contiguous [frame][row][column] destination using independent worker scratch.

Frame order and floating-point results within each frame are unchanged by worker_count. Values larger than the frame count are capped; zero is rejected. A worker count of one executes directly without spawning.

pub fn forward_intensity_stack( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, worker_count: usize, ) -> Result<Vec<f64>>

Allocation-owning convenience wrapper for Self::forward_intensity_stack_into.

pub fn residual( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, frame: usize, measured: &[f64], ) -> Result<Vec<f64>>

pub fn amplitude_projection( &self, field: &[Complex64], measured_intensity: &[f64], frame: usize, epsilon: f64, ) -> Result<Vec<Complex64>>

pub fn frame_loss( &self, object_spectrum: &Array2<Complex64>, pupil: &Pupil, frame: usize, measured: &[f64], loss_type: LossType, ) -> Result<f64>