26 namespace seqan3::detail
51 template <
typename score_t>
55 class score_matrix_single_column
59 using physical_column_t =
std::vector<score_t, aligned_allocator<score_t,
alignof(score_t)>>;
63 class matrix_iterator;
66 physical_column_t optimal_column{};
68 physical_column_t horizontal_column{};
70 virtual_column_t vertical_column{};
72 size_t number_of_columns{};
78 score_matrix_single_column() =
default;
79 score_matrix_single_column(score_matrix_single_column
const &) =
default;
80 score_matrix_single_column(score_matrix_single_column &&) =
default;
81 score_matrix_single_column & operator=(score_matrix_single_column
const &) =
default;
82 score_matrix_single_column & operator=(score_matrix_single_column &&) =
default;
83 ~score_matrix_single_column() =
default;
113 template <std::
integral column_index_t, std::
integral row_index_t>
114 void resize(column_index_type<column_index_t>
const number_of_columns,
115 row_index_type<row_index_t>
const number_of_rows,
116 score_t
const initial_value = score_t{})
118 this->number_of_columns = number_of_columns.get();
119 optimal_column.clear();
120 horizontal_column.clear();
121 optimal_column.resize(number_of_rows.get(), initial_value);
122 horizontal_column.resize(number_of_rows.get(), initial_value);
123 vertical_column =
views::repeat_n(initial_value, number_of_rows.get());
130 matrix_iterator
begin()
132 return matrix_iterator{*
this, 0u};
136 matrix_iterator
begin()
const =
delete;
139 matrix_iterator
end()
141 return matrix_iterator{*
this, number_of_columns};
145 matrix_iterator
end()
const =
delete;
160 template <
typename score_t>
161 class score_matrix_single_column<score_t>::matrix_iterator
166 using matrix_column_t = decltype(
views::zip(std::declval<physical_column_t &>(),
167 std::declval<physical_column_t &>(),
168 std::declval<virtual_column_t &>()));
174 using fwd_tuple_t = decltype(tpl);
175 return affine_cell_proxy<std::remove_cvref_t<fwd_tuple_t>>{std::forward<fwd_tuple_t>(tpl)};
179 score_matrix_single_column * host_ptr{
nullptr};
181 size_t current_column_id{};
188 using value_type = decltype(std::declval<matrix_column_t>() | transform_to_affine_cell);
190 using reference = value_type;
192 using pointer = void;
202 matrix_iterator() noexcept = default;
203 matrix_iterator(matrix_iterator const &) noexcept = default;
204 matrix_iterator(matrix_iterator &&) noexcept = default;
205 matrix_iterator & operator=(matrix_iterator const &) noexcept = default;
206 matrix_iterator & operator=(matrix_iterator &&) noexcept = default;
207 ~matrix_iterator() = default;
214 explicit matrix_iterator(score_matrix_single_column & host_matrix,
size_t const initial_column_id) noexcept :
216 current_column_id{initial_column_id}
224 reference operator*()
const
226 return views::zip(host_ptr->optimal_column, host_ptr->horizontal_column, host_ptr->vertical_column)
227 | transform_to_affine_cell;
235 matrix_iterator & operator++()
252 friend bool operator==(matrix_iterator
const & lhs, matrix_iterator
const & rhs) noexcept
254 return lhs.current_column_id == rhs.current_column_id;
258 friend bool operator!=(matrix_iterator
const & lhs, matrix_iterator
const & rhs) noexcept
260 return !(lhs == rhs);
Provides seqan3::detail::affine_cell_proxy.
Provides concepts for core language types and relations that don't have concepts in C++20 (yet).
decltype(detail::transform< trait_t >(list_t{})) transform
Apply a transformation trait to every type in the list and return a seqan3::type_list of the results.
Definition: traits.hpp:471
constexpr auto zip
A zip view.
Definition: zip.hpp:29
constexpr auto repeat_n
A view factory that repeats a given value n times.
Definition: repeat_n.hpp:95
A type that satisfies std::is_arithmetic_v<t>.
Provides seqan3::detail::matrix_index, seqan3::detail::matrix_coordinate and associated strong types.
Adaptations of concepts from the Ranges TS.
Provides seqan3::aligned_allocator.
Provides seqan3::simd::simd_concept.
Provides seqan3::views::repeat_n.
Provides seqan3::views::zip.