main.cpp
3.51 KB
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
#include <iostream>
#include <string>
#include <vector>
#include <unordered_map>
// import
#include "Importer/fasta.h"
#include "multithread/broadcaster.h"
// multithread
#include "multithread/buffer.h"
// Feature
#include "feature.h"
#include "feature_maker.h"
// callable
#include "feature_callable.h"
// Storage
#include "data/storage_buffer_flow.h"
#include "data/data_store_map.h"
// Data
#include "Importer/entry.h"
using namespace std;
void create_feature_maker(feature::Feature_callable *c,
string id,
data::Data_Store *ds,
multithread::Broadcaster *re);
std::vector<multithread::Buffer<data::Data_basic*>*> data_buffers;
std::vector<feature::Feature_maker*> feature_makers;
std::vector<feature::Feature_callable*> callable_entry;
std::unordered_map<std::string,data::Data_Store*> ds;
int main(int argc, char* argv[])
{
string input_file = string(argv[1]);
cout << input_file << endl;
string output = string(argv[2]);
cout << output << endl;
multithread::Buffer<data::Data_basic*> *fasta_out = new multithread::Buffer<data::Data_basic*>(THREAD_NB*BUFFER_SIZE);
fasta_out->close();
data_buffers.push_back(fasta_out);
multithread::Broadcaster re(fasta_out);
data::Data_Store_Map *tmp_ds;
/* ************** KMER ******************/
feature::Feature_callable *kmer = static_cast<feature::Feature_callable*>(feature::Kmer::get_callable(3));
tmp_ds = new data::Data_Store_Map(output+"KMER3.txt");
ds["KMER3"] = tmp_ds;
create_feature_maker(kmer,feature::Kmer::TYPE+std::to_string(6),tmp_ds,&re);
// kmer = static_cast<feature::Feature_callable*>(feature::Kmer::get_callable(6));
// tmp_ds = new data::Data_Store_Map(output+"KMER6.txt");
// ds["KMER6"] = tmp_ds;
// create_feature_maker(kmer,feature::Kmer::TYPE+std::to_string(6),tmp_ds,&re);
/* *********** ORF *********************/
/* ORF callable */
feature::Feature_callable *orf = static_cast<feature::Feature_callable*>(feature::ORF::get_callable());
tmp_ds = new data::Data_Store_Map(output+"ORF.txt");
ds["ORF"] = tmp_ds;
create_feature_maker(orf,feature::ORF::TYPE,tmp_ds,&re);
/* *********** Codon Biases ***************/
feature::Feature_callable *cb = static_cast<feature::Feature_callable*>(feature::CodonBiases::get_callable());
tmp_ds = new data::Data_Store_Map(output+"CB.txt");
ds["CB"] = tmp_ds;
create_feature_maker(cb,feature::CodonBiases::TYPE,tmp_ds,&re);
/* *********** Fasta reader *************/
import::Fasta reader(fasta_out,input_file);
re.start();
reader.start();
/* **** join process ******/
reader.join();
re.join();
// shutdown pointer
while(!feature_makers.empty()){
feature::Feature_maker *fm = feature_makers.back();
feature_makers.pop_back();
delete fm;
}
while(!data_buffers.empty()){
multithread::Buffer<data::Data_basic*> *tmp_dk = data_buffers.back();
data_buffers.pop_back();
delete tmp_dk;
}
while(!callable_entry.empty()){
feature::Feature_callable *tmp_dk = callable_entry.back();
callable_entry.pop_back();
delete tmp_dk;
}
auto it = ds.begin();
while(it != ds.end()){
delete it->second;
++it;
}
return 0;
}
void create_feature_maker(feature::Feature_callable *c,
string id,
data::Data_Store *ds,
multithread::Broadcaster *re){
callable_entry.push_back(c);
multithread::Buffer<data::Data_basic*> *buf = new multithread::Buffer<data::Data_basic*>(THREAD_NB*BUFFER_SIZE);
data_buffers.push_back(buf);
feature::Feature_maker *fm = new feature::Feature_maker(c,buf,ds,id);
feature_makers.push_back(fm);
re->add_worker(buf);
}