25namespace seqan3::detail
65template <simd_concept simd_score_t, semialphabet alphabet_t,
typename alignment_t>
66 requires (seqan3::alphabet_size<alphabet_t> > 1)
67 && (std::same_as<alignment_t, align_cfg::method_local> || std::same_as<alignment_t, align_cfg::method_global>)
68class simd_match_mismatch_scoring_scheme
72 using scalar_type =
typename simd_traits<simd_score_t>::scalar_type;
74 using alphabet_ranks_type = simd_score_t;
78 static constexpr scalar_type padding_symbol =
static_cast<scalar_type
>(1u << (bits_of<scalar_type> - 1));
84 constexpr simd_match_mismatch_scoring_scheme() =
default;
86 constexpr simd_match_mismatch_scoring_scheme(simd_match_mismatch_scoring_scheme
const &) =
default;
88 constexpr simd_match_mismatch_scoring_scheme(simd_match_mismatch_scoring_scheme &&) =
default;
90 constexpr simd_match_mismatch_scoring_scheme & operator=(simd_match_mismatch_scoring_scheme
const &) =
default;
92 constexpr simd_match_mismatch_scoring_scheme & operator=(simd_match_mismatch_scoring_scheme &&) =
default;
94 ~simd_match_mismatch_scoring_scheme() =
default;
97 template <
typename scoring_scheme_t>
99 constexpr explicit simd_match_mismatch_scoring_scheme(scoring_scheme_t
const & scoring_scheme)
101 initialise_from_scalar_scoring_scheme(scoring_scheme);
105 template <
typename scoring_scheme_t>
107 constexpr simd_match_mismatch_scoring_scheme & operator=(scoring_scheme_t
const & scoring_scheme)
109 initialise_from_scalar_scoring_scheme(scoring_scheme);
143 constexpr simd_score_t score(alphabet_ranks_type
const & ranks1, alphabet_ranks_type
const & ranks2)
const noexcept
145 typename simd_traits<simd_score_t>::mask_type mask;
148 if constexpr (std::same_as<alignment_t, align_cfg::method_global>)
149 mask = (ranks1 ^ ranks2) <= simd::fill<simd_score_t>(0);
151 mask = (ranks1 ^ ranks2) == simd::fill<simd_score_t>(0);
153 return mask ? match_score : mismatch_score;
158 constexpr auto padding_match_score() noexcept
160 return match_score[0];
164 template <
typename alphabet_ranks_t>
165 requires std::same_as<std::remove_cvref_t<alphabet_ranks_t>, alphabet_ranks_type>
166 constexpr alphabet_ranks_t make_score_profile(alphabet_ranks_t && ranks)
const noexcept
168 return std::forward<alphabet_ranks_t>(ranks);
187 template <
typename scoring_scheme_t>
188 constexpr void initialise_from_scalar_scoring_scheme(scoring_scheme_t
const & scoring_scheme)
190 using score_t =
decltype(std::declval<scoring_scheme_t const &>().score(alphabet_t{}, alphabet_t{}));
191 using simd_scalar_t =
typename simd_traits<simd_score_t>::scalar_type;
193 score_t scalar_match_score =
195 score_t scalar_mismatch_score =
199 if constexpr (
sizeof(simd_scalar_t) <
sizeof(score_t))
205 "for the selected scalar type of the simd type."};
209 match_score = simd::fill<simd_score_t>(
static_cast<simd_scalar_t
>(scalar_match_score));
210 mismatch_score = simd::fill<simd_score_t>(
static_cast<simd_scalar_t
>(scalar_mismatch_score));
213 simd_score_t match_score;
214 simd_score_t mismatch_score;
Provides algorithms to modify seqan3::simd::simd_type.
Provides global and local alignment configurations.
Core alphabet concept and free function/type trait wrappers.
Provides utility functions for bit twiddling.
Adaptions of concepts from the Cereal library.
constexpr auto assign_rank_to
Assign a rank to an alphabet object.
Definition: concept.hpp:293
A concept that requires that type be able to score two letters.
Provides seqan3::scoring_scheme_for.
Provides seqan3::simd::simd_concept.