Re: [Help-gsl] Re: exponential of a matrix

[Liam Healy]

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 <address@hidden> 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 <address@hidden> 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
> address@hidden
> http://lists.gnu.org/mailman/listinfo/help-gsl
>