Backward-Euler dualtime implementation
This will require several things:
-
Calculation of Jacobians. We should have a mechanism to either generate analytic jacobians or generate them numerically. The numeric generation may be a bit tricky. This may require independent graphs that can be executed to generate each column in the jacobian? For now, we should just mandate analytic Jacobians probably. See issue #20 (closed). -
Solving the pointwise linear system. See also James_Research_Group/SpatialOps#12, which provides support for doing this through SpatialOps.