ただのバグだった件

昨日ごちゃごちゃと書きましたが、Gaussの消去法がうまくいかないのは、絶対値を取り忘れるという初歩的なバグが原因だったことが判明。このぶんだと、まだまだバグが仕込んでありそうで・・・
おまけに何も考えずに解かせると25分くらいかかるので、高速化したほうがよさそうな。Poisson方程式をtime stepごとに解かないといけないのですが、もし次のstep以降も反復法が失敗するようだと、25分 x step数とか、ちょっと大変・・・