PS: I remember Mr.Jens Kuska having replied to a similar post long ago
with a Runge-Kutta integrator that evaluates the soln at just one
point. I downloaded that but that does not seem to work with Matrix
DEs since it may have been written for an older Mathematica version!
To integrate systems of differential algebraic equations or polynomial
equations look at the IMTEK Mathematica Supplement:
or a BDF solver:
http://www.imtek.uni-freiburg.de/simulation/mathematica/IMSweb/imsTOC/Differential%20Equation%20Systems/System%20Theory/BDFDocu.html
(see example 4)
you might want to try a sub-space projection. this projects your original
system of 1920 degrees of freedom to say 100 degrees of freedom, which of
course are much faster to time integrate and use less memory.
hth,
Oliver
Oliver Ruebenkoenig, <ruebenko AT uni-freiburg.de>
NDSolve[eqs, vars, {t, tf, tf}]
instead of
NDSolve[eqs, vars, {t, ti, tf}]
The result is still returned as an InterpolatingFunction but it only contains
the data at the end of the integration tf.
If you make sure to use the same variable for the solution in the loop
then previous references will be freed when the new solution is allocated.
Mark Sofroniou,
Wolfram Research