求教一个关于运算符重载和GC的问题！

如下代码中，我尝试重构Matrix之间和其与f64的所有矩阵运算法则。我想求问一下，所有步骤我都返回一个新的Matrix算是一个比较适合的方法吗？比如let ci = &a + c1 * c2 * &b;相当于let ci = &Matrix1 + Matrix2，在我这个代码里，Matrix2不会被消耗掉，但也不与任何变量绑定，它会在什么时候drop？是赋值给ci的瞬间，它就被drop了吗？还是在main函数结束时？

第二种办法是只允许&Matrix参与运算，得到的答案也是&Matrix。此时，按我的理解是，内存里有一个Matrix但不与任何变量绑定，而我们手头上只能存在有它的引用&Matrix。是这样吗？

哪种办法是比较合理，或者说，计算效率最高的呢？我的理解是计算效率是一致的。但是第一种方法，是否会消耗额外的内存。而第二种方法，写起来方便（只用写关于&Matrix那一套），但手上只有&Matrix引用会比较不方便吗？

此外，f64 + Matrix和Matrix + f64需要写两套，很不方便。比如这里我就忘了写Matrix + &Matrix。我看ndarray貌似是用宏解决这个问题，有更简单的方法吗？

fn main() {
    let a = Matrix::Array(array![[0.0, 1.0], [2.0, 3.0]]);
    let b = Matrix::Array(array![[4.0, 5.0], [6.0, 7.0]]);
    let c1 = 0.001f64;
    let c2 = 0.5f64;
    let ci = &a + c1 * c2 * &b;
}

#[derive(Debug, PartialEq, Clone)]
pub enum Matrix {
    Raw(Vec<Vec<f64>>),
    RawPureDiagonal(Vec<f64>),
    RawSparse(Vec<(usize, usize, f64)>),
    Array(Array2<f64>),
    PureDiagonal(Array1<f64>),
}

impl<'a, 'b> Add<&'b Matrix> for &'a Matrix {
    type Output = Matrix;
    // --snip--
}

impl<'a> Add<Matrix> for &'a Matrix {
    type Output = Matrix;

    fn add(self, other: Matrix) -> Matrix {
        match (self, &other) {
            // ndarray库运算规则, &Array2 + &Array2 不消耗任何变量，生成一个新的Array2
            (Matrix::Array(x), Matrix::Array(y)) => {Matrix::Array(x + y)},
            _ => {panic!("Unsupported operations!")},
        }
    }
}

impl<'a> Mul<f64> for &'a Matrix {
    type Output = Matrix;
    // --snip--
}

impl Mul<f64> for Matrix {
    type Output = Matrix;
    // --snip--
}

impl<'a> Mul<&'a Matrix> for f64 {
    type Output = Matrix;

    fn mul(self, other: &'a Matrix) -> Matrix {
        match (self, other) {
            // ndarray库运算规则, &Array2 * f64 只消耗f64，生成一个新的Array2
            (x, Matrix::Array(y)) => {Matrix::Array(x * y)},
            _ => {panic!("Unsupported operations!")},
        }
    }   
}

impl Mul<Matrix> for f64 {
    type Output = Matrix;
    // --snip--
}