int rp_operator_newton::reduce(rp_box b) { bool proof = true; for (int i=0; i<_arity; ++i) { int v = _v[i]; rp_interval aux, save; rp_interval_add(aux,rp_ivector_elem(_unknown,i),rp_box_elem(_midpoint,v)); rp_interval_copy(save,rp_box_elem(b,v)); if (!(rp_interval_strictly_included(aux,save))) proof = false; rp_interval_inter(rp_box_elem(b,v),save,aux); if (rp_interval_empty(rp_box_elem(b,v))) { return( 0 ); } } if (proof) rp_box_set_safe(b); return( 1 ); }
void rp_operator_newton::compute_midpoint(rp_box b) { for (int i=0; i<_arity; ++i) { rp_interval_set_point(rp_box_elem(_midpoint,_v[i]), rp_interval_midpoint(rp_box_elem(b,_v[i]))); } }
void rp_operator_newton::compute_unknown(rp_box b) { for (int i=0; i<_arity; ++i) { rp_interval_sub(rp_ivector_elem(_unknown,i), rp_box_elem(b,_v[i]), rp_box_elem(_midpoint,_v[i])); } }
// Application function // b(x) := hull (b(x) inter initial_domain) int rp_operator_domain::apply(rp_box b) { DREAL_LOG_DEBUG << "rp_operator_domain::apply"; /* Rounding for discrete variables */ if (rp_variable_integer(_x)) { rp_interval_trunc(rp_box_elem(b,_id)); if (rp_interval_empty(rp_box_elem(b,_id))) { return( 0 ); } } /* Intersection with initial domain */ return( rp_union_inter_iu(rp_box_elem(b,_id),rp_variable_domain(_x)) ); }
// Application of operator to reduce the box b int rp_operator_piecewise::apply(rp_box b) { DREAL_LOG_DEBUG << "rp_operator_piecewise::apply"; // Check each piece Ij:Cj, if Cj violated then the elements of Ij // are removed from the domain of the main variable of _c for (int i=0; i<rp_ctr_piecewise_arity(_c); ++i) { int violated = 0, j = 0; while ((!violated) && (j<rp_ctr_piecewise_elem_size(_c,i))) { if (rp_ctr_num_unfeasible(rp_ctr_piecewise_elem_ctrnum(_c,i,j),b)) { violated = 1; } else ++j; } if (violated) { // domain restriction dom(var) := dom(var) \ Ij rp_interval aux; rp_interval_copy(aux,rp_box_elem(b,rp_ctr_piecewise_var(_c))); rp_interval_setminus(rp_box_elem(b,rp_ctr_piecewise_var(_c)), aux, rp_ctr_piecewise_elem_dom(_c,i)); if (rp_interval_empty(rp_box_elem(b,rp_ctr_piecewise_var(_c)))) { return( 0 ); } } } // Check whether the domain of the main variable of _c // intersects at least one Ij int intersect = 0, i = 0; while ((!intersect) && (i<rp_ctr_piecewise_arity(_c))) { if (!rp_interval_disjoint(rp_box_elem(b,rp_ctr_piecewise_var(_c)), rp_ctr_piecewise_elem_dom(_c,i))) { intersect = 1; } else ++i; } return( intersect ); }
void pprint_vars(FILE* out, rp_problem p, rp_box b) { for(int i = 0; i < rp_problem_nvar(p); i++) { fprintf(out, "%s", rp_variable_name(rp_problem_var(p, i))); fprintf(out, " is in: "); interval_cout_local(rp_box_elem(b,i), 6, RP_INTERVAL_MODE_BOUND); if (i != rp_problem_nvar(p) - 1) fprintf(out, ";"); fprintf(out, "\n"); } }
// Application of operator to reduce the box b int rp_operator_condvar::apply(rp_box b) { DREAL_LOG_DEBUG << "rp_operator_condvar::apply"; // check whether the guard is certainly satisfied if (rp_interval_included(rp_box_elem(b,_v),_dom)) { return( _o->apply(b) ); } else { return( 1 ); } }
// Application of operator to reduce the box b int rp_operator_3b::apply(rp_box b) { DREAL_LOG_DEBUG << "rp_operator_3b::apply"; // domain to be reduced double u = rp_binf(rp_box_elem(b,_v)); double v = rp_bsup(rp_box_elem(b,_v)); double x; // Reduction of left bound [a,x] RP_ROUND_UPWARD(); x = u + _improve * (v-u); if (x>=v) { // reduction of whole domain return _o->apply(b); } else { rp_box_copy(_baux,b); rp_bsup(rp_box_elem(_baux,_v)) = x; if (!_o->apply(_baux)) { rp_binf(rp_box_elem(b,_v)) = x; } } // Reduction of right bound [x,b] RP_ROUND_DOWNWARD(); x = v - _improve * (v-u); rp_box_copy(_baux,b); rp_binf(rp_box_elem(_baux,_v)) = x; if (!_o->apply(_baux)) { rp_bsup(rp_box_elem(b,_v)) = x; } return( !(rp_interval_empty(rp_box_elem(b,_v))) ); }
void rp_splitter_time::apply(rp_box_set& bs, int var) { rp_interval i1, i2; rp_box b1 = bs.remove_insert(); this->observe(b1, var); rp_box b2 = bs.insert(b1); if (quick_split_vars[&b1] == NULL){ quick_split_vars[&b1] = new std::set<int>(); } quick_split_vars[&b2] = new std::set<int>(*quick_split_vars[&b1]); if (is_time_variable(var) && !did_quick_split(&b1, var)){ if (this->real_hole(rp_box_elem(b1, var), rp_variable_domain(rp_problem_var(*_problem, var)), i1, i2)){ rp_interval_copy(rp_box_elem(b1, var), i1); rp_interval_copy(rp_box_elem(b2, var), i2); } else { // std::cout << "TIME SPLIT" << std::endl; rp_interval &vi = rp_box_elem(b1, var); double split = std::min(rp_binf(vi) + precision, rp_split_point(rp_binf(vi), rp_bsup(vi), 1000, 1)); // Real variable: [a,b] --> [center,b] and [a,center] rp_binf(rp_box_elem(b1, var)) = rp_bsup(rp_box_elem(b2, var)) = split; quick_split_vars[&b1]->insert(var); quick_split_vars[&b2]->insert(var); // rp_interval_midpoint(rp_box_elem(b1,var)); } } else if (rp_variable_integer(rp_problem_var(*_problem, var))){ if (this->integer_hole(rp_box_elem(b1, var), rp_variable_domain(rp_problem_var(*_problem, var)), i1, i2)) { rp_interval_copy(rp_box_elem(b1, var), i1); rp_interval_copy(rp_box_elem(b2, var), i2); } else { // no hole found // Integer variable: [a,b] --> [a+1,b] and [a,a] ++rp_binf(rp_box_elem(b1, var)); rp_bsup(rp_box_elem(b2, var)) = rp_binf(rp_box_elem(b2, var)); } } else { if (this->real_hole(rp_box_elem(b1, var), rp_variable_domain(rp_problem_var(*_problem, var)), i1, i2)) { rp_interval_copy(rp_box_elem(b1, var), i1); rp_interval_copy(rp_box_elem(b2, var), i2); } else { // Real variable: [a,b] --> [center,b] and [a,center] rp_binf(rp_box_elem(b1, var)) = rp_bsup(rp_box_elem(b2, var)) = rp_interval_midpoint(rp_box_elem(b1, var)); } } }
void rp_splitter_mixed_hybrid::apply(rp_box_set& bs, int var) { rp_interval i1, i2; rp_box b1 = bs.remove_insert(); this->observe(b1, var); rp_box b2 = bs.insert(b1); rp_box b1_copy; rp_box_clone(&b1_copy, b1); rp_box suggestion = m_ode_sim_heuristic->sim(b1_copy, var); DREAL_LOG_DEBUG << "rp_splitter_mixed_hybrid::apply() " << "suggestion = " << var << " [" << rp_binf(rp_box_elem(suggestion, var)) << ", " << rp_bsup(rp_box_elem(suggestion, var)) << "]"; if (rp_variable_integer(rp_problem_var(*_problem, var))) { if (this->integer_hole(rp_box_elem(b1, var), rp_variable_domain(rp_problem_var(*_problem, var)), i1, i2)){ if (rp_interval_included(i1, rp_box_elem(suggestion, var))){ rp_interval_copy(rp_box_elem(b1, var), i1); rp_interval_copy(rp_box_elem(b2, var), i2); } else{ rp_interval_copy(rp_box_elem(b1, var), i2); rp_interval_copy(rp_box_elem(b2, var), i1); } } else { // no hole found ++rp_binf(rp_box_elem(b1, var)); if (rp_interval_included(rp_box_elem(b1, var), rp_box_elem(suggestion, var))){ // Integer variable: [a,b] --> [a+1,b] and [a,a] rp_bsup(rp_box_elem(b2, var)) = rp_binf(rp_box_elem(b2, var)); } else { --rp_binf(rp_box_elem(b1, var)); ++rp_binf(rp_box_elem(b2, var)); rp_bsup(rp_box_elem(b1, var)) = rp_binf(rp_box_elem(b1, var)); } } } else { if (this->real_hole(rp_box_elem(b1, var), rp_variable_domain(rp_problem_var(*_problem, var)), i1, i2)) { if (rp_interval_included(i1, rp_box_elem(suggestion, var))){ rp_interval_copy(rp_box_elem(b1, var), i1); rp_interval_copy(rp_box_elem(b2, var), i2); } else { rp_interval_copy(rp_box_elem(b1, var), i2); rp_interval_copy(rp_box_elem(b2, var), i1); } } else { // Real variable: [a,b] --> [center,b] and [a,center] double mid = rp_interval_midpoint(rp_box_elem(b1, var)); rp_binf(rp_box_elem(b1, var)) = rp_bsup(rp_box_elem(b2, var)) = mid; double b1Intersection = std::min(rp_bsup(rp_box_elem(b1, var)), rp_bsup(rp_box_elem(suggestion, var))) - std::max(rp_binf(rp_box_elem(b1, var)), rp_binf(rp_box_elem(suggestion, var))); double b2Intersection = std::min(rp_bsup(rp_box_elem(b2, var)), rp_bsup(rp_box_elem(suggestion, var))) - std::max(rp_binf(rp_box_elem(b2, var)), rp_binf(rp_box_elem(suggestion, var))); if (b2Intersection > b1Intersection // rp_interval_included(rp_box_elem(suggestion, var), rp_box_elem(b2, var)) ){ // okay DREAL_LOG_DEBUG << "rp_splitter_mixed_hybrid::apply() " << "*[" << rp_binf(rp_box_elem(b2, var)) << ", " << rp_bsup(rp_box_elem(b2, var)) << "]* [" << rp_binf(rp_box_elem(b1, var)) << ", " << rp_bsup(rp_box_elem(b1, var)) << "]"; // cout << "*[" << rp_binf(rp_box_elem(b2, var)) << ", " << rp_bsup(rp_box_elem(b2, var)) // << "]* [" << rp_binf(rp_box_elem(b1, var)) << ", " << rp_bsup(rp_box_elem(b1, var)) // << "] \t"; } else if (b2Intersection < b1Intersection){ // rp_interval_included(rp_box_elem(suggestion, var), rp_box_elem(b1, var))){ // reverse double b2_sup = rp_bsup(rp_box_elem(b1, var)); rp_binf(rp_box_elem(b1, var)) = rp_binf(rp_box_elem(b2, var)); rp_bsup(rp_box_elem(b1, var)) = mid; rp_binf(rp_box_elem(b2, var)) = mid; rp_bsup(rp_box_elem(b2, var)) = b2_sup; DREAL_LOG_DEBUG << "rp_splitter_mixed_hybrid::apply() " << "[" << rp_binf(rp_box_elem(b1, var)) << ", " << rp_bsup(rp_box_elem(b1, var)) << "] *[" << rp_binf(rp_box_elem(b2, var)) << ", " << rp_bsup(rp_box_elem(b2, var)) << "]*"; // cout << "[" << rp_binf(rp_box_elem(b1, var)) << ", " << rp_bsup(rp_box_elem(b1, var)) // << "] *[" << rp_binf(rp_box_elem(b2, var)) << ", " << rp_bsup(rp_box_elem(b2, var)) // << "]* \t"; } else { DREAL_LOG_DEBUG << "rp_splitter_mixed_hybrid::apply() suggestion not found"; DREAL_LOG_DEBUG << "rp_splitter_mixed_hybrid::apply() " << "-[" << rp_binf(rp_box_elem(b2, var)) << ", " << rp_bsup(rp_box_elem(b2, var)) << "]- [" << rp_binf(rp_box_elem(b1, var)) << ", " << rp_bsup(rp_box_elem(b1, var)) << "]"; // cout << "-[" << rp_binf(rp_box_elem(b2, var)) << ", " << rp_bsup(rp_box_elem(b2, var)) // << "]- [" << rp_binf(rp_box_elem(b1, var)) << ", " << rp_bsup(rp_box_elem(b1, var)) // << "] \t"; } } } rp_box_destroy(&b1_copy); }