線型方程式 続き

id:yshl:20121209:1355059627 の続きですが、OCaml で Array ではなく List を使ってみたら、48 秒かかりました。まあ、明らかに Array 向きですよね。
あと、Array の方は unsafe オプションを付けると 0.9 秒になりました。
追記: Java で 0.9 秒。あと、Lapack で完全ピボット選択で解く場合、dgetc2 で LU 分解してから dgetsc2 で解を出せばいいのか? Cで自分で解くと 1.7 秒程度なのに、Lapack で 5.5秒 かかってしまう。
追記: boost の ublas を使ったら 34 秒。
追記: Scheme (Stalin-0.11) 追加。0.77 秒程度で計算できるけど、計算結果の精度がちょっと悪い。他が 1e-13 以下なのに、Stalin は 1e-8 程度の誤差になった。