fn internal_iteration(
rhs: &DMatrix<Complex>,
apply: impl Fn(MatRef<'_, Complex>) -> DMatrix<Complex>,
) -> Result<Option<DMatrix<Complex>>>Expand description
The unrestarted Krylov solve of the operator, or of its adjoint, for a thin batch, or
None when it fails numerically and the caller should factor the operator instead.