問題タブ [lapack]

$A \times B$ の SVD を見つける必要があります。ここで、$A$ は上三角行列、B は下三角行列です。LAPACKでこれを行う効率的な方法とルーチンは何ですか?

lapack lib から F77_NAME(dgeqrf) 関数を呼び出そうとしていますが、プログラムがクラッシュします。何が間違っていると思いますか？

ac プログラムで Lapack lib の F77_NAME(dgeqrf) 関数を使用して qr 分解を計算したいと考えています。

行列 3x3 の場合:

出力 3x3 (R 行列と Q の構築に使用されるいくつかのベクトルの組み合わせ)(線形形式) を取得します。

次に、Lapack の F77_NAME(dorgqr) を使用して Q 行列を抽出し、出力 3x3 (線形形式) を取得します。

これはウィキペディアから取った例で、私の Q はウィキペディア Q とは異なるようです: http://en.wikipedia.org/wiki/QR_decomposition#Example_2

Fortran と C の配列表現の違いが原因でしょうか? 初期行列の転置は問題を解決しますか?

Blasいくつかのルーティンを使用するためにライブラリを使用したいのですがLapack、でそれらを使用する方法がわかりませんVisual C++ 2010。

このコンテキストでそれらを使用する方法は？

Visual C/C++ 2010 で LAPACK と BLAS を使用したい (Visual C/C++ 2010 に LAPACK と BLAS ライブラリを追加する) ため、ここで非常によく似た質問を読んでいました(同じですが、Visual C/C++ 2008 と同じです)、彼は主張します彼はそれを次のように解決しました：

私は自分の問題に対する答えを見つけました。問題は、LAPACK .lib ファイルを含むディレクトリをライブラリ ディレクトリのリスト ([ツール] -> [オプション] -> [プロジェクトとソリューション] -> [VC++ ディレクトリ] -> [ライブラリ ファイル]) に配置する必要があることを知っていたのに、そうしなければならないことを忘れていたことでした。 Project -> (Project Name) Properties -> Linker -> Input の「Additional Dependencies」に追加します。

Tools -> Options -> Projects and Solutions -> VC++ Directories Visualを取得した後、彼のソリューションに従おうとしたとき、C/C++ 2010 は VC++ ディレクトリが使用されていないと言います...

では、Visual C/C++ 2010 が認識できるように .lib を追加するにはどうすればよいでしょうか。

例のように、マトリックスとして見ることができる結果テーブルを取得したとします (マトリックスは、のテーブルに適用されたStore Procedureまたはの結果です) 。UDFSQL Server 2008

したがって、次のようなものが得られたとします。

この種の結果 (ストア プロシージャまたは UDF から取得) を操作したいのですがoutside the DBMS、次の 2 つの方法で行いたいと思います。

1.-Send the complete matrixある種のバイナリファイルとして、つまりavoid passing through a text file

MATLAB（それから私はorLAPACKを使いたいですread this binary file and manipulate it）。

2.- この完全な行列を、または行列の情報を持つdirect memory pointerサイズのメモリのチャンクとして送信し、次に、またはこのポインタまたはメモリのチャンクを送信し、行列を操作できるようにします。1 MBMATLABLAPACK would access

• この手順のいずれかが実行可能ですか?, もしそうなら
• 手順は何ですか？
• このすべてを で実行できますStore procedureか? Visual C/C++ を使用する必要がありますか?

LAPACK bindingsそれらを使用してメモリのチャンクを読み取るにはどうすればよいですか？

• この行列「テーブル」を含むメモリチャンクを作成するにはどうすればよいSQL Server 2008ですか？
• それが不可能な場合、マトリックスを含むメモリポインタは実行可能ですか？

R で lapack lib を使用して行列 A の QR 分解を取得するための C コードを作成しています。結果は、R 言語コマンドを使用して取得した結果とは異なります。

これはlapackの問題ですか、それとも私のコードのバグですか?

行列の場合 (行優先) :

r 結果は :

私の結果：

この 1 週間、独立したサードパーティによって生成されたファイルを使用して、大規模な静的ライブラリ (1000 以上の C ファイル) を構築してきました。

現在、私はいくつかの激しい計算を行うアプリをプログラミングしており、CLAPACK ライブラリを追加する必要があるようです。私は加速フレームワークが進むべき道だと信じていますが、私はそれを動作させることができないようです:

フレームワークがないと、cblas_zgemm と dgetrf が不足していることを示すリンカー エラーが多数表示されます (合計で 12 以上のエラーが 30 個のファイルに分散しています)。

フレームワークをプロジェクトに追加した後、コードを 1 イオタ変更せずに、cblas_zgemm リンカー エラーは消えますが、dgetrf は残ります (フレームワークに属しているにもかかわらず)。

を追加する#include (or #import) <Accelerate/Accelerate.h>と、1000 をはるかに超えるリンカ エラーが発生し、一連の列挙子 (CblasTrans など) を再宣言していて、xyz の型が競合していることがわかります... エラー メッセージが繰り返されているようですが、dgetrf は未定義のままです。 .

私が間違っていること/何が欠けていますか?

ありがとう

編集：追加後の完全なエラーメッセージは次のとおりです。

編集 2:フレームワークを追加した後の元のリンカー エラー:

Octave と R を使用して、単純な行列を使用して SVD を計算し、2 つの異なる答えを得ています。コードを以下に示します。

R

オクターブ

私は Octave と R の両方に慣れていないので、最初の質問は、これを正しく行っているかということです。もしそうなら、どれが「正しい」ですか？どちらも正しいですか？また、numpy でこれを試し、LAPACK 関数 dgesdd および dgesvd を直接呼び出しました。Numpy は Octave に似た答えをくれ、LAPACK 関数を呼び出すと R に似た答えが返ってきます。

ありがとう！