PerlでNaNを使う．

Perlで数値解析をする際，ゼロ除算を適切に扱う必要がある．Perlはv5.6.1以上であればNaNを使える．NaNはどんな演算をしてもNaNとなるが，実数NaN共にゼロ除算ではエラーを吐いて止まってしまう．ゼロ除算時にNaNを発生し伝搬させるためには，分母が0であればNaNとなる様な関数を書いてやればよい(下の例ではsub DivNaN)．
ちなみに無限も'inf'として使える．

BACK全選択

#!/bin/perl
$value = 10; 
$value = &DivNaN($value , 20);
printf("10/20= $value\n");
 
$value = &DivNaN('inf' , 20);
printf("inf/20= $value\n");
 
$value = &DivNaN('nan' , 20);
printf("nan/20= $value\n");
 
$value = &DivNaN(10 , 0); 
printf("10/0= $value\n");
 
$value = &DivNaN('nan' , 'nan');
printf("nan/nan= $value\n");
 
$value = &DivNaN('nan' , 0); 
printf("nan/0= $value\n");
 
sub DivNaN {
  my ($numerator, $denominator) = @_; 
  if ($denominator eq 0){ 
    return 'nan';
  } else {
    return $numerator / $denominator;
  }  
}

#!/bin/perl $value = 10; $value = &DivNaN($value , 20); printf("10/20= $value\n"); $value = &DivNaN('inf' , 20); printf("inf/20= $value\n"); $value = &DivNaN('nan' , 20); printf("nan/20= $value\n"); $value = &DivNaN(10 , 0); printf("10/0= $value\n"); $value = &DivNaN('nan' , 'nan'); printf("nan/nan= $value\n"); $value = &DivNaN('nan' , 0); printf("nan/0= $value\n"); sub DivNaN { my ($numerator, $denominator) = @_; if ($denominator eq 0){ return 'nan'; } else { return $numerator / $denominator; } }

BACK全選択

% perl test.pl
10/20= 0.5
inf/20= inf
nan/20= nan
10/0= nan
nan/nan= nan
nan/0= nan

% perl test.pl 10/20= 0.5 inf/20= inf nan/20= nan 10/0= nan nan/nan= nan nan/0= nan

octaveで関数プロット時に"error: expecting integer index..."と表示されエラーになる

octave では変数や関数を自作する事ができるが，変数と関数に同じ名前をつけることが出来ない．関数を変数で上書きした後に関数を呼び出すと，"error: expecting integer index..."とエラーになる．

octave:1> sin(pi)
ans =  1.2246e-16
octave:2> sin=10
sin =  10
octave:3> sin(pi)
error: expecting integer index, found 3.141593 ← エラーになる
octave:3> clear sin
octave:4> sin(pi)
ans =  1.2246e-16
octave:5> 

はまったのでメモ．

参考リンク．ハンガリー語のスライド．
http://xymarkus.web.elte.hu/Jegyzet/Fiznum_1_ea-10-11_csabai/lecture02.pdf