Function solve_with_pullback 

pub fn solve_with_pullback(
    a: MatrixView<'_>,
    b: MatrixView<'_>,
    tol: Tolerance,
) -> Result<(Matrix, SolvePullback), SolverError>