使用faer crate从向量创建对角矩阵时触发断言失败问题排查

错误根源 #

你调用的column_vector_as_diagonal方法要求输入必须是列矩阵（行数为n、列数为1的矩阵），但你把向量转换成了1行5列的行矩阵，导致矩阵列数为5，触发了断言检查self.ncols() == 1失败。

正确实现方式 #

方式一：直接创建对角矩阵（推荐） #

faer提供了更直接的diagonal函数，可直接从向量生成对角矩阵，无需先转成列矩阵：

use faer::{mat, DiagonalMatrix};

fn main() {
    let vec = vec![1.0, 2.0, 3.0, 4.0, 5.0];
    let diag_mat = DiagonalMatrix::from_vec(vec);
    // 如需转为普通矩阵：
    let mat = diag_mat.to_mat();
}

方式二：先转成列矩阵再调用方法 #

如果一定要用column_vector_as_diagonal，需把向量转换成n行1列的列矩阵：

use faer::{mat, Mat};

fn main() {
    let vec = vec![1.0, 2.0, 3.0, 4.0, 5.0];
    // 创建n行1列的列矩阵
    let col_mat = Mat::from_column_vec(vec);
    // 调用方法生成对角矩阵
    let diag_mat = col_mat.column_vector_as_diagonal();
}

补充说明 #

• Mat::from_vec默认按行优先生成矩阵，传入长度为5的向量时会创建1行5列的矩阵，这是你之前出错的核心原因。
• Mat::from_column_vec会直接生成n行1列的列矩阵，满足column_vector_as_diagonal的参数要求。

内容的提问来源于stack exchange，提问作者Alberson Miranda