comparison.cuh

Go to the documentation of this file.

#pragma once
 
#include <cuda/std/cstdint>
 
#include <cuddl/detail/register.cuh>
#include <cuddl/pairwise_counts.cuh>
 
namespace cuddl::detail {
 
__device__ inline void
classify(pairwise_counts& target, uint32_t this_score, uint32_t other_score) noexcept {
    auto const this_zero = static_cast<uint32_t>(this_score == 0U);
    auto const lt = static_cast<uint32_t>(this_score < other_score);
    auto const gt = static_cast<uint32_t>(this_score > other_score);
    auto const eq = 1U - lt - gt;
    target.both_empty += eq & this_zero;
    target.lower += lt;
    target.higher += gt;
    target.equal += eq & (1U - this_zero);
}
 
}  // namespace cuddl::detail