前回うまくいったと思ったプログラムは思いっきり間違っていました。で、直してみると…全然速度がでない。この版の解はあってるんだけど…。
それで、微妙に表現を変えてみると…解が変わる。どうやら、微妙な最適化の具合で、私のプログラムは動いているらしい…。もしかしたら、最適化の具合が悪いせいで丸め誤差が発生するから、今までまともに計算できた人がいないのではないか?とさえ思ってきた…
今までは、ダラッと数式をプログラム上に書いていたので、コンパイラによる最適化の自由度があったんではないかなぁ…