Skip to content
Snippets Groups Projects
Commit bf440c2b authored by Vandamme Léa's avatar Vandamme Léa
Browse files

add memory usage for 1 kmer

parent 95232450
No related branches found
No related tags found
No related merge requests found
# Results
## Indexing with k-mers of size 15
### Using phage_5kb_10percent.fa
| | Time | Resource usage | Resource usage (1 k-mer) |
|--- |:-: |:-: |:-: |
| Naive version | 4 seconds | 842 Mo | 0.17 Ko |
| Robin Hood Hashtable <br/>& K_mer uint64_t | FACE WITH TEARS OF JOY | ```&#x1F602;``` | ```&#128514;``` |
| Delta encoding ||||
\ No newline at end of file
......@@ -43,6 +43,7 @@ void Index::index_fasta(const string& read_file, int k){
umap hmap;
uint32_t num_read = 0;
uint32_t num_kmers = 0;
ifstream fichier(read_file, ios::in);
if(fichier)
{
......@@ -51,6 +52,7 @@ void Index::index_fasta(const string& read_file, int k){
getline(fichier,ligne);
if (ligne[0] != '>' && !ligne.empty()){ // ONLY GET SEQUENCES
for (long unsigned int i=0; i<=ligne.length()-k; i++){
num_kmers ++;
string kmer = ligne.substr(i, k);
hmap[kmer].push_back(num_read); // CHECK IF THE KMER IS ALREADY IN THE HASHMAP AND INSERT
}
......@@ -62,7 +64,7 @@ void Index::index_fasta(const string& read_file, int k){
else{
cerr << "Error opening the file." << endl;
}
set_nb_kmers(num_read);
set_nb_kmers(num_kmers);
set_index(hmap);
}
......
......@@ -84,7 +84,9 @@ int main(int argc, char *argv[])
}
}
cout << "Resource usage : " << to_string(getMemorySelfMaxUsed()) << "Ko" << endl;
uint64_t memory_used = getMemorySelfMaxUsed();
cout << "Resource usage : " << memory_used << "Ko" << endl;
cout << "Resource usage (for 1 k-mer) : " << (double)memory_used/index.get_nb_kmers() << "Ko" << endl;
}
return 0;
}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment