Function right_environment 

pub fn right_environment<T: SVDScalar + EinsumScalar>(
    mpo_a: &MPO<T>,
    mpo_b: &MPO<T>,
    site: usize,
    cache: &mut Vec<Option<Matrix2<T>>>,
) -> Result<Matrix2<T>>
where
    <T as ComplexFloat>::Real: Into<f64>,

Compute the right environment at site i for MPO contraction

The right environment R[i] represents the contraction of all sites i+1..L for the product of two MPOs A and B.

R[i] has shape (right_a_i, right_b_i) representing the accumulated contraction from the right.