Skip to content

Commit

Permalink
Set split 1 at the start of Pref Index Reader so it writes correct meta
Browse files Browse the repository at this point in the history
  • Loading branch information
milot-mirdita committed Jan 13, 2025
1 parent d9823e8 commit 0898eb9
Showing 1 changed file with 4 additions and 5 deletions.
9 changes: 4 additions & 5 deletions src/prefiltering/PrefilteringIndexReader.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -58,6 +58,10 @@ void PrefilteringIndexReader::createIndexFile(const std::string &outDB,
bool compBiasCorrection, int alphabetSize, int kmerSize, int maskMode,
int maskLowerCase, float maskProb, int maskNrepeats, int kmerThr, int targetSearchMode, int splits,
int indexSubset) {
const bool noKmerIndex = (indexSubset & Parameters::INDEX_SUBSET_NO_PREFILTER) != 0;
if (noKmerIndex) {
splits = 1;
}

const int SPLIT_META = splits > 1 ? 0 : 0;
const int SPLIT_SEQS = splits > 1 ? 1 : 0;
Expand Down Expand Up @@ -190,11 +194,6 @@ void PrefilteringIndexReader::createIndexFile(const std::string &outDB,
(Parameters::isEqualDbtype(seqType, Parameters::DBTYPE_NUCLEOTIDES) || Parameters::isEqualDbtype(seqType, Parameters::DBTYPE_AMINO_ACIDS))
? alphabetSize -1: alphabetSize;

const bool noKmerIndex = (indexSubset & Parameters::INDEX_SUBSET_NO_PREFILTER) != 0;
if (noKmerIndex) {
splits = 1;
}

ScoreMatrix s3;
ScoreMatrix s2;
if (Parameters::isEqualDbtype(seqType, Parameters::DBTYPE_HMM_PROFILE) == false && noKmerIndex == false) {
Expand Down

0 comments on commit 0898eb9

Please sign in to comment.