Eclipseでlapackを用いたコード(C++)のコンパイル[Mac OS X ver. 10.10.3]

対象

EclipseでMac OS標準で装備されているlapackを用いたコード(C++)のコンパイル方法を知りたい方
ネットでの情報が少なかったため、メモを残す。

Mac OSのバージョン

$ sw_vers
ProductName:    Mac OS X
ProductVersion: 10.10.3
BuildVersion:   14D136

Eclipseのバージョン

Version: Luna Service Release 2 (4.4.2)
「Eclipse」 -> 「Eclipseについて」で確認

ヘッダファイルのシンボリックリンク作成

clapack.h等のヘッダファイルは以下のパスに格納されている。

$ ls -l /System/Library/Frameworks/Accelerate.framework/Frameworks/vecLib.framework/Versions/Current/Headers/
LinearAlgebra/    clapack.h         vBigNum.h         vDSP_translate.h  vecLib.h          vectorOps.h
cblas.h           vBasicOps.h       vDSP.h            vForce.h          vecLibTypes.h     vfp.h

Eclipseでビルドするため、以下のようにファイルのシンボリックリンクを作成する。

$ sudo ln -s /System/Library/Frameworks/Accelerate.framework/Frameworks/vecLib.framework/Versions/Current/Headers/* /usr/include/

Eclipseでの設定

1. パッケージエクスプローラにあるビルド対象プロジェクトを右クリックし、「プロパティ」をクリックする。
2. C/C++ビルドの「設定」をクリックする。
3. 「ツール設定」タブの「MacOS X C++ Linker」->「その他」-> リンカー・フラグに以下の文字列を設定する。 -Wl,-framework -Wl,Accelerate

以上で設定は完了です。

コンパイル

ビルド対象プロジェクトを右クリックし、「プロジェクトのビルド」を選択する。

以上でコンパイルは完了です。

以下の実対称行列対角化のサンプルコードを例に試してみて下さい。

#include <iostream>
#include <cmath>
#include <clapack.h>
using namespace std;

int main() {
    // 行列サイズ
    int M = 3;

    // 対角化する行列(固有ベクトルが格納される)
    double* A = new double[M * M];
    for (int i = 0; i < M; i++) {
        for (int j = 0; j < M; j++) {
            A[i * M + j] = 1.0;
        }
    }

    // 固有値
    double w[M];
    int info;
    int lwork = -1;
    double* work = new double[1];

    char* V = const_cast<char*> (((string)"V").c_str());
    char* U = const_cast<char*> (((string)"U").c_str());
    dsyev_(V, U, &M, A, &M, w, work, &lwork, &info);
    lwork = work[0];
    delete[] work;

    work = new double[lwork];
    dsyev_(V, U, &M, A, &M, w, work, &lwork, &info);
    delete[] work;

    if (info != 0) {
        cout << "error!" << endl;
        return EXIT_FAILURE;
    } else {
        cout << "eigen values:" << endl;
        for (int i = 0; i < M; i++) {
            printf("%f \t", w[i]);
        }
        cout << endl;

        cout << "eigen vector:" << endl;
        for (int i = 0; i < M; i++) {
            for (int j = 0; j < M; j++) {
                cout << A[i * M + j] << " \t";
            }
        }
    }
    delete[] A;
    return EXIT_SUCCESS;
}