Option to not save previous q for Strang splitting in 3D PyClaw

[Grady Lemoine]

Hi all,

A notion just occurred to me regarding reducing memory usage.  With Strang splitting, all variants of Clawpack currently save the previous time step's q array before advancing by the first half time step of the source term, in case the hyperbolic step has to be rejected.  For some problems (I'm thinking of my poroelasticity code here, but the situation is the same for any linear problem and maybe well-behaved nonlinear ones) time steps will never rejected, so I was thinking that it might be worthwhile to add a run time option to not save the previous time step.

I think I'm going to go ahead and add this into my local copy of PyClaw, since it will give me a little more grid size to play with, but is the rest of the PyClaw community interested in this?

--Grady

[Kyle Mandli]

Seems like an easy and worthwhile thing to include.  Regardless, it will be easy to add it to the Python options without everyone have to change their setup.py files.

--
You received this message because you are subscribed to the Google Groups "claw-dev" group.
To unsubscribe from this group and stop receiving emails from it, send an email to claw-dev+u...@googlegroups.com.
To post to this group, send email to claw...@googlegroups.com.
Visit this group at http://groups.google.com/group/claw-dev?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.