def old_integer_coeff(float_list):
relations = []
possibilities = []
for fl1 in float_list:
if almost_equal(fl1, 0.0, tolerance=TOLERANCE):
# No es un buen candidato para buscar MCM
continue
aux = []
fl1 = '%s'%abs(fl1)
rat1 = Fraction(fl1)
for fl2 in float_list:
if almost_equal(fl2, 0.0, tolerance=TOLERANCE):
aux.append(0.0)
continue
modif2 = -1 if fl2 < 0 else 1
fl2 = '%s'%(abs(fl2))
rat2 = Fraction(fl2)
rel = rat1/rat2
rel = format(modif2*float(rel), '.%dg'%TRUNCATE)
rel = Fraction(rel)
aux.append(rel)
relations.append(aux)
mcm = float(lcm(*[abs(x) for x in aux if x!= 0]))
integerified = [x and int(mcm/x) or 0 for x in aux]
possibilities.append(integerified)
#if all([x.is_integer() for x in integerified]):
# return [int(x) for x in integerified]
ret = None
sup = None
# Busco la "mejor" inecuación
# (en el sentido de menor +:abs() de los coef y el TI)
# de entre las posibles
for pos in possibilities:
max_val = 0
for coef in pos:
max_val += abs(coef)
if sup is None or max_val < sup:
sup = max_val
ret = pos
return ret
def _act(self, sensor_value):
super(PController, self)._act(sensor_value)
error = (self.target_temperature - self.temperature)
self.print(" (error: {})".format(round(error, 5)))
if almost_equal(error, 0.0, decimals=2):
self.steady_state_ticks += 1
if self.steady_state_ticks >= self.STEADY_STATE_THRESHOLD:
self.running = False
self.ticks -= self.STEADY_STATE_THRESHOLD # hack to make run() return the right val
self.print("Controller reached steady state in {} ticks".format(self.ticks))
else:
self.steady_state_ticks = 0
output = error * self.PROPORTIONAL_GAIN
self.temperature_delta = output
def integer_coeff(float_list):
minimum = None
abs_float = []
modif = []
for fl in float_list:
if almost_equal(fl, 0.0, tolerance=TOLERANCE):
fl = 0.0
if fl >= 0:
modif.append(1)
else:
fl = abs(fl)
modif.append(-1)
abs_float.append(fl)
if fl and (minimum is None or fl <= minimum):
minimum = fl
bigger_one = []
for fl in float_list:
fl = fl / minimum
bigger_one.append(int(floor(fl)))
return bigger_one
