且构网

分享程序员开发的那些事...
且构网 - 分享程序员编程开发的那些事

使用CLAPACK未定义的引用错误

更新时间:2022-10-25 20:21:53

问题解决..我不得不说初学者的文档是可怕的...



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


I am trying to get clapack working on ubuntu 14... But it keeps telling me undefined reference to "ssyev_"(actually it happens to any of those c functions in "clapack.h")

I just used the cmd "g++ lineartest3.cpp -o linear3"

I am pretty sure i followed the instructions shown below.

  • cd CLAPACK-3.2.1
  • cp make.inc.example make.inc
  • make f2clib
  • make blaslib
  • make (this takes a while)
  • 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**

And my little codes.

#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 help ..

Problem solved.. I have to say the documents for beginners are terrible...

copy blas_LINUX.a lapack_LINUX.a F2CLIBS/libf2c.a to /usr/local/lib and rename to libblas.a liblapack.a libf2c.a separately.

change the permission "chmod 777 libblas.a liblapack.a libf2c.a"

write a test.cpp and call "g++ test.cpp -llapack -lblas -lf2c -o test"