Venho explorando algoritmos que exigem algum trabalho em matrizes e obtive algum código simples trabalhando na minha máquina Linux. Aqui está um trecho:

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);
}

tudo o que eu tinha que fazer erasudo apt-get install liblapack e vincule-o à biblioteca.

gora estou tentando obter esse código funcionando no MinGW usando as dll de 32 bits deaqu mas estou vendo segfaults e saída inválida. Continuarei com o gdb para determinar a localização do erro, mas suspeito que haja uma maneira melhor, mais limpa e mais portátil de fazer isso.

que eu fiz para compilar foi instalar o fortran for mingw mingw-get install fortran) e link para as dlls de 32 bits BLAS e LAPACK do link anterio

Não tenho certeza do quanto estou perdendo aqui ... Como todos os outros usam o LAPACK ao codificar com gcc para win32?

O que estou procurando é uma interface C fácil de usar. Não quero aulas de wrapper por todo o lad

Tentei encontrar um download para o Intel MKL ... Nem sequer é um software gratuito!?