gsl 의 사용 (우분투)

먼저, 우분투 소프트웨어 센터에서

gsl binary package,

gsl development package,

gsl library package 

를 선택해서 깐다

(필요 없는 것도 있을 것이나 일단 3개를 선택해서 깔아봤다)

여기 까지 깔면 컴파일해서 inclue 문제는 해결된다.

다만 링크 문제가 남는데

gcc example.c -lgsl -lgslcblas -lm

처럼 링크를 해주면 에러 없이 컴파일 된다.

다음은 예제이다.

http://www.gnu.org/software/gsl/manual/html_node/Eigenvalue-and-Eigenvector-Examples.html

#include <stdio.h>

     #include <gsl/gsl_math.h>

     #include <gsl/gsl_eigen.h>

     

     int

     main (void)

     {

       double data[] = { 1.0  , 1/2.0, 1/3.0, 1/4.0,

                         1/2.0, 1/3.0, 1/4.0, 1/5.0,

                         1/3.0, 1/4.0, 1/5.0, 1/6.0,

                         1/4.0, 1/5.0, 1/6.0, 1/7.0 };

     

       gsl_matrix_view m 

         = gsl_matrix_view_array (data, 4, 4);

     

       gsl_vector *eval = gsl_vector_alloc (4);

       gsl_matrix *evec = gsl_matrix_alloc (4, 4);

     

       gsl_eigen_symmv_workspace * w = 

         gsl_eigen_symmv_alloc (4);

       

       gsl_eigen_symmv (&m.matrix, eval, evec, w);

     

       gsl_eigen_symmv_free (w);

     

       gsl_eigen_symmv_sort (eval, evec, 

                             GSL_EIGEN_SORT_ABS_ASC);

       

       {

         int i;

     

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

           {

             double eval_i 

                = gsl_vector_get (eval, i);

             gsl_vector_view evec_i 

                = gsl_matrix_column (evec, i);

     

             printf ("eigenvalue = %g\n", eval_i);

             printf ("eigenvector = \n");

             gsl_vector_fprintf (stdout, 

                                 &evec_i.vector, "%g");

           }

       }

     

       gsl_vector_free (eval);

       gsl_matrix_free (evec);

     

       return 0;

     }

Here is the beginning of the output from the program,

     $ ./a.out

     eigenvalue = 9.67023e-05

     eigenvector =

     -0.0291933

     0.328712

     -0.791411

     0.514553

     ...