結局,陽的RK,陰的RK+単純反復, +Newton, +簡易Newton, +簡易Newton+Tridiagonalization(SPARK流実装)と,陽的・陰的Euler法まで作ってしまった。アホかワシ。SPARK流実装についてはエラーが出なくなったという段階なので,明日以降デバッグ作業。今週中にまともに動くといいなぁ(希望的観測)。
来週は紀要原稿兼JSIAM発表準備のため,ちと進みが遅くなる。とはいえ,とりあえず何としても固有値計算とcuBLASとMAMGMAの解説は終わらせたい。しかしsparse BLASの話のおさまりが悪くてどうしようか思案中也。密行列だけにしておこうかなぁとも思いつつ,MAGMAに疎行列の反復法まで入っちゃったからなぁ。この辺使ってみたいんだよなぁ。悩ましい也