lapack

1熱

1答えて

C++、LAPACK：ライブラリがリンクされた `dsyev 'への未定義参照

私はLAPACKのライブラリサブルーチンdsyevを使って行列を対角化する必要があるプログラムを書いています。まず、BLASライブラリとLAPACKライブラリを正しくリンクしていることを確認するためにテストコードを使用しました。このコードは次のとおりです。 // LAPACK test code //compile with: g++ -Wall -L/lib64 -llapack -lbl

0熱

1答えて

win7x64マシンでC、C++、Fortranのコンパイラを見つける方法は？

私はLAPACK 3.7.0のバイナリファイルをwhat is told hereに基づいてCMAKEを使ってビルドしようとしています。私のVisual Studio 2013 C++ projectでそれらを使用するために。しかし、私は本当にここでどのオプションを選択すべきか分かりません。また、次のコンパイラを選択する必要がありますか？私はProgram FilesまたはProgram Fi

1熱

1答えて

C++（LAPACK、sgels）とPython（Numpy、lstsq）の結果の違い

C++とPython計算の数値結果を比較しています。 C++では、私はLAPACKのsgels関数を使って、線形回帰問題の係数を計算します。 Pythonでは、同様の作業にNumpyのlinalg.lstsq関数を使用します。 sgelsとlinalg.lstsqで使用される方法の数学的な違いは何ですか？結果（すなわち、回帰係数）を数値的に比較すると、予想される誤差（6桁など）はどのくらいですか

1熱

1答えて

Lapackルーチンのセグメンテーションフォールトdstev

現在、C++の有限差分法を使って微分方程式を解く小さなプログラムを書いています。問題は私の三重対角行列を解くためにlapackを使用していることですが、ルーチンが呼び出された後にSegmentation Faultが発生しています。ルーチンはinfo = 0で終了しますので、問題の内容がわかりません。ここにプログラムがあります。 #include <iostream> #include <cma

2熱

1答えて

タイプA * X = Bの方程式を解くdgtsv_またはsgtsv_を使用

SWIFTでA * X = B型の線形方程式の系を解こうとしています。これは、O（N^2）メモリーを消費するLU分解に基づくアルゴリズムを使用してこれを実行できました。私の配列は一般的に大きく（10000サンプル以上）、LAPACKはO（N）メモリ空間しか消費しないトリジアゴナル行列に特有のいくつかの機能を持っています。&がより効率的です。 http://www.netlib.org/lapa

0熱

1答えて

lapsackのcgesvの "info> 0"エラーはどうなりますか？

は、私は、設定された線形方程式を解くためにLAPACKを使用しようとしているが、私はプログラムを実行するとき、私はそれは要素U（10011001）が正確にゼロであることを述べインフォメーション> 0でエラーが発生します。因数分解は完了しましたが、因数Uは正確に特異なので、解を計算できませんでした。これはどういう意味ですか？私はcgesvサブルーチンを使用しています。私はAとBの1001×1

0熱

1答えて

lapackがlibgcc_s_sjlj-1.dllについて文句を言う

私は、Visual Studioでarmadilloを使用するプログラムを作成しようとしています。私はhttp://icl.cs.utk.edu/lapack-for-windows/lapack/からlapackとblas用のPrebuiltライブラリをダウンロードしました。 MinGWをダウンロードし、C：/ MinGW/binをシステムPATHに追加しました。 C：/ MinGW/binには

0熱

3答えて

デフォルトのGCCコンパイラをLinux上でMPIで使用するために変更する方法CentOS

Linux（CentOS）マシンに2つのGCCコンパイラがインストールされています。古いバージョンのGCC（4.4.7）はデフォルトのフォルダ（CentOSに付属）にあり、新しいものは/usr/local/gcc/4.9.3/にあります。私のコードでは、MPIとLAPACK/LAPACKE/BLASライブラリを使用しています。古いGCCでは、次のようにソースをコンパイルしています（例えば "mai

0熱

2答えて

SpectraのSymEigsShiftSolverは固有ベクトルを返しませんか？

大規模な疎行列の固有値を計算するためにSymEigsShiftSolverを実装しましたが、固有ベクトルは返しません。まだ実装されていないのですか？ void Eigens::computeEigenvectors(Matrices m) { SparseSymShiftSolve<double> op(m.Lpl); SymEigsShiftSolver< double, SMALLEST