カテゴリー: 日記

　結局，陽的RK，陰的RK+単純反復, +Newton, +簡易Newton, +簡易Newton+Tridiagonalization(SPARK流実装)と，陽的・陰的Euler法まで作ってしまった。アホかワシ。SPARK流実装についてはエラーが出なくなったという段階なので，明日以降デバッグ作業。今週中にまともに動くといいなぁ（希望的観測）。

　来週は紀要原稿兼JSIAM発表準備のため，ちと進みが遅くなる。とはいえ，とりあえず何としても固有値計算とcuBLASとMAMGMAの解説は終わらせたい。しかしsparse BLASの話のおさまりが悪くてどうしようか思案中也。密行列だけにしておこうかなぁとも思いつつ，MAGMAに疎行列の反復法まで入っちゃったからなぁ。この辺使ってみたいんだよなぁ。悩ましい也

SPARK流の実装をするために必要な記述を埋めつつ，最終目標の実装開始。頭がごちゃごちゃになってきたので明日もう一度仕切り直し。

　週末は案の定何にもできなかったので，ブロック三重対角行列のLU分解解説から始める。ついでにRow majorとColumn majorの解説も追加。思い出しながら書いていたので前進・後退代入部分まで手が回んなかった。残りは明日，プログラムも完成させたいぞと。

　午前中非常勤講義，午後は模擬講義と実験室の大掃除（の指示と後始末）をしてグッタリ。PCの移設も加わって更に面倒なことになった。係数の埋め込みのみ行って終いにする。明日だ明日！

　一応自分のルーチンとしては作ってあるんだが，今作っているLAPACK/BLASベースの実装の方が少なくとも倍精度に関しては良さげなので，係数もあることだし，SPARK流のRuductionを組み込んだルーチンも作っておくことに。

　卒研発表後にやったのはソースの切り分けと係数の埋め込みのみ。後は明日以降。今週中にケリがつけばよし。この辺で止めておかないと何の本だかわからなくなるしなぁ。