#include "div.hxx"
namespace Arith {
template <typename T> constexpr bool test_div(T x, T y) {
  T d = div(x, y);
  T m = mod(x, y);
  return d * y + m == x && m >= 0 && m < abs(y);
}
// This function is compiled, but not executed. The tests are "run" at
// compile time. If this function compiles, the tests pass.
void TestDiv() {
  static_assert(test_div(0, 1));
  static_assert(test_div(0, 2));
  static_assert(test_div(0, 3));
  static_assert(test_div(1, 1));
  static_assert(test_div(1, 2));
  static_assert(test_div(1, 3));
  static_assert(test_div(2, 1));
  static_assert(test_div(2, 2));
  static_assert(test_div(2, 3));
  static_assert(test_div(3, 1));
  static_assert(test_div(3, 2));
  static_assert(test_div(3, 3));
  static_assert(test_div(-1, 1));
  static_assert(test_div(-1, 2));
  static_assert(test_div(-1, 3));
  static_assert(test_div(-2, 1));
  static_assert(test_div(-2, 2));
  static_assert(test_div(-2, 3));
  static_assert(test_div(-3, 1));
  static_assert(test_div(-3, 2));
  static_assert(test_div(-3, 3));
}
} // namespace Arith

#ifndef DIV_HXX
#define DIV_HXX
#include <cassert>
#include <cmath>
#include <cctk.h>
#ifdef CCTK_DEBUG
#define ARITH_INLINE
#else
#define ARITH_INLINE CCTK_ATTRIBUTE_ALWAYS_INLINE
#endif
namespace Arith {
using namespace std;
template <typename T> constexpr T div(T x, T y) ARITH_INLINE {
  // C++ division truncates; we want to round towards -infinity instead
  if (y < 0) {
    x = -x;
    y = -y;
  }
  return x >= 0 ? x / y : (x - y + 1) / y;
}
template <typename T> constexpr T mod(const T &x, const T &y) ARITH_INLINE {
  return x - div(x, y) * y;
}
} // namespace Arith
#undef ARITH_INLINE
#endif // #ifndef DIV_HXX

SRCS = dual.cxx mat.cxx rational.cxx sum.cxx vec.cxx

SRCS = div.cxx dual.cxx mat.cxx rational.cxx sum.cxx vec.cxx