Unitary partitioning in general constraint preserving DAE integrators

A number of numerical algorithms have been developed for various special classes of DAEs. This paper describes a new variable step size, constraint preserving integrator for general nonlinear fully implicit higher index DAEs. Numerical implementation issues are discussed. Numerical examples illustrate the effectiveness of the new method.