Skip to main content

lsqr_complex

qsm_core::inversion::ilsqr

Function lsqr_complex 

Source 
pub fn lsqr_complex<F, G>(
    apply_a: F,
    apply_ah: G,
    b: &[Complex64],
    tol: f64,
    max_iter: usize,
    verbose: bool,
) -> Vec<Complex64> ⓘ
where
    F: Fn(&[Complex64]) -> Vec<Complex64>,
    G: Fn(&[Complex64]) -> Vec<Complex64>,
Expand description

LSQR iterative solver for Ax = b (complex version)

Solves the least squares problem min ||Ax - b||² using the LSQR algorithm. Based on Paige & Saunders (1982), with convergence tests matching MATLAB’s lsqr.

Convergence tests (matching MATLAB):

  1. ||r|| / ||b|| <= btol + atol * ||A|| * ||x|| / ||b|| (residual test)
  2. ||A’r|| / (||A|| * ||r||) <= atol (normal equations test)