Calculate entropy using modern C++

entropy.cpp

#include <cassert>
#include <cmath>
#include <concepts>

namespace phd {

  template<typename T>
  concept numeric = std::integral<T> || std::floating_point<T>;

  template<numeric Arg, std::integral Base>
  auto log(Arg argument, Base base = 10)
  {
      assert(base > 0);

      return std::log(argument) / std::log(base);
  }
  
}

template<std::integral Choice, std::integral... Uncertainty>
auto get_entropy(Choice choices, Uncertainty... uncertainties)
{
    assert(choices > 0);
    assert((... && (uncertainties >= 1)));

    return (... + (((long double) uncertainties / choices) *
                    phd::log((long double) choices / uncertainties, 2)));
}