lapack

8熱

1答えて

LapackをPHP（http://php.net/manual/en/book.lapack.php）と一緒に使用するには、特異値の計算が含まれている唯一のPHPライブラリです。私はApache2を搭載したUbuntu 14.04サーバを使用しています。私はgfortranとliblapack-devを他のウェブサイトの指示どおりにインストールしましたが、それらはすべてどのように記述するかに

1熱

1答えて

make中にファイルがない（cholrl.a）

私はLuback-3.6.1をUbutunにインストールしました。make中に問題があり、ファイルがない（cholrl.a）。エラーメッセージはこれです。 [1]作る：終了ディレクトリ '/usr/local/src/lapack-3.6.1/BLAS/SRC'（CD SRC /変異体;作る）[1]作る：入力ディレクトリ「/ usr/local/src/lapack-3.6.1/SRC/VARI

1熱

2答えて

LAPACKE_zheevx（）が収束に失敗しました - C++で2 * DLAMCH（ 'S'）を使用してABSTOLを増加させる方法は？

これはC++のLAPACKE_zheevx（）関数による固有値計算の収束のための適切な公差（ "abstol"）の設定に関する質問です。デフォルト値「abstol」（すなわちabstol = -1）で固有値/固有ベクトルを計算するときにLAPACKE_zheev（）が収束しないとき、LAPACKマニュアルにabstol = 2 * DLAMCH（ 'S'）を設定するように指示されます。しかし、D

0熱

1答えて

多次元配列の製品では、軸の総和の有無をどのように整列させるのですか？

いくつかの繰り返しインデックスが合算され、それ以外のインデックスが合算されていない場合、配列演算を行う最良の方法は何ですか？これらの操作にはeinsumを使用する必要があるようですが、整列されているが合計されていないディメンションのフラグを持つtensordotの代替があった方が良いでしょう。誰かが、テンポルドのように動作する高速の数値ルーチン（lapackでおそらく？）を知っていますか？ ==

1熱

2答えて

lapackとblasをC++で呼び出す

私のC++コードでlapackとblasを使う必要があります.MKL、ACML、またはデフォルトのlapackとblasをLinuxシステムにリンクしたいと思っています。残念なことに、彼らはcで異なる慣習を持っています。例えば、MKL（mkl_blas.h）でZDOTCは zdotc(&result, &n, x, &incx, y, &incy); と、デフォルトのLAPACKおよびBLAS

1熱

1答えて

LAPACK dgesv VS MATLAB mldivide

私は有限差分法プログラムを書いています。私はIntel Math Kernel Libraryを使用しています。たとえば、1000x1000の行列Aと1000x1000の行列B. インテルMKLでは、cblas_dgemm（）関数を使用するA * Bは約600ミリ秒かかりました。 MATLABでは、A * Bに約800ミリ秒かかりました。私はMKLが非常に速いと思った。しかし、1000×1

0熱

1答えて

LAPACKとBLASサブルーチンへの未定義の参照

BLASとLAPACKがFortranでどのように動作しているかを理解しようとしているので、行列を生成して反転するコードを作った。は、ここで私が呼んでいるファイル内の行列Aがあるコード program test Implicit none external ZGETRF external ZGETRI integer ::M complex*16,allocatable,dimensi

1熱

1答えて

ラエルのzheevの下限と上限の入力に対して異なる固有ベクトルが得られる

次のプログラムで、zheevの 'L'（行列の下三角部分を読む）を 'U'（上三角入力）出力固有ベクトルが異なることが分かりました。誰でも助けますか？ program zheev_test INTEGER N PARAMETER (N = 4) INTEGER LDA PARAMETER (LDA = N) INTEGER LWMA

0熱

1答えて

固有ライブラリ関数はlapackと同じdsyev_

lapack/blasを使用していくつかのopensourceをダウンロードしました。これをSIMDコード生成のEigenベースのソースに変更したいと思います。 Eigenライブラリにはdsyevと同じ機能があり、LAPACKにありますか？ dsyveは、いくつかの目的でinfo値を返します。私が知る限り、eigensolverはEigenにライブラリeigenvalueまたはeigenvect

0熱

1答えて

Matlab mexでDGESVを使用する際のエラー

私はmexファイルでDGESVを使用して線形システムを解こうとしています。 2x2システムの場合、mexファイルは正常に動作し、エラーは発生しませんでしたが、システムが2より大きい場合、MATLABシステムエラーダイアログボックスのapperasが発生し、内部問題が発生し、終了する必要があります。 Imはコンパイルラインは、64ビットのウィンドウ10とインテル作曲XE 2013 にMatlabのr