boost-react-native-bundle

<html> <head> <meta http-equiv="Content-Type" content="text/html; charset=US-ASCII"> <title>Continued Fraction Evaluation</title> <link rel="stylesheet" href="../../math.css" type="text/css"> <meta name="generator" content="DocBook XSL Stylesheets V1.78.1"> <link rel="home" href="../../index.html" title="Math Toolkit 2.1.0"> <link rel="up" href="../internals1.html" title="Utilities & internals"> <link rel="prev" href="series_evaluation.html" title="Series Evaluation"> <link rel="next" href="rational.html" title="Polynomial and Rational Function Evaluation"> </head> <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF"> <table cellpadding="2" width="100%"><tr> <td valign="top"><img alt="Boost C++ Libraries" width="277" height="86" src="../../../../../../boost.png"></td> <td align="center"><a href="../../../../../../index.html">Home</a></td> <td align="center"><a href="../../../../../../libs/libraries.htm">Libraries</a></td> <td align="center"><a href="http://www.boost.org/users/people.html">People</a></td> <td align="center"><a href="http://www.boost.org/users/faq.html">FAQ</a></td> <td align="center"><a href="../../../../../../more/index.htm">More</a></td> </tr></table> <hr> <div class="spirit-nav"> <a accesskey="p" href="series_evaluation.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../internals1.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="rational.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a> </div> <div class="section"> <div class="titlepage"><div><div><h3 class="title"> <a name="math_toolkit.internals1.cf"></a><a class="link" href="cf.html" title="Continued Fraction Evaluation">Continued Fraction Evaluation</a> </h3></div></div></div> <h5> <a name="math_toolkit.internals1.cf.h0"></a> <a name="math_toolkit.internals1.cf.synopsis"></a><a class="link" href="cf.html#math_toolkit.internals1.cf.synopsis">Synopsis</a> </h5> <pre class="programlisting">#include <boost/math/tools/fraction.hpp> </pre> <pre class="programlisting">namespace boost{ namespace math{ namespace tools{ template <class Gen, class U> typename detail::fraction_traits<Gen>::result_type continued_fraction_b(Gen& g, const U& tolerance, boost::uintmax_t& max_terms) template <class Gen, class U> typename detail::fraction_traits<Gen>::result_type continued_fraction_b(Gen& g, const U& tolerance) template <class Gen, class U> typename detail::fraction_traits<Gen>::result_type continued_fraction_a(Gen& g, const U& tolerance, boost::uintmax_t& max_terms) template <class Gen, class U> typename detail::fraction_traits<Gen>::result_type continued_fraction_a(Gen& g, const U& tolerance) // // These interfaces are present for legacy reasons, and are now deprecated: // template <class Gen> typename detail::fraction_traits<Gen>::result_type continued_fraction_b(Gen& g, int bits); template <class Gen> typename detail::fraction_traits<Gen>::result_type continued_fraction_b(Gen& g, int bits, boost::uintmax_t& max_terms); template <class Gen> typename detail::fraction_traits<Gen>::result_type continued_fraction_a(Gen& g, int bits); template <class Gen> typename detail::fraction_traits<Gen>::result_type continued_fraction_a(Gen& g, int bits, boost::uintmax_t& max_terms); }}} // namespaces </pre> <h5> <a name="math_toolkit.internals1.cf.h1"></a> <a name="math_toolkit.internals1.cf.description"></a><a class="link" href="cf.html#math_toolkit.internals1.cf.description">Description</a> </h5> <a href="http://en.wikipedia.org/wiki/Continued_fraction" target="_top">Continued fractions are a common method of approximation. </a> These functions all evaluate the continued fraction described by the generator type argument. The functions with an "_a" suffix evaluate the fraction: <img src="../../../equations/fraction2.png"> and those with a "_b" suffix evaluate the fraction: <img src="../../../equations/fraction1.png"> This latter form is somewhat more natural in that it corresponds with the usual definition of a continued fraction, but note that the first a value returned by the generator is discarded. Further, often the first a and b values in a continued fraction have different defining equations to the remaining terms, which may make the "_a" suffixed form more appropriate. The generator type should be a function object which supports the following operations: <div class="informaltable"><table class="table"> <colgroup> <col> <col> </colgroup> <thead><tr> <th> Expression </th> <th> Description </th> </tr></thead> <tbody> <tr> <td> Gen::result_type </td> <td> The type that is the result of invoking operator(). This can be either an arithmetic type, or a std::pair<> of arithmetic types. </td> </tr> <tr> <td> g() </td> <td> Returns an object of type Gen::result_type. Each time this operator is called then the next pair of a and b values is returned. Or, if result_type is an arithmetic type, then the next b value is returned and all the a values are assumed to 1. </td> </tr> </tbody> </table></div> In all the continued fraction evaluation functions the tolerance parameter is the precision desired in the result, evaluation of the fraction will continue until the last term evaluated leaves the relative error in the result less than tolerance. The deprecated interfaces take a number of digits precision here, internally they just convert this to a tolerance and forward call. If the optional max_terms parameter is specified then no more than max_terms calls to the generator will be made, and on output, max_terms will be set to actual number of calls made. This facility is particularly useful when profiling a continued fraction for convergence. <h5> <a name="math_toolkit.internals1.cf.h2"></a> <a name="math_toolkit.internals1.cf.implementation"></a><a class="link" href="cf.html#math_toolkit.internals1.cf.implementation">Implementation</a> </h5> Internally these algorithms all use the modified Lentz algorithm: refer to Numeric Recipes in C++, W. H. Press et all, chapter 5, (especially 5.2 Evaluation of continued fractions, p 175 - 179) for more information, also Lentz, W.J. 1976, Applied Optics, vol. 15, pp. 668-671. <h5> <a name="math_toolkit.internals1.cf.h3"></a> <a name="math_toolkit.internals1.cf.examples"></a><a class="link" href="cf.html#math_toolkit.internals1.cf.examples">Examples</a> </h5> The <a href="http://en.wikipedia.org/wiki/Golden_ratio" target="_top">golden ratio phi = 1.618033989...</a> can be computed from the simplest continued fraction of all: <img src="../../../equations/fraction3.png"> We begin by defining a generator function: <pre class="programlisting">template <class T> struct golden_ratio_fraction { typedef T result_type; result_type operator() { return 1; } }; </pre> The golden ratio can then be computed to double precision using: <pre class="programlisting">continued_fraction_a( golden_ratio_fraction<double>(), std::numeric_limits<double>::epsilon()); </pre> It's more usual though to have to define both the a's and the b's when evaluating special functions by continued fractions, for example the tan function is defined by: <img src="../../../equations/fraction4.png"> So its generator object would look like: <pre class="programlisting">template <class T> struct tan_fraction { private: T a, b; public: tan_fraction(T v) : a(-v*v), b(-1) {} typedef std::pair<T,T> result_type; std::pair<T,T> operator()() { b += 2; return std::make_pair(a, b); } }; </pre> Notice that if the continuant is subtracted from the b terms, as is the case here, then all the a terms returned by the generator will be negative. The tangent function can now be evaluated using: <pre class="programlisting">template <class T> T tan(T a) { tan_fraction<T> fract(a); return a / continued_fraction_b(fract, std::numeric_limits<T>::epsilon()); } </pre> Notice that this time we're using the "_b" suffixed version to evaluate the fraction: we're removing the leading a term during fraction evaluation as it's different from all the others. </div> <table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr> <td align="left"></td> <td align="right"><div class="copyright-footer">Copyright © 2006-2010, 2012-2014 Nikhar Agrawal, Anton Bikineev, Paul A. Bristow, Marco Guazzone, Christopher Kormanyos, Hubert Holin, Bruno Lalande, John Maddock, Johan Råde, Gautam Sewani, Benjamin Sobotta, Thijs van den Berg, Daryle Walker and Xiaogang Zhang Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at <a href="http://www.boost.org/LICENSE_1_0.txt" target="_top">http://www.boost.org/LICENSE_1_0.txt</a>) </div></td> </tr></table> <hr> <div class="spirit-nav"> <a accesskey="p" href="series_evaluation.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../internals1.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="rational.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a> </div> </body> </html>