📁 File Manager Pro
v10.0.3 | PHP: 8.1.34
Server: Apache
2026-06-21 21:11:36
📂
/ (Root)
/
opt
/
alt
/
ruby19
/
share
/
doc
/
ruby
/
html
/
d7
/
deb
📍 /opt/alt/ruby19/share/doc/ruby/html/d7/deb
🔄 Refresh
✏️
Editing: acosh_8c_source.html
Read Only
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/> <meta http-equiv="X-UA-Compatible" content="IE=9"/> <meta name="generator" content="Doxygen 1.8.14"/> <meta name="viewport" content="width=device-width, initial-scale=1"/> <title>Ruby: missing/acosh.c Source File</title> <link href="../../tabs.css" rel="stylesheet" type="text/css"/> <script type="text/javascript" src="../../jquery.js"></script> <script type="text/javascript" src="../../dynsections.js"></script> <link href="../../doxygen.css" rel="stylesheet" type="text/css" /> </head> <body> <div id="top"><!-- do not remove this div, it is closed by doxygen! --> <div id="titlearea"> <table cellspacing="0" cellpadding="0"> <tbody> <tr style="height: 56px;"> <td id="projectalign" style="padding-left: 0.5em;"> <div id="projectname">Ruby  <span id="projectnumber">1.9.3p551(2014-11-13revision48407)</span> </div> </td> </tr> </tbody> </table> </div> <!-- end header part --> <!-- Generated by Doxygen 1.8.14 --> <script type="text/javascript" src="../../menudata.js"></script> <script type="text/javascript" src="../../menu.js"></script> <script type="text/javascript"> /* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */ $(function() { initMenu('../../',false,false,'search.php','Search'); }); /* @license-end */</script> <div id="main-nav"></div> <div id="nav-path" class="navpath"> <ul> <li class="navelem"><a class="el" href="../../dir_f3bfeebb553c3f6ecfb19202628b4493.html">missing</a></li> </ul> </div> </div><!-- top --> <div class="header"> <div class="headertitle"> <div class="title">acosh.c</div> </div> </div><!--header--> <div class="contents"> <a href="../../d7/deb/acosh_8c.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/**********************************************************************</span></div><div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"></span></div><div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment"> acosh.c -</span></div><div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"></span></div><div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> $Author: akr $</span></div><div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"> created at: Fri Apr 12 00:34:17 JST 2002</span></div><div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"></span></div><div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"> public domain rewrite of acosh(3), asinh(3) and atanh(3)</span></div><div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"></span></div><div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment">**********************************************************************/</span></div><div class="line"><a name="l00011"></a><span class="lineno"> 11</span> </div><div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="preprocessor">#include <errno.h></span></div><div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="preprocessor">#include <float.h></span></div><div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="preprocessor">#include <math.h></span></div><div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="preprocessor">#include "ruby.h"</span></div><div class="line"><a name="l00016"></a><span class="lineno"> 16</span> </div><div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment">/* DBL_MANT_DIG must be less than 4 times of bits of int */</span></div><div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="preprocessor">#ifndef DBL_MANT_DIG</span></div><div class="line"><a name="l00019"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#aeb343417fe0051b641fc7921022c7b0d"> 19</a></span> <span class="preprocessor">#define DBL_MANT_DIG 53 </span><span class="comment">/* in this case, at least 12 digit precision */</span><span class="preprocessor"></span></div><div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="preprocessor">#endif</span></div><div class="line"><a name="l00021"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#add8667ccee1e991eee871e95598a1317"> 21</a></span> <span class="preprocessor">#define BIG_CRITERIA_BIT (1<<DBL_MANT_DIG/2)</span></div><div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="preprocessor">#if BIG_CRITERIA_BIT > 0</span></div><div class="line"><a name="l00023"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#abe45bab9f28b4e54546e3294cd3d002a"> 23</a></span> <span class="preprocessor">#define BIG_CRITERIA (1.0*BIG_CRITERIA_BIT)</span></div><div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="preprocessor">#else</span></div><div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#define BIG_CRITERIA (1.0*(1<<DBL_MANT_DIG/4)*(1<<(DBL_MANT_DIG/2+1-DBL_MANT_DIG/4)))</span></div><div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="preprocessor">#endif</span></div><div class="line"><a name="l00027"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#a2557f6f1399be0eb3f96d7bfcbf64f55"> 27</a></span> <span class="preprocessor">#define SMALL_CRITERIA_BIT (1<<(DBL_MANT_DIG/3))</span></div><div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="preprocessor">#if SMALL_CRITERIA_BIT > 0</span></div><div class="line"><a name="l00029"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#aa45bb2d04070dad96f9c67a33fc59c01"> 29</a></span> <span class="preprocessor">#define SMALL_CRITERIA (1.0/SMALL_CRITERIA_BIT)</span></div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="preprocessor">#else</span></div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="preprocessor">#define SMALL_CRITERIA (1.0*(1<<DBL_MANT_DIG/4)*(1<<(DBL_MANT_DIG/3+1-DBL_MANT_DIG/4)))</span></div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="preprocessor">#endif</span></div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span> </div><div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="preprocessor">#ifndef HAVE_ACOSH</span></div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="keywordtype">double</span></div><div class="line"><a name="l00036"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#ac862445bcb46911b1782c284829d5e46"> 36</a></span> <a class="code" href="../../d7/deb/acosh_8c.html#ac862445bcb46911b1782c284829d5e46">acosh</a>(<span class="keywordtype">double</span> x)</div><div class="line"><a name="l00037"></a><span class="lineno"> 37</span> {</div><div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  <span class="keywordflow">if</span> (x < 1)</div><div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  x = -1; <span class="comment">/* NaN */</span></div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span> (x == 1)</div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  <span class="keywordflow">return</span> 0;</div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span> (x > <a class="code" href="../../d7/deb/acosh_8c.html#abe45bab9f28b4e54546e3294cd3d002a">BIG_CRITERIA</a>)</div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  x += x;</div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span>  <span class="keywordflow">else</span></div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  x += sqrt((x + 1) * (x - 1));</div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <span class="keywordflow">return</span> log(x);</div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span> }</div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span> <span class="preprocessor">#endif</span></div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span> </div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="preprocessor">#ifndef HAVE_ASINH</span></div><div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="keywordtype">double</span></div><div class="line"><a name="l00052"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#a907297107e93ea9c8a0c8409a8e65079"> 52</a></span> <a class="code" href="../../d7/deb/acosh_8c.html#a907297107e93ea9c8a0c8409a8e65079">asinh</a>(<span class="keywordtype">double</span> x)</div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span> {</div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  <span class="keywordtype">int</span> <a class="code" href="../../df/d73/time_8c.html#a83f11f6624f6648aa502330ab8408073">neg</a> = x < 0;</div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <span class="keywordtype">double</span> z = fabs(x);</div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span> </div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  <span class="keywordflow">if</span> (z < <a class="code" href="../../d7/deb/acosh_8c.html#aa45bb2d04070dad96f9c67a33fc59c01">SMALL_CRITERIA</a>) <span class="keywordflow">return</span> x;</div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  <span class="keywordflow">if</span> (z < (1.0/(1<<<a class="code" href="../../d7/deb/acosh_8c.html#aeb343417fe0051b641fc7921022c7b0d">DBL_MANT_DIG</a>/5))) {</div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="keywordtype">double</span> x2 = z * z;</div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  z *= 1 + x2 * (-1.0/6.0 + x2 * 3.0/40.0);</div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  }</div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span> (z > <a class="code" href="../../d7/deb/acosh_8c.html#abe45bab9f28b4e54546e3294cd3d002a">BIG_CRITERIA</a>) {</div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  z = log(z + z);</div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  }</div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  <span class="keywordflow">else</span> {</div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  z = log(z + sqrt(z * z + 1));</div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  }</div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  <span class="keywordflow">if</span> (<a class="code" href="../../df/d73/time_8c.html#a83f11f6624f6648aa502330ab8408073">neg</a>) z = -z;</div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  <span class="keywordflow">return</span> z;</div><div class="line"><a name="l00070"></a><span class="lineno"> 70</span> }</div><div class="line"><a name="l00071"></a><span class="lineno"> 71</span> <span class="preprocessor">#endif</span></div><div class="line"><a name="l00072"></a><span class="lineno"> 72</span> </div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span> <span class="preprocessor">#ifndef HAVE_ATANH</span></div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span> <span class="keywordtype">double</span></div><div class="line"><a name="l00075"></a><span class="lineno"><a class="line" href="../../d7/deb/acosh_8c.html#ac07fd7b77fa341563abccf9efcc1a711"> 75</a></span> <a class="code" href="../../d7/deb/acosh_8c.html#ac07fd7b77fa341563abccf9efcc1a711">atanh</a>(<span class="keywordtype">double</span> x)</div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span> {</div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  <span class="keywordtype">int</span> <a class="code" href="../../df/d73/time_8c.html#a83f11f6624f6648aa502330ab8408073">neg</a> = x < 0;</div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  <span class="keywordtype">double</span> z = fabs(x);</div><div class="line"><a name="l00079"></a><span class="lineno"> 79</span> </div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <span class="keywordflow">if</span> (z < <a class="code" href="../../d7/deb/acosh_8c.html#aa45bb2d04070dad96f9c67a33fc59c01">SMALL_CRITERIA</a>) <span class="keywordflow">return</span> x;</div><div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  z = log(z > 1 ? -1 : (1 + z) / (1 - z)) / 2;</div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  <span class="keywordflow">if</span> (<a class="code" href="../../df/d73/time_8c.html#a83f11f6624f6648aa502330ab8408073">neg</a>) z = -z;</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="keywordflow">if</span> (<a class="code" href="../../d3/d90/missing_8h.html#ae7a5ae88899707625ed6ab2ec7045583">isinf</a>(z))</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span> <span class="preprocessor">#if defined(ERANGE)</span></div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  <a class="code" href="../../d4/db0/__sdbm_8c.html#ad65a8842cc674e3ddf69355898c0ecbf">errno</a> = ERANGE;</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span> <span class="preprocessor">#elif defined(EDOM)</span></div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  <a class="code" href="../../d4/db0/__sdbm_8c.html#ad65a8842cc674e3ddf69355898c0ecbf">errno</a> = EDOM;</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span> <span class="preprocessor">#else</span></div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  ;</div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span> <span class="preprocessor">#endif</span></div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  <span class="keywordflow">return</span> z;</div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span> }</div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span> <span class="preprocessor">#endif</span></div><div class="ttc" id="time_8c_html_a83f11f6624f6648aa502330ab8408073"><div class="ttname"><a href="../../df/d73/time_8c.html#a83f11f6624f6648aa502330ab8408073">neg</a></div><div class="ttdeci">#define neg(x)</div><div class="ttdef"><b>Definition:</b> <a href="../../df/d73/time_8c_source.html#l00171">time.c:171</a></div></div> <div class="ttc" id="acosh_8c_html_a907297107e93ea9c8a0c8409a8e65079"><div class="ttname"><a href="../../d7/deb/acosh_8c.html#a907297107e93ea9c8a0c8409a8e65079">asinh</a></div><div class="ttdeci">double asinh(double x)</div><div class="ttdef"><b>Definition:</b> <a href="../../d7/deb/acosh_8c_source.html#l00052">acosh.c:52</a></div></div> <div class="ttc" id="missing_8h_html_ae7a5ae88899707625ed6ab2ec7045583"><div class="ttname"><a href="../../d3/d90/missing_8h.html#ae7a5ae88899707625ed6ab2ec7045583">isinf</a></div><div class="ttdeci">RUBY_EXTERN int isinf(double)</div><div class="ttdef"><b>Definition:</b> <a href="../../d9/d89/isinf_8c_source.html#l00056">isinf.c:56</a></div></div> <div class="ttc" id="__sdbm_8c_html_ad65a8842cc674e3ddf69355898c0ecbf"><div class="ttname"><a href="../../d4/db0/__sdbm_8c.html#ad65a8842cc674e3ddf69355898c0ecbf">errno</a></div><div class="ttdeci">int errno</div></div> <div class="ttc" id="acosh_8c_html_aeb343417fe0051b641fc7921022c7b0d"><div class="ttname"><a href="../../d7/deb/acosh_8c.html#aeb343417fe0051b641fc7921022c7b0d">DBL_MANT_DIG</a></div><div class="ttdeci">#define DBL_MANT_DIG</div><div class="ttdef"><b>Definition:</b> <a href="../../d7/deb/acosh_8c_source.html#l00019">acosh.c:19</a></div></div> <div class="ttc" id="acosh_8c_html_abe45bab9f28b4e54546e3294cd3d002a"><div class="ttname"><a href="../../d7/deb/acosh_8c.html#abe45bab9f28b4e54546e3294cd3d002a">BIG_CRITERIA</a></div><div class="ttdeci">#define BIG_CRITERIA</div><div class="ttdef"><b>Definition:</b> <a href="../../d7/deb/acosh_8c_source.html#l00023">acosh.c:23</a></div></div> <div class="ttc" id="acosh_8c_html_ac862445bcb46911b1782c284829d5e46"><div class="ttname"><a href="../../d7/deb/acosh_8c.html#ac862445bcb46911b1782c284829d5e46">acosh</a></div><div class="ttdeci">double acosh(double x)</div><div class="ttdef"><b>Definition:</b> <a href="../../d7/deb/acosh_8c_source.html#l00036">acosh.c:36</a></div></div> <div class="ttc" id="acosh_8c_html_aa45bb2d04070dad96f9c67a33fc59c01"><div class="ttname"><a href="../../d7/deb/acosh_8c.html#aa45bb2d04070dad96f9c67a33fc59c01">SMALL_CRITERIA</a></div><div class="ttdeci">#define SMALL_CRITERIA</div><div class="ttdef"><b>Definition:</b> <a href="../../d7/deb/acosh_8c_source.html#l00029">acosh.c:29</a></div></div> <div class="ttc" id="acosh_8c_html_ac07fd7b77fa341563abccf9efcc1a711"><div class="ttname"><a href="../../d7/deb/acosh_8c.html#ac07fd7b77fa341563abccf9efcc1a711">atanh</a></div><div class="ttdeci">double atanh(double x)</div><div class="ttdef"><b>Definition:</b> <a href="../../d7/deb/acosh_8c_source.html#l00075">acosh.c:75</a></div></div> </div><!-- fragment --></div><!-- contents --> <!-- start footer part --> <hr class="footer"/><address class="footer"><small> Generated by  <a href="http://www.doxygen.org/index.html"> <img class="footer" src="../../doxygen.png" alt="doxygen"/> </a> 1.8.14 </small></address> </body> </html>
💾 Save Changes
❌ Cancel