LAPACK在Win32上

本文介绍了LAPACK在Win32上的处理方法，对大家解决问题具有一定的参考价值，需要的朋友们下面随着小编来一起学习吧！

问题描述

我一直在探索需要对矩阵进行一些工作的算法，我已经在我的Linux机器上得到了一些简单的代码。下面是一段摘录：

I have been exploring algorithms that require some work on matrices, and I have gotten some straightforward code working on my Linux machine. Here is an excerpt:

extern "C" {
    // link w/ LAPACK
    extern void dpptrf_(const char *uplo, const int *n, double *ap, int *info);
    extern void dpptri_(const char *uplo, const int *n, double *ap, int *info);
    // BLAS todo: get sse2 up in here (ATLAS?)
    extern void dgemm_(const char *transa, const char *transb, const int *m,
            const int *n, const int *k, const double *alpha, const double *a,
            const int *lda, const double *b, const int *ldb, const double *beta,
            double *c, const int *ldc);
}

// in-place: be sure that (N*(N+1)/2) doubles have been initialized
inline void invert_mat_sym_packed(double *vd, int n) {
    int out = 0;
    dpptrf_("U",&n,vd,&out);
    ASSERT(!out);
    dpptri_("U",&n,vd,&out);
    ASSERT(!out);
}

// use with col-major ordering!!!
inline void mult_cm(double *a, double *b, double alpha, int m, int k, int n, double *c) {
    int lda = m, ldb = k, ldc = m; double beta = 1.0;
    dgemm_("N","N",&m,&n,&k,&alpha,a,&lda,b,&ldb,&beta,c,&ldc);
}

所有我需要做的是 sudo apt-get安装liblapack ，并链接库。

all I had to do was sudo apt-get install liblapack, and link against the library.

我现在试图让这个代码从MinGW使用32位的dll从，但我看到segfaults和无效的输出。我将继续使用gdb来确定错误的位置，但我怀疑有一个更好，更清洁，更便携的方式来完成这个。

I am now trying to get this code working from MinGW using the 32-bit dll's from here but I am seeing segfaults and invalid output. I will proceed with gdb to determine the location of the error but I suspect there's a better, cleaner, more portable way to get this done.

我做了什么，编译是安装fortran mingw（ mingw-get install fortran ）并链接到32位BLAS和LAPACK DLL从早期的链接。

What I did to get it to compile was install fortran for mingw (mingw-get install fortran) and link to the 32bit BLAS and LAPACK dll's from the earlier link.

我不知道我在这里丢失了多少...当用gcc编译win32编码时，其他人怎么会得到他们的LAPACK？

I'm not sure how much I'm missing here... How does everybody else get their LAPACK going when coding with gcc for win32?

我要找的是一个易于使用的C接口。我不想在所有的地方包装类。

What I'm looking for is an easy-to-use C interface. I don't want wrapper classes all over the place.

我试图找到一个英特尔MKL的下载...甚至是免费软件吗？

I tried to find a download for Intel MKL... Ain't even free software!?

LAPACK在Win32上

问题描述

推荐答案