Orz日記 by Akio Morita

_ [SAD]Re/Im/Conjugate実装の怪

Listに対しては、Map的に作用する1引数数学関数として、tfefunrefから tfcmplxf@src/tfeexpr.f経由で実装しており、mode引数は Re/Im/Conjugateに対して1/2/3が割り当てられているのだが、List再帰条件がtflistq(it,ia) .or. mode .eq. 3で設定されているので、

Re[{1, 2+I, I}] -> {1, 2, 0}
Re[a[1, 2+I, I]] -> Re[a[1, (2+I), I]]
Im[{1, 2+I, I}] -> {0, 1, 1}
Im[a[1, 2+I, I]] -> Im[a[1, (2+I), I]]
Conjugate[{1, 2+I, I}] -> {1, (2-I), (-I)}
Conjugate[a[1, 2+I, I]] -> {1, (2-I), (-I)}

のような謎動作となる。 ちぐはぐな動作なので、実装バクっぽぃ

仮に、List以外のコンテナ型に対して作用させるなら

Conjugate[a[1, 2+I, I]] -> a[1, (2-I), (-I)]

を返すべきだと思うが、Complex型もコンテナ型なので Conjugate[Complex[a, b]]は、Complex[Conjugate[a], Conjugate[b]]のように評価すべきという話となり一貫性を欠く

他の数学関数と同様のframeworkで綺麗に実装出来るのだが、この変な挙動に依存したコードが存在したらどうしよう…

_ [SAD]OddQバグ

内部的にOddQ[x_Real] := EvenQ[x+1]で実装しており、内部表現がIEEE754倍精度浮動小数点なので、OddQ[2^53]がTrueになってしまう

2^53は、仮数部の最小単位が2になっているので、正の方向の次の浮動小数点は2^53 + 2となり隣の整数値が表現出来ない