NDSolve-memory
- To: mathgroup at smc.vnet.net
- Subject: [mg32024] NDSolve-memory
- From: Pieter-Jan.DeSmet at fys.kuleuven.ac.be
- Date: Wed, 19 Dec 2001 04:29:00 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
Dear readers, I have a system of differential equations F( y(t) ,y'(t) ,t )=0, with initial condition y(0) = y0. Here y(t) is a vector-valued function, let us denote the dimension of this vector by n. When I solve this system in Mathematica with NDSolve, the program uses too much memory and crashes if n is too big. (In my case this happens when n = 6). I am however only interested in the value of y at t=1 ( y(1) ). So, I think that Mathematica can be more memory efficient if it throws away all previous values of y(t) it does not need anymore in its difference scheme. Is there a method to impose this? Of course, the result of NDSolve would not give an InterpolatingFunction object any more. Other suggestions to reduce the needed memory are also appreciated. Lowering PrecisionGoal and AccuracyGoal does not seem have much effect on the memory. Pieter-Jan De Smet