Skip to main content

initial_solve_single

qsm_core::inversion::ilsqr

Function initial_solve_single 

Source 
fn initial_solve_single(
    n: [usize; 3],
    p: &[f64],
    w: &[f64],
    dk: &[f64],
    max_iter: usize,
) -> (Vec<f64>, LsqrOutcome)
Expand description

Step 3 of ilsqr as STI computes it (IlsqrPrecision::Single): the right-hand side D W φ in double, rounded to single; the operator D W D evaluated in single precision (single 3-D FFTs, single(D), single(W)) and MATLAB’s lsqr on single data. Returns the real part of the solution in double, as STI’s double(x1).