我正在使用n代数,我想通过将一个矩阵设置为另一个具有兼容维度的列来修改它,如下所示:
let zero: T = convert(0.0);
let mut basis = DMatrix::from_element(rows, rank, zero);
// matrix is an input with dimensions rows x cols
for i in 0..location.len()
{
basis.column_mut(i) = matrix.column(location[i]);
}
我还try 了取消对任务双方的引用,并寻找了某种"分配"方法,但没有成功.
set_column
不起作用,因为DMatrix
不实现DimName
我目前的工作是这样的,但我一点也不喜欢:
// Construct basis vectors initialized to 0.
let zero: T = convert(0.0);
let mut basis = DMatrix::from_element(rows, rank, zero);
for i in 0..location.len()
{
// Copy the pivot column from the matrix.
let col = location[i];
for r in 0..matrix.nrows()
{
basis[(r, i)] = matrix[(r, col)];
}
}