Memory problem using NDSolve in For Loop
- To: mathgroup at smc.vnet.net
- Subject: [mg83337] Memory problem using NDSolve in For Loop
- From: DBK <boydkramer at gmail.com>
- Date: Sat, 17 Nov 2007 05:36:36 -0500 (EST)
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]}]];