math.hpp

Go to the documentation of this file.

#pragma once
 
#include "math/constants.hpp"
#include "math/differences.hpp"
#include "math/la.hpp"
 
namespace precice::math {

inline int sign(double number)
{
  if (greater(number, 0.0)) {
    return 1;
  } else if (greater(0.0, number)) {
    return -1;
  }
  return 0;
}

template <int iexp, typename T>
inline constexpr T pow_int(const T x)
{
  static_assert(iexp >= 0, "Exponent must be an integer greater or equal to zero.");
 
  if (iexp == 0)
    return static_cast<T>(1.);
  else
    // exponentiation by squaring
    return ((iexp % 2 == 1) ? x * pow_int<iexp / 2>(x * x) : pow_int<iexp / 2>(x * x));
}
 
} // namespace precice::math