25 namespace seqan3::detail
52 template <
typename trace_t,
bool coordinate_only = false>
53 class alignment_trace_matrix_full_banded :
54 protected alignment_trace_matrix_base<trace_t>,
55 public alignment_matrix_column_major_range_base<alignment_trace_matrix_full_banded<trace_t, coordinate_only>>
58 static_assert(std::same_as<trace_t, trace_directions> || simd_concept<trace_t>,
59 "Value type must either be a trace_directions object or a simd vector.");
62 using matrix_base_t = alignment_trace_matrix_base<trace_t>;
64 using range_base_t = alignment_matrix_column_major_range_base<alignment_trace_matrix_full_banded<trace_t,
70 using typename matrix_base_t::element_type;
71 using typename matrix_base_t::coordinate_type;
72 using typename range_base_t::alignment_column_type;
75 decltype(std::views::iota(coordinate_type{}, coordinate_type{})),
78 std::views::iota(coordinate_type{}, coordinate_type{})))>;
85 using value_type = alignment_trace_matrix_proxy<coordinate_type,
87 detail::ignore_t
const,
90 using reference = value_type;
92 using iterator =
typename range_base_t::iterator;
94 using sentinel =
typename range_base_t::sentinel;
95 using typename matrix_base_t::size_type;
102 constexpr alignment_trace_matrix_full_banded() =
default;
104 constexpr alignment_trace_matrix_full_banded(alignment_trace_matrix_full_banded
const &) =
default;
106 constexpr alignment_trace_matrix_full_banded(alignment_trace_matrix_full_banded &&) =
default;
108 constexpr alignment_trace_matrix_full_banded & operator=(alignment_trace_matrix_full_banded
const &) =
default;
110 constexpr alignment_trace_matrix_full_banded & operator=(alignment_trace_matrix_full_banded &&) =
default;
112 ~alignment_trace_matrix_full_banded() =
default;
129 template <std::ranges::forward_range first_sequence_t, std::ranges::forward_range second_sequence_t>
130 constexpr alignment_trace_matrix_full_banded(first_sequence_t && first,
131 second_sequence_t && second,
132 align_cfg::band_fixed_size
const & band,
133 [[maybe_unused]] trace_t
const initial_value = trace_t{})
135 matrix_base_t::num_cols =
static_cast<size_type
>(std::ranges::distance(first) + 1);
136 matrix_base_t::num_rows =
static_cast<size_type
>(std::ranges::distance(second) + 1);
138 band_col_index = std::min<int32_t>(std::max<int32_t>(band.upper_diagonal, 0),
139 matrix_base_t::num_cols - 1);
140 band_row_index = std::min<int32_t>(std::abs(std::min<int32_t>(band.lower_diagonal, 0)),
141 matrix_base_t::num_rows - 1);
142 band_size = band_col_index + band_row_index + 1;
145 if constexpr (!coordinate_only)
147 matrix_base_t::data =
typename matrix_base_t::pool_type{number_rows{
static_cast<size_type
>(band_size)},
148 number_cols{matrix_base_t::num_cols}};
149 matrix_base_t::cache_left.resize(band_size + 1, initial_value);
155 auto trace_path(matrix_coordinate
const & trace_begin)
157 static_assert(!coordinate_only,
"Requested trace but storing the trace was disabled!");
159 using matrix_iter_t = std::ranges::iterator_t<typename matrix_base_t::pool_type>;
160 using trace_iterator_t = trace_iterator_banded<matrix_iter_t>;
161 using path_t = std::ranges::subrange<trace_iterator_t, std::default_sentinel_t>;
163 if (trace_begin.row >=
static_cast<size_t>(band_size) || trace_begin.col >= matrix_base_t::num_cols)
166 return path_t{trace_iterator_t{matrix_base_t::data.begin() + matrix_offset{trace_begin},
167 column_index_type{band_col_index}},
168 std::default_sentinel};
172 int32_t band_col_index{};
174 int32_t band_row_index{};
180 constexpr alignment_column_type initialise_column(size_type
const column_index) noexcept
182 int32_t slice_begin = std::max<int32_t>(0, band_col_index - column_index);
183 int32_t row_end_index = column_index - band_col_index + band_size;
184 int32_t slice_end = band_size - std::max<int32_t>(row_end_index - matrix_base_t::num_rows, 0);
186 assert(row_end_index >= 0);
187 assert(slice_begin >= 0);
188 assert(slice_end > 0);
189 assert(slice_begin < slice_end);
191 coordinate_type row_begin{column_index_type{column_index}, row_index_type{
static_cast<size_type
>(slice_begin)}};
192 coordinate_type row_end{column_index_type{column_index}, row_index_type{
static_cast<size_type
>(slice_end)}};
193 if constexpr (coordinate_only)
195 return alignment_column_type{*
this,
196 column_data_view_type{std::views::iota(
std::move(row_begin),
201 matrix_coordinate band_begin{row_index_type{
static_cast<size_type
>(slice_begin)},
202 column_index_type{column_index}};
203 size_type slice_size = slice_end - slice_begin;
209 return alignment_column_type{*
this, column_data_view_type{
std::move(col)}};
214 template <std::random_access_iterator iter_t>
215 constexpr value_type make_proxy(iter_t host_iter) noexcept
217 if constexpr (coordinate_only)
219 return {*host_iter, std::ignore, std::ignore, std::ignore, std::ignore};
223 return {std::get<2>(*host_iter),
224 std::get<0>(*host_iter),
225 std::get<1>(*(host_iter + 1)),
226 std::get<1>(*host_iter),
227 matrix_base_t::cache_up
Provides seqan3::detail::align_config_band.
Provides seqan3::detail::alignment_matrix_column_major_range_base.
Provides seqan3::detail::alignment_trace_matrix_base.
Provides seqan3::detail::alignment_trace_matrix_proxy.
constexpr auto zip
A zip view.
Definition: zip.hpp:29
auto const move
A view that turns lvalue-references into rvalue-references.
Definition: move.hpp:74
Provides C++20 additions to the <iterator> header.
Adaptations of concepts from the Ranges TS.
Provides seqan3::detail::trace_iterator_banded.
Provides seqan3::views::zip.