使用CLAPACK未定义的引用错误

时间:2014-06-26 12:50:43

标签: c++ linux

我试图让clapack在ubuntu 14上工作......但它一直告诉我对“ssyev _”的未定义引用(实际上它发生在“clapack.h”中的任何c函数中)

我刚使用了cmd“g ++ lineartest3.cpp -o linear3”

我很确定我按照下面的说明进行操作。

  • cd CLAPACK-3.2.1
  • cp make.inc.example make.inc
  • 制作f2clib
  • make blaslib
  • make(这需要一段时间)
  • sudo cp INCLUDE / f2c.h / usr / include
  • sudo cp INCLUDE / clapack.h / usr / include
  • sudo chmod 644 /usr/include/f2c.h
  • sudo chmod 644 /usr/include/clapack.h **

我的小代码。

#include <iostream>
using namespace std;
extern "C" {               
#include "f2c.h"           
#include "clapack.h"       
}


int main()
{
int i,j;
float *A = new float[16];
A[0] = 1; A[1] = 1; A[2] = 1;  A[3] = 1;
A[4] = 1; A[5] = 2; A[6] = 3;  A[7] = 4;
A[8] = 1; A[9] = 3; A[10]= 6;  A[11]= 10;
A[12]= 1; A[13]= 4; A[14]= 10; A[15]= 20;
cout << "A = (";
for (i=0; i<15; i++) cout << A[i] << ",";
cout << A[15] << ")" << endl;

// use call-by-reference parameters, in another adaptation to Fortran
int n = 4;
char jobz = 'V';
char uplo = 'U';
integer N = n;
integer LDA=n;
float* eigvalue=new float[n]; 
int worksize;    // must have lwork>=3n-1, so n^2 is not big enough for n=2!
if (n==2) worksize = 5; else worksize = n*n;
float *work; 
work=new float[worksize];
integer lwork = worksize;
integer info;    // exit status

ssyev_(&jobz, &uplo, &N, A, &LDA, eigvalue, work, &lwork, &info);

cout << "info = " << info << endl;
if (info>=0)

{

    cout << "Eigenvalues = (" 

     << eigvalue[0] << "," << eigvalue[1] << "," 

     << eigvalue[2] << "," << eigvalue[3] << ")" << endl;

    cout << "Correct answer in GvL: (.0380, .4538, 2.2034 (typo?), 26.3047)" 

     << endl;

    // eigenvectors are stored in A (overwriting the original matrix):

    // ith eigenvector = ith column of A

    for (i=0; i<4; i++)

      {

    cout << "Eigenvector " << i << " = (";
    for (j=0; j<3; j++)
      cout << A[j+i*4] << ","; // jth elt of ith col=jth elt of ith evector
    cout << A[3+i*4] << ")" << endl;
      }
}
delete [] A;
delete [] eigvalue;  // be a good boy scout
delete [] work;
return 0;

}   

Plz帮助..

1 个答案:

答案 0 :(得分:1)

问题解决了......我不得不说初学者的文件很糟糕......

将blas_LINUX.a lapack_LINUX.a F2CLIBS / libf2c.a复制到/ usr / local / lib并分别重命名为libblas.a liblapack.a libf2c.a。

更改权限“chmod 777 libblas.a liblapack.a libf2c.a”

编写一个test.cpp并调用“g ++ test.cpp -llapack -lblas -lf2c -o test”