Enjoy our world with numerical computations!
SPARK流の実装，割とあっさり完成。まぁ一度作ってあるしな。ちょろっと動かしてみたが，やっぱり早いなぁ。しかしステップ幅制御がテキトーすぎるのでrejectが割と頻繁。この辺は動かしながら調整しよう。
固有値問題もコード整理して，非対称行列と対称行列の固有値・固有ベクトル計算完了。とはいえ，Intel Math KernelのC99のdouble complex, MKL_Complex16, magmaDoubleComplexの非互換性についての解説はどうしたもんやらと思案中。も～複素数ったら扱いが杜撰なのかみんなテキトーに実装しすぎだよ全く。大体，magmaDoubleComplexって実部と虚部がx, yってふざけてんのか。
まぁこの辺で第一稿の完成が見えてきたのは喜ばしい。今後は
あたりを埋めたら全体を整理して3月中に関係者にチェックを希う予定。
結局，陽的RK，陰的RK+単純反復, +Newton, +簡易Newton, +簡易Newton+Tridiagonalization(SPARK流実装)と，陽的・陰的Euler法まで作ってしまった。アホかワシ。SPARK流実装についてはエラーが出なくなったという段階なので，明日以降デバッグ作業。今週中にまともに動くといいなぁ（希望的観測）。
来週は紀要原稿兼JSIAM発表準備のため，ちと進みが遅くなる。とはいえ，とりあえず何としても固有値計算とcuBLASとMAMGMAの解説は終わらせたい。しかしsparse BLASの話のおさまりが悪くてどうしようか思案中也。密行列だけにしておこうかなぁとも思いつつ，MAGMAに疎行列の反復法まで入っちゃったからなぁ。この辺使ってみたいんだよなぁ。悩ましい也
SPARK流の実装をするために必要な記述を埋めつつ，最終目標の実装開始。頭がごちゃごちゃになってきたので明日もう一度仕切り直し。
週末は案の定何にもできなかったので，ブロック三重対角行列のLU分解解説から始める。ついでにRow majorとColumn majorの解説も追加。思い出しながら書いていたので前進・後退代入部分まで手が回んなかった。残りは明日，プログラムも完成させたいぞと。
午前中非常勤講義，午後は模擬講義と実験室の大掃除（の指示と後始末）をしてグッタリ。PCの移設も加わって更に面倒なことになった。係数の埋め込みのみ行って終いにする。明日だ明日！