c++ - ublas：将 ublas::vector 包装为 ublas::matrix_expression

Question

我是 Boost::uBLAS 的菜鸟。

我有一个以 aublas::matrix_expression<double>作为输入的函数：

namespace ublas = boost::numeric::ublas;

void Func(const ublas::matrix_expression<double>& in,
                ublas::matrix_expression<double>& out);

调用者将行向量保存为ublas::vector<double>，我希望将其传递给Func。

到目前为止，我还没有找到任何方法来做到这一点。
最好的方法是什么，最好没有任何临时分配？

谢谢。

Answer 1

好吧，有一个选项可以将内存的连续区域的只读适配器创建为只读矩阵。看看例子 3。使用起来非常简单：

#include "storage_adaptors.hpp"
#include <boost/numeric/ublas/assignment.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>

ublas::vector<double> v(6);
v <<= 1, 2, 3, 4, 5, 6;
ublas::matrix<double> m = ublas::make_matrix_from_pointer(2, 3, &v(0));
std::cout << m << std::endl;

可能您可以进一步调整以适应您的需求/示例。

Answer 2

如果你准备牺牲一些乘法，你可以避免分配，使用

outer_prod(scalar_vector<double>(1, 1), vec)

将向量转换为矩阵表达式。此外，您的功能可能应该是

template<class C>
void Func(const matrix_expression<C>& in...

matrix_expression 本身并不模拟矩阵表达式概念，它只是基类。