math.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-3-Clause
// Copyright (C) 2026 Vallés Puig, Ramon
 
#pragma once
 
#include <cmath>
#include <type_traits>
 
#include "units/angular.hpp"
#include "units/dimensionless.hpp"
 
namespace qtty {
 
template <typename U> Ratio ratio_to(const Quantity<U> &a, const Quantity<U> &b) {
  return Ratio(a.value() / b.value());
}
 
// ── Transcendental helpers on dimensionless quantities ──────────────────────
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Ratio exp(const Quantity<U> &q) {
  return Ratio(std::exp(q.value()));
}
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Ratio ln(const Quantity<U> &q) {
  return Ratio(std::log(q.value()));
}
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Ratio powi(const Quantity<U> &q, int n) {
  return Ratio(std::pow(q.value(), static_cast<double>(n)));
}
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Ratio powf(const Quantity<U> &q, const Ratio &exponent) {
  return Ratio(std::pow(q.value(), exponent.value()));
}
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Radian asin_angle(const Quantity<U> &q) {
  return Radian(std::asin(q.value()));
}
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Radian acos_angle(const Quantity<U> &q) {
  return Radian(std::acos(q.value()));
}
 
template <typename U, std::enable_if_t<is_dimensionless_v<U>, int> = 0>
Radian atan_angle(const Quantity<U> &q) {
  return Radian(std::atan(q.value()));
}
 
} // namespace qtty