Implementation and test for #194
Use sloppy bounds on true/false to improve big-M stability. They will be presolved anyway in solver