[ruby-list:37740] Re: missing/erf.c

まつもと ゆきひろです

In message "[ruby-list:37739] Re: missing/erf.c"
    on 03/05/31, NISHIMATSU Takeshi <[e-mail]> writes:

|[ruby-list:37738] でまつもとさん:
|（数学関数について）
|> 厳密に言うとPOSIXに入ってないものに関してはちょっと慎重にな
|> ろうと思っています。ですから、 j0, j1, jn, y0, y1, yn の方が
|> 取り込みやすいですね。
|
|どのような点について慎重になろうとされているのでしょうか.
|（１）取り込みそのものについて. やたら関数が増えるのはよくない.
|（２）関数名.
|（３）（不正な引数が与えられたときの）エラーの処理について.
|（４）missing/*.c の精度.
|でしょうか.
|
|今回の erf, erfc については
|（１）SVID3, XPG4.2, BSD4.3, C99 で採用されているので, そんなに
|      特殊ではない. 個人的にはないと不便.
|（２）erf, erfc 以外になることはおそらくない.
|（３）全領域にわたってなめらかな関数なので, エラーは返さない.
|      （画像を http://www.vni.com/products/jmsl/api/com/imsl/math/Sfun.html
|      に見つけました. Googleの画像検索は便利ですね.）
|（４）じつはあまり注意していません.
|      また, inexact, underflow はセットしなくてかまわないのですよね.
|と思います.

「きりがないから」というのがその理由です。で、一応POSIXをひ
とつの目安にしてます。が、C99に入ってるなら問題ないでしょう。

|まずその行は実行されないハズで, 実行されたとしても最小
|の桁が浮動小数点数演算のため振動した結果だと思うので,
|［１］その行は削除.
|［２］printf("missing/erf.c:q_gamma(): could not converge.");
|      などにしておく.
|もしくは,
|［３］http://ftp.rutgers.edu/pub/FreeBSD/FreeBSD-current/src/lib/libm/common_source/erf.c
|      のような近似曲線を使ったものにする.（これをその
|      まま使うのはライセンス的にまずいのですよね.）
|という策があります.

どれが良いのかは、私がerfを使わないので分かりません。

少なくとも[3]のものにライセンス的な問題はありません、BSDライ
センスなら。FreeBSDの一部なのできっとそうでしょう。