夜行録

最小点探索

例によって gsl だが、１変数の minimization をやらせてみたら、これが探索範囲の上限下限の他に候補値を入れてやらないといけない。
ところがこれが、たとえば上限下限が a, b で候補値が c とすると、 f(c) が f(a), f(b) の両方より小さくないと気に入らないようで困る。たとえば最小値がレンジの中には確かにあるんだがレンジ境界ギリギリにあったりすると、何も考えずに候補値として (a+b)/2 なんかを突っ込むと f(a) < f(c) < f(b) とかなってくれてあっさりハネられる。
全域嘗めて候補値を探してもいいんだが、そこまでやるくらいならそのついでに最低点を発見できそうだし。っていうか、ライブラリの側でやれ。

で、１番の問題は、こうして候補値が不適当だった時、セットする関数が返り値で notice するのではなく、あっさりライブラリの中でエラーを吐いて abort してくれやがる点だ。

なんでエラーになるのかと二日くらい悩んで debug して、こういうことだということがやっと判ってドキュメント見ると確かに

The value of the function at x must be less than the value of the function at the ends of the interval

とか書いてある…。うーむ、ライブラリの仕様の話だったのか。騙された気分でいっぱいだ。