用C++类替换Numerical配方的dmatrix

Question

我正在修改一个旧的应用程序，它广泛地使用了Numerical Recipes的dmatrix。因为我在这个应用程序上工作的原因之一是它的代码即将被打开，所以我想用可以免费分发的代码替换所有的Numerical Recipes代码。

dmatrix是一个返回双精度矩阵的函数。被调用提供了每个索引的下限和上限，如下所示：

double **mat = dmatrix(1,3,1,3);

mat现在有3行(从1到3 )和3列(从1到3)，所以mat[1][1]是第一个元素，mat[3][3]是最后一个元素。

我研究了各种C++矩阵实现，它们都不允许我指定每个维度的下限。有没有什么我可以使用的，或者我必须为此编写另一个矩阵类？

Answer 1

我相信你可以很容易地对其他一些矩阵实现进行包装，以添加下限特性。示例(未测试)：

class Matrix {
    OtherMatrix m;
    int lowerX, lowerY;
public:

    Matrix(int lx, int hx, int ly, int hy) :
        m(hx-lx, hy-ly),
        lowerX(lx), lowerY(ly) { }

    MatrixCol operator[] (int x) {
        return {this, x};
    } 
};

class MatrixCol {
    friend class Matrix;
    Matrix* mm;
    int x;
public:
    double& operator[] (int y) {
        return mm->m[x - mm->lowerX, y - mm->lowerY];
    } 
};

这可能需要更健壮的实现，这取决于您的用例。但这是基本的想法，从它扩展。