マリー=アントワネット様の素晴らしいお言葉に勇気付けられながら，IA-32 + MSVC++でのlog1pの1実装． 以下の実装はIA-32 + MSVC++な環境でのみ有効なコードです． // log(1.0 + x)を計算する．特に|x| << 1.0で高精度であることが期待できる． double log1p…

VC7.1のリリースビルドだとexp(-∞)が無効な値になるようなアセンブリ吐くじゃんか，ぶ〜ぶ〜ぶ〜． // exp(x)のリリースビルドが吐くコード // st(0) <- x fldl2e; fmulp st(1), st; fld st(0); frndint; fxch st(1); fsub st, st(1); // x = -∞だとここが(-…