void writer::Wstats::print_all_statistics(int stage, const Statistics& info, const ProblemInstance& cfg, const MBGraph& graph) {
if (stage == 0) {
println("Initial graph:");
} else {
println("After Stage " + toString(stage) + " graph:");
}
print_complete_edges(graph);
print_connected_components(graph);
print_rear_characters(info.get_compl_stat(graph));
#ifndef VERSION2
print_estimated_dist(stage, cfg, graph);
#endif
print_fair_edges(graph);
print_not_compl_characters(info.get_no_compl_stat(graph));
}
int main(int argc, char const *argv[])
{
char* s;
std::srand(std::time(0)); //use current time as seed for random generator
int r = rand() % 1000;
for(int i = 0; i < r; i++)
{
rand();
}
if(argc < 3)
{
int forestSize = strtol(argv[1], &s, 10);
for(int i = 0 ; i < forestSize ; i++)
{
printf("%lf\n",fRand(1,std::sqrt(10)));
}
return 1;
}
int forestSize = strtol(argv[1], &s, 10);
int iterations = strtol(argv[2], &s, 10);
double SIDE = std::sqrt(forestSize);
SIDE = fRand(std::sqrt(SIDE),std::sqrt(2)*SIDE);
double R = 1;
double begin, end;
std::vector<int> empty;
std::vector<Tree*> Forest;
std::vector< std::vector<int> > neighbors(forestSize,empty);
std::vector< std::vector<double> > metrics(iterations,std::vector<double>(forestSize,0.0));
//Parallel variables
int num_threads;
std::vector<int> order;
begin = omp_get_wtime();
#pragma omp parallel shared(Forest,neighbors,metrics,forestSize,iterations,order)
{
#pragma omp master
{
// INIT VARIABLES
num_threads = omp_get_num_threads();
std::vector<Point> positions;
std::cout << "Running " << forestSize << " trees for " << iterations << " iterations on " << num_threads << " processors" << std::endl;
printf("SIDE = %lf, R = %lf\n",SIDE,R);
for(int i = 0; i < forestSize; i++)
{
// double x = std::fabs((SIDE-1)*std::sin(i));
// double y = std::fabs(SIDE*std::cos(i*i));
double x = fRand(0,SIDE);
double y = fRand(0,SIDE);
Point p = {x,y};
Tree *T = new MonopodialTree();
Forest.push_back(T);
positions.push_back(p);
for(int j = 0 ; j < i ; j++)
{
Point q = positions[j];
if(pointDistance(p,q) < R)
{
neighbors[j].push_back(i);
neighbors[i].push_back(j);
}
}
}
order = get_order(neighbors);
for(int i = 0; i < order.size(); i++)
std::cout << order[i] << " ";
std::cout << std::endl;
}
#pragma omp barrier
int thread_num = omp_get_thread_num();
// ITERATE
int N = forestSize;
int T = iterations;
int P = omp_get_num_threads();
int x = thread_num;
int y = 0;
while( x+N*y < N*T)
{
int i = order[x];
// printf("%d (%d, %d)\n",thread_num,y,i );
while(Forest[i]->iteration < y);
bool ready = false;
while(!ready)
{
ready = true;
for(int k = 0; k < neighbors[i].size() ; k++)
{
if( Forest[ neighbors[i][k] ]->iteration < y)
{
ready = false;
break;
}
}
}
if(y > 0)
{
Forest[i]->updateMetric(metrics[y],neighbors[i]);
}
Forest[i]->next();
double metric = Forest[i]->calculateMetric();
#pragma omp critical(metrics)
{
metrics[y][i] = metric;
}
x+=P;
if(x >= N)
{
x -= N;
y++;
}
}
}
end = omp_get_wtime();
print_forest(Forest, neighbors, metrics[iterations-1]);
std::vector< std::vector<int> > connected_components = get_connected_components(neighbors);
print_connected_components( connected_components);
char buffer[80];
FILE *f = fopen("Results_lookahead.txt", "a");
if(f != NULL)
{
fprintf(f, "%s\n", gettime(buffer));
fprintf(f,"%d threads\n",num_threads);
fprintf(f,"%d trees\n",forestSize);
fprintf(f,"%d iterations\n",iterations);
for(int i = 0; i < connected_components.size(); i++)
{
fprintf(f, "%d ", connected_components[i].size());
}
fprintf(f, "\n");
fprintf(f,"Time : %f seconds\n", end-begin);
fprintf(f,"\n=====================\n");
}
for(int i = 0; i < Forest.size() ; i++)
{
delete Forest[i];
}
return 0;
}