Re: [Help-gsl] Help-gsl Eigenvectors and eigenvalues

[Juan Pablo Amorocho D.]

Hi Luke,

GSL 1.8 is rather old. The current version is 1.16. I would give it a try.
Even so, I would recommend you run the test battery after installation. I
have observed some issues in Mac OS X that don't appear on Linux with the
same  GSL version.

Hope this helps.


-- Juan Pablo

2015-06-18 9:14 GMT+02:00 Luke <address@hidden>:

> Hello
>
> I am using this example taken from gsl-ref.pdf:
>
>
> #include <iostream>
> #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");
>         }
>     }
>
>     return 0;
> }
>
>
> When I run the program I see:
>
> eigenvalue = 6.94149e-308
> eigenvector =
>
> Process returned -1073741819 (0xC0000005)   execution time : 2.658 s
> Press any key to continue.
>
>
>
>
> For other example codes like:
>
> computing the value of the Bessel function J0(x)
> or
> FFT routines for complex data
>
> taken from gsl-ref.pdf
> everything is fine.
>
> I am using gsl-1.8 on Windows 8 64-bit.
>
>
> I would be grateful for any help.
> Best regards,
> Luke
>
>
>
>
>