int witness_set_test(Workspace& workspace_cpu, CPUInstHom& cpu_inst_hom, Parameter path_parameter, CT* x_cpu, CT t, int device_option) {
std::cout << "--------- Witness Set Test ----------" << std::endl;
bool cpu_test = 0;
bool gpu_test = 0;
if(device_option == 1){
cpu_test = 1;
gpu_test = 0;
}
else if(device_option == 2){
cpu_test = 0;
gpu_test = 1;
}
else{
std::cout << "Device_option Invalid. Choose from the following:" << std::endl
<< " 1. CPU" << std::endl
<< " 2. GPU" << std::endl
<< "Your device_option = " << device_option << std::endl;
}
CT* x_new = NULL;
PolySolSet witness_set_start(cpu_inst_hom.dim);
PolySolSet witness_set_target(cpu_inst_hom.dim);
witness_set_start.add_diff_sol(x_cpu);
std::cout << "Witness Set1" << std::endl;
witness_set_start.print_short();
double timeSec_Predict = 0;
double timeSec_Eval = 0;
double timeSec_MGS = 0;
double timeSec = 0;
for(int i=0; i<10; i++){
clock_t begin = clock();
bool success = 0;
int alpha_round = 0;
while(success==0){
cpu_inst_hom.path_data.clear();
t = CT(0.0,0.0);
workspace_cpu.init_x_t_idx();
workspace_cpu.update_x_t(x_cpu, t);
cpu_inst_hom.update_alpha();
if(cpu_test == 1){
success = path_tracker(workspace_cpu, cpu_inst_hom, path_parameter,\
timeSec_Predict, timeSec_Eval, timeSec_MGS);
}
if(gpu_test == 1){
success = GPU_Path(cpu_inst_hom, path_parameter, x_cpu, t, x_new);
}
alpha_round++;
if(success == 0 && alpha_round > 5){
break;
}
}
if(success == 0 && alpha_round > 5){
std::cout << "Path Fail! Start -> Target" << std::endl \
<< "alpha round = " << alpha_round << std::endl \
<< "path round = " << i << std::endl;
break;
}
cpu_inst_hom.path_data.print_phc();
clock_t end = clock();
timeSec += (end - begin) / static_cast<double>( CLOCKS_PER_SEC );
if(cpu_test == 1){
workspace_cpu.copy_x_last(x_cpu);
}
if(gpu_test == 1){
delete[] x_cpu;
x_cpu = x_new;
}
witness_set_target.add_diff_sol(x_cpu);
std::cout << "Witness Set Target" << std::endl;
witness_set_target.print_short();
if(witness_set_target.n_sol == 4 and witness_set_start.n_sol == 4){
break;
}
// Path Tracking Reverse Test
success = 0;
alpha_round = 0;
while(success==0){
cpu_inst_hom.path_data.clear();
t = CT(0.0,0.0);
workspace_cpu.init_x_t_idx();
workspace_cpu.update_x_t(x_cpu, t);
cpu_inst_hom.update_alpha();
//success = path_tracker(workspace_cpu, cpu_inst_hom, path_parameter,\
// timeSec_Predict, timeSec_Eval, timeSec_MGS, 1);
if(cpu_test == 1){
success = path_tracker(workspace_cpu, cpu_inst_hom, path_parameter,\
timeSec_Predict, timeSec_Eval, timeSec_MGS, 1);
}
if(gpu_test == 1){
success = GPU_Path(cpu_inst_hom, path_parameter, x_cpu, t, x_new, 1);
}
alpha_round++;
if(success == 0 && alpha_round > 5){
break;
}
}
cpu_inst_hom.path_data.print_phc();
if(success == 0 && alpha_round > 5){
std::cout << "Path Fail! Target -> Start" << std::endl \
<< "alpha round = " << alpha_round << std::endl \
<< "path round = " << i << std::endl;
break;
}
if(cpu_test == 1){
workspace_cpu.copy_x_last(x_cpu);
}
if(gpu_test == 1){
delete[] x_cpu;
x_cpu = x_new;
}
witness_set_start.add_diff_sol(x_cpu);
std::cout << "Witness Set Start" << std::endl;
witness_set_start.print_short();
if(witness_set_target.n_sol == 4 and witness_set_start.n_sol == 4){
break;
}
}
std::cout << "Witness Set Start" << std::endl;
witness_set_start.print_short();
std::cout << "Witness Set Target" << std::endl;
witness_set_target.print_short();
/*for(int i=0; i<cpu_inst_hom.dim; i++) {
std::cout << i << " " << x_cpu[i];
}*/
/*CT* x_gpu = GPU_Path(cpu_inst_hom, sol0, t, n_predictor, max_it, err_max_delta_x, max_step);
cout << "Path CPU Test Time: "<< timeSec << endl;
std::cout << "--------- Path Tracking Error CPU vs GPU ----------" << std::endl;
T1 err = err_check_workspace(x_cpu, x_gpu, cpu_inst_hom.dim);
std::cout << " x_cpu[0] = " << x_cpu[0];
std::cout << " x_gpu[0] = " << x_gpu[0];
free(x_gpu);*/
//T1 err = err_check_workspace(sol0, x_cpu, cpu_inst_hom.dim);
cout << "Path CPU Predict Time: "<< timeSec_Predict << endl;
cout << "Path CPU Eval Time: "<< timeSec_Eval << endl;
cout << "Path CPU MGS Time: "<< timeSec_MGS << endl;
return witness_set_start.n_sol;
}
