a********r
发帖数: 92 | 1 in function amotry() for downhill simplex minimization,
why fac2=fac2=fac1-fac; instead of fac2=-fac;
then this may not be a reflection or extrapolation, because
ytry[j] = center * (1-fac) + p[ihi][j]*face - (1-fac)/ndim * p[ihi][j]
the first two terms is the extrapolation, what's the point for the 3rd term?
It will be incorrect in some case, reflection will not be reflection anymore.
float amotry(float **p, float y[], float psum[], int ndim,
float (*funk)(float []), int ihi, float fac)
Extrapo | m****n
发帖数: 51 | 2 Both Numerical Recipes in C and Numerical Recvipes in Fortran should be used.
anymore.
【在 a********r 的大作中提到】
: in function amotry() for downhill simplex minimization,
: why fac2=fac2=fac1-fac; instead of fac2=-fac;
: then this may not be a reflection or extrapolation, because
: ytry[j] = center * (1-fac) + p[ihi][j]*face - (1-fac)/ndim * p[ihi][j]
: the first two terms is the extrapolation, what's the point for the 3rd term?
: It will be incorrect in some case, reflection will not be reflection anymore.
: float amotry(float **p, float y[], float psum[], int ndim,
: float (*funk)(float []), int ihi, float fac)
: Extrapo
|
|
