math.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2024 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
#include <concepts>
 
namespace Ikarus {
namespace Impl {
  template <typename T>
  constexpr T sqrt_helper(T x, T lo, T hi) {
    if (lo == hi)
      return lo;
 
    const T mid = (lo + hi + 1) / 2;
 
    if (x / mid < mid)
      return sqrt_helper<T>(x, lo, mid - 1);
    else
      return sqrt_helper(x, mid, hi);
  }
 
} // namespace Impl
 
template <typename T>
requires std::integral<T>
constexpr T ct_sqrt(T x) {
  return Impl::sqrt_helper<T>(x, 0, x / 2 + 1);
}
} // namespace Ikarus