ということで、log(0) は EDOM でも ERANGE でもいいっぽいですねえ
違うよ、それは違うよ。Sun コンパイラが言っているのは、要するに
- SVID と X/Open 両規格で要求する結果が一致しないから、規格を明示してくれたらその期待通りの出力してあげる
であって、SVID 指定なら EDOM だし、X/Open 指定なら ERANGE、です
(デフォルトは X/Open)。
どっちも原規格にあたっていないので、shall なのか should なのか知りませんけど。
log(0.0) の結果を printf で表示するようにして、
指定によって変わる様を示しておきます。
% cat log.c
#include <math.h>
#include <errno.h>
#include <stdio.h>
int main(int ac, char **av)
{
double d;
d=log(0.0);
perror("");
printf("%d\n", errno);
printf("%f\n", d);
}
% cc -lm -Xt log.c
% a.out
log: SING error
Argument out of domain
33
-340282346638528859811704183484516925440.000000
% cc -lm -Xa log.c
% a.out
Result too large
34
-Inf
他に関係する規格の規定はこんな感じ、
- C99:「値域エラーになることもある」(定義域エラーにはならない)
- IEEE754: 定義域内であり(log の定義域は [0,+∞]) 結果値は -Inf 、とするのが推奨 (recommended)。
したがって、環境がどの規格に準拠しているのかで、定義域エラーになるか値域エラーになるかが違うワケです。が、定義域エラーになるのはめずらしいような。
なお、どれでも、log(-1) は定義域エラーです。
_ [雑記] つづき。
python の元の問題を試してみたのさ。
import math
for x in range(0, 10):
z = 10 ** int(math.log10(x))
if x > z and x % z:
continue
print(x)
さっきと同じ環境だと
% python log.py
Traceback (most recent call last):
File "log.py", line 3, in <module>
z = 10 ** int(math.log10(x))
OverflowError: math range error
ふむ。さっきの log(0.0) の結果 (-Xa) と同じ。
別の環境
% gcc -lm log.c
% a.out
Numerical result out of range
29
-inf
おっと、errno が上書きされてる。perror() と errno 表示をひっくりかえすと
% a.out
34
Numerical result out of range
-inf
この環境で、
% python log.py
Traceback (most recent call last):
File "log.py", line 3, in <module>
z = 10 ** int(math.log10(x))
ValueError: math domain error
Huh? お前は何をいっているのだ？ ちなみに ruby だとおかしなことにはなっていない。
% ruby -e 'print Math.log(0.0)'
-e:1:in `log': Numerical result out of range - log (Errno::ERANGE)
from -e:1
_ うんの [算数.org の方にもメールしましたが、明日は休講になりました。 ゆっくり休んでくださいませ。]