dna.hpp

Go to the documentation of this file.

#pragma once
 
#include <cuda/std/cstdint>
 
 
namespace cuddl::detail {
 
#ifdef __CUDACC__
__host__ __device__
#endif
inline constexpr uint8_t encode_base(char base) noexcept {
    auto const byte = static_cast<uint8_t>(base);
    auto const upper = static_cast<uint8_t>(byte & 0xDFu);
    auto const x = (byte >> 1u) & 3u;
    auto const valid =
        static_cast<uint8_t>((upper == 'A') | (upper == 'C') | (upper == 'G') | (upper == 'T'));
    auto const mask = static_cast<uint8_t>(0u - valid);
    return static_cast<uint8_t>((x & mask) | (0xFFu & ~mask));
}
 
inline constexpr uint32_t missing_cell = 0x00FFFFFFU;
 
}  // namespace cuddl::detail