23namespace seqan3::detail
40template <
typename traits_t,
typename alignment_matrix_t>
41 requires (is_type_specialisation_of_v<traits_t, alignment_configuration_traits> &&
42 requires (alignment_matrix_t & matrix,
typename traits_t::score_type
const initial_score)
44 { matrix.resize(column_index_type{
size_t{}}, row_index_type{
size_t{}}, initial_score) };
46class policy_alignment_matrix
50 using score_type =
typename traits_t::score_type;
52 using matrix_index_type =
typename traits_t::matrix_index_type;
55 int32_t lower_diagonal{};
57 int32_t upper_diagonal{};
59 bool last_column_is_free{};
61 bool last_row_is_free{};
66 policy_alignment_matrix() =
default;
67 policy_alignment_matrix(policy_alignment_matrix
const &) =
default;
68 policy_alignment_matrix(policy_alignment_matrix &&) =
default;
69 policy_alignment_matrix & operator=(policy_alignment_matrix
const &) =
default;
70 policy_alignment_matrix & operator=(policy_alignment_matrix &&) =
default;
71 ~policy_alignment_matrix() =
default;
86 template <
typename alignment_configuration_t>
87 requires (is_type_specialisation_of_v<alignment_configuration_t, configuration>)
88 policy_alignment_matrix(alignment_configuration_t
const & config)
94 lower_diagonal = band.lower_diagonal;
95 upper_diagonal = band.upper_diagonal;
97 bool invalid_band = upper_diagonal < lower_diagonal;
98 std::string error_cause = (invalid_band) ?
" The upper diagonal is smaller than the lower diagonal." :
"";
100 if constexpr (traits_t::is_global)
102 auto method_global_config = get<seqan3::align_cfg::method_global>(config);
104 bool first_row_is_free = method_global_config.free_end_gaps_sequence1_leading;
105 bool first_column_is_free = method_global_config.free_end_gaps_sequence2_leading;
107 last_row_is_free = method_global_config.free_end_gaps_sequence1_trailing;
108 last_column_is_free = method_global_config.free_end_gaps_sequence2_trailing;
110 invalid_band |= (upper_diagonal < 0 && !first_column_is_free) || (lower_diagonal > 0 && !first_row_is_free);
111 error_cause +=
" The band starts in a region without free gaps.";
115 throw invalid_alignment_configuration{
"The selected band [" +
std::to_string(lower_diagonal) +
":"
117 +
"] cannot be used with the current "
118 "alignment configuration:"
145 auto acquire_matrices(
size_t const sequence1_size,
146 size_t const sequence2_size,
147 score_type initial_score = score_type{})
const
152 if constexpr (traits_t::is_banded)
153 check_valid_band_configuration(sequence1_size, sequence2_size);
155 static thread_local alignment_matrix_t alignment_matrix{};
156 static thread_local coordinate_matrix<matrix_index_type> index_matrix{};
159 size_t const column_count = sequence1_size + 1;
160 size_t row_count = sequence2_size + 1;
162 index_matrix.resize(column_index_type{column_count}, row_index_type{row_count});
164 if constexpr (traits_t::is_banded)
166 assert(upper_diagonal - lower_diagonal + 1 > 0);
168 row_count = std::min<int64_t>(upper_diagonal - lower_diagonal + 2, row_count);
171 alignment_matrix.resize(column_index_type{column_count}, row_index_type{row_count}, initial_score);
173 return std::tie(alignment_matrix, index_matrix);
184 void check_valid_band_configuration(
size_t const sequence1_size,
size_t const sequence2_size)
const
186 bool const upper_diagonal_ends_before_last_cell = (upper_diagonal + sequence2_size) < sequence1_size;
187 bool const lower_diagonal_ends_behind_last_cell = (-lower_diagonal + sequence1_size) < sequence2_size;
189 bool invalid_band =
false;
192 if constexpr (traits_t::is_global)
195 invalid_band |= (lower_diagonal_ends_behind_last_cell && !last_column_is_free)
196 || (upper_diagonal_ends_before_last_cell && !last_row_is_free);
197 error_cause =
"The band ends in a region without free gaps.";
201 throw invalid_alignment_configuration{
"The selected band [" +
std::to_string(lower_diagonal) +
":"
203 +
"] cannot be used with the current "
204 "alignment configuration: "
Includes customized exception types for the alignment module .
Provides helper type traits for the configuration and execution of the alignment algorithm.
Configuration element for setting a fixed size band.
Definition: align_config_band.hpp:63
Provides seqan3::configuration and utility functions.
Provides seqan3::detail::coordinate_matrix.
constexpr auto const & get(configuration< configs_t... > const &config) noexcept
This is an overloaded member function, provided for convenience. It differs from the above function o...
Definition: configuration.hpp:415
Provides type traits for working with templates.