Ошибка при использовании dgesv_

Question

Ошибка при использовании dgesv_

Я пытаюсь решить простую проблему Ax = b с помощью функции dgesv_. Однако я столкнулся с проблемой, от которой я не могу избавиться. Мой код:

#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

У кого-нибудь есть идеи, как это решить?

Благодарю.

3

c++lapack linear

Решение

Другие решения

Других решений пока нет …

Источник

Accepted Answer

Ваша ошибка может быть связана с тем, что 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;
}

1