为什么Blitz++代码不能编译?

我看过Blitz++的源代码。 尽管看起来令人惊讶,但有 不 TinyMatrix的模板构造函数,但有一个template=运算符。 这意味着你不能做你想做的事情。所以我建议忘记你的矩阵 常量 。或者找到另一种方法,比如创建一个非常量矩阵,将其作为常量引用作为函数的参数传递。

现在,只是为了好玩,事实是你可以复制模板运算符=()代码,并在blitz/tinymat.h中创建一个新的模板构造函数,但我不会!如果你 真的 想要,以下是方法:

// Turn the following code...

template<typename T_expr>
TinyMatrix<T_numtype, N_rows, N_columns>&
operator=(_bz_tinyMatExpr<T_expr> expr)
{
    // USEFUL CODE
    return *this;
}

// ...into this :
template<typename T_expr>
TinyMatrix(_bz_tinyMatExpr<T_expr> expr)
{
    // USEFUL CODE
}

你的代码编译(并且可能工作)。

//const m33 d=blitz::product(a,b);  // Fails to compile.
const m33 d;
const_cast<m33>(d) = blitz::product(a,b);

抱歉,我不知道闪电战,但有些想法可能会有所帮助。

闪电战:TinyMatrix有合适的构造函数吗?
也许可以尝试其他形式的复制构造函数

const m33 d( blitz::product(a,b) );

使用标准的const关键字真的合适吗?也许你需要

typedef blitz::TinyMatrix<const double,3,3> const_m33;
typedef blitz::ConstTinyMatrix<double,3,3> const_m33;