我正在改造一个旧的应用程序,该应用程序广泛使用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];
}
};
这可能需要一个更健壮的实现,具体取决于您的用例。但这是基本的想法,从中扩展。