void twoOpt() {
if (done_twoOpt) return;
cur_energy = cur.getDist();
int start, finish, temp;
finish = -1;
start = 0;
while (finish < start) {
start = cur.getDist();
for (int r = 0; r < NUM_COLORS - 1; r++) {
for (int c = r + 1; c < NUM_COLORS; c++) {
if (!cur.swapIsBeneficial(r, c)) continue;
cur.twoOptSwap(r, c);
cur_energy = cur.getUpdatedDistWhenTwoOptSwap(cur_energy, r, c);
if (ANIMATE) {
cur.draw();
update();
}
}
}
finish = cur.getDist();
if (finish < start)
best_tour.updateColors(cur.colors);
}
cout << "Final distance after 2-opt: " << best_tour.getDist() << endl;
done_twoOpt = true;
}
Tour * TourGA::crossover (Tour * parent1, Tour * parent2)
{
Tour * child = new Tour();
int startPos = (int)rand() % parent1->tourLength();
int endPos = (int)rand() % parent1->tourLength();
for (int index = 0; index < child->tourLength(); index++)
{
if ((startPos < endPos) && (index > startPos) && (index < endPos))
{
child->setCity(index, parent1->getCity(index));
}
else if (startPos > endPos)
{
if (!((index < startPos) && (index > endPos)))
{
child->setCity (index, parent1->getCity(index));
}
}
}
for (int index = 0; index < parent2->tourLength(); index++)
{
if (!child->containsCity(parent2->getCity(index)))
{
for (int innerLoop = 0; innerLoop < child->tourLength(); innerLoop++)
{
if (child->getCity(innerLoop) == NULL)
{
child->setCity(innerLoop, parent2->getCity(index));
break;
}
}
}
}
return child;
}
void MainWindow::drawPath(QItemSelection item){
QCPCurve *curve = qobject_cast<QCPCurve*>(ui->citiesPlot->plottable(0));
curve->clearData();
curve->setLineStyle(QCPCurve::lsLine);
ui->citiesPlot->replot();
QModelIndex index = item.indexes().first();
Tour tour = qvariant_cast<Tour>(ui->listView->model()->data(index, Qt::UserRole));
for(int i=0; i<tour.tourSize(); i++){
int x = tour.getCity(i).getX();
int y = tour.getCity(i).getY();
curve->addData(x, y);
}
curve->addData(tour.getCity(0).getX(), tour.getCity(0).getY());
ui->citiesPlot->replot();
}
int main(int argc, char *argv[]) {
QApplication a(argc, argv);
// define 4 points forming a square
Point p(100.0, 100.0);
Point q(500.0, 100.0);
Point r(500.0, 500.0);
Point s(100.0, 500.0);
// Set up a Tour with those four points
// The constructor should link p->q->r->s->p
Tour squareTour(p, q, r, s);
squareTour.show();
string filename = "tsp10.txt";
ifstream input;
input.open(filename);
// get dimensions
int width;
int height;
input >> width;
input >> height;
// setup graphical window
QGraphicsView *view = new QGraphicsView();
QGraphicsScene *scene = new QGraphicsScene();
view->setScene(scene);
view->scale(1, -1); //screen y-axis is inverted
view->setSceneRect(0, 0, width, height);
view->show();
// run insertion heuristic
Tour tour;
double x;
double y;
while (input >> x >> y) {
Point p(x, y);
//tour.insertNearest(p);
tour.insertSmallest(p);
//uncomment the 4 lines below to animate
tour.draw(scene);
std::chrono::milliseconds dura(50);
std::this_thread::sleep_for(dura);
a.processEvents();
}
input.close();
// print tour to standard output
cout << "Tour distance: " << std::fixed << std::setprecision(4)
<< std::showpoint << tour.distance() << endl;
cout << "Number of points: " << tour.size() << endl;
tour.show();
// draw tour
tour.draw(scene);
return a.exec(); // start Qt event loop
}
int main(int argc, char * argv[])
{
srand ( time(NULL) );
CSDLManagerLite::getInstance()->initializeSDL(WIDTH, HEIGHT, TITLE);
cout << "Initial distance: " << cur.getDist() << endl;
cout << "Colors: " << NUM_COLORS << endl;
int a, b;
double cur_energy, new_energy;
while (!CSDLInputManagerLite::getInstance() -> isExit())
{
CSDLManagerLite::getInstance()->delay(1000);
twoOpt();
best_tour.draw();
update();
}
CSDLManagerLite::getInstance()->clean();
}
/* Do a depth-first search to find the best tour
* starting with the partially complete tour t
*/
void DFS(Tour& t) {
if (!t.better_than(best)) return; // if we've already seen something better
// If we've added all the cities and we're still here, this must be
// the best tour
if (t.get_num_cities() == num_cities) {
best = t;
#ifdef DEBUG
cout << "New best:" << best;
#endif
}
// Otherwise, try to add the next city and recursively call DFS
// if successful
else {
for (unsigned int i=0; i < num_cities; i++) {
if (t.add_city(i, best)) {
DFS(t);
t.remove_last();
}
}
}
}
//Handle controls
void key_callback(int key, int action){
if(!tour) {
if ((key == 'W' || key == 'w') && action == GLFW_PRESS){
viewer->setForwardVelocity(1);
} else if ((key == 'W' || key == 'w') && action == GLFW_RELEASE){
viewer->setForwardVelocity(0);
} else if ((key == 'S' || key == 's') && action == GLFW_PRESS){
viewer->setForwardVelocity(-1);
} else if ((key == 'S' || key == 's') && action == GLFW_RELEASE){
viewer->setForwardVelocity(0);
} else if ((key == 'A' || key == 'a') && action == GLFW_PRESS){
viewer->setStrafeVelocity(-1);
} else if ((key == 'A' || key == 'a') && action == GLFW_RELEASE){
viewer->setStrafeVelocity(0);
} else if ((key == 'D' || key == 'd') && action == GLFW_PRESS){
viewer->setStrafeVelocity(1);
} else if ((key == 'D' || key == 'd') && action == GLFW_RELEASE){
viewer->setStrafeVelocity(0);
} else if ((key == GLFW_KEY_UP) && action == GLFW_PRESS){
viewer->changeVelocity(0,0,0.2);
} else if ((key == GLFW_KEY_DOWN) && action == GLFW_PRESS){
viewer->changeVelocity(0,0,-0.2);
} else if ((key == GLFW_KEY_PAGEUP) && action == GLFW_PRESS){
viewer->setElevationVelocity(1);
} else if ((key == GLFW_KEY_PAGEUP) && action == GLFW_RELEASE){
viewer->setElevationVelocity(0);
} else if ((key == GLFW_KEY_PAGEDOWN) && action == GLFW_PRESS){
viewer->setElevationVelocity(-1);
} else if ((key == GLFW_KEY_PAGEDOWN) && action == GLFW_RELEASE){
viewer->setElevationVelocity(0);
} else if ((key == GLFW_KEY_LEFT) && action == GLFW_PRESS){
viewer->setCameraRotationVelocity(2);
} else if ((key == GLFW_KEY_RIGHT) && action == GLFW_PRESS){
viewer->setCameraRotationVelocity(-2);
} else if ((key == GLFW_KEY_LEFT) && action == GLFW_RELEASE){
viewer->setCameraRotationVelocity(0);
} else if ((key == GLFW_KEY_RIGHT) && action == GLFW_RELEASE){
viewer->setCameraRotationVelocity(0);
} else if ((key == GLFW_KEY_SPACE) && action == GLFW_PRESS){
viewer->setCameraRotationVelocity(0);
viewer->setVelocity(0,0,0);
} else if ((key == 'P' || key == 'p') && action == GLFW_PRESS){
viewer->setCameraRotationVelocity(0);
viewer->setVelocity(0,0,0);
viewer->resetElevation();
viewer->gotoLocation(glm::vec3(13.8057, 0.3, 0.657926), glm::vec3(-6.42059,0.25837,9.794014) - glm::vec3(13.8057, 0.3, 0.657926) );
} else if ((key == 'T' || key == 't') && action == GLFW_PRESS){
tour = 1;
//Stop any motion
viewer->setCameraRotationVelocity(0);
viewer->setVelocity(0,0,0);
t.restart();
}
} else {
if((key == 'E' || key == 'e') && action == GLFW_PRESS){
tour = 0;
}
}
if ((key == 'H' || key == 'h') && action == GLFW_PRESS){
show_help = !show_help;
}
}
int main(int argc, char *argv[]){
//Initialise GLFW
if (!glfwInit()) {
fprintf(stderr, "There was a problem initializing glfw");
exit(EXIT_FAILURE);
}
//glfw hints
glfwOpenWindowHint(GLFW_FSAA_SAMPLES, 4);
glfwOpenWindowHint(GLFW_OPENGL_VERSION_MAJOR, 3);
glfwOpenWindowHint(GLFW_OPENGL_VERSION_MINOR, 3);
glfwOpenWindowHint(GLFW_OPENGL_PROFILE, GLFW_OPENGL_CORE_PROFILE);
//Open a window
if (!glfwOpenWindow(WINDOW_WIDTH, WINDOW_HEIGHT, 0, 0, 0, 0, 0, 0, GLFW_WINDOW)) {
fprintf(stderr, "There was a problem opening a window");
exit(EXIT_FAILURE);
}
//Turn on experimental features to prevent seg fault
glewExperimental = GL_TRUE;
//Intitialize GLEW
if (glewInit() != GLEW_OK) {
fprintf(stderr, "There was a problem initialising GLEW");
exit(EXIT_FAILURE);
}
//various glfw settings
glfwEnable(GLFW_STICKY_KEYS);
glfwSetKeyCallback(&key_callback);
glfwSetWindowTitle(TITLE);
glEnable(GL_CULL_FACE);
glCullFace(GL_BACK);
glEnable(GL_DEPTH_TEST);
glEnable(GL_DEPTH_CLAMP);
//Create text generator
GLuint textShaderID = setupShaders("shaders/text/vert.gls", "shaders/text/frag.gls");
TextGenerator tg((char*)"textures/font.png", textShaderID, ' ', '~', 16, 8, WINDOW_WIDTH, WINDOW_HEIGHT);
loading(&tg);
//Load in and set program (Shaders)
// GLuint simpleShaderID = setupShaders("shaders/simple/vert.gls", "shaders/simple/frag.gls");
// GLuint phongShaderID = setupShaders("shaders/phong/vert.gls", "shaders/phong/frag.gls");
GLuint perFragmentShaderID = setupShaders("shaders/perFragment/vert.gls", "shaders/perFragment/frag.gls");
//Load in terrain and sky box
HeightMapLoader groundObj = HeightMapLoader("textures", "img2.png");
ObjLoader skyboxObj = ObjLoader("models","skybox.obj");
Object ground(&groundObj, perFragmentShaderID, viewer, GL_FILL);
Skybox skybox(&skyboxObj, perFragmentShaderID, viewer, GL_FILL);
loading(&tg);
//Load in thunder birds
ObjLoader thunderBird1Obj = ObjLoader("models", "thunderbird1.obj");
ObjLoader t2ContainerObj = ObjLoader("models","t2container.obj");
ObjLoader thunderBird2Obj = ObjLoader("models","thunderbird2.obj");
ObjLoader thunderBird3Obj = ObjLoader("models","thunderbird3.obj");
ObjLoader eggObj = ObjLoader("models", "egg.obj");
Object thunderBird1(&thunderBird1Obj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird1B(&thunderBird1Obj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird1C(&thunderBird1Obj, perFragmentShaderID, viewer, GL_FILL);
Object t2container(&t2ContainerObj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird2(&thunderBird2Obj, perFragmentShaderID, viewer, GL_FILL);
Object t2containerB(&t2ContainerObj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird2B(&thunderBird2Obj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird3(&thunderBird3Obj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird3B(&thunderBird3Obj, perFragmentShaderID, viewer, GL_FILL);
Object thunderBird3C(&thunderBird3Obj, perFragmentShaderID, viewer, GL_FILL);
loading(&tg);
//Set up collisions
viewer->addTerrain(&ground);
loading(&tg);
thunderBird1.setPosition(glm::vec3(-2,0.1,0));
t2container.setPosition(glm::vec3(-7.5,0.15,5));
t2container.setScale(glm::vec3(0.2,0.2,0.2));
t2container.setRotation(90, glm::vec3(0,1,0));
thunderBird2.setPosition(glm::vec3(-7.5,0.15,5));
thunderBird2.setScale(glm::vec3(0.2,0.2,0.2));
thunderBird2.setRotation(90, glm::vec3(0,1,0));
t2containerB.setPosition(glm::vec3(-13.7941, 0.0, -11.716));
t2containerB.setScale(glm::vec3(0.2,0.2,0.2));
t2containerB.setRotation(120, glm::vec3(0,1,0));
thunderBird2B.setPosition(glm::vec3(-13.7941, 0.0, -11.716));
thunderBird2B.setScale(glm::vec3(0.2,0.2,0.2));
thunderBird2B.setRotation(120, glm::vec3(0,1,0));
thunderBird3.setPosition(glm::vec3(2,2,12.5));
thunderBird3.setScale(glm::vec3(0.4,0.4,0.4));
thunderBird3.setRotation(-90, glm::vec3(1,0,0));
thunderBird3B.setPosition(glm::vec3(22.3989, 0.541667, -5.00961));
thunderBird3B.setScale(glm::vec3(0.4,0.4,0.4));
thunderBird3B.setRotation(-90, glm::vec3(1,0,0));
thunderBird3B.setPosition(glm::vec3(2,2,12.5));
thunderBird3B.setScale(glm::vec3(0.4,0.4,0.4));
thunderBird3B.setRotation(-90, glm::vec3(1,0,0));
ground.setScale(glm::vec3(0.35,0.025,0.35));
ground.setPosition(glm::vec3(-groundObj.width/2 * 0.35,0,-groundObj.height/2 * 0.35));
std::vector<glm::vec3> positions;
positions.push_back(glm::vec3(-6.42059,0.25837,4.494014));
positions.push_back(glm::vec3(-3.85007,0.1,7.12377));
positions.push_back(glm::vec3(-0.601204,0.1,3.592329));
positions.push_back(glm::vec3(8.2538,0.1,2.5428));
positions.push_back(glm::vec3(10.6809,0.38337,-3.84069));
positions.push_back(glm::vec3(-5.13897,0.1,-11.4451));
positions.push_back(glm::vec3(-14.3374,0.1,-8.11471));
positions.push_back(glm::vec3(-20.2777, 0.1, 8.26625));
positions.push_back(glm::vec3(-22.3929, 1.35833,19.6362));
positions.push_back(glm::vec3(1.64004, 0.1, 33.6688));
positions.push_back(glm::vec3(30.9513, 0.683333, 18.214));
positions.push_back(glm::vec3(32.9145, 1.56667, -11.0445));
positions.push_back(glm::vec3(31.17666, 0.808333, -22.8479));
positions.push_back(glm::vec3(-29.7966, 0.35, -27.9126));
positions.push_back(glm::vec3(-33.6018, 0.433333, -25.1177));
positions.push_back(glm::vec3(-29.1175, 2.39167, -17.2737));
positions.push_back(glm::vec3(-27.6434, 2.1833, -14.5711));
positions.push_back(glm::vec3(-33.0654, 3.5, 2.62864));
positions.push_back(glm::vec3(-20.6909, 4, 3.49496));
positions.push_back(glm::vec3(-3.85007, 3.5 , 7.12377));
positions.push_back(glm::vec3(1.42417, 2.116667, 0.276567));
WayPoint w1 = WayPoint(positions[0], positions[1] - positions[0]);
WayPoint w2 = WayPoint(positions[1], positions[2] - positions[1]);
WayPoint w3 = WayPoint(positions[2], positions[3] - positions[2]);
WayPoint w4 = WayPoint(positions[3], positions[4] - positions[3]);
WayPoint w5 = WayPoint(positions[4], positions[5] - positions[4]);
WayPoint w6 = WayPoint(positions[5], positions[6] - positions[5]);
WayPoint w7 = WayPoint(positions[6], positions[7] - positions[6]);
WayPoint w8 = WayPoint(positions[7], positions[8] - positions[7]);
WayPoint w9 = WayPoint(positions[8], positions[9] - positions[8]);
WayPoint w10 = WayPoint(positions[9], positions[10] - positions[9]);
WayPoint w11 = WayPoint(positions[10], positions[11] - positions[10]);
WayPoint w12 = WayPoint(positions[11], positions[12] - positions[11]);
WayPoint w13 = WayPoint(positions[12], positions[13] - positions[12]);
WayPoint w14 = WayPoint(positions[13], positions[14] - positions[13]);
WayPoint w15 = WayPoint(positions[14], positions[15] - positions[14]);
WayPoint w16 = WayPoint(positions[15], positions[16] - positions[15]);
WayPoint w17 = WayPoint(positions[16], positions[17] - positions[16]);
WayPoint w18 = WayPoint(positions[17], positions[18] - positions[17]);
WayPoint w19 = WayPoint(positions[18], positions[19] - positions[18]);
WayPoint w20 = WayPoint(positions[19], positions[20] - positions[19]);
WayPoint w21 = WayPoint(positions[20], glm::vec3(0,0,1));
t.addWayPoint(0.0, &w1, 1);
t.addWayPoint(2.0, &w2, 1);
t.addWayPoint(4.0, &w3, 1);
t.addWayPoint(6.0, &w4, 1);
t.addWayPoint(8.0, &w5, 1);
t.addWayPoint(10.0, &w6, 1);
t.addWayPoint(12.0, &w7, 1);
t.addWayPoint(14.0, &w8, 1);
t.addWayPoint(16.0, &w9, 1);
t.addWayPoint(18.0, &w10, 1);
t.addWayPoint(20.0, &w11, 1);
t.addWayPoint(22.0, &w12, 1);
t.addWayPoint(24.0, &w13, 1);
t.addWayPoint(26.0, &w14, 1);
t.addWayPoint(28.0, &w15, 1);
t.addWayPoint(30.0, &w16, 1);
t.addWayPoint(32.0, &w17, 1);
t.addWayPoint(34.0, &w18, 1);
t.addWayPoint(36.0, &w19, 1);
t.addWayPoint(38.0, &w20, 1);
t.addWayPoint(40.0, &w21, 1);
Animutator* tb1Anim = new Animutator();
KeyFrame tb1KF1 = KeyFrame(glm::vec3(0,1,1), 90, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1KF2 = KeyFrame(glm::vec3(10,1,1), 90, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1KF3 = KeyFrame(glm::vec3(10,1,1), -90, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1KF4 = KeyFrame(glm::vec3(0,1,1), -90, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
tb1Anim->addKeyFrame(0.0, &tb1KF1);
tb1Anim->addKeyFrame(10.0, &tb1KF2);
tb1Anim->addKeyFrame(12.5, &tb1KF3);
tb1Anim->addKeyFrame(22.5, &tb1KF4);
tb1Anim->addKeyFrame(25.0, &tb1KF1);
thunderBird1.addAnimutator(tb1Anim);
Animutator* tb1BAnim = new Animutator();
KeyFrame tb1BKF1 = KeyFrame(glm::vec3(-23.027, 1.35, 16.35409), -110.328, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF2 = KeyFrame(glm::vec3(-24.5907, 2.29167, 15.7748), -110.328, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF3 = KeyFrame(glm::vec3(-24.5907, 2.29167, 15.7748), -38.604, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF4 = KeyFrame(glm::vec3(-32.497, 1, 25.6773), -38.604, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF5 = KeyFrame(glm::vec3(-32.497, 1, 25.6773), 5.98309, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF6 = KeyFrame(glm::vec3(-31.8602, 1.55833, 31.7533), 5.98309, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF7 = KeyFrame(glm::vec3(-31.8602, 1.55833, 31.7533), 38.384, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF8 = KeyFrame(glm::vec3(-28.5349, 1, 35.9512), 38.384, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF9 = KeyFrame(glm::vec3(-28.5349, 1, 35.9512), 99.0798, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF10 = KeyFrame(glm::vec3(-13.2714, 1, 33.5119), 99.0798, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF11 = KeyFrame(glm::vec3(-13.2714, 1, 33.5119), -150.378, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1BKF12 = KeyFrame(glm::vec3(-23.027, 1.35, 16.35409), -150.378, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
tb1BAnim->addKeyFrame(0.0, &tb1BKF1);
tb1BAnim->addKeyFrame(5.0, &tb1BKF2);
tb1BAnim->addKeyFrame(5.5, &tb1BKF3);
tb1BAnim->addKeyFrame(10.5, &tb1BKF4);
tb1BAnim->addKeyFrame(11.0, &tb1BKF5);
tb1BAnim->addKeyFrame(16.0, &tb1BKF6);
tb1BAnim->addKeyFrame(16.5, &tb1BKF7);
tb1BAnim->addKeyFrame(21.5, &tb1BKF8);
tb1BAnim->addKeyFrame(22.0, &tb1BKF9);
tb1BAnim->addKeyFrame(27.0, &tb1BKF10);
tb1BAnim->addKeyFrame(27.5, &tb1BKF11);
tb1BAnim->addKeyFrame(32.5, &tb1BKF12);
tb1BAnim->addKeyFrame(33.0, &tb1BKF1);
thunderBird1B.addAnimutator(tb1BAnim);
Animutator* tb1CAnim = new Animutator();
KeyFrame tb1CKF1 = KeyFrame(glm::vec3(25.7527, 1.5, -29.4807), -77.7031, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1CKF2 = KeyFrame(glm::vec3(-24.549, 2.76667, -18.516), -77.7031, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1CKF3 = KeyFrame(glm::vec3(-24.549, 2.76667, -18.516), 102.297, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
KeyFrame tb1CKF4 = KeyFrame(glm::vec3(25.7527, 1.5, -29.4807), 102.297, glm::vec3(0,1,0), glm::vec3(0.5,0.5,0.5));
tb1CAnim->addKeyFrame(0.0, &tb1CKF1);
tb1CAnim->addKeyFrame(7.5, &tb1CKF2);
tb1CAnim->addKeyFrame(9, &tb1CKF3);
tb1CAnim->addKeyFrame(16.5, &tb1CKF4);
tb1CAnim->addKeyFrame(18, &tb1CKF1);
thunderBird1C.addAnimutator(tb1CAnim);
Animutator* tb2cAnim = new Animutator();
KeyFrame tbc2KF1 = KeyFrame(glm::vec3(-7.5, 0.15, 5), 90, glm::vec3(0,1,0), glm::vec3(0.2,0.2,0.2));
KeyFrame tbc2KF2 = KeyFrame(glm::vec3(-7.5, -0.3, 5), 90, glm::vec3(0,1,0), glm::vec3(0.2,0.2,0.2));
tb2cAnim->addKeyFrame(0.0, &tbc2KF1);
tb2cAnim->addKeyFrame(5.0, &tbc2KF2);
tb2cAnim->addKeyFrame(5.5, &tbc2KF2);
tb2cAnim->addKeyFrame(10.0, &tbc2KF1);
tb2cAnim->addKeyFrame(10.5, &tbc2KF1);
t2container.addAnimutator(tb2cAnim);
Animutator* tb2cBAnim = new Animutator();
KeyFrame tbc2BKF1 = KeyFrame(glm::vec3(-13.7941, 0.0, -11.716), 120, glm::vec3(0,1,0), glm::vec3(0.2,0.2,0.2));
KeyFrame tbc2BKF2 = KeyFrame(glm::vec3(-13.7941, -0.4, -11.716), 120, glm::vec3(0,1,0), glm::vec3(0.2,0.2,0.2));
tb2cBAnim->addKeyFrame(0.0, &tbc2BKF1);
tb2cBAnim->addKeyFrame(5.0, &tbc2BKF2);
tb2cBAnim->addKeyFrame(5.5, &tbc2BKF2);
tb2cBAnim->addKeyFrame(10.0, &tbc2BKF1);
tb2cBAnim->addKeyFrame(10.5, &tbc2BKF1);
t2containerB.addAnimutator(tb2cBAnim);
Animutator* tb3cAnim = new Animutator();
KeyFrame tb3KF1 = KeyFrame(glm::vec3(2, 1.90, 12.5), -90, glm::vec3(1,0,0), glm::vec3(0.4,0.4,0.4));
KeyFrame tb3KF2 = KeyFrame(glm::vec3(2, 4, 12.5), -90, glm::vec3(1,0,0), glm::vec3(0.4,0.4,0.4));
tb3cAnim->addKeyFrame(0.0, &tb3KF1);
tb3cAnim->addKeyFrame(5.0, &tb3KF2);
tb3cAnim->addKeyFrame(10.0, &tb3KF1);
thunderBird3.addAnimutator(tb3cAnim);
Animutator* tb3BcAnim = new Animutator();
KeyFrame tb3BKF1 = KeyFrame(glm::vec3(22.3989, 2, -5.00961), -90, glm::vec3(1,0,0), glm::vec3(0.4,0.4,0.4));
KeyFrame tb3BKF2 = KeyFrame(glm::vec3(22.3989, 4.5, -5.00961), -90, glm::vec3(1,0,0), glm::vec3(0.4,0.4,0.4));
tb3BcAnim->addKeyFrame(0.0, &tb3BKF1);
tb3BcAnim->addKeyFrame(5.0, &tb3BKF2);
tb3BcAnim->addKeyFrame(10.0, &tb3BKF1);
thunderBird3B.addAnimutator(tb3BcAnim);
Animutator* tb3CcAnim = new Animutator();
KeyFrame tb3CKF1 = KeyFrame(glm::vec3(-30.3314, 1.9, 2.7756), -90, glm::vec3(1,0,0), glm::vec3(0.4,0.4,0.4));
KeyFrame tb3CKF2 = KeyFrame(glm::vec3(-30.3314, 4.6, 2.7756), -90, glm::vec3(1,0,0), glm::vec3(0.4,0.4,0.4));
tb3CcAnim->addKeyFrame(0.0, &tb3CKF1);
tb3CcAnim->addKeyFrame(5.0, &tb3CKF2);
tb3CcAnim->addKeyFrame(10.0, &tb3CKF1);
thunderBird3C.addAnimutator(tb3CcAnim);
int frame_count = 0;
int last_fps = 0;
GLfloat lastTime = glfwGetTime();
loading(&tg);
do {
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT | GL_STENCIL_BUFFER_BIT);
if(tour){
t.update();
} else {
viewer->update();
}
thunderBird1.draw();
thunderBird1B.draw();
thunderBird1C.draw();
thunderBird2.draw();
t2container.draw();
thunderBird2B.draw();
t2containerB.draw();
thunderBird3.draw();
thunderBird3B.draw();
thunderBird3C.draw();
ground.draw();
skybox.draw();
if(glfwGetTime() - lastTime < 1) {
frame_count++;
} else {
lastTime = glfwGetTime();
last_fps = frame_count;
frame_count = 0;
}
if(show_help){
showHelp(tg, last_fps);
}
glfwSwapBuffers();
} while ( (glfwGetKey(GLFW_KEY_ESC) != GLFW_PRESS) &&
(glfwGetKey('q') != GLFW_PRESS) &&
(glfwGetKey('Q') != GLFW_PRESS) &&
(glfwGetWindowParam( GLFW_OPENED )) );
glfwTerminate();
//Everything was okay, return success
exit(EXIT_SUCCESS);
}
void no_recursive_search_reduce(dgrafo & g, Tour & t, Tour & best_tour, bool pt = 0) {
tour_pila_t pila;
pila.push(t);
while (!pila.empty()){
Tour actual = pila.top();
pila.pop();
if (g.size() == actual.size()) {
myVector vecinos = g.get_edges(*(actual.firt_city().first));
int p = g.peso(*(actual.last_city().first), *(actual.firt_city().first));
city temp(actual.firt_city().first, p);
actual.insert(temp);
if (pt){ actual.print();}
if (best_tour.ranking() > actual.ranking()){
best_tour = actual;}
//actual.remove_last_city();
}else{
myVector vecinos = g.get_edges(*(actual.last_city().first));
for (int i=vecinos.size()-1;i>=1;i--){
if (actual.feasible(*(vecinos[i].first))) {
actual.insert(vecinos[i]);
pila.push(actual);
actual.remove_last_city();
}
}
}
}
}
void no_recursive_search_parallel(dgrafo & g, Tour & t, Tour & best_tour) {
tour_pila_t pila;
pila.push(t);
omp_lock_t writelock;
omp_init_lock(&writelock);
int count = 0;
#pragma omp parallel
{
#pragma omp for
for (int u = 0; u<5;u++){
Tour actual;
if(!pila.empty()){
omp_set_lock(&writelock);
count++;
omp_unset_lock(&writelock);
//cout <<"Hola"<<endl;
actual = pila.top();
pila.pop();
if (g.size() == actual.size()) {
myVector vecinos = g.get_edges(*(actual.firt_city().first));
int p = g.peso(*(actual.last_city().first), *(actual.firt_city().first));
city temp(actual.firt_city().first, p);
actual.insert(temp);
omp_set_lock(&writelock);
if (best_tour.ranking() > actual.ranking()){
best_tour = actual;}
omp_unset_lock(&writelock);
//actual.remove_last_city();
}else{
myVector vecinos = g.get_edges(*(actual.last_city().first));
for (int i=vecinos.size()-1;i>=1;i--){
omp_set_lock(&writelock);
if (actual.feasible(*(vecinos[i].first))) {
actual.insert(vecinos[i]);
pila.push(actual);
actual.remove_last_city();
}
omp_unset_lock(&writelock);
}
}
}
}
}
cout<<"Contador: "<<count<<endl;
omp_destroy_lock(&writelock);
}
int main(int argc, char* argv[]) {
Tour* tour;
// If the last argument was -c or -a
if (argv[argc - 1] == NUM_CITY_DELIMITER || argv[argc - 1] == ALGORITHM_DELIMITER)
{
std::cerr << "Insufficient arguments. " + correctUsage << std::endl;
return 1;
}
// If the NO_GUI flag was found
if (std::find(argv, argv + argc, NO_GUI) != argv + argc)
{
std::cout << "GUI disabled" << std::endl;
Tour::disableGUI = true;
}
#ifndef SFML_FOUND
Tour::disableGUI = true;
#endif
// If the NUM_CITY flag was given
if (std::find(argv, argv + argc, NUM_CITY_DELIMITER) != argv + argc)
{
// If the NUM_CITY flag was given in conjunction with the city list argument
if (argv[1] == ALGORITHM_DELIMITER || argv[1] == NO_GUI)
{
std::cerr << "Incorrect usage. " + correctUsage << std::endl;
return 1;
}
// Try to parse the argument after the -c flag
try {
int numCities = std::stoi(*(std::find(argv, argv + argc, NUM_CITY_DELIMITER) + 1));
tour = new Tour(numCities);
} catch (...) {
std::cerr << "City number not given. " + correctUsage << std::endl;
return 1;
}
}
else
{
// If a city list was provided
if (argc > 1 && argv[1] != ALGORITHM_DELIMITER && argv[1] != NO_GUI)
{
std::cout << "Reading from file " << argv[1] << std::endl;
tour = new Tour(argv[1]);
}
else
{
std::cout << "Defaulting to 100 cities" << std::endl;
tour = new Tour(DEFAULT_CITY_NUM);
}
}
path * solution;
if (std::find(argv, argv + argc, ALGORITHM_DELIMITER) != argv + argc)
{
std::string algorithm;
try {
algorithm = std::string(*(std::find(argv, argv + argc, ALGORITHM_DELIMITER) + 1));
} catch (...) {
std::cerr << "Algorithm not given. " + correctUsage << std::endl;
return 1;
}
if (algorithm == ALGORITHM_RANDOM)
{
std::cout << "Using the random algorithm" << std::endl;
solution = new path(tour->solveRandom());
} else if (algorithm == ALGORITHM_TWO_OPT)
{
std::cout << "Using the 2-opt algorithm" << std::endl;
solution = new path(tour->solveRandomWithSwitches(3500.0, 2000));
} else if (algorithm == ALGORITHM_SIMULATED_ANNEALING)
{
std::cout << "Using simulated annealing" << std::endl;
solution = new path(tour->solveSimulatedAnnealing(20000.0, 20, 0.1));
}
else
{
std::cerr << "Invalid algorithm. " + correctUsage << std::endl;
return 1;
}
}
else
solution = new path(tour->solveSimulatedAnnealing(20000.0, 20, 0.1));
tour->printGraph();
std::cout << std::endl;
Tour::printPath(*solution);
std::cout << std::endl << "Tour Length: " << tour->getTourLength(*solution) << std::endl;
if (!Tour::disableGUI)
tour->getDisplayThread()->join();
delete tour;
return 0;
}
int main()
{
Tour tour;
tour.run();
return 0;
}
int main(int argc, char** argv)
{
vector<string> args(0);
vector<string> values(0);
char* test = "node1=13&nodet2=34&node3=u&node4=u&edge1-2=1&edge2-1=1&edge2-3=1&edge3-2=1&edge3-4=1&edge4-3=1&edge4-1=1&edge1-4=1&edge1-3=2&edge3-1=2&edge4-2=2&edge2-4=2&start=3&goal=3";
//strcpy(test,"http://localhost/cgi-bin/little?node1=u&node2=u&node3=u&node4=u&edge1-2=1&edge2-1=1&edge2-3=1&edge3-2=1&edge3-4=1&edge4-3=1&edge4-1=1&edge1-4=1&edge1-3=2&edge3-1=2&edge4-2=2&edge2-4=2&start=2&goal=2");
char* tachatte = getenv("QUERY_STRING");
string query_string(test);
parseQueryString(query_string,args,values);
/*
for(int i=0;i<args.size();i++)
{
cout << args[i] << ":" << values[i] << endl;
}
*/
ofstream file("zobite.txt",ios::out | ios::trunc);
Graph g(false,true);
bool boule = createGraph(g,args,values);
Tour response;
int start = atoi(values[values.size()-1].c_str());
if(boule)
response = g.little(start);
boule = boule && response.isCircuit() && ((int)g.getOrder()==(int)response.size());
cout << "Content-Type:text/xml";
cout << "<little>";
string status;
(boule) ? status="OK" : status="Aucune tournee trouvee";
cout << "<status>" << status << "</status>";
cout << "<cost>" << response.getCost() << "</cost>";
cout << "<tour>";
int* tab = response.toTab();
for(int i=0;i<response.length();i++)
{
cout << "<id>" << tab[i] << "</id>";
}
cout << "</tour>";
cout << "</little>";
file.close();
return 0;
}
