MathGroup Archive 2007

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

Search the Archive

Re: Memory problem using NDSolve in For Loop

  • To: mathgroup at smc.vnet.net
  • Subject: [mg83361] Re: Memory problem using NDSolve in For Loop
  • From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
  • Date: Sun, 18 Nov 2007 04:50:10 -0500 (EST)
  • Organization: The Open University, Milton Keynes, UK
  • References: <fhmhrn$atl$1@smc.vnet.net>

DBK wrote:

> I am running a number of simulations (sim) with different initial
> conditions (ic). The code I have written runs fine until I increase
> the number of simulations. In my "For" loop below, I thought that
> previous solutions (sol[sim]) were being cleared but this obviously is
> not happening. How do I clear sol and reduce my memory usage? Sorry if
> this is pretty basic stuff.
> 
> For[sim = 1, sim < 2500, sim = sim + 1,
>   sol[sim] =
>    NDSolve[{equations, a[0] == 10, b[0] == .45, c[0] == .30,
>      d[0] == .40, e[0] == ic[[sim, 1]], f[0] == ic[[sim, 2]]}, {a, b,
>      c, d, e, f}, {t, 0, time}, MaxSteps -> Infinity,
>     WorkingPrecision -> 5, MaxStepSize -> 0.10];
>   AppendTo[outlist, {ic[[sim, 1]], ic[[sim, 2]], e[time] /. sol[sim],
>     f[time] /. sol[sim]}]];

As written, the for loop creates 2500 variables sol[1], sol[2], ..., 
sol[2500]. Apparently, you are just holding the result of NDSolve, 
sol[n], during the current cycle, say n, and then never touch it again. 
So why, in first instance, would you create so many symbols? You can 
just use one symbol, say sol, as in the following snippet:

For[sim = 1, sim < 2500, sim = sim + 1,
   sol = NDSolve[{equations, a[0] == 10, b[0] == .45, c[0] == .30,
      d[0] == .40, e[0] == ic[[sim, 1]], f[0] == ic[[sim, 2]]}, {a, b,
      c, d, e, f}, {t, 0, time}, MaxSteps -> Infinity,
     WorkingPrecision -> 5, MaxStepSize -> 0.10];
   AppendTo[
    outlist, {ic[[sim, 1]], ic[[sim, 2]], e[time] /. sol,
     f[time] /. sol}]];

Regards,
-- 
Jean-Marc


  • Prev by Date: Re: Minimizing Mathematica file size for source code control?
  • Next by Date: Re: memory release problem in mathematica6.0
  • Previous by thread: Memory problem using NDSolve in For Loop
  • Next by thread: Re: Memory problem using NDSolve in For Loop