数值食谱中显然未定义的变量

时间:2016-01-27 23:38:31

标签: c++ undefined

我正在尝试使用Numerical Recipes,3rd Edition,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的理解。

1 个答案:

答案 0 :(得分:3)

在这一行:

 MatDoub pp(ndim+1,ndim);

pp未被使用正在声明和构建它。在此行pp之后是一个初始化对象,后续用法很好。您可以找到有关constructor initialization here.

的更多信息