Skip to content

Commit d65afaf

Browse files
committed
fix case where number of kmer < max_entries AND no formatting :)
1 parent 3c6e38f commit d65afaf

File tree

1 file changed

+4
-2
lines changed

1 file changed

+4
-2
lines changed

modiscolite/gapped_kmer.py

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -93,8 +93,10 @@ def _extract_gkmers(X, min_k, max_k, max_gap, max_len, max_entries):
9393

9494
idxs = np.argsort(-np.abs(scores_), kind='mergesort')[:max_entries]
9595

96-
keys[xi] = keys_[idxs]
97-
scores[xi] = scores_[idxs]
96+
# in case there are less kmers than max_entries
97+
n_kmer = len(idxs)
98+
keys[xi, :n_kmer] = keys_[idxs]
99+
scores[xi, :n_kmer] = scores_[idxs]
98100

99101
return keys, scores
100102

0 commit comments

Comments
 (0)