問題タブ [cusp-library]

cudaにcuspライブラリを使用するのは初めてです。CUDA の改訂されたシンプレックス アルゴリズムを実装しようとしています。そのためには、基本行列を更新するために 2 つの疎行列を乗算する必要があります。

問題は、cuspライブラリを使用して2つのスパース行列を乗算するにはどうすればよいですか(どの形式でも問題ありません)。また、結果行列に含まれる非ゼロ要素の数を知る方法はありますか(メモリ割り当ての目的で)?

私は試した：

cusp:multiply() は、私が見つけた乗算の唯一の関数でした。

カスプ v.0.4 CUDA v.5.5

CUSP ホスト機能も openMP と並行して実装されていますか? githubの議論では「openMPに機能が追加された」と書かれているため、少し混乱していますが、openMPで完全に実装されているという明確な声明はカスプにありません。

私はCUSPで次のことをしようとしています:

これは

一方、CUSP からは奇妙な回答が返ってきます。

結果は次のようになります

最後のものはあまり正しくないようです。私が間違っていることを知っている人はいますか？コードの使い方が間違っているのでしょうか、それとも非常に優れた推測ソリューションと前提条件を使用する必要があるのでしょうか?

CUSP を使用して、GMRES メソッドを使用して複雑な行列を解こうとしています。

コンパイルすると、「"cusp::complex" から "float" への適切な変換関数が存在しません」というエラーが表示されます。

エラーの原因を確認すると、gmres.inl の 143 行目が表示されます。

resid と s type はどこにありますか

両方の型が複雑なのに、なぜ float と表示されるのですか?

code:` // 空の疎行列構造を作成します (CSR 形式) cusp::csr_matrix,cusp::device_memory>A;

Fortran から呼び出した CUSP ソルバーを使用して、CG/BicGSTAB を繰り返し解きたいと考えています。転送を避けるために、Fortran データを CUSP に直接渡しています。コードはコンパイルされますが、実行時にフラグが立てられます。

コードのコアは言うまでもなく、印刷ストリームも発生していません。もちろんコードは準備段階ですが、どこが悪いのか気になります。

これが不可能な場合は、別のアプローチに移行できますが、正確性について確信が持てないため、決定できません。どんな助けでも大歓迎です。

CUSP と CUDA を使用して疎行列を csr 形式に変換しようとしています。CUSP ドキュメントに示されているコードを少し変更しました。

しかし、得られる結果は、インデックスが 1 シフトされた coo 形式の行列です。

助けはありますか？

ありがとう。