float128.html 19 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190
  1. <html>
  2. <head>
  3. <meta http-equiv="Content-Type" content="text/html; charset=US-ASCII">
  4. <title>float128</title>
  5. <link rel="stylesheet" href="../../../multiprecision.css" type="text/css">
  6. <meta name="generator" content="DocBook XSL Stylesheets V1.79.1">
  7. <link rel="home" href="../../../index.html" title="Chapter&#160;1.&#160;Boost.Multiprecision">
  8. <link rel="up" href="../floats.html" title="floating-point Numbers">
  9. <link rel="prev" href="mpfr_float.html" title="mpfr_float">
  10. <link rel="next" href="fp_eg.html" title="Examples">
  11. </head>
  12. <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
  13. <table cellpadding="2" width="100%"><tr>
  14. <td valign="top"><img alt="Boost C++ Libraries" width="277" height="86" src="../../../../../../../boost.png"></td>
  15. <td align="center"><a href="../../../../../../../index.html">Home</a></td>
  16. <td align="center"><a href="../../../../../../../libs/libraries.htm">Libraries</a></td>
  17. <td align="center"><a href="http://www.boost.org/users/people.html">People</a></td>
  18. <td align="center"><a href="http://www.boost.org/users/faq.html">FAQ</a></td>
  19. <td align="center"><a href="../../../../../../../more/index.htm">More</a></td>
  20. </tr></table>
  21. <hr>
  22. <div class="spirit-nav">
  23. <a accesskey="p" href="mpfr_float.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../floats.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="fp_eg.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
  24. </div>
  25. <div class="section">
  26. <div class="titlepage"><div><div><h4 class="title">
  27. <a name="boost_multiprecision.tut.floats.float128"></a><a class="link" href="float128.html" title="float128">float128</a>
  28. </h4></div></div></div>
  29. <p>
  30. <code class="computeroutput"><span class="preprocessor">#include</span> <span class="special">&lt;</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">multiprecision</span><span class="special">/</span><span class="identifier">float128</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">&gt;</span></code>
  31. </p>
  32. <pre class="programlisting"><span class="keyword">namespace</span> <span class="identifier">boost</span><span class="special">{</span> <span class="keyword">namespace</span> <span class="identifier">multiprecision</span><span class="special">{</span>
  33. <span class="keyword">class</span> <span class="identifier">float128_backend</span><span class="special">;</span>
  34. <span class="keyword">typedef</span> <span class="identifier">number</span><span class="special">&lt;</span><span class="identifier">float128_backend</span><span class="special">,</span> <span class="identifier">et_off</span><span class="special">&gt;</span> <span class="identifier">float128</span><span class="special">;</span>
  35. <span class="special">}}</span> <span class="comment">// namespaces</span>
  36. </pre>
  37. <p>
  38. The <code class="computeroutput"><span class="identifier">float128</span></code> number type
  39. is a very thin wrapper around GCC's <code class="computeroutput"><a class="link" href="float128.html" title="float128">float128</a></code>
  40. or Intel's <code class="computeroutput"><span class="identifier">_Quad</span></code> data types
  41. and provides an real-number type that is a drop-in replacement for the
  42. native C++ floating-point types, but with a 113 bit mantissa, and compatible
  43. with FORTRAN's 128-bit QUAD real.
  44. </p>
  45. <p>
  46. All the usual standard library and <code class="computeroutput"><span class="identifier">numeric_limits</span></code>
  47. support are available, performance should be equivalent to the underlying
  48. native types: for example the LINPACK benchmarks for GCC's <code class="computeroutput"><a class="link" href="float128.html" title="float128">float128</a></code>
  49. and <code class="computeroutput"><span class="identifier">boost</span><span class="special">::</span><span class="identifier">multiprecision</span><span class="special">::</span><span class="identifier">float128</span></code> both achieved 5.6 MFLOPS<a href="#ftn.boost_multiprecision.tut.floats.float128.f0" class="footnote" name="boost_multiprecision.tut.floats.float128.f0"><sup class="footnote">[3]</sup></a>.
  50. </p>
  51. <p>
  52. As well as the usual conversions from arithmetic and string types, instances
  53. of <code class="computeroutput"><span class="identifier">float128</span></code> are copy constructible
  54. and assignable from GCC's <code class="computeroutput"><a class="link" href="float128.html" title="float128">float128</a></code>
  55. and Intel's <code class="computeroutput"><span class="identifier">_Quad</span></code> data
  56. types.
  57. </p>
  58. <p>
  59. It's also possible to access the underlying <code class="computeroutput"><a class="link" href="float128.html" title="float128">float128</a></code>
  60. or <code class="computeroutput"><span class="identifier">_Quad</span></code> type via the
  61. <code class="computeroutput"><span class="identifier">data</span><span class="special">()</span></code>
  62. member function of <code class="computeroutput"><span class="identifier">float128_backend</span></code>.
  63. </p>
  64. <p>
  65. Things you should know when using this type:
  66. </p>
  67. <div class="itemizedlist"><ul class="itemizedlist" style="list-style-type: disc; ">
  68. <li class="listitem">
  69. Default constructed <code class="computeroutput"><span class="identifier">float128</span></code>s
  70. have the value zero.
  71. </li>
  72. <li class="listitem">
  73. This backend supports rvalue-references and is move-aware, making instantiations
  74. of <code class="computeroutput"><span class="identifier">number</span></code> on this backend
  75. move aware.
  76. </li>
  77. <li class="listitem">
  78. This type is fully <code class="computeroutput"><span class="keyword">constexpr</span></code>
  79. aware - basic constexpr arithmetic is supported from C++14 and onwards,
  80. comparisons, plus the functions <code class="computeroutput"><span class="identifier">fabs</span></code>,
  81. <code class="computeroutput"><span class="identifier">abs</span></code>, <code class="computeroutput"><span class="identifier">fpclassify</span></code>, <code class="computeroutput"><span class="identifier">isnormal</span></code>,
  82. <code class="computeroutput"><span class="identifier">isfinite</span></code>, <code class="computeroutput"><span class="identifier">isinf</span></code> and <code class="computeroutput"><span class="identifier">isnan</span></code>
  83. are also supported if either the compiler implements C++20's <code class="computeroutput"><span class="identifier">std</span><span class="special">::</span><span class="identifier">is_constant_evaluated</span><span class="special">()</span></code>,
  84. or if the compiler is GCC.
  85. </li>
  86. <li class="listitem">
  87. It is not possible to round-trip objects of this type to and from a
  88. string and get back exactly the same value when compiled with Intel's
  89. C++ compiler and using <code class="computeroutput"><span class="identifier">_Quad</span></code>
  90. as the underlying type: this is a current limitation of our code. Round
  91. tripping when using <code class="computeroutput"><a class="link" href="float128.html" title="float128">float128</a></code>
  92. as the underlying type is possible (both for GCC and Intel).
  93. </li>
  94. <li class="listitem">
  95. Conversion from a string results in a <code class="computeroutput"><span class="identifier">std</span><span class="special">::</span><span class="identifier">runtime_error</span></code>
  96. being thrown if the string can not be interpreted as a valid floating-point
  97. number.
  98. </li>
  99. <li class="listitem">
  100. Division by zero results in an infinity being produced.
  101. </li>
  102. <li class="listitem">
  103. Type <code class="computeroutput"><span class="identifier">float128</span></code> can be
  104. used as a literal type (constexpr support).
  105. </li>
  106. <li class="listitem">
  107. Type <code class="computeroutput"><span class="identifier">float128</span></code> can be
  108. used for full <code class="computeroutput"><span class="keyword">constexpr</span></code>
  109. arithmetic from C++14 and later with GCC. The functions <code class="computeroutput"><span class="identifier">abs</span></code>, <code class="computeroutput"><span class="identifier">fabs</span></code>,
  110. <code class="computeroutput"><span class="identifier">fpclassify</span></code>, <code class="computeroutput"><span class="identifier">isnan</span></code>, <code class="computeroutput"><span class="identifier">isinf</span></code>,
  111. <code class="computeroutput"><span class="identifier">isfinite</span></code> and <code class="computeroutput"><span class="identifier">isnormal</span></code> are also <code class="computeroutput"><span class="keyword">constexpr</span></code>,
  112. but the transcendental functions are not.
  113. </li>
  114. <li class="listitem">
  115. When using the Intel compiler, the underlying type defaults to <code class="computeroutput"><a class="link" href="float128.html" title="float128">float128</a></code>
  116. if it's available and <code class="computeroutput"><span class="identifier">_Quad</span></code>
  117. if not. You can override the default by defining either <code class="computeroutput"><span class="identifier">BOOST_MP_USE_FLOAT128</span></code> or <code class="computeroutput"><span class="identifier">BOOST_MP_USE_QUAD</span></code>.
  118. </li>
  119. <li class="listitem">
  120. When the underlying type is Intel's <code class="computeroutput"><span class="identifier">_Quad</span></code>
  121. type, the code must be compiled with the compiler option <code class="computeroutput"><span class="special">-</span><span class="identifier">Qoption</span><span class="special">,</span><span class="identifier">cpp</span><span class="special">,--</span><span class="identifier">extended_float_type</span></code>.
  122. </li>
  123. <li class="listitem">
  124. When compiling with <code class="computeroutput"><span class="identifier">gcc</span></code>,
  125. you need to use the flag <code class="computeroutput"><span class="special">--</span><span class="identifier">std</span><span class="special">=</span><span class="identifier">gnu</span><span class="special">++</span><span class="number">11</span><span class="special">/</span><span class="number">14</span><span class="special">/</span><span class="number">17</span></code>,
  126. as the suffix 'Q' is a GNU extension. Compilation fails with the flag
  127. <code class="computeroutput"><span class="special">--</span><span class="identifier">std</span><span class="special">=</span><span class="identifier">c</span><span class="special">++</span><span class="number">11</span><span class="special">/</span><span class="number">14</span><span class="special">/</span><span class="number">17</span></code> unless you also use <code class="computeroutput"><span class="special">-</span><span class="identifier">fext</span><span class="special">-</span><span class="identifier">numeric</span><span class="special">-</span><span class="identifier">literals</span></code>.
  128. </li>
  129. </ul></div>
  130. <h6>
  131. <a name="boost_multiprecision.tut.floats.float128.h0"></a>
  132. <span class="phrase"><a name="boost_multiprecision.tut.floats.float128.float128_example"></a></span><a class="link" href="float128.html#boost_multiprecision.tut.floats.float128.float128_example">float128
  133. example:</a>
  134. </h6>
  135. <pre class="programlisting"><span class="preprocessor">#include</span> <span class="special">&lt;</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">multiprecision</span><span class="special">/</span><span class="identifier">float128</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">&gt;</span>
  136. <span class="preprocessor">#include</span> <span class="special">&lt;</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">math</span><span class="special">/</span><span class="identifier">special_functions</span><span class="special">/</span><span class="identifier">gamma</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">&gt;</span>
  137. <span class="preprocessor">#include</span> <span class="special">&lt;</span><span class="identifier">iostream</span><span class="special">&gt;</span>
  138. <span class="keyword">int</span> <span class="identifier">main</span><span class="special">()</span>
  139. <span class="special">{</span>
  140. <span class="keyword">using</span> <span class="keyword">namespace</span> <span class="identifier">boost</span><span class="special">::</span><span class="identifier">multiprecision</span><span class="special">;</span>
  141. <span class="comment">// Operations at 128-bit precision and full numeric_limits support:</span>
  142. <span class="identifier">float128</span> <span class="identifier">b</span> <span class="special">=</span> <span class="number">2</span><span class="special">;</span>
  143. <span class="comment">// There are 113-bits of precision:</span>
  144. <span class="identifier">std</span><span class="special">::</span><span class="identifier">cout</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">numeric_limits</span><span class="special">&lt;</span><span class="identifier">float128</span><span class="special">&gt;::</span><span class="identifier">digits</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">endl</span><span class="special">;</span>
  145. <span class="comment">// Or 34 decimal places:</span>
  146. <span class="identifier">std</span><span class="special">::</span><span class="identifier">cout</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">numeric_limits</span><span class="special">&lt;</span><span class="identifier">float128</span><span class="special">&gt;::</span><span class="identifier">digits10</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">endl</span><span class="special">;</span>
  147. <span class="comment">// We can use any C++ std lib function, lets print all the digits as well:</span>
  148. <span class="identifier">std</span><span class="special">::</span><span class="identifier">cout</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">setprecision</span><span class="special">(</span><span class="identifier">std</span><span class="special">::</span><span class="identifier">numeric_limits</span><span class="special">&lt;</span><span class="identifier">float128</span><span class="special">&gt;::</span><span class="identifier">max_digits10</span><span class="special">)</span>
  149. <span class="special">&lt;&lt;</span> <span class="identifier">log</span><span class="special">(</span><span class="identifier">b</span><span class="special">)</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">endl</span><span class="special">;</span> <span class="comment">// print log(2) = 0.693147180559945309417232121458176575</span>
  150. <span class="comment">// We can also use any function from Boost.Math:</span>
  151. <span class="identifier">std</span><span class="special">::</span><span class="identifier">cout</span> <span class="special">&lt;&lt;</span> <span class="identifier">boost</span><span class="special">::</span><span class="identifier">math</span><span class="special">::</span><span class="identifier">tgamma</span><span class="special">(</span><span class="identifier">b</span><span class="special">)</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">endl</span><span class="special">;</span>
  152. <span class="comment">// And since we have an extended exponent range we can generate some really large </span>
  153. <span class="comment">// numbers here (4.02387260077093773543702433923004111e+2564):</span>
  154. <span class="identifier">std</span><span class="special">::</span><span class="identifier">cout</span> <span class="special">&lt;&lt;</span> <span class="identifier">boost</span><span class="special">::</span><span class="identifier">math</span><span class="special">::</span><span class="identifier">tgamma</span><span class="special">(</span><span class="identifier">float128</span><span class="special">(</span><span class="number">1000</span><span class="special">))</span> <span class="special">&lt;&lt;</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">endl</span><span class="special">;</span>
  155. <span class="comment">//</span>
  156. <span class="comment">// We can declare constants using GCC or Intel's native types, and the Q suffix,</span>
  157. <span class="comment">// these can be declared constexpr if required:</span>
  158. <span class="keyword">constexpr</span> <span class="identifier">float128</span> <span class="identifier">pi</span> <span class="special">=</span> <span class="number">3.1415926535897932384626433832795028841971693993751058</span><span class="identifier">Q</span><span class="special">;</span>
  159. <span class="keyword">return</span> <span class="number">0</span><span class="special">;</span>
  160. <span class="special">}</span>
  161. </pre>
  162. <div class="footnotes">
  163. <br><hr style="width:100; text-align:left;margin-left: 0">
  164. <div id="ftn.boost_multiprecision.tut.floats.float128.f0" class="footnote"><p><a href="#boost_multiprecision.tut.floats.float128.f0" class="para"><sup class="para">[3] </sup></a>
  165. On 64-bit Ubuntu 11.10, GCC-4.8.0, Intel Core 2 Duo T5800.
  166. </p></div>
  167. </div>
  168. </div>
  169. <table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
  170. <td align="left"></td>
  171. <td align="right"><div class="copyright-footer">Copyright &#169; 2002-2019 John Maddock
  172. and Christopher Kormanyos<p>
  173. Distributed under the Boost Software License, Version 1.0. (See accompanying
  174. 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>)
  175. </p>
  176. </div></td>
  177. </tr></table>
  178. <hr>
  179. <div class="spirit-nav">
  180. <a accesskey="p" href="mpfr_float.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../floats.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="fp_eg.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
  181. </div>
  182. </body>
  183. </html>