20namespace seqan3::detail
37template <
typename traits_t,
typename alignment_matrix_t>
38 requires (is_type_specialisation_of_v<traits_t, alignment_configuration_traits> &&
39 requires (alignment_matrix_t & matrix,
typename traits_t::score_type
const initial_score)
41 { matrix.resize(column_index_type{
size_t{}}, row_index_type{
size_t{}}, initial_score) };
43class policy_alignment_matrix
47 using score_type =
typename traits_t::score_type;
49 using matrix_index_type =
typename traits_t::matrix_index_type;
52 int32_t lower_diagonal{};
54 int32_t upper_diagonal{};
56 bool last_column_is_free{};
58 bool last_row_is_free{};
63 policy_alignment_matrix() =
default;
64 policy_alignment_matrix(policy_alignment_matrix
const &) =
default;
65 policy_alignment_matrix(policy_alignment_matrix &&) =
default;
66 policy_alignment_matrix & operator=(policy_alignment_matrix
const &) =
default;
67 policy_alignment_matrix & operator=(policy_alignment_matrix &&) =
default;
68 ~policy_alignment_matrix() =
default;
83 template <
typename alignment_configuration_t>
84 requires (is_type_specialisation_of_v<alignment_configuration_t, configuration>)
85 policy_alignment_matrix(alignment_configuration_t
const & config)
91 lower_diagonal = band.lower_diagonal;
92 upper_diagonal = band.upper_diagonal;
94 bool invalid_band = upper_diagonal < lower_diagonal;
95 std::string error_cause = (invalid_band) ?
" The upper diagonal is smaller than the lower diagonal." :
"";
97 if constexpr (traits_t::is_global)
99 auto method_global_config = get<seqan3::align_cfg::method_global>(config);
101 bool first_row_is_free = method_global_config.free_end_gaps_sequence1_leading;
102 bool first_column_is_free = method_global_config.free_end_gaps_sequence2_leading;
104 last_row_is_free = method_global_config.free_end_gaps_sequence1_trailing;
105 last_column_is_free = method_global_config.free_end_gaps_sequence2_trailing;
107 invalid_band |= (upper_diagonal < 0 && !first_column_is_free) || (lower_diagonal > 0 && !first_row_is_free);
108 error_cause +=
" The band starts in a region without free gaps.";
112 throw invalid_alignment_configuration{
"The selected band [" +
std::to_string(lower_diagonal) +
":"
114 +
"] cannot be used with the current "
115 "alignment configuration:"
142 auto acquire_matrices(
size_t const sequence1_size,
143 size_t const sequence2_size,
144 score_type initial_score = score_type{})
const
149 if constexpr (traits_t::is_banded)
150 check_valid_band_configuration(sequence1_size, sequence2_size);
152 static thread_local alignment_matrix_t alignment_matrix{};
153 static thread_local coordinate_matrix<matrix_index_type> index_matrix{};
156 size_t const column_count = sequence1_size + 1;
157 size_t row_count = sequence2_size + 1;
159 index_matrix.resize(column_index_type{column_count}, row_index_type{row_count});
161 if constexpr (traits_t::is_banded)
163 assert(upper_diagonal - lower_diagonal + 1 > 0);
165 row_count = std::min<int64_t>(upper_diagonal - lower_diagonal + 2, row_count);
168 alignment_matrix.resize(column_index_type{column_count}, row_index_type{row_count}, initial_score);
170 return std::tie(alignment_matrix, index_matrix);
181 void check_valid_band_configuration(
size_t const sequence1_size,
size_t const sequence2_size)
const
183 bool const upper_diagonal_ends_before_last_cell = (upper_diagonal + sequence2_size) < sequence1_size;
184 bool const lower_diagonal_ends_behind_last_cell = (-lower_diagonal + sequence1_size) < sequence2_size;
186 bool invalid_band =
false;
189 if constexpr (traits_t::is_global)
192 invalid_band |= (lower_diagonal_ends_behind_last_cell && !last_column_is_free)
193 || (upper_diagonal_ends_before_last_cell && !last_row_is_free);
194 error_cause =
"The band ends in a region without free gaps.";
198 throw invalid_alignment_configuration{
"The selected band [" +
std::to_string(lower_diagonal) +
":"
200 +
"] cannot be used with the current "
201 "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:60
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:412
Provides type traits for working with templates.