Matrix is singular to working precision problem

The problem derives from the internal representation of the model. The inclusion of an auxiliary variable and corresponding equation can cause numerical problems as it may introduce spurious terminal conditions.