Skip to main content

lsmr_diag

qsm_core::inversion::rts

Function lsmr_diag 

Source 
fn lsmr_diag(
    d: &[f64],
    b: &[Complex64],
    nx: usize,
    maxit: usize,
) -> Vec<Complex64> ⓘ
Expand description

LSMR for the diagonal operator diag(d) on a spectrum, from x = 0, with lambda = 0, atol = btol = 0 and conlim = 1e8, as QSM.jl’s rts calls its copy of IterativeSolvers.jl’s lsmr! (v0.9.2).