https://gcc.gnu.org/bugzilla/show_bug.cgi?id=102971

--- Comment #2 from Andrew Pinski <pinskia at gcc dot gnu.org> ---
#include <cstdint>

std::uint64_t umul_naive(std::uint64_t a,std::uint64_t b,std::uint64_t& high)
noexcept
{
        std::uint32_t a0(static_cast<std::uint32_t>(a));
        std::uint32_t a1(static_cast<std::uint32_t>(a>>32));
        std::uint32_t b0(static_cast<std::uint32_t>(b));
        std::uint32_t b1(static_cast<std::uint32_t>(b>>32));
        std::uint64_t c00(static_cast<std::uint64_t>(a0)*b0);
        std::uint64_t c01(static_cast<std::uint64_t>(a0)*b1);
        std::uint64_t c10(static_cast<std::uint64_t>(a1)*b0);
        std::uint64_t c11(static_cast<std::uint64_t>(a1)*b1);

        std::uint64_t d0{static_cast<std::uint32_t>(c00)};
        std::uint64_t c00_high{c00>>32};
        std::uint64_t c01_low{static_cast<std::uint32_t>(c01)};
        std::uint64_t c01_high{c01>>32};
        std::uint64_t c10_low{static_cast<std::uint32_t>(c10)};
        std::uint64_t c10_high{c10>>32};

        std::uint64_t d1{c00_high+c01_low+c10_low};
        std::uint64_t d2{(d1>>32)+c10_high+c01_high};
        std::uint64_t d3{(d2>>32)+c11};
        high=d2|(d3<<32);
        return d0|(d1<<32);
}

Reply via email to