dynamic_generation/main.cpp at master · varunk122/dynamic_generation · 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
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
#include<generator.hpp>

std::vector<Particle> particles; //this list will store all the particles


int main(int argc, char *argv[]) {

    if (argc !=9) {
        std::cout << "wrong number of input arguments\n";
        std::cout << "total number of input arguments provided " << argc - 1 << "\n";
        return 0;
    }

    RandomGenerator generator;

    generator.radius_ = std::atof (argv[1]);
    generator.insertion_pf_ = std::atof (argv[2]);

    int x0 = std::atoi(argv[3]);
    int y0 = std::atoi(argv[4]);
    int z0 = std::atoi(argv[5]);

    int x1 = std::atoi(argv[6]);
    int y1 = std::atoi(argv[7]);
    int z1 = std::atoi(argv[8]);


    Domain simDomain (Point3d(x0,y0,z0),Point3d(x1,y1,z1));
    generator.randomParticleGenerator(particles, simDomain);
    std::cout <<"Packing fraction after inserting a lot of particles: "<< generator.calculatePackingFraction(particles,simDomain) << std::endl;
    // deleteParticles(particles);
    generator.addParticlesToGrid(particles, simDomain);
    // printGrid();

    generator.newDeleteOverlappingParticles(100);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 100 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(60);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 60 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(50);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 50 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(40);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 40 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.newDeleteOverlappingParticles(30);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 30 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(25);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 25 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(20);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 20 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(15);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 15 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(13);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 13 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(12);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 12 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(10);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 10 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(8);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 8 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(6);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 6 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(5);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 5 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(4);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 4 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);

    // generator.deleteInGrid(1e-4);


    generator.newDeleteOverlappingParticles(3);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 3 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);
    generator.newDeleteOverlappingParticles(2);
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction 2 overlap: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToList(particles);
    // std::cout << "Packing fraction after deleting overlap particles: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;

    generator.overlapAnalysis(particles);

    generator.addParticlesToGrid(particles, simDomain);

    generator.newDeleteOverlappingParticles(1);
    generator.addParticlesToList(particles);
    std::cout << "Final Packing Fraction: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.addParticlesToGrid(particles, simDomain);

    std::cout << "Total cells: " << generator.totalCellsInGrid() << std::endl;
    std::cout << "Total number of empty cells: " << generator.emptyCellsCount() << std::endl;
    // generator.deleteInGrid(1e-4);
    // std::cout << "Total number of empty cells: " << generator.emptyCellsCount() << std::endl;

    // generator.addParticlesInEmptyCells();
    // std::cout << "Total number of empty cells: " << generator.emptyCellsCount() << std::endl;
    // generator.deleteInGrid();
    // std::cout << "Total number of empty cells: " << generator.emptyCellsCount() << std::endl;
    generator.addParticlesToList(particles);
    std::cout << "Packing fraction after deleting overlap particles: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    generator.randomDeleting(particles, simDomain, 0.6);
    std::cout << "Packing fraction after deleting particles randomly to achieve desired packing fraction: " << generator.calculatePackingFraction(particles,simDomain) << std::endl;
    // printParticlesList(particles);
    generator.saveToCSV(particles);

    return 0;
}