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