c++ - 处理 boost ublas 稀疏或密集矩阵的 API

Question

我对 Boost ublas 文档有点困惑。我似乎不清楚稀疏矩阵类和密集矩阵类共享一个共同的父类——我相信这是设计使然。但是，我如何设计一个可以接受稀疏或密集矩阵的 API，假设它只需要使用 operator() 访问器对矩阵的条目进行操作，例如。例如，像这样：

float doMatrixMath(matrix_base<float> m) 
{
  return m(1,1)+m(2,2);
}

也许我对这个的想法是错误的。任何有关如何考虑 ublas 类的对象建模的指导将不胜感激！

score 1 · Accepted Answer

不幸的是模板。您可以使用非常通用的类型，或者深入研究并找到更具体的内容，但总体思路是：

template< typename MatrixType >
float doMatrixMath(MatrixType m) 
{
  return m(1,1)+m(2,2);
}

当然，这可以通过更具体的类型和返回值检测来增强......

1 回答 1