Mpfr_value.hpp

Go to the documentation of this file.

/*******************************************************************************
 Causal Dynamical Triangulations in C++ using CGAL
 
 Copyright © 2026 Adam Getchell
 ******************************************************************************/

 
#ifndef CDT_PLUSPLUS_MPFR_VALUE_HPP
#define CDT_PLUSPLUS_MPFR_VALUE_HPP
 
#include <CGAL/Gmpfr.h>
 
#include <stdexcept>
#include <type_traits>
 
#include "Settings.hpp"
 
namespace cdt::mpfr_values
{
  using Value                         = CGAL::Gmpfr;

  inline constexpr auto rounding_mode = MPFR_RNDN;
 
  static_assert(std::is_nothrow_destructible_v<Value>,
                "MPFR values must release their resources without throwing.");

  inline constexpr auto precision =
      static_cast<Value::Precision_type>(PRECISION);

  [[nodiscard]] inline auto zero() -> Value { return Value{0, precision}; }

  [[nodiscard]] inline auto from_integer(long const value) -> Value
  {
    auto result = zero();
    mpfr_set_si(result.fr(), value, rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto from_long_double(long double const value) -> Value
  {
    auto result = zero();
    mpfr_set_ld(result.fr(), value, rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto from_decimal(char const* value) -> Value
  {
    if (value == nullptr)
    {
      throw std::invalid_argument("MPFR decimal value must not be null.");
    }
    auto result = zero();
    if (mpfr_set_str(result.fr(), value, 10, rounding_mode) != 0)
    {
      throw std::invalid_argument("Invalid decimal MPFR value.");
    }
    return result;
  }

  [[nodiscard]] inline auto pi() -> Value
  {
    auto result = zero();
    mpfr_const_pi(result.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto add(Value const& left, Value const& right) -> Value
  {
    auto result = zero();
    mpfr_add(result.fr(), left.fr(), right.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto subtract(Value const& left, Value const& right)
      -> Value
  {
    auto result = zero();
    mpfr_sub(result.fr(), left.fr(), right.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto multiply(Value const& left, Value const& right)
      -> Value
  {
    auto result = zero();
    mpfr_mul(result.fr(), left.fr(), right.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto divide(Value const& numerator,
                                   Value const& denominator) -> Value
  {
    auto result = zero();
    mpfr_div(result.fr(), numerator.fr(), denominator.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto square_root(Value const& value) -> Value
  {
    auto result = zero();
    mpfr_sqrt(result.fr(), value.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto inverse_hyperbolic_sine(Value const& value) -> Value
  {
    auto result = zero();
    mpfr_asinh(result.fr(), value.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto arc_cosine(Value const& value) -> Value
  {
    auto result = zero();
    mpfr_acos(result.fr(), value.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto negate(Value const& value) -> Value
  {
    auto result = zero();
    mpfr_neg(result.fr(), value.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto exponential(Value const& value) -> Value
  {
    auto result = zero();
    mpfr_exp(result.fr(), value.fr(), rounding_mode);
    return result;
  }

  [[nodiscard]] inline auto to_double(Value const& value) -> double
  { return mpfr_get_d(value.fr(), rounding_mode); }

  [[nodiscard]] inline auto to_long_double(Value const& value) -> long double
  { return mpfr_get_ld(value.fr(), rounding_mode); }
}  // namespace cdt::mpfr_values
 
#endif  // CDT_PLUSPLUS_MPFR_VALUE_HPP