INCLUDES HEADER: fixmath.hxx IN fixmath.hxx

USES INCLUDE HEADER: fixmath.hxx

constexpr bool isnan(const Arith::dual<T, U> &x);

constexpr bool isnan1(const Arith::dual<T, U> &x);

constexpr bool isnan(const Arith::dual<T, U> &x) {
  return isnan(x.val) || x.eps.isnan().any();

constexpr bool isnan1(const Arith::dual<T, U> &x) {
  return isnan1(x.val) || x.eps.isnan1().any();

#ifndef FIXMATH_HXX
#define FIXMATH_HXX
// Include this file BEFORE including <cctk.h>, best before including any other
// include files
#include <cmath>
namespace std {
template <typename T> constexpr T copysign1(T x, T y) {
  return std::copysign(x, y);
}
template <typename T> constexpr int fpclassify1(T x) {
  return std::fpclassify(x);
}
template <typename T> constexpr bool isfinite1(T x) { return std::isfinite(x); }
template <typename T> constexpr bool isinf1(T x) { return std::isinf(x); }
template <typename T> constexpr bool isnan1(T x) { return std::isnan(x); }
template <typename T> constexpr bool isnormal1(T x) { return std::isnormal(x); }
template <typename T> constexpr bool signbit1(T x) { return std::signbit(x); }
} // namespace std
#include <cctk.h>
#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit
#endif // #ifndef FIXMATH_HXX

#include <fixmath.hxx> // include this before <cctk.h>

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

isnan(const Arith::vect<T, D> &x) {

isnan1(const Arith::vect<T, D> &x) {

    r.elts[d] = isnan(x.elts[d]);

    r.elts[d] = isnan1(x.elts[d]);

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

template <typename T> T fmax1(T x, T y) {
  if (CCTK_isnan(CCTK_REAL(x)))

template <typename T> constexpr T fmax1(T x, T y) noexcept {
  if (isnan1(x))

  if (CCTK_isnan(CCTK_REAL(y)))

  if (isnan1(y))

template <typename T> T fmin1(T x, T y) {
  if (CCTK_isnan(CCTK_REAL(x)))

template <typename T> constexpr T fmin1(T x, T y) noexcept {
  if (isnan1(x))

  if (CCTK_isnan(CCTK_REAL(y)))

  if (isnan1(y))

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

#undef copysign
#undef fpclassify
#undef isfinite
#undef isinf
#undef isnan
#undef isnormal
#undef signbit

  assert(!isnan(etlen) && etlen >= 1.0e-12);

  assert(!isnan1(etlen) && etlen >= 1.0e-12);

      assert(!isnan(det(a)(b)));

      assert(!isnan1(det(a)(b)));

      assert(!isnan(dephi(a)(b)));

      assert(!isnan1(dephi(a)(b)));

      assert(!isnan(detheta(a)(b)));

      assert(!isnan1(detheta(a)(b)));

      assert(!isnan(der(a)(b)));

      assert(!isnan1(der(a)(b)));

  constexpr vec4<bool, dnup> isnan() const {
    return vec4<bool, dnup>(std::isnan(elts));

  constexpr vec4<bool, dnup> isnan1() const {
    return vec4<bool, dnup>(std::isnan1(elts));