aminoacid_base.hpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2006-2025 Knut Reinert & Freie Universität Berlin
// SPDX-FileCopyrightText: 2016-2025 Knut Reinert & MPI für molekulare Genetik
// SPDX-License-Identifier: BSD-3-Clause

 
#pragma once
 
#include <seqan3/alphabet/alphabet_base.hpp>
#include <seqan3/alphabet/aminoacid/concept.hpp>
#include <seqan3/alphabet/detail/concept.hpp>
#include <seqan3/alphabet/detail/convert.hpp>
#include <seqan3/utility/char_operations/transform.hpp>
 
namespace seqan3
{

template <typename derived_type, auto size>
class aminoacid_base : public aminoacid_empty_base, public alphabet_base<derived_type, size, char>
{
private:
    using base_t = alphabet_base<derived_type, size, char>;

    friend base_t;

    constexpr aminoacid_base() noexcept = default;                                   
    constexpr aminoacid_base(aminoacid_base const &) noexcept = default;             
    constexpr aminoacid_base(aminoacid_base &&) noexcept = default;                  
    constexpr aminoacid_base & operator=(aminoacid_base const &) noexcept = default; 
    constexpr aminoacid_base & operator=(aminoacid_base &&) noexcept = default;      
    ~aminoacid_base() noexcept = default;                                            


    friend derived_type;
 
protected:
    // Import from base:
    using typename base_t::char_type;
    using typename base_t::rank_type;
 
public:
    using base_t::alphabet_size;
    using base_t::to_rank;

    // This constructor needs to be public, because constructor templates are not inherited otherwise
    template <typename other_aa_type>
        requires (!std::same_as<aminoacid_base, other_aa_type>) && (!std::same_as<derived_type, other_aa_type>)
              && aminoacid_alphabet<other_aa_type>
              && detail::convertable_to_through_char_representation<other_aa_type, derived_type>
    explicit constexpr aminoacid_base(other_aa_type const other) noexcept
    {
        if constexpr (is_constexpr_default_constructible_v<other_aa_type>
                      && detail::writable_constexpr_alphabet<other_aa_type>)
        {
            static_cast<derived_type &>(*this) =
                detail::convert_through_char_representation<other_aa_type, derived_type>[seqan3::to_rank(other)];
        }
        else
        {
            seqan3::assign_char_to(seqan3::to_char(other), static_cast<derived_type &>(*this));
        }
    }


    static constexpr bool char_is_valid(char_type const c) noexcept
    {
        return valid_char_table[static_cast<uint8_t>(c)];
    }
 
private:
    static constexpr std::array<bool, 256> valid_char_table{
        []() constexpr
        {
            std::array<bool, 256> ret{};
 
            ret.fill(false); // Default constructor does not initialise!
 
            // the original valid chars and their lower cases
            for (size_t rank = 0u; rank < derived_type::alphabet_size; ++rank)
            {
                uint8_t c = derived_type::rank_to_char(rank);
                ret[c] = true;
                ret[to_lower(c)] = true;
            }
 
            return ret;
        }()};
};
 
} // namespace seqan3