Background: I am working on a energy system model (RMIP) where it minimizes the total system cost, while accounting for energy demand balancing constraints on a certain spatial and temporal resolution. I also have constraints on the upper limit of the possible energy production by the generation technologies, at the time slides of the temporal resolution. My solver is CPLEX.

Issue: Sometimes while CPLEX shows finding an optimal solution, the inequality constraint of production upper limit would get relaxed (and thus not truly an optimal solution) and lines of lines of “**slack removed by singularity**” are produced before the “optimal” solution was being found. Any advice on how to avoid of this happening and making this constraints more binding?

Thank you for any advice!

Also: I cannot share my code.