nonlinear continuity conditions

[ghorn]

Basically right now g(x_k) == x_{k+1}, (where g is like an integration or extrapolation), but you might want to say something like c( g(x_k), x_{k+1} ) == 0. This could be useful for re-normalizing a quaternion or DCM, for example