乙嫁語り 5巻 (ビームコミックス)

DGETC2 っていうのは、LAPACK の補助サブルーチンで、完全ピボット選択で LU 分解を行うやつです。 Fortran の 2次元配列は、メモリ上では列優先で A(1,1), A(2,1), ... A(1,2), A(2,2), ... と並んでいるため、2次元配列全体にメモリアクセスする場合、同じ…

聖☆おにいさん(8) (モーニング KC)

id:yshl:20121209:1355059627 の続きですが、OCaml で Array ではなく List を使ってみたら、48 秒かかりました。まあ、明らかに Array 向きですよね。 あと、Array の方は unsafe オプションを付けると 0.9 秒になりました。 追記: Java で 0.9 秒。あと、L…