昨日作ったMatrix Market (mtx)ファイル(Coordinateタイプ)を読み込んで係数行列とするJacobi反復法のプログラムを作成。まずはSparse BLASで動くようにした。
この後はIMKL→MAGMA→cuSPARSEと3バージョン作る予定なのだが，COO(rdinate)形式をサポートしてないのにmtxファイルを直接読み込むAPIのあるMAGMAの説明は変えなきゃいかん。cuSPARSEはCUDA 7で登場したcuSOLVERとの絡みも説明しようと思ったが，API見てたら行列分解ぐらいしかないのでどうしようかなと思案中。Thrustは完全C++で，直接cuSPARSEを扱えるかどうか調べている暇がないので，余裕があれば使うことにして今回はスルー予定。