SPARK流の実装，割とあっさり完成。まぁ一度作ってあるしな。ちょろっと動かしてみたが，やっぱり早いなぁ。しかしステップ幅制御がテキトーすぎるのでrejectが割と頻繁。この辺は動かしながら調整しよう。
固有値問題もコード整理して，非対称行列と対称行列の固有値・固有ベクトル計算完了。とはいえ，Intel Math KernelのC99のdouble complex, MKL_Complex16, magmaDoubleComplexの非互換性についての解説はどうしたもんやらと思案中。も～複素数ったら扱いが杜撰なのかみんなテキトーに実装しすぎだよ全く。大体，magmaDoubleComplexって実部と虚部がx, yってふざけてんのか。
まぁこの辺で第一稿の完成が見えてきたのは喜ばしい。今後は
- 固有値問題，リダクションとQR分解の解説
- 疎行列の解説，BiCG, BiCGSTABの解説
- 斉次線型ODEの解析解計算
- 並列計算(PthreadとOpenMP)の解説
あたりを埋めたら全体を整理して3月中に関係者にチェックを希う予定。