On Mon, Mar 21, 2011 at 6:11 PM, Jed Brown <[email protected]> wrote:

On Tue, Mar 22, 2011 at 00:06, Matthew Knepley <[email protected]

...

wrote:

...

3) It should go into TS

Really? It depends on a specific DM implementation. If you offered up a generic algebraic interface, then it could just be a normal explicit TS. But it seems to me that DM is a better choice for particle tracking, at least all the spatial queries and such.

Rationale: The function of characteristic is to handle the material derivative in a timestepping method. Caveats: I agree that the current implementation is dependent on DMDA. This should be remedied. Matt

On Tue, Mar 22, 2011 at 00:06, Matthew Knepley <[email protected]>wrote:

...

3) It should go into TS

Really? It depends on a specific DM implementation. If you offered up a generic algebraic interface, then it could just be a normal explicit TS. But it seems to me that DM is a better choice for particle tracking, at least all the spatial queries and such.

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener