問題タブ [diagonal]

境界ブロック システムを解決する必要があります。ここには 2 つのバージョンがあります。正常に動作するシリアル バージョンと動作しないパラレル バージョン (MPI を使用) で、理由がわかりません... 誰かが私を理解するのを手伝ってくれるかもしれませんなぜ？

================================================== =============================== = アプリケーションプロセスの 1 つの不適切な終了 = 終了コード: 11 = クリーンアップ中残りのプロセス = 以下のクリーンアップ メッセージは無視できます ======================================= ========================================= アプリケーションが終了文字列で終了しました: セグメンテーション違反 (シグナル 11) これは通常、アプリケーションの問題を示しています。デバッグの提案については、FAQ ページを参照してください。

シリアル版:

パラレルバージョン:

これは res.h http://bpaste.net/show/yAQCCwFsZlQCCiw7GKC9/で、先生から提供されたものです。

下手な英語といくつかのイタリア語のコメントで申し訳ありません。コードは明確である必要があります。私の目標は、コードのシリアル部分を MPI で並列化することです。プロセッサの数は常に 4 の倍数です。よろしくお願いします。

R を使用して行列 L から対角要素 (diagL) を削除するにはどうすればよいですか? 私は以下を使用してみました：

しかし、私は0の数字を取得します...

LAPACK++、Armadillo、IT++、Eigen、またはその他のライブラリを選択する必要があります (Windows の場合)。

私がする必要があるのは、大きな (約 10,000*10,000) 行列が対角化可能かどうかを確認することだけです。もしそうなら、D=(P^(-1))*A*P. これは、できるだけ早く実行する必要があります。どのライブラリを使用すればよいかわかりません。

また、これらの各ライブラリの長所と短所を一般的に教えていただければ幸いです。

簡単な質問です。疎行列 A とインデックス リストb = [2, 3, 4 ...]があり、 を設定したいと思い A(2, 2), A(3, 3), A(4, 4)... = 0ます。A(b, b) = 0 を試してみました。結果は、A(2, 3), A(2, 4), A(3, 2), A(3, 4), A(4, 2), A(4, 3)この行列の一部の非対角要素もゼロに設定されています。これをどのように処理しますか?

これは、現在、配列の対角線に要素を追加する問題を解決する方法です。他の方法をいくつか試しましたが、満足できません。if二重ループとステートメントを使用するよりスマートな代替手段があるはずです。どんな提案でも大歓迎です。