void UbermasterProcess::critreal(WitnessSet & W, vec_mp *pi, VertexSet & V)
{
SystemRandomizer randomizer;
randomizer.setup(W.num_variables()-1-W.dimension(),solve_options.PPD.num_funcs);
print_point_to_screen_matlab(pi[0],"pi");
FILE *OUT;
OUT = safe_fopen_write("most_recent_pi");
print_vec_out_mp(OUT, pi[0]); /* Print vector to file */
fclose(OUT);
NullspaceConfiguration ns_config; // this is set up in the nullspace call.
SolverOutput solve_out;
compute_crit_nullspace_left(solve_out, // the returned value
W, // input the original witness set
std::make_shared<SystemRandomizer>(randomizer),
pi,
W.dimension(), // dimension of ambient complex object
W.dimension(), // target dimension to find
W.dimension(), // COdimension of the critical set to find.
program_options,
solve_options,
&ns_config);
WitnessSet W_crit, W_nonsing, W_sing;
solve_out.get_noninfinite_w_mult_full(W_crit);
solve_out.get_nonsing_finite_multone(W_nonsing);
solve_out.get_sing_finite(W_sing);
W_crit.write_dehomogenized_coordinates("crit_data");
W_nonsing.write_dehomogenized_coordinates("nonsingular");
W_sing.write_dehomogenized_coordinates("singular");
return;
}
int UbermasterProcess::main_loop()
{
boost::timer::auto_cpu_timer t;
program_options.splash_screen();
//parse the options
program_options.startup(); // tests for existence of necessary files, etc.
//if desired, display the options
if (program_options.verbose_level()>=3)
program_options.display_current_options();
parse_preproc_data("preproc_data", &solve_options.PPD);
int num_vars = get_num_vars_PPD(solve_options.PPD);
NumericalIrreducibleDecomposition witness_data;
witness_data.populate(&solve_options.T);
WitnessSet W = witness_data.choose(program_options);
if (W.num_points()==0) {
std::cout << "no witness points, cannot decompose anything..." << std::endl;
return 1;
}
W.get_variable_names(num_vars);
W.set_input_filename(program_options.input_filename());
if (program_options.verbose_level()>=1) {
W.print_to_screen();
}
VertexSet V(num_vars);
V.set_tracker_config(&solve_options.T);
V.set_same_point_tolerance(1e1*solve_options.T.real_threshold);
vec_mp *pi = (vec_mp *) br_malloc(W.dimension()*sizeof(vec_mp ));
for (int ii=0; ii<W.dimension(); ii++) {
init_vec_mp2(pi[ii],W.num_variables(), solve_options.T.AMP_max_prec);
pi[ii]->size = W.num_variables();
}
get_projection(pi, program_options, W.num_variables(), W.dimension());
for (int ii=0; ii<W.dimension(); ii++) {
V.add_projection(pi[ii]);
}
if (program_options.primary_mode()==BERTINIREAL) {
bertini_real(W,pi,V);
}
else if(program_options.primary_mode()==CRIT)
{
critreal(W,pi,V);
}
for (int ii=0; ii<W.dimension(); ii++)
clear_vec_mp(pi[ii]);
free(pi);
// dismiss the workers
int sendme = TERMINATE;
MPI_Bcast(&sendme, 1, MPI_INT, 0, MPI_COMM_WORLD);
return SUCCESSFUL;
}
void UbermasterProcess::bertini_real(WitnessSet & W, vec_mp *pi, VertexSet & V)
{
W.set_incidence_number(get_incidence_number( W.point(0), program_options, program_options.input_filename()));
boost::filesystem::path temp_name = program_options.output_dir();
std::stringstream converter;
converter << "_dim_" << W.dimension() << "_comp_" << W.component_number();
temp_name += converter.str();
program_options.output_dir(temp_name);
switch (W.dimension()) {
case 1:
{
Curve C;
// curve
C.main(V, W, pi, program_options, solve_options);
if (program_options.verbose_level()>=2)
printf("outputting data\n");
C.output_main(program_options.output_dir());
V.print(program_options.output_dir()/ "V.vertex");
}
break;
case 2:
{
Surface S;
// surface
S.main(V, W, pi, program_options, solve_options);
S.output_main(program_options.output_dir());
V.print(program_options.output_dir()/ "V.vertex");
}
break;
default:
{
std::cout << "bertini_real not programmed for components of dimension " << W.dimension() << std::endl;
}
break;
}
}
