为什么这次闪电战不发生？ 代码编译？

Question

我是 blitz++ 新手。 到目前为止，一切都很好，但我有点困惑为什么下面代码中的注释行无法

error: conversion from ‘blitz::_bz_tinyMatExpr<blitz::_bz_tinyMatrixMatrixProduct<double, double, 3, 3, 3, 3, 1, 3, 1> >’ to non-scalar type ‘const m33’ requested

在 Debian/Lenny（g++ 4.3.2，Blitz 0.9）上编译。 这是代码：

#include <blitz/blitz.h>
#include <blitz/array.h>
#include <blitz/matrix.h>
#include <blitz/matexpr.h>
#include <blitz/tinymat.h>
#include <blitz/tinymatexpr.h>
#include <blitz/tinyvec.h>
#include <blitz/tinyvec-et.h>
#include <blitz/vector.h>
#include <blitz/vector-et.h>

typedef blitz::TinyVector<double,3> v3;
typedef blitz::TinyMatrix<double,3,3> m33;

int main(int,char**)
{
  const m33 a;
  const m33 b;
  m33 c;c=blitz::product(a,b);  // Compiles OK
  //const m33 d=blitz::product(a,b);  // Fails to compile.
  return 0;
}

我确实喜欢完全 const-ed，所以很高兴知道是否有办法让它工作（到目前为止 Blitz++ 的经验表明这可能只是选择正确的问题）包括...）。

Answer 1

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

现在，只是为了好玩，事实是您可以复制模板operator=() 代码并在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
}

以及您的代码的编译（并且可能有效）。

Answer 2

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

Answer 3

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

blitz::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;