問題タブ [lapack]

いくつかの三角ソルバーを探していますが、2つのソルバーに出くわしました。1つはBLAS：dtrsm で、もう1つはLAPACK：dtrtrsです。見た目からすると、両方とも共通の機能を持っているように見えますが、dtrsmにはもう少し機能があります（システムを解決する前に右側をスケーリングします）。

知りたいのですが
1）これらの機能は他にどのように違うのですか？
2）同じ操作を実行する場合、どちらが速いですか？
3）（2）の答えが明らかでない場合、dtrsmがdtrtrsよりも提案されるのはいつですか。その逆も同様です。

LAPACK を使用する目的は、行列のコレスキーを計算することです。Fedora で C/C++ でプログラミングしていますが、どちらの lapack をインストールするかについて混乱しています - LAPACK と lapacke または clpack?

Lapack++2.5.4とATLAS3.8.4をFedora12にインストールしました。

Eclipseを使用してlapack++をテストする簡単なプログラムを作成しました。

私はこれらのパスをEclipseに設定しました：

インクルードパス：/trunk/lapack/lapackpp-2.5.4/include

ライブラリ：lapackpp

ライブラリパス：/ usr / local / lib

エラーなしでビルドしますが、実行しようとすると吐き出されます

私は解決策を探していましたが、解決策が見つからないようです。-lg2cを含めてみましたが、コンパイラが見つかりません。助けてください。

サイズの異なる行列の固有値と固有ベクトルを計算しようとしています。私は非常に単純な Fortran90 コードを使用しており、Ubuntu で動作する私のマシンで利用可能な Intel MKL パッケージに含まれる適切な Lapack ライブラリにリンクしてコンパイルしています。メッセージの末尾に「matrix_diag_01.f90」というコードが添付されています。「random」モジュールには、Numerical Recipes の「ran」乱数ジェネレーターが含まれています。コードは次を使用して適切にコンパイルされます

小さい行列が与えられた場合、実行可能ファイルはうまく機能します。ただし、サイズが 3000x3000 の行列の場合、奇妙な動作が発生します。最初にこのエラーが発生します

ただし、SSYEVD の呼び出しには 3 つのパラメーターしかありません。次に、固有ベクトルを返しますが、固有値は返しません。より大きなメモリを搭載した別のマシンでコンパイルして確認しましたが、結果は同じでした。

誰でも助けてもらえますか？

ありがとう！

lapack ライブラリclapckを使用したいのですが、Visual Studio 2008 の C++ プロジェクトでどのように使用できますか...

また、私はテンプレート数値ツールキットを見直していました

同様のことをしたことがありますか?Visual Studio C++ プロジェクトで何を使用することをお勧めしますか? 含める必要があるヘッダーを含む簡単な例を投稿できますか?

正方形で密な行列Bと、長方形で疎な行列Aがあります。

製品を効率的に計算する方法はありますB^-1 * Aか?

これまでのところ、私は（numpyで）使用しています

これが A のスパース性を生み出していると私は信じています。スパース メソッドを使用することを考えて numpy.sparse.linalg.spsolveいましたが、これには A ではなく B がスパースである必要があります。

物事をスピードアップする別の方法はありますか？

Accelerate (CLAPACK) で圧縮された列ストレージ マトリックス (ccs) を解決するために使用する関数/メソッドを知っている人はいますか?サンプル コード、サンプルへのリンク、マトリックスの解決方法に関するヒントを探しています。私はドキュメントを調べてきましたが、ほとんどすべてがより複雑なグラフィカルシステムに関係しており、一見似た方法が何百もあります。少し早いですがお礼を

プロジェクトArmadillo'sでlapackを含む「C++線形代数ライブラリ」コードを使用する方法はありますか?C#

簡単な例を投稿していただけませんか...

.dllを取得してプロジェクトにインポートすることを考えていましたがC#、これは実現可能ですか、それとも別の方法がありますか?

私は、並列計算とオプションとトリックの検索が初めてです。

^{私の現在の問題は、タイプ B = AX の 2 32 個}の独立した線形方程式を見つけることです。A[8x8]、B[8x1]、および X[8x1] の寸法。A は密行列です。

私は CLAPACK - zgelss 関数を Intel コア 7i プロセッサ、4 コア、8 スレッド、および Mac マシンで順番に使用しています。

さらに、私の行列サイズは、内部最適化を使用するのに十分小さいようです。

このループに並列処理を適用したい:

それを行うための私のオプションは何ですか？

数人が OpenMP の使用を提案しましたが、OpenMP について読んだ後、プリコンパイルされたコードの部分に OpenMP ディレクティブを適用できないことがわかりました。

そうすれば動作しますが、プリコンパイルされた関数内の一時変数をプライベート化できないため、安全ではない可能性があります。{正しい/間違っている/その他をコメントしてください}

可能であれば、どのようにすればよいかお聞きしたいです。いくつかのトリック....何でも...

私はいくつかのソフトウェアのリリースを行おうとしていますが、現在ビルドプロセスのスクリプトを使用しています。x86_64 linuxでLAPACKを静的にリンクするという、思いもよらなかった何かに固執しています。構成中AC_SEARCH_LIB([main],[lapack])は機能しますが、lapackユニットのコンパイルは機能しません。たとえば、undefiend reference to 'dsyev_'lapack/blasルーチンが見過ごされることはありません。

ライブラリがインストールされていることを確認し、同じ結果で静的にするための適切なオプションを使用してライブラリを自分でコンパイルしました。

数年前にLAPACKを初めて使用したときに使用した例を次に示します。これは動的に機能しますが、静的には機能しません。http: //pastebin.com/cMm3wcwF

私がコンパイルに使用している2つの方法は、次のとおりです。