Re: ODE systems
- To: mathgroup at smc.vnet.net
- Subject: [mg40629] Re: [mg40521] ODE systems
- From: Selwyn Hollis <selwynh at earthlink.net>
- Date: Fri, 11 Apr 2003 02:04:33 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
On Wednesday, April 9, 2003, at 01:32 AM, Kyriakos Chourdakis wrote: > Dear All, > > I have a few questions related to ODE systems: > > 1. Is there a special way for Mathematica to handle large > (e.g. 10000x10000) sparse matrices? Something similar > to the sparse-type commands of other systems perhaps. See Developer`LinearExpressionToSparseMatrix` in the Help Browser. > > 2. Having an ODE system of the form x'=A.x can Mathematica be > forced to use the Krylov-type approximations to > compute the Matrix exponential? If anyone has had any > experience with this approach, I would welcome your > comments on the speed/accuracy. I know the ExpoKit > package for another system and Fortran, but I have not seen > anything similar for Mathematica yet. > > 3. And finally a question that must be silly: does the > Krylov subspace approximation solve complex valued > ODEs as well, or it restricted to a real valued matrix > A? Beats me... Can you give me a reference on the topic? ----- Selwyn Hollis http://www.math.armstrong.edu/faculty/hollis