将数值配方的 dmatrix 替换为 C++ 类

Replace Numerical Recipe's dmatrix with a C++ class

本文关键字：替换 C++ dmatrix      更新时间：2023-10-16

我正在改造一个旧的应用程序，该应用程序广泛使用Numerical Recipes的dmatrix。由于我开发该应用程序的原因之一是它的代码即将打开，所以我想用可以自由分发的代码替换所有的Numerical Recipes代码。

CCD_ 2是一个返回一个二重矩阵的函数。called为每个索引提供下限和上限，如下所示：

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

mat现在有从1到3的3行和从1到三的3列，因此mat[1][1]是第一个元素，mat[3][3]是最后一个元素。

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

我相信您可以很容易地制作其他矩阵实现的包装器来添加下界功能。示例（未经测试）：

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

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

相关文章：