1

我正在修改一个dmatrix非常广泛地使用数字食谱的旧应用程序。由于我正在开发该应用程序的原因之一是它的代码即将被打开，因此我想用可以自由分发的代码替换所有数字食谱代码。

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

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

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

我查看了各种 C++ 矩阵实现，没有一个允许我指定每个维度的下限。有什么我可以使用的，还是我必须为此编写另一个矩阵类？

4

1 回答 1

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];
    } 
};

根据您的用例，这可能需要更健壮的实现。但这是基本思想，从它展开。

于 2014-09-23T11:41:59.950 回答