我正在尝试使用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的理解。
答案 0 :(得分:3)
在这一行:
MatDoub pp(ndim+1,ndim);
pp
未被使用正在声明和构建它。在此行pp
之后是一个初始化对象,后续用法很好。您可以找到有关constructor initialization
here.