Question

我想用这样的boost定义一个常量3x3矩阵，它在执行过程中永远不会改变：

[1 2 3
 4 5 6
 7 8 9]

此矩阵将成为类的成员。那么，我可以像原始类型一样定义和初始化一个常量矩阵变量作为类成员吗？当我尝试为someMatrix变量键入const时，我无法在构造函数中分配矩阵数据并得到此错误：

error: assignment of read-only location '((Test*)this)->Test::someMatrix.boost::numeric::ublas::matrix<double>::operator()(0, 0)'

以下是代码：

Test.h

#ifndef TEST_H_
#define TEST_H_

#include <boost/numeric/ublas/matrix.hpp>

namespace bnu = boost::numeric::ublas;

class Test {
private:
    const double a = 1;
    const double b = 2;
    const double c = 3;
    const double d = 4;
    const double e = 5;
    const double f = 6;
    const double g = 7;
    const double h = 8;
    const double i = 9;
    const bnu::matrix<double> someMatrix;

public:
    Test();
    virtual ~Test();
};

#endif /* TEST_H_ */

Test.cpp的

Test::Test(){
    someMatrix(0,0) = a;
}

Main.cpp的

include "Test.h"

int main() {
    Test * t = new Test();

}

我真正想要的是找到一种方法来定义someMatrix，如下所示：

const bnu::matrix<double> someMatrix(3,3) = {a,b,c,d,e,f,g,h,i};

Answer 1

您可以编写辅助函数来执行此操作

class Test {
private:
    const bnu::matrix<double> someMatrix;
    static bnu::matrix<double> initSomeMatrix();
public:
    Test();
    virtual ~Test();
}

Test::Test() : someMatrix(initSomeMatrix()) {
}

bnu::matrix<double> Test::initSomeMatrix() {
    bnu::matrix<double> temp(3, 3);
    temp(0,0) = 1;
    ...
    return temp;
}

RVO应该使这个效率合理。

Answer 2

使用<boost/numeric/ublas/assignment.hpp>，您可以使用ublas::matrix将值插入ublas::vector或<<=，这样您就可以像这样设置矩阵：

bnu::matrix<double> a(3,3); a <<=  0, 1, 2,
                                   3, 4, 5,
                                   6, 7, 8;

要使它保持不变只需复制它：

const bnu::matrix<double> b = a;

HERE是从here

复制的有效最小示例

具有增强ublas的恒定矩阵

2 个答案: