Do you do everything mentioned here?
Note dependent mass matrices are less stably handled by Rosenbrock methods, so did you try DAEProblem with IDA?