使用dgesv_时发生错误

Question

我正在尝试使用dgesv_函数解决一个简单问题Ax = b。但是我遇到了一个我无法克服的问题。我的代码是：

#include <cstdio> 
#include <f2c.h> 
#include <clapack.h> 

void main(void) 
{ 
    /* 3x3 matrix A 
    * 76 25 11 
    * 27 89 51 
    * 18 60 32 
    */ 
    double A[9] = {76, 27, 18, 25, 89, 60, 11, 51, 32}; 
    double b[3] = {10, 7, 43}; 

    int N = 3; 
    int nrhs = 1; 
    int lda = 3; 
    int ipiv[3]; 
    int ldb = 3; 
    int info; 

    dgesv_(&N, &nrhs, A, &lda, ipiv, b, &ldb, &info); 

} 

我认为的代码是正确的，但是，每当我运行它，我得到以下错误：

LINK : C:\...\Eye Tracker.exe not found or not built by the last incremental link; performing full link 
1>  Creating library C:\...\Eye Tracker\Debug\Eye Tracker.lib and object C:\Users\Daniel\documents\visual studio 2010\Projects\Eye Tracker\Debug\Eye Tracker.exp 
1>ellipse_fit.obj : error LNK2019: unresolved external symbol "void __cdecl dgesv_(int const *,int const *,double *,int const *,int *,double *,int const *,int *)" (?dgesv_@@YAXPBH0PAN0PAH102@Z) referenced in function "void __cdecl ttt(void)" (?ttt@@YAXXZ) 
1>C:\Users\Daniel\documents\visual studio 2010\Projects\Eye Tracker\Debug\Eye Tracker.exe : fatal error LNK1120: 1 unresolved externals 

任何人有任何想法如何解决这个问题？

谢谢。

Answer 1

您的错误可能来自LAPACK未链接到您的程序。 CLAPACK很容易链接到C程序，但链接到C++需要添加几行。据http://wwwx.cs.unc.edu/~cquammen/wp/2010/08/12/calling-clapack-code-from-c/以下行应该做的伎俩：

extern "C" { 
    #include <f2c.h> 
    #include <clapack.h> 
} 

如果没有足够的，这里是一个代码的工作件，它可以通过g++ main.cpp -o main -llapack -lblas -lm编译：

#include <iostream> 

using namespace std; 
extern "C" 
{ 
void dgesv_(int* n,int* nrhs,double* a,int* lda,int* ipiv, double* b,int* ldb,int* info); 
} 

int main(void) 
{ 
    /* 3x3 matrix A 
    * 76 25 11 
    * 27 89 51 
    * 18 60 32 
    */ 
    double A[9] = {1, 0, 0, 0, 2, 0, 0, 0,4}; 
    double b[3] = {42, 84, 168}; 

    int N = 3; 
    int nrhs = 1; 
    int lda = 3; 
    int ipiv[3]; 
    int ldb = 3; 
    int info; 

    dgesv_(&N, &nrhs, A, &lda, ipiv, b, &ldb, &info); 

    cout<<"solution : "<<b[0]<<" "<<b[1]<<" "<<b[2]<<endl; 
}