This is implemented in a new method of apply_simple_constraint!, which is now called throughout the code.

This makes the IPOPT optimiser more robust

This is to support Pulseq, where this behaviour is allowed.