MathGroup Archive 1992

[Date Index] [Thread Index] [Author Index]

Search the Archive

solving large (53000 equations) systems of linear ODEs with Mma


I have a systems of 53000 linear ODEs:
                                       P'(t) = A P(t)
where P(t) is a vector of length 53000, and A is a very large matrix
(constant) with only 5-diagonal of non-zero elements.
I've solved this problem using standard Runge-Kutta (in C and Fortran)
and using the Krylov approximation method for the matrix exponential operator.
I'd like to do the same thing (i.e., solve the system above) using Mma, if possible.
So I ask you fellow Mma users: Has anyone used Mma for this purpose ? Can Mma
solve (numerically, of course) such a large number of linear ODEs ?

Thanks for any information.

   Roque
   oliveria at engin.umich.edu
   Nuclear Eng Dept, The U. of Michigan





  • Prev by Date: time consumption FindRoot
  • Previous by thread: time consumption FindRoot