Function restart_gmres_with_truncation 

pub fn restart_gmres_with_truncation<T, F, Tr>(
    apply_a: F,
    b: &T,
    x0: Option<&T>,
    options: &RestartGmresOptions,
    truncate: Tr,
) -> Result<RestartGmresResult<T>>
where
    T: TensorLike,
    F: Fn(&T) -> Result<T>,
    Tr: Fn(&mut T) -> Result<()>,

Solve A x = b using restarted GMRES with truncation.

This wraps gmres_with_truncation with an outer loop that recomputes the true residual at each restart. This is particularly useful for MPS/MPO computations where truncation can cause the inner GMRES residual to be inaccurate.

Algorithm

for outer_iter in 0..max_outer_iters:
    r = b - A*x0          // Compute true residual
    r = truncate(r)
    if ||r|| / ||b|| < rtol:
        return x0         // Converged
    x' = gmres_with_truncation(A, r, 0, inner_options, truncate)
    x0 = truncate(x0 + x')

Type Parameters

• T - A tensor type implementing TensorLike
• F - A function that applies the linear operator: F(x) = A x
• Tr - A function that truncates a tensor in-place: Tr(&mut x)

Arguments

• apply_a - Function that applies the linear operator A to a tensor
• b - Right-hand side tensor
• x0 - Initial guess (if None, starts from zero)
• options - Solver options
• truncate - Function that truncates a tensor to control bond dimension

Returns

A RestartGmresResult containing the solution and convergence information.