Function lsmr 

pub fn lsmr<T: ComplexField>(
    out: MatMut<'_, T>,
    right_precond: impl BiPrecond<T>,
    mat: impl BiLinOp<T>,
    rhs: MatRef<'_, T>,
    params: LsmrParams<T::Real>,
    callback: impl FnMut(MatRef<'_, T>),
    par: Par,
    stack: &mut MemStack,
) -> Result<LsmrInfo<T::Real>, LsmrError<T::Real>>

executes lsmr using the provided preconditioner

note

this function is also optimized for a rhs with multiple columns