数值配方中显然未定义的变量

Apparently undefined variable in Numerical Recipes

本文关键字:未定义 变量 方中显      更新时间:2023-10-16

我正在尝试使用Numerical Recipes,第3版,Amoeba.h,但找不到定义变量的位置。这是代码的开头:

struct Amoeba {
    const Doub ftol;
    Int nfunc;
    Int mpts;
    Int ndim;
    Doub fmin;
    VecDoub y;
    MatDoub p;
    Amoeba(const Doub ftoll) : ftol(ftoll) {}
    template <class T>
    VecDoub minimize(VecDoub_I &point, const Doub del, T &func)
    {
        VecDoub dels(point.size(),del);
        return minimize(point,dels,func);
    }
    template <class T>
    VecDoub minimize(VecDoub_I &point, VecDoub_I &dels, T &func)
    {
        Int ndim=point.size();
        MatDoub pp(ndim+1,ndim);           <--------- XXXXXXX
        for (Int i=0;i<ndim+1;i++) {
            for (Int j=0;j<ndim;j++)
                pp[i][j]=point[j];
            if (i !=0 ) pp[i][i-1] += dels[i-1];
        }
        return minimize(pp,func);
    }
    template <class T>
    VecDoub minimize(MatDoub_I &pp, T &func)
    {
        const Int NMAX=5000;
        const Doub TINY=1.0e-10;
        Int ihi,ilo,inhi;
        mpts=pp.nrows();
        ndim=pp.ncols();
        VecDoub psum(ndim),pmin(ndim),x(ndim);
        p=pp;
.....

在标有X的行上,使用了变量pp,但似乎没有定义。在下一个方法中,pp 用作定义为 MatDoub_I &pp 的参数变量(MatDoub 是一个实现矩阵的类),但在第一个函数中它似乎是未定义的。这只是作者的错误,还是我对这里的C++有什么不了解的?

请注意,只有 p 被定义为上面的MatDoub。这是否意味着pp与(MatDoub)(MatDoub)相同?这违反了我对C的理解。

在这一行中:

 MatDoub pp(ndim+1,ndim);

pp没有被使用,它正在被声明和建造。在此行之后pp是一个初始化的对象,后续使用就可以了。您可以在此处找到有关constructor initialization的更多信息。