53template <simd_concept simd_score_t, semialphabet alphabet_t,
typename alignment_t>
54 requires (std::same_as<alignment_t, align_cfg::method_local> || std::same_as<alignment_t, align_cfg::method_global>)
66 "The selected simd scalar type is not large enough to represent the given alphabet including an "
67 "additional padding symbol!");
70 "The selected simd scalar type is not large enough to represent the given alphabet including an "
71 "additional padding symbol!");
74 static constexpr bool is_global = std::same_as<alignment_t, align_cfg::method_global>;
76 static constexpr size_t index_offset = seqan3::alphabet_size<alphabet_t> + 1;
78 static constexpr scalar_type score_for_padding_symbol = (is_global) ? 1 : -1;
98 template <
typename scoring_scheme_t>
101 initialise_from_scalar_scoring_scheme(scoring_scheme);
105 template <
typename scoring_scheme_t>
108 initialise_from_scalar_scoring_scheme(scoring_scheme);
140 simd_score_t
const matrix_index = score_profile + ranks;
141 simd_score_t result{};
143 for (
size_t idx = 0; idx < simd_traits<simd_score_t>::length; ++idx)
144 result[idx] = scoring_scheme_data.data()[
matrix_index[idx]];
153 return score_for_padding_symbol;
168 return ranks * simd::fill<simd_score_t>(index_offset);
180 template <
typename scoring_scheme_t>
184 using score_t =
decltype(std::declval<scoring_scheme_t const &>().score(alphabet_t{}, alphabet_t{}));
187 [[maybe_unused]]
auto check_score_range = [&]([[maybe_unused]] score_t score)
192 if constexpr (
sizeof(
scalar_type) <
sizeof(score_t))
197 if (score > max_score_value || score < min_score_value)
199 "for the selected scalar type of the simd type."};
204 scoring_scheme_data.resize(index_offset * index_offset, score_for_padding_symbol);
208 auto data_it = scoring_scheme_data.begin();
209 for (alphabet_size_t lhs_rank = 0; lhs_rank < seqan3::alphabet_size<alphabet_t>; ++lhs_rank)
211 for (alphabet_size_t rhs_rank = 0; rhs_rank < seqan3::alphabet_size<alphabet_t>; ++rhs_rank, ++data_it)
216 check_score_range(tmp_score);
217 *data_it = tmp_score;
Provides algorithms to modify seqan3::simd::simd_type.
Provides global and local alignment configurations.
Core alphabet concept and free function/type trait wrappers.
Adaptions of concepts from the Cereal library.
A vectorised scoring scheme to handle scoring matrices using gather strategy.
Definition simd_matrix_scoring_scheme.hpp:56
constexpr simd_matrix_scoring_scheme & operator=(scoring_scheme_t const &scoring_scheme)
Store the given scoring scheme matrix into a private member variable.
Definition simd_matrix_scoring_scheme.hpp:106
constexpr simd_score_t score(simd_score_profile_type const &score_profile, simd_alphabet_ranks_type const &ranks) const noexcept
Given two simd vectors, compute an element-wise score and return another simd vector.
Definition simd_matrix_scoring_scheme.hpp:137
simd_score_t simd_alphabet_ranks_type
The type of the simd vector representing the alphabet ranks of one sequence batch.
Definition simd_matrix_scoring_scheme.hpp:63
~simd_matrix_scoring_scheme()=default
Defaulted.
constexpr simd_matrix_scoring_scheme & operator=(simd_matrix_scoring_scheme const &)=default
Defaulted.
constexpr simd_matrix_scoring_scheme()=default
Defaulted.
constexpr simd_matrix_scoring_scheme(scoring_scheme_t const &scoring_scheme)
Store the given scoring scheme matrix into a private member variable.
Definition simd_matrix_scoring_scheme.hpp:99
constexpr scalar_type padding_match_score() const noexcept
Returns the score used when aligning a padding symbol.
Definition simd_matrix_scoring_scheme.hpp:151
constexpr void initialise_from_scalar_scoring_scheme(scoring_scheme_t const &scoring_scheme)
Store the given scoring scheme matrix into a private member variable.
Definition simd_matrix_scoring_scheme.hpp:182
constexpr simd_matrix_scoring_scheme(simd_matrix_scoring_scheme &&)=default
Defaulted.
constexpr simd_matrix_scoring_scheme & operator=(simd_matrix_scoring_scheme &&)=default
Defaulted.
constexpr simd_matrix_scoring_scheme(simd_matrix_scoring_scheme const &)=default
Defaulted.
simd_score_t simd_score_profile_type
The score profile type used for this scoring scheme, which is the same as the simd score type.
Definition simd_matrix_scoring_scheme.hpp:61
typename simd_traits< simd_score_t >::scalar_type scalar_type
The underlying scalar type of the simd vector.
Definition simd_matrix_scoring_scheme.hpp:59
constexpr simd_score_profile_type make_score_profile(simd_alphabet_ranks_type const &ranks) const noexcept
Converts the simd alphabet ranks into a score profile used for scoring it later with the alphabet ran...
Definition simd_matrix_scoring_scheme.hpp:166
constexpr auto assign_rank_to
Assign a rank to an alphabet object.
Definition alphabet/concept.hpp:290
A concept that requires that type be able to score two letters.
The internal SeqAn3 namespace.
Definition aligned_sequence_concept.hpp:26
Provides seqan3::scoring_scheme_for.
A representation of a location or offset within a two-dimensional matrix.
Definition matrix_coordinate.hpp:87
seqan3::simd::simd_traits is the trait class that provides uniform interface to the properties of sim...
Definition simd_traits.hpp:38
Provides seqan3::simd::simd_concept.