52nlp上无约术最优化的文章:http://www.52nlp.cn/unconstrained-optimization-one ,后面提到了一个无约束最优化的C++实现http://www.chokkan.org/software/liblbfgs/。
研究这个简单的实现,算是做了一些小解释。可以参考52npl上的文章看这个程序。程序需要VC 2010编译。
sample
lbfgs.c
{
lm;
计算原函数和梯度fx ,g;
d = - g;
计算xnorm,gnorm;
判定gnorm/xnorm 是否满足条件;