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 <[email protected]>: > 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 > > > > >
