modindices() now work if linear equality constraints (using the "==") are included in the model