Question

我正在修改一个使用Numerical Recipes＆＃39;的旧应用程序。 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];
    } 
};

根据您的使用情况，这可能需要更强大的实施。但这是基本的想法，从中扩展。

用C ++类替换Numerical Recipe的dmatrix

1 个答案: