lapack

    6

    3答えて

    誰かが私にライブラリ/コードを教えて、コレスキー分解で低ランクの更新を実行できるようにすることができますか? Matlabはこの機能を「cholupdate」という機能として提供しています。 LINPACKにもこの機能がありますが、LAPACKにまだ移植されていないため、(私の知る限り)LAPACKには移植されていません。サイフィ。 私はscikits.sparseがCHOLMODに基づいて同様の

    1

    1答えて

    私は行列に関するいくつかの作業を必要とするアルゴリズムを探求しています。私はLinuxマシンで作業する簡単なコードを手に入れました。以下はその抜粋です: extern "C" { // link w/ LAPACK extern void dpptrf_(const char *uplo, const int *n, double *ap, int *info); e

    2

    1答えて

    状況は単純ですが、奇妙です。 LinearAlgebra.oソース(LAPACKへのリンクが必要)を使わずにプログラムをコンパイルすると、C++の例外がキャッチされて処理されます。コンパイルユニットを含めずにライブラリ(-llapack -lblas)にリンクしている場合は、例外がキャッチされて処理されます。しかし、いったんそれを取得すると(コードは正常に実行されます)、C++の例外は正しく処理さ

    0

    1答えて

    私は多くのノルム、ドットプロダクト、そして最も重要なのはマトリックスベクトル乗算を実行する必要があるアプリケーションを持っています。 マトリックスおよびベクトルは、巨大です。 while(condition) /* usually iterations=dimension of matrix, so around 1 million iterations are *at least* requi

    3

    1答えて

    私は現在、たくさんのblasルーチンの呼び出しを伴うアプリケーションに取り組んでいます。私が発見したメモリリークを定期的にチェックして、dgemmコールでバイトが失われていることを確認しました。コールは次のようになります。 // I want to multiply 2 nxn matrices and put the result into C - an nxn matrix double z

    1

    2答えて

    私はより速く実行しようとしている線形代数コードを持っています。これまでのところ、私はMATMUL(Fortran Lib。)、DGEMVを使用しましたが、私自身のMVコードをOpenMPで作成しようとしましたが、アルゴリズムはスケーラビリティの面で優れていません。スピードアップは、私が割り当てているプロセッサの数に関係なく、わずか3.5 - 4です(私は64プロセッサを試しました)。 プロファイリ

    0

    1答えて

    MATLABに行列の逆行列をとるなどの単純な線形演算を実行できるMEXファイルを書くのに問題があります。私はVisual Studio 2010を使用して行列の逆行列を取り除くことに成功し、MEXファイルを正常に作成したので、問題を抱えている唯一のものはこれらの2つの概念を一緒にしています。私はMathWorksサイトから入手したMEXサンプルコードをコンパイルしようとしましたが、運がありません。

    2

    1答えて

    Linux gccにLAPACKをセットアップするには、助けが必要です。私はLAPACKには新しく、Fortranを使用する知識はありません。 私はlapack-3.4.0をダウンロードし、ライブラリを liblapack.aとlibrefblas.aにするようにしました。その後 、私は私のプログラムにこれらのライブラリをリンク: -llapack -lrefblas 私はdgetriなどdpot

    3

    1答えて

    を計算するdstevを使用してゼロ値を取得しています、私はインテルのmkl_lapack.hライブラリがインストールされています。 プログラムにはNxNの三重対角行列がありますので、行列の値を格納するために2つのベクトルを使用します。 "d"ベクトルが主対角線である場合、副対角線の値は "e"に格納されます。私はマトリックスが16×16であるため、その後、私は2つのベクトル(私はすべてのために一度

    9

    1答えて

    FortranとLAPACKを使って実際の対称行列を三重対角化したい。 LAPACKは基本的に2つのルーチンを提供します.1つはフル・マトリックスで動作し、もう1つはパックド・ストレージのマトリックス上で動作します。後者は確かに少ないメモリを使用していますが、スピードの違いについて何か言えることがありますか?