2

我正在尝试使用 Boost 的 ublas 部分，但由于某种原因，我无法将矩阵相乘并将结果分配给其他矩阵。

这有效：

#include <boost/numeric/ublas/symmetric.hpp>
#include <boost/numeric/ublas/io.hpp>

using namespace boost::numeric::ublas;

typedef symmetric_matrix<int,lower> symatrix;

int main() {
  int N = 10;
  symatrix foo(N,N);
  for (int i = 0; i < N; i++)
    for(int j = 0; j <= i; j++) {
      foo(i,j) = i - j + 1;
    }
  symatrix goo(foo);
  //goo = prod(foo,foo);
  std::cout << prod(foo,foo)<< std::endl;

}

但是，如果我取消注释该行goo = prod(foo,foo);或尝试类似：

symatrix goo = prod(foo,foo);

我收到一个我无法破译的运行时错误。

Check failed in file /usr/include/boost/numeric/ublas/detail/matrix_assign.hpp at line 761:
detail::expression_type_check (m, cm)
terminate called after throwing an instance of 'boost::numeric::ublas::external_logic'
  what():  external logic
Aborted

如何将矩阵相乘并分配结果？

4

1 回答 1

2

当您将两个对称矩阵相乘时，不能保证总是得到一个对称矩阵。因此，此错误可能与此有关，尽管我不知道为什么当我将您的类型更改为 symmetric_matrix 类型为双精度时代码会起作用。

于 2011-05-02T18:38:12.577 回答