-
Notifications
You must be signed in to change notification settings - Fork 80
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
6 changed files
with
259 additions
and
6 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,46 @@ | ||
#include <seqan3/alphabet/nucleotide/dna4.hpp> | ||
#include <seqan3/core/debug_stream.hpp> | ||
#include <seqan3/range/views/kmer_hash.hpp> | ||
#include <seqan3/search/dream_index/interleaved_bloom_filter.hpp> | ||
|
||
using seqan3::operator""_dna4; | ||
|
||
int main() | ||
{ | ||
seqan3::interleaved_bloom_filter ibf{seqan3::bin_count{8u}, | ||
seqan3::bin_size{8192u}, | ||
seqan3::hash_function_count{2u}}; | ||
|
||
auto const sequence1 = "ACTGACTGACTGATC"_dna4; | ||
auto const sequence2 = "GTGACTGACTGACTCG"_dna4; | ||
auto const sequence3 = "AAAAAAACGATCGACA"_dna4; | ||
auto hash_adaptor = seqan3::views::kmer_hash(seqan3::ungapped{5u}); | ||
|
||
// Insert all 5-mers of sequence1 into bin 0 | ||
for (auto && value : sequence1 | hash_adaptor) | ||
ibf.emplace(value, seqan3::bin_index{0u}); | ||
|
||
// Insert all 5-mers of sequence2 into bin 4 | ||
for (auto && value : sequence2 | hash_adaptor) | ||
ibf.emplace(value, seqan3::bin_index{4u}); | ||
|
||
// Insert all 5-mers of sequence3 into bin 7 | ||
for (auto && value : sequence3 | hash_adaptor) | ||
ibf.emplace(value, seqan3::bin_index{7u}); | ||
|
||
auto agent = ibf.counting_agent(); | ||
|
||
// Count all 5-mers of sequence1 for all bins | ||
seqan3::debug_stream << agent.bulk_count(sequence1 | hash_adaptor) << '\n'; // [11,0,0,0,9,0,0,0] | ||
|
||
// Search for specific values | ||
std::vector<size_t> const values{92, 1238, 812, 81273}; | ||
seqan3::debug_stream << agent.bulk_count(values) << '\n'; // [0,0,0,0,0,0,0,0] | ||
seqan3::debug_stream << agent.bulk_count(std::views::iota(0u, 1024u)) << '\n'; // [6,0,0,0,7,0,0,10] | ||
|
||
// The default counters are 16 bit unsigned integer. | ||
// An optional template parameter can be used to specify the counter type | ||
auto agent2 = ibf.counting_agent<uint8_t>(); | ||
// The returned counts are now 8 bit unsigned integers. | ||
seqan3::debug_stream << agent.bulk_count(sequence1 | hash_adaptor) << '\n'; // [11,0,0,0,9,0,0,0] | ||
} |
18 changes: 18 additions & 0 deletions
18
test/snippet/search/dream_index/counting_agent_construction.cpp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,18 @@ | ||
#include <seqan3/search/dream_index/interleaved_bloom_filter.hpp> | ||
|
||
int main() | ||
{ | ||
// Construct an Interleaved Bloom Filter to be used with the counting_agent. | ||
seqan3::interleaved_bloom_filter ibf{seqan3::bin_count{43u}, | ||
seqan3::bin_size{8192u}, | ||
seqan3::hash_function_count{3}}; | ||
|
||
// The counting_agent can now be constructed by calling `counting_agent` on the Interleaved Bloom Filter. | ||
auto agent = ibf.counting_agent(); | ||
|
||
// Calling `increase_bin_number_to` invalidates the agent. | ||
ibf.increase_bin_number_to(seqan3::bin_count{60u}); | ||
|
||
// So make sure to construct a new counting_agent. | ||
agent = ibf.counting_agent(); | ||
} |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters