How to solve PDAEs (incorporate mass matrix into PDE solver)

Is your DAE index 1? I.e., are all of your constraints functions of the algebraic variables?