kl_arbitrary_precision_root.h File Reference

#include <boost/math/tools/roots.hpp>
#include <limits>

Include dependency graph for kl_arbitrary_precision_root.h:

Go to the source code of this file.

Data Structures
struct	cbrt_functor< T >

Namespaces
	boost
	boost::math
	boost::math::tools

Functions
template<class F , class T >
T	boost::math::tools::newton_raphson_iterate (F f, T guess, T min, T max, int digits)
template<class F , class T >
T	boost::math::tools::newton_raphson_iterate (F f, T guess, T min, T max, int digits, boost::uintmax_t &max_iter)
template<class F , class T >
T	boost::math::tools::halley_iterate (F f, T guess, T min, T max, int digits)
template<class F , class T >
T	boost::math::tools::halley_iterate (F f, T guess, T min, T max, int digits, boost::uintmax_t &max_iter)
template<class F , class T >
T	boost::math::tools::schroeder_iterate (F f, T guess, T min, T max, int digits)
template<class F , class T >
T	boost::math::tools::schroeder_iterate (F f, T guess, T min, T max, int digits, boost::uintmax_t &max_iter)
template<class T >
T	kl_arbitrary_precision_root (T z, unsigned int precision)

Function Documentation

template<class T >

T kl_arbitrary_precision_root	(	T	z,
		unsigned int	precision
	)

Definition at line 183 of file kl_arbitrary_precision_root.h.

References boost::math::tools::halley_iterate().

{
   using namespace std;
   int exp;
   frexp(z, &exp);
   T min = ldexp(0.5, exp/3);
   T max = ldexp(2.0, exp/3);
   T guess = ldexp(1.0, exp/3);
   int digits = std::numeric_limits<T>::digits / 2;

    cout<<"The typeid for this datatype is : "<<typeid(z).name();
    
    string zType =typeid(z).name();
    int usingArb = zType.compare(typeid(boost::math::ntl::RR).name());
    if (usingArb==0)
    {
        cout<<"We're using arb precision"<<endl;
        digits = precision /2;
    }

   T ans =boost::math::tools::halley_iterate(cbrt_functor<T>(z), guess, min, max, digits);
   return ans;
   //return tools::halley_iterate(detail::cbrt_functor<T>(z), guess, min, max, digits);
}

Here is the call graph for this function: