21 namespace seqan3::detail
26 struct policy_max_error
30 search_cfg::max_error_total total{};
32 search_cfg::max_error_substitution substitution{};
34 search_cfg::max_error_insertion insertion{};
36 search_cfg::max_error_deletion deletion{};
39 bool only_max_error_total{
false};
41 bool has_max_error_total{
false};
46 policy_max_error() =
default;
47 policy_max_error(policy_max_error
const &) =
default;
48 policy_max_error(policy_max_error &&) =
default;
49 policy_max_error & operator=(policy_max_error
const &) =
default;
50 policy_max_error & operator=(policy_max_error &&) =
default;
51 ~policy_max_error() =
default;
62 template <
typename configuration_t>
64 requires is_type_specialisation_of_v<configuration_t, seqan3::configuration>
66 explicit policy_max_error(configuration_t
const & config)
68 using search_traits_t = search_traits<configuration_t>;
69 only_max_error_total = search_traits_t::only_max_error_total;
70 has_max_error_total = search_traits_t::has_max_error_total;
72 total = config.get_or(search_cfg::max_error_total{search_cfg::error_count{0}});
73 substitution = config.get_or(search_cfg::max_error_substitution{search_cfg::error_count{0}});
74 insertion = config.get_or(search_cfg::max_error_insertion{search_cfg::error_count{0}});
75 deletion = config.get_or(search_cfg::max_error_deletion{search_cfg::error_count{0}});
87 template <std::ranges::forward_range query_t>
88 auto max_error_counts(query_t && query)
90 detail::search_param errors{0, 0, 0, 0};
94 errors.total = to_error_count(total.error, query_size);
95 errors.substitution = to_error_count(substitution.error, query_size);
96 errors.insertion = to_error_count(insertion.error, query_size);
97 errors.deletion = to_error_count(deletion.error, query_size);
100 if (only_max_error_total)
101 errors.substitution = errors.insertion = errors.deletion = errors.total;
103 else if (!has_max_error_total)
104 errors.total = std::min<uint32_t>(255, errors.substitution + errors.insertion + errors.deletion);
108 if (errors.substitution > errors.total)
109 throw std::invalid_argument{
"The substitution error threshold is higher than the total error threshold."};
110 if (errors.insertion > errors.total)
111 throw std::invalid_argument{
"The insertion error threshold is higher than the total error threshold."};
112 if (errors.deletion > errors.total)
113 throw std::invalid_argument{
"The deletion error threshold is higher than the total error threshold."};
135 [[maybe_unused]]
size_t const query_size)
139 if constexpr (std::same_as<decltype(error), search_cfg::error_count>)
144 if (0.0 > error.get() || error.get() > 1.0)
148 uint8_t
const calculated_error_count =
std::clamp(error.get() * query_size, 0.0, 255.0);
149 return calculated_error_count;