si_math.hpp

Go to the documentation of this file.

#pragma once
#include "si_internal.hpp"
#include "si_prefix.hpp"
#include <cmath>
namespace Si {
namespace Internal {
template <typename T> struct std_pow {
  T operator()(const T &base, const T &power) const {
    return std::pow(base, power);
  }
};
} // namespace Internal
namespace Math {
template <typename Power, typename T, typename Function = Internal::std_pow<T>,
          typename... Powers>
Internal::Unit<T, std::ratio_multiply<Power, Powers>...>
pow(const Internal::Unit<T, Powers...> &base) {
  return Internal::Unit<T, std::ratio_multiply<Power, Powers>...>(
      Function()(static_cast<T>(base),
                 static_cast<T>(Power::num) / static_cast<T>(Power::den)));
}
template <
    typename Power, typename Unit,
    typename Function = Internal::std_pow<decltype(+std::declval<Unit>())>,
    typename PrefixRatio>
Internal::GenericPrefix<
    typename std::ratio<static_cast<std::intmax_t>(std::pow(PrefixRatio::num,
                                                            Power::num)),
                        static_cast<std::intmax_t>(std::pow(PrefixRatio::den,
                                                            Power::num))>::type,
    decltype(pow<Power>(std::declval<Unit>()))>
pow(const Internal::GenericPrefix<PrefixRatio, Unit> &base) {
  static_assert(
      Power::den == 1,
      "std::ratio for power in case of unit with prefix has to be integer");
  return Internal::GenericPrefix<
      typename std::ratio<static_cast<std::intmax_t>(
                              std::pow(PrefixRatio::num, Power::num)),
                          static_cast<std::intmax_t>(
                              std::pow(PrefixRatio::den, Power::num))>::type,
      decltype(pow<Power>(std::declval<Unit>()))>(
      +pow<Power, decltype(+std::declval<Unit>()), Function>(base.raw()));
}
} // namespace Math
} // namespace Si