kv-0.4.52 - kashiの日記

2021/08/05(木)kv-0.4.52

kvライブラリを0.4.52にアップデートしました。

変更点は以下の通りです。

test-rounding.ccで、非正規化数を強制的に0にする"flush to zero"モードになっていないかどうかをチェックするようにした。
iccでコンパイルできるように修正。
optimize.hppで、「端にない」区間で傾きが単調ならばそこに最小値は存在しないというチェックを入れる改良を行った。
defint-newtoncotes.hpp, double-newtoncotes.hppで、最適な分割数の推定がうまくできていなかったバグを修正。
defint.hpp, defint-singular.hppで、急減少関数などの特定の被積分関数に対してstep sizeがうまく計算できていなかった問題を修正。
psa.hppの中のinv関数を、係数のオーバーフローが発生しづらいように改良。
gamma.hppの中のdigamma関数の計算アルゴリズムをより高精度なものに変更。
その他細かい修正。

iccはデフォルトで非正規化数を強制的に0にする"flush to zero"になっているので、そのままだと区間演算の結果が異常になってしまいます。例えば、上向き丸めモードで2^{-1022}/2を計算すると、0になってしまいます。すなわち、IEEE754に従っていません。オプション-no-ftz (windowsなら/Qftz-)を付けるとこの問題は回避できるようです。iccを使うときは、test-rounding.ccをコンパイル、実行してチェックしてみて下さい。

また、digammaを高精度化したのは、ある問題を解くのに必要になったためです。忘れないうちにその話も記事にするつもりです。