BLAS 2用Fortran vs

Question

私は多くのノルム、ドットプロダクト、そして最も重要なのはマトリックスベクトル乗算を実行する必要があるアプリケーションを持っています。

マトリックスおよびベクトルは、巨大です。

while(condition) 
/* usually iterations=dimension of matrix, so around 1 million iterations are *at least* required (if not more) */ 
matrix-vector multiplication 
3 dot prods 
2 norms 

私は現在、インテルMKLとインテル®Fortranを使用しています：行列の次元は、ループ構造がある100000x100000

なる傾向にあります。インテル®MKLのインテル®Cでコードを書き直すのに役立ちますか？ 誰かが（特にDGEMVの）あらゆる種類のベンチマークを実行しましたか？ コードの書き換えは大きな痛みですが、理由がわかったら書き直しても大丈夫です。

編集：私は間違っています：マトリックスの寸法は百万ではなく100000です。かなり重大なエラー：|

はい、行列は、密度が高く、密度が高い必要があります。です。 さらに、それは対称ではなく、正定でさえありません。 私のアルゴリズムはQMRの修正版です。

Answer 1

パフォーマンスは、CやFortranで完全に同じになります。これは、ライブラリ呼び出しをバッキングする実際の実装が同じであり、本質的にすべての時間がこれらのライブラリ呼び出しに費やされるためです。