interval_c[0] = c_start_point if fabs(delta_2) < fabs(delta_1): delta_1 = delta_2 if delta_1 > 0: interval_a[1] = a_start_point else: interval_a[0] = a_start_point return None b_beg = uniform(0, 1) b_end = uniform(4, 4) a = uniform(0.8, 4.2) c = uniform(0.4, 6.6) same_b_beg_a = formula.compare(b_beg, a)[1] same_b_beg_c = formula.compare(b_beg, c)[1] same_b_end_a = formula.compare(b_end, a)[1] same_b_end_c = formula.compare(b_end, c)[1] print a.args, b_beg.args, b_end.args, c.args print '-------------' print formula.compare(b_beg, a) print formula.compare(b_beg, c) print formula.compare(b_end, a) print formula.compare(b_end, c) print formula.compare(a, c) print '-------------' a, b_beg, b_end, c = find_config(same_b_beg_a, same_b_beg_c, same_b_end_a, same_b_end_c) print a.args, b_beg.args, b_end.args, c.args print '-------------'
from math import fabs from scipy.stats import uniform from spatiotemporal.temporal_events import RelationFormulaGeometricMean from spatiotemporal.temporal_events.trapezium import RelationFormulaTrapezium __author__ = 'keyvan' b_beg = uniform(10, 3) b_end = uniform(1001, 3) c_actual = uniform(11, 991) f = RelationFormulaTrapezium() same_beg = f.compare(b_beg, c_actual)[1] same_end = f.compare(b_end, c_actual)[1] interval = [10, 13.0] x = [] for i in xrange(40): c_start_point = (interval[0] + interval[1]) / 2 length = ((13.0 - c_start_point) / same_beg)**2 / 3.0 c_heuristic = uniform(c_start_point, length) same_heuristic = f.compare(b_end, c_heuristic)[1] same_heuristic_beg = f.compare(b_beg, c_heuristic)[1] delta = same_end - same_heuristic x.append(delta) if fabs(delta) < 0.000000001: print same_beg print same_heuristic_beg print same_end print same_heuristic
interval_c[0] = c_start_point if fabs(delta_2) < fabs(delta_1): delta_1 = delta_2 if delta_1 > 0: interval_a[1] = a_start_point else: interval_a[0] = a_start_point return None b_beg = uniform(0, 1) b_end = uniform(4, 4) a = uniform(0.8, 4.2) c = uniform(0.4, 6.6) same_b_beg_a = formula.compare(b_beg, a)[1] same_b_beg_c = formula.compare(b_beg, c)[1] same_b_end_a = formula.compare(b_end, a)[1] same_b_end_c = formula.compare(b_end, c)[1] print a.args, b_beg.args, b_end.args, c.args print '-------------' print formula.compare(b_beg, a) print formula.compare(b_beg, c) print formula.compare(b_end, a) print formula.compare(b_end, c) print formula.compare(a, c) print '-------------' a, b_beg, b_end, c = find_config(same_b_beg_a, same_b_beg_c, same_b_end_a, same_b_end_c) print a.args, b_beg.args, b_end.args, c.args