When I inquired about this before I was told it was unsupported because it lacked tests.
http://lists.gnu.org/archive/html/help-gsl/2008-08/msg00004.html A superficial test I tried was correct. It would nice if we get got some thorough tests so that it could be put officially into GSL. Liam On Sun, Dec 13, 2009 at 4:56 PM, Carl Boettiger <[email protected]> wrote: > Hi again, > > The command seems to work fine in my example below. Output in my example > seems to agree with octave's expm function when the third argument to > exponential_ss is small enough (type gsl_mode_t, I gather this a precision > control ). Anyone use this function or is some other method preferable? > > -Carl > > #include <stdio.h> > #include <gsl/gsl_linalg.h> > > int > main (void) > { > double a[9] = {1,2,3,4,5,6,7,8,9}; > double b[9] = {0}; > > gsl_matrix_view m > = gsl_matrix_view_array(a, 3, 3); > gsl_matrix_view em > = gsl_matrix_view_array(b, 3, 3); > > gsl_linalg_exponential_ss(&m.matrix, &em.matrix, .01); > > gsl_matrix_fprintf(stdout, &em.matrix, "%g"); > > return 0; > } > > > On Sun, Dec 13, 2009 at 1:14 PM, Carl Boettiger <[email protected]> wrote: > >> Hi, >> >> I'm looking to exponentiate a matrix. It seems there are functions for >> this in exponential.c under linalg, but I don't see anything about them in >> the documentation. How do I call the exponential of a matrix using >> gsl_linalg_exponential_ss() ? >> >> I understand that the matmult functions in linalg can be replaced by blas, >> but not sure why the exponential is undocumented. An example of how to >> compute this would be great. Thanks! >> >> -Carl >> >> -- >> Carl Boettiger >> Population Biology, UC Davis >> http://two.ucdavis.edu/~cboettig <http://two.ucdavis.edu/%7Ecboettig> >> > > > > -- > Carl Boettiger > Population Biology, UC Davis > http://two.ucdavis.edu/~cboettig > _______________________________________________ > Help-gsl mailing list > [email protected] > http://lists.gnu.org/mailman/listinfo/help-gsl > _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
