added saturation calculation onto IP and OP codes.

OP.update_results(sat = 0.5)
IP.update_results(IPsat = 0.5)