This was missing and rendered impossible simulations with lagrangian multipliers. Fixed.