align_result_selector.hpp

Go to the documentation of this file.

// -----------------------------------------------------------------------------------------------------
// Copyright (c) 2006-2020, Knut Reinert & Freie Universität Berlin
// Copyright (c) 2016-2020, Knut Reinert & MPI für molekulare Genetik
// This file may be used, modified and/or redistributed under the terms of the 3-clause BSD-License
// shipped with this file and also available at: https://github.com/seqan/seqan3/blob/master/LICENSE.md
// -----------------------------------------------------------------------------------------------------
 
#pragma once
 
#include <optional>
#include <type_traits>
 
#include <seqan3/alignment/configuration/align_config_debug.hpp>
#include <seqan3/alignment/configuration/align_config_result.hpp>
#include <seqan3/alignment/matrix/alignment_coordinate.hpp>
#include <seqan3/alignment/matrix/detail/two_dimensional_matrix.hpp>
#include <seqan3/alignment/matrix/trace_directions.hpp>
#include <seqan3/alignment/pairwise/detail/type_traits.hpp>
#include <seqan3/alphabet/gap/gapped.hpp>
#include <seqan3/core/algorithm/configuration.hpp>
#include <seqan3/core/type_traits/basic.hpp>
#include <seqan3/core/type_traits/range.hpp>
#include <seqan3/core/type_traits/template_inspection.hpp>
#include <seqan3/core/type_traits/transformation_trait_or.hpp>
#include <seqan3/core/type_list/traits.hpp>
#include <seqan3/core/type_list/type_list.hpp>
#include <seqan3/range/decorator/gap_decorator.hpp>
#include <seqan3/range/views/type_reduce.hpp>
#include <seqan3/std/ranges>
 
namespace seqan3::detail
{
 
template <typename first_t, typename second_t>
struct alignment_type;
 
template <typename first_t, typename second_t>
    requires std::ranges::random_access_range<first_t> &&
             std::ranges::sized_range<first_t> &&
             std::ranges::random_access_range<second_t> &&
             std::ranges::sized_range<second_t>
struct alignment_type<first_t, second_t>
{
    using type = std::tuple<gap_decorator<type_reduce_view<first_t &>>, gap_decorator<type_reduce_view<second_t &>>>;
};
 
template <std::ranges::forward_range first_range_t,
          std::ranges::forward_range second_range_t,
          typename configuration_t>
    requires is_type_specialisation_of_v<remove_cvref_t<configuration_t>, configuration>
struct align_result_selector
{
private:
    using score_type = typename alignment_configuration_traits<configuration_t>::original_score_t;
 
    static constexpr auto select()
    {
        static_assert(configuration_t::template exists<align_cfg::result>());
 
        if constexpr (configuration_t::template exists<align_cfg::result<with_back_coordinate_type, score_type>>())
        {
            return alignment_result_value_type<uint32_t,
                                               score_type,
                                               alignment_coordinate>{};
        }
        else if constexpr (configuration_t::template exists<align_cfg::result<with_front_coordinate_type, score_type>>())
        {
            return alignment_result_value_type<uint32_t,
                                               score_type,
                                               alignment_coordinate,
                                               alignment_coordinate>{};
        }
        else if constexpr (configuration_t::template exists<align_cfg::result<with_alignment_type, score_type>>())
        {
            // Due to an error with gcc8 we define these types beforehand.
            using first_gapped_seq_type = gapped<value_type_t<first_range_t>>;
            using second_gapped_seq_type = gapped<value_type_t<second_range_t>>;
 
            // We use vectors of gapped sequence if the gap decorator cannot be used.
            using fallback_t = std::tuple<std::vector<first_gapped_seq_type>, std::vector<second_gapped_seq_type>>;
 
            // If the ranges are RandomAccess and Sized we can use the Gap Decorator, otherwise fallback_t.
            using decorator_t = alignment_type<first_range_t, second_range_t>;
 
            return alignment_result_value_type<uint32_t,
                                               score_type,
                                               alignment_coordinate,
                                               alignment_coordinate,
                                               detail::transformation_trait_or_t<decorator_t, fallback_t>>{};
        }
        else
        {
            return alignment_result_value_type<uint32_t, score_type>{};
        }
    }
 
    template <typename alignment_result_value_t>
    static constexpr auto augment_if_debug(alignment_result_value_t)
    {
        if constexpr (configuration_t::template exists<detail::debug_mode>())
        {
            using as_type_list = transfer_template_args_onto_t<alignment_result_value_t, type_list>;
            using score_matrix_t = two_dimensional_matrix<std::optional<score_type>,
                                                          std::allocator<std::optional<score_type>>,
                                                          matrix_major_order::column>;
            using trace_matrix_t = two_dimensional_matrix<std::optional<trace_directions>,
                                                          std::allocator<std::optional<trace_directions>>,
                                                          matrix_major_order::column>;
 
            if constexpr (configuration_t::template exists<align_cfg::result<with_alignment_type>>())
            {
                using with_score_t = list_traits::replace_at<score_matrix_t, 5, as_type_list>;
                return transfer_template_args_onto_t<list_traits::replace_at<trace_matrix_t, 6, with_score_t>,
                                                     alignment_result_value_type>{};
            }
            else
            {
                return transfer_template_args_onto_t<list_traits::replace_at<score_matrix_t, 5, as_type_list>,
                                                     alignment_result_value_type>{};
            }
        }
        else  // Return as is.
        {
            return alignment_result_value_t{};
        }
    }
 
public:
    using type = decltype(augment_if_debug(select()));
};
 
} // namespace seqan3::detail