rational_horner1_2.hpp 1.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142
  1. // (C) Copyright John Maddock 2007.
  2. // Use, modification and distribution are subject to the
  3. // Boost Software License, Version 1.0. (See accompanying file
  4. // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
  5. //
  6. // This file is machine generated, do not edit by hand
  7. // Polynomial evaluation using Horners rule
  8. #ifndef BOOST_MATH_TOOLS_POLY_RAT_2_HPP
  9. #define BOOST_MATH_TOOLS_POLY_RAT_2_HPP
  10. namespace boost{ namespace math{ namespace tools{ namespace detail{
  11. template <class T, class U, class V>
  12. inline V evaluate_rational_c_imp(const T*, const U*, const V&, const mpl::int_<0>*) BOOST_MATH_NOEXCEPT(V)
  13. {
  14. return static_cast<V>(0);
  15. }
  16. template <class T, class U, class V>
  17. inline V evaluate_rational_c_imp(const T* a, const U* b, const V&, const mpl::int_<1>*) BOOST_MATH_NOEXCEPT(V)
  18. {
  19. return static_cast<V>(a[0]) / static_cast<V>(b[0]);
  20. }
  21. template <class T, class U, class V>
  22. inline V evaluate_rational_c_imp(const T* a, const U* b, const V& x, const mpl::int_<2>*) BOOST_MATH_NOEXCEPT(V)
  23. {
  24. if(x <= 1)
  25. return static_cast<V>((a[1] * x + a[0]) / (b[1] * x + b[0]));
  26. else
  27. {
  28. V z = 1 / x;
  29. return static_cast<V>((a[0] * z + a[1]) / (b[0] * z + b[1]));
  30. }
  31. }
  32. }}}} // namespaces
  33. #endif // include guard