-
Notifications
You must be signed in to change notification settings - Fork 0
/
Experiment.cpp
53 lines (46 loc) · 2.19 KB
/
Experiment.cpp
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
//
// Created by osmangokalp on 10.06.2020.
//
#include "Experiment.h"
#include "GRASP.h"
#include <fstream>
#include <iostream>
void Experiment::saveToAFile(std::string outputFileName, std::string line) const {
std::ofstream outFile;
outFile.open(experimentsFolder + std::string("/") + outputFileName, std::ios_base::app);
if (outFile.fail()) {
std::cerr << "unable to open file << outputFileName << for writing" << std::endl;
exit(1);
}
outFile << line << std::endl;
outFile.close();
}
Experiment::Experiment(const std::string &experimentsFolder) : experimentsFolder(experimentsFolder) {}
void Experiment::doGRASPExperiment(int numTry, int maxEval, int motifLength, std::string problemInstanceName) const {
GRASP grasp;
double alphaArr[] = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0};
double candidateRatioArr[] = {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0};
bool firstImpArr[] = {true, false};
Problem problem(problemInstanceName);
for (int i = 0; i < 2; ++i) {
bool firstImp = firstImpArr[i];
for (int j = 0; j < 11; ++j) {
double alpha = alphaArr[j];
for (int k = 0; k < 10; ++k) {
double candidateRatio = candidateRatioArr[k];
for (int l = 0; l < numTry; ++l) {
std::default_random_engine generator(100 + l);
Solution *solution = grasp.GRASPMotifSearch(problem, motifLength, alpha, candidateRatio, firstImp,
generator, maxEval);
std::string line = problemInstanceName + " " +
std::to_string(l) + " " + std::to_string(firstImp) + " " +
std::to_string(alpha) + " " +
std::to_string(candidateRatio) + " " + std::to_string(maxEval) + " " +
std::to_string(motifLength) + " " + solution->consensusString + " " +
std::to_string(solution->similarityScore);
saveToAFile("allResults.txt", line);
}
}
}
}
}