主页 / user-7761783

Quantumania's questions

Martin Hope
Quantumania
Asked: 2023-12-31 20:44:28 +0800 CST
  • 7

我是 Rust 新手，我正在尝试使用 Rust 对多个矩阵执行顺序克罗内克积，即我想做

在此输入图像描述

在Python中，我知道我能做到

import numpy as np
from functools import reduce

X = np.array([[0, 1], [1, 0])
matrices = [X for _ in range(8)]
product = reduce(np.kron, matrices)

并得到想要的结果。如何在 Rust 中做同样的事情？

我的 Rust 代码目前如下所示：

use ndarray::{array, ArrayBase, OwnedRepr, Dim};
use ndarray::linalg::kron;
use num::complex::Complex64 as Complex;



fn X() -> ArrayBase<OwnedRepr<Complex>, Dim<[usize; 2]>> {
    array![
        [Complex::new(0.0, 0.0), Complex::new(1.0, 0.0)],
        [Complex::new(1.0, 0.0), Complex::new(0.0, 0.0)]
    ]
}

fn main() {    
    let mut matrices = Vec::new();
    for _ in 0..8 {
        matrices.push(X());
    }
    let product = matrices
        .iter()
        .reduce(|g1, g2| kron(g1, g2));

}

这是我得到的误差矩阵：

错误

这是我不明白的事情：

  • OwnedRepr：我不确定这到底是做什么用的，但这就是锈迹分析仪建议我做的
  • 我尝试根据错误消息的建议修改我的代码，但我 新的错误消息 没有得到有关如何处理的进一步建议。

我想知道如何解决这个问题，如果有任何关于如何使用 Rust 进行科学计算的初学者级示例文件的建议，我将不胜感激。

math