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 non-linear fully implicit higher index DAEs. Numerical implementation issues are discussed. Numerical examples illustrate the effectiveness of the new method.