26   typedef gsl_matrix* Matrix;  typedef
28   static Matrix createMatrix(int rows, int cols) { return gsl_matrix_alloc(rows,cols); }  in createMatrix()
30   static void free(Matrix& m) { gsl_matrix_free(m); m=0; }  in free()
32 …static void prod(const Matrix& m, const Vector& v, Vector& x) { gsl_blas_dgemv(CblasNoTrans,1,m,v,…  in prod()
33   static void cholesky(Matrix& m) { gsl_linalg_cholesky_decomp(m); }  in cholesky()
34 …static void cholesky_solve(const Matrix& m, const Vector& b, Vector& x) { gsl_linalg_cholesky_solv…  in cholesky_solve()
35   static void eigen_symm(const Matrix& m, Vector& eval, Matrix& evec)  in eigen_symm()
38     Matrix a = createMatrix(m->size1, m->size2);  in eigen_symm()
45   static void eigen_symm_gen(const Matrix& m, const Matrix& _b, Vector& eval, Matrix& evec)  in eigen_symm_gen()
48     Matrix a = createMatrix(m->size1, m->size2);  in eigen_symm_gen()
49     Matrix b = createMatrix(_b->size1, _b->size2);  in eigen_symm_gen()
61   typedef gsl_matrix_complex* Matrix;
63   static Matrix createMatrix(int rows, int cols) { return gsl_matrix_complex_alloc(rows,cols); }
65   static void free(Matrix& m) { gsl_matrix_complex_free(m); m=0; }
67   static void cholesky(Matrix& m) { gsl_linalg_complex_cholesky_decomp(m); }
68 …static void cholesky_solve(const Matrix& m, const Vector& b, Vector& x) { gsl_linalg_complex_chole…
69   static void prod(const Matrix& m, const Vector& v, Vector& x)
71   static void eigen_symm(const Matrix& m, gsl_vector* &eval, Matrix& evec)
74     Matrix a = createMatrix(m->size1, m->size2);
81   static void eigen_symm_gen(const Matrix& m, const Matrix& _b, gsl_vector* &eval, Matrix& evec)
84     Matrix a = createMatrix(m->size1, m->size2);
85     Matrix b = createMatrix(_b->size1, _b->size2);