Re: [R] Krylov subspace computations of matrix exponentials

From: Mario Valle <mvalle_at_cscs.ch>
Date: Thu, 14 Apr 2011 10:27:53 +0200

Look at ExpoKit (http://www.maths.uq.edu.au/expokit/) It implement Krylov methods and is really fast. Unfortunately no R wrapper that I know yet.

If you plan to implement it, I can provide my testing code. Ciao!

                     mario

On 14-Apr-11 10:16, Niels Richard Hansen wrote:
> I use the very nice expm functions available from the expm and Matrix
> packages. My understanding is that for large sparse matrices the
> currently best methods available are Krylov subspace methods, but
> they are as far as I can tell not implemented in either of the packages
> mentioned, nor in any other R package I have found.
>
> Does anybody know if Krylov subspace methods are available from
> any R package?
>
> If not, is there anybody working on this or planning to do so?
>
> If not, I might be interested in making such an implementation, because
> I find myself in the need of fast methods for very large sparse matrices.
>
> A third question. There exists a Fortran implementation called Expokit
> with an ad hoc license stating that
>
> "Permission to use, copy, modify, and distribute EXPOKIT and its
> supporting documentation for non-commercial purposes, is hereby
> granted without fee, provided that this permission message and
> copyright notice appear in all copies."
>
> If I choose to use Expokit and just write an interface to R, would
> there be any problems with having such a copyright notice in the package?
>
> Thanks, Niels
>

-- 
Ing. Mario Valle
Data Analysis and Visualization Group            | http://www.cscs.ch/~mvalle
Swiss National Supercomputing Centre (CSCS)      | Tel:  +41 (91) 610.82.60
v. Cantonale Galleria 2, 6928 Manno, Switzerland | Fax:  +41 (91) 610.82.82

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Thu 14 Apr 2011 - 08:30:38 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 14 Apr 2011 - 08:50:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive