一応陰的RK法の３バージョン(単純反復，Newton, 簡易Newton)はでけた。係数はRow MajorなのにJacobi行列はColumn Majorってのは変態かしらん？ でもTestsetのFortranルーチンをそのまんま使おうとするとColumn Majorの方がいいのよね。cuBLASにも対応楽だし。
つーことで，次は並列化→GPU化という段階。さて明日にはできるかしらん・・・。
ところでRow majorは「行方向格納」，Column majorは「列方向格納」って訳語を作ったけど，もうちっとマシな言葉にならないかしらん。