void create_sliced_system(boost::filesystem::path input_file, boost::filesystem::path output_file,
vec_mp * linears, int num_to_add,
const WitnessSet & W)
{
#ifdef functionentry_output
std::cout << "create_sliced_system" << std::endl;
#endif
if (W.num_var_names()==0) {
std::cout << "trying to create a sliced system, but witness set does not have the variable names." << std::endl;
deliberate_segfault();
}
int *declarations = NULL;
partition_parse(&declarations, input_file, "func_input", "config", 0); // the 0 means not self conjugate.
free(declarations);
FILE *OUT = safe_fopen_write(output_file);
FILE *IN = safe_fopen_read("func_input");
fprintf(OUT,"INPUT\n\n");
copyfile(IN,OUT);
fclose(IN);
std::vector< int > indicators;
for (int ii=0; ii<num_to_add; ii++) {
indicators.push_back(rand());
}
for (int ii=0; ii<num_to_add; ii++) {
std::stringstream linname;
linname << "supp_lin_" << indicators[ii];
write_vector_as_constants(linears[ii], linname.str(), OUT);
linname.clear(); linname.str("");
}
for (int ii=0; ii<num_to_add; ii++) {
fprintf(OUT,"function supp_lin_%d;\n",indicators[ii]);
}
for (int ii=0; ii<num_to_add; ii++) {
fprintf(OUT,"supp_lin_%d = supp_lin_%d_1",indicators[ii],indicators[ii]);
for (int jj=1; jj<W.num_variables(); jj++) {
fprintf(OUT," + %s*supp_lin_%d_%d",W.name(jj).c_str(), indicators[ii],jj+1);
}
fprintf(OUT, ";\n\n");
}
fprintf(OUT,"END;\n\n\n\n");
for (unsigned int ii=0; ii<W.num_patches(); ii++) {
std::stringstream linname;
linname << "patch_" << ii;
write_vector_as_constants(W.patch(ii), linname.str(), OUT);
linname.clear(); linname.str("");
}
fclose(OUT);
}
void WitnessSet::merge(const WitnessSet & W_in, tracker_config_t * T)
{
//error checking first
if ( (num_vars_==0) && (W_in.num_vars_!=0) && (num_points()==0)) {
num_vars_ = W_in.num_vars_;
num_natty_vars_ = W_in.num_natty_vars_;
}
if (W_in.num_natty_vars_ != this->num_natty_vars_) {
std::stringstream ss;
ss << "merging two witness sets with differing numbers of natural variables. "<< W_in.num_natural_variables() <<" merging set, "<< this->num_natural_variables() << " existing\n",
throw std::logic_error(ss.str());
}
//just mindlessly add the linears. up to user to ensure linears get merged correctly. no way to know what they want...
for (unsigned int ii = 0; ii<W_in.num_linears(); ii++) {
WitnessSet::add_linear(W_in.linear(ii));
}
for (unsigned int ii = 0; ii<W_in.num_points(); ii++) {
int is_new = 1;
vec_mp & in_point = W_in.point(ii);
for (unsigned int jj = 0; jj<num_points(); jj++){
vec_mp & curr_point = this->point(jj);
//cache the sizes
int in_size = in_point->size;
int curr_size = curr_point->size;
in_point->size = this->num_natural_variables();
curr_point->size = this->num_natural_variables();
if (isSamePoint_homogeneous_input(curr_point, in_point, T->final_tol_times_mult)) {
is_new = 0;
in_point->size = in_size;
curr_point->size = curr_size;
break;
}
// restore the sizes
in_point->size = in_size;
curr_point->size = curr_size;
}
if (is_new==1)
WitnessSet::add_point( (in_point) );
}
for (unsigned int ii = 0; ii<W_in.num_patches(); ii++) {
int is_new = 1;
for (unsigned int jj = 0; jj<this->num_patches(); jj++){
if ( this->patch(jj)->size == W_in.patch(ii)->size) {
if (isSamePoint_inhomogeneous_input(patch(jj), W_in.patch(ii), T->final_tol_times_mult)) {
is_new = 0;
break;
}
}
}
if (is_new==1)
WitnessSet::add_patch(W_in.patch(ii));
}
return;
}//re: merge_witness_sets
