MathGroup Archive 2001

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

Search the Archive

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



  • Prev by Date: Re: Simplify
  • Next by Date: Searching over real values only
  • Previous by thread: no math font in exporting plots to .eps files
  • Next by thread: Re: NDSolve-memory