numbers.hpp

Go to the documentation of this file.

#pragma once
 
#include <cmath>
#include <numbers>
 
namespace MeltPoolDG::numbers
{
  // @note We did not use std::numeric_limits<double>::lowest() on purpose.
  // @p invalid_double should serve as a default value for an optional parameter.
  // During parsing, the ParameterHandler converts the default value (e.g.
  // set to invalid_double) to a string via Patterns::Tools::to_string<double>()
  // and then back to a double value. This entails a loss in the number of decimal
  // places and thus round-off errors.
  constexpr double invalid_double = -1e100;
 
  // Check if a @p number is invalid, if @p number is smaller or equal to
  // numbers::invalid_double.
  inline bool
  is_invalid(const double &number)
  {
    return number <= invalid_double;
  }
 
  template <typename number>
  number
  compute_angle_in_radians(const number angle_deg)
  {
    return angle_deg * number(std::numbers::pi_v<number> / 180.0);
  }
 
  constexpr double sqrt_pi = std::numbers::pi_v<double> * std::numbers::inv_sqrtpi_v<double>;
} // namespace MeltPoolDG::numbers