SpherePacking/main.cpp at master · t-fi/SpherePacking · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
#include "simulator.h"
#include "renderer.h"
#include <iostream>
#include <limits>
#include <iomanip>
#include <cmath>

#define PI 3.14159265359

int main(int argc, char * argv[])
{
	// argument parsing
	int numPoints, seed;
	long int steps;
	double lambda, sigma;

	// reading arguments
	for(int i=0; i<argc; i++){
		if(i==1) numPoints = atoi(argv[i]);
		if(i==2) steps     = atol(argv[i]);
		if(i==3) lambda    = atof(argv[i]);
		if(i==4) sigma     = atof(argv[i]);
		if(i==5) seed      = atof(argv[i]);
	}

	// setting defaults for unspecified args
	for(int i=5; i>=argc; i--){
		if(i==1) numPoints = 12;
		if(i==2) steps     = 1000000;
		if(i==3) lambda    = 1.;
		if(i==4) sigma     = 1.;
		if(i==5) seed      = 0;
	}

	// error handling
	if(numPoints <= 0  ||
	   steps     <= 0  ||
     lambda    <= 0. ||
	   sigma     <= 0.)
	{
		std::cout << "An argument was invalid! Program terminated." << std::endl;
		return 1;
	}

	// printing arguments
	std::cout << '#' << " numPoints = " << numPoints << std::endl;
	std::cout << '#' << " steps     = " << steps     << std::endl;
	std::cout << '#' << " lambda    = " << lambda    << std::endl;
	std::cout << '#' << " sigma     = " << sigma     << std::endl;
	std::cout << '#' << " seed      = " << seed      << std::endl;

	// starting simulator
	Simulator simulator(numPoints, lambda, sigma, seed);

	//succes rate (SR) counters
	double percSRradius = 0.0;
	double percSRposition = 0.0;
	long int SRradius = 0;
	long int SRposition = 0;
	long int trials = 0;

	// main loop
	for(long int i=0; i<steps; i++){

		// log every 1000th step
		if(i%(steps/1000)==0){
			percSRradius = (double)SRradius/(double)(trials*numPoints);
			percSRposition = (double)SRposition/(double)trials;
			trials = 0;
			SRradius = 0;
			SRposition = 0;
			std::cerr << simulator.elapsedTime() << " " << i << " " << std::fixed << std::setprecision(19) << 2*simulator.radius << " " << std::fixed << std::setprecision(10) << percSRradius << " " << percSRposition << std::endl;
			//simulator.saveFiles(i);
		}

		trials++;
		for(long int j = 0; j<numPoints;j++){
			SRradius += simulator.increaseRadius();
		}
		SRposition += simulator.movePoint();
	}

	simulator.saveFiles(0);

	return 0;
}