-
Notifications
You must be signed in to change notification settings - Fork 1
v2 createMQF() entered infinite loop #96
Copy link
Copy link
Open
Description
When I used this FASTA file with the following script, createMQF() and createBMQF() functions entered an infinite loop.
#include <iostream>
#include "kDataFrame.hpp"
#include <string>
#include <vector>
#include "algorithms.hpp"
#include <math.h>
#include "kDataframes/kDataFrameMQF.hpp"
#include "kDataframes/kDataFrameSTL.hpp"
int main(int argc, char** argv) {
string file_name = argv[1];
int chunk_size = 10;
int kSize = 25;
readingModes RM = KMERS;
hashingModes HM = integer_hasher;
auto* KF = kDataFrameFactory::createPHMAP(kSize, HM);
kmerDecoder* KD = kmerDecoder::getInstance(file_name, chunk_size, RM, HM, { {"kSize", kSize} });
while (!KD->end()) {
KD->next_chunk();
for (const auto& seq : *KD->getKmers())
for (const auto& kmer : seq.second)
KF->insert(kmer.hash);
}
cout << "converting to MQF..." << endl;
kDataFrame* kf_mqf = kDataFrameFactory::createMQF(KF);
kf_mqf->save("test");
}Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels