commons-commits mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From er...@apache.org
Subject [42/51] [partial] commons-complex git commit: Removed "site-content" directory.
Date Fri, 06 Jan 2017 17:10:22 GMT
http://git-wip-us.apache.org/repos/asf/commons-complex/blob/b3576eeb/site-content/.svn/pristine/01/01537ca4c7defdef6d3b73f3ae9ddd8477d5e619.svn-base
----------------------------------------------------------------------
diff --git a/site-content/.svn/pristine/01/01537ca4c7defdef6d3b73f3ae9ddd8477d5e619.svn-base b/site-content/.svn/pristine/01/01537ca4c7defdef6d3b73f3ae9ddd8477d5e619.svn-base
deleted file mode 100644
index b2b3e9c..0000000
--- a/site-content/.svn/pristine/01/01537ca4c7defdef6d3b73f3ae9ddd8477d5e619.svn-base
+++ /dev/null
@@ -1,671 +0,0 @@
-<!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" xml:lang="en" lang="en">
-<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-<title>FastMathCalc xref</title>
-<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" />
-</head>
-<body>
-<div id="overview"><a href="../../../../../../apidocs/org/apache/commons/math3/util/FastMathCalc.html">View Javadoc</a></div><pre>
-<a class="jxr_linenumber" name="L1" href="#L1">1</a>   <em class="jxr_comment">/*</em>
-<a class="jxr_linenumber" name="L2" href="#L2">2</a>   <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em>
-<a class="jxr_linenumber" name="L3" href="#L3">3</a>   <em class="jxr_comment"> * contributor license agreements.  See the NOTICE file distributed with</em>
-<a class="jxr_linenumber" name="L4" href="#L4">4</a>   <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em>
-<a class="jxr_linenumber" name="L5" href="#L5">5</a>   <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em>
-<a class="jxr_linenumber" name="L6" href="#L6">6</a>   <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em>
-<a class="jxr_linenumber" name="L7" href="#L7">7</a>   <em class="jxr_comment"> * the License.  You may obtain a copy of the License at</em>
-<a class="jxr_linenumber" name="L8" href="#L8">8</a>   <em class="jxr_comment"> *</em>
-<a class="jxr_linenumber" name="L9" href="#L9">9</a>   <em class="jxr_comment"> *      <a href="http://www.apache.org/licenses/LICENSE-2." target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.</a>0</em>
-<a class="jxr_linenumber" name="L10" href="#L10">10</a>  <em class="jxr_comment"> *</em>
-<a class="jxr_linenumber" name="L11" href="#L11">11</a>  <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em>
-<a class="jxr_linenumber" name="L12" href="#L12">12</a>  <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em>
-<a class="jxr_linenumber" name="L13" href="#L13">13</a>  <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em>
-<a class="jxr_linenumber" name="L14" href="#L14">14</a>  <em class="jxr_comment"> * See the License for the specific language governing permissions and</em>
-<a class="jxr_linenumber" name="L15" href="#L15">15</a>  <em class="jxr_comment"> * limitations under the License.</em>
-<a class="jxr_linenumber" name="L16" href="#L16">16</a>  <em class="jxr_comment"> */</em>
-<a class="jxr_linenumber" name="L17" href="#L17">17</a>  <strong class="jxr_keyword">package</strong> org.apache.commons.math3.util;
-<a class="jxr_linenumber" name="L18" href="#L18">18</a>  
-<a class="jxr_linenumber" name="L19" href="#L19">19</a>  <strong class="jxr_keyword">import</strong> java.io.PrintStream;
-<a class="jxr_linenumber" name="L20" href="#L20">20</a>  
-<a class="jxr_linenumber" name="L21" href="#L21">21</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.exception.DimensionMismatchException;
-<a class="jxr_linenumber" name="L22" href="#L22">22</a>  
-<a class="jxr_linenumber" name="L23" href="#L23">23</a>  <em class="jxr_javadoccomment">/** Class used to compute the classical functions tables.</em>
-<a class="jxr_linenumber" name="L24" href="#L24">24</a>  <em class="jxr_javadoccomment"> * @since 3.0</em>
-<a class="jxr_linenumber" name="L25" href="#L25">25</a>  <em class="jxr_javadoccomment"> */</em>
-<a class="jxr_linenumber" name="L26" href="#L26">26</a>  <strong class="jxr_keyword">class</strong> <a href="../../../../../org/apache/commons/math3/util/FastMathCalc.html">FastMathCalc</a> {
-<a class="jxr_linenumber" name="L27" href="#L27">27</a>  
-<a class="jxr_linenumber" name="L28" href="#L28">28</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L29" href="#L29">29</a>  <em class="jxr_javadoccomment">     * 0x40000000 - used to split a double into two parts, both with the low order bits cleared.</em>
-<a class="jxr_linenumber" name="L30" href="#L30">30</a>  <em class="jxr_javadoccomment">     * Equivalent to 2^30.</em>
-<a class="jxr_linenumber" name="L31" href="#L31">31</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L32" href="#L32">32</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> HEX_40000000 = 0x40000000L; <em class="jxr_comment">// 1073741824L</em>
-<a class="jxr_linenumber" name="L33" href="#L33">33</a>  
-<a class="jxr_linenumber" name="L34" href="#L34">34</a>      <em class="jxr_javadoccomment">/** Factorial table, for Taylor series expansions. 0!, 1!, 2!, ... 19! */</em>
-<a class="jxr_linenumber" name="L35" href="#L35">35</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> FACT[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[]
-<a class="jxr_linenumber" name="L36" href="#L36">36</a>          {
-<a class="jxr_linenumber" name="L37" href="#L37">37</a>          +1.0d,                        <em class="jxr_comment">// 0</em>
-<a class="jxr_linenumber" name="L38" href="#L38">38</a>          +1.0d,                        <em class="jxr_comment">// 1</em>
-<a class="jxr_linenumber" name="L39" href="#L39">39</a>          +2.0d,                        <em class="jxr_comment">// 2</em>
-<a class="jxr_linenumber" name="L40" href="#L40">40</a>          +6.0d,                        <em class="jxr_comment">// 3</em>
-<a class="jxr_linenumber" name="L41" href="#L41">41</a>          +24.0d,                       <em class="jxr_comment">// 4</em>
-<a class="jxr_linenumber" name="L42" href="#L42">42</a>          +120.0d,                      <em class="jxr_comment">// 5</em>
-<a class="jxr_linenumber" name="L43" href="#L43">43</a>          +720.0d,                      <em class="jxr_comment">// 6</em>
-<a class="jxr_linenumber" name="L44" href="#L44">44</a>          +5040.0d,                     <em class="jxr_comment">// 7</em>
-<a class="jxr_linenumber" name="L45" href="#L45">45</a>          +40320.0d,                    <em class="jxr_comment">// 8</em>
-<a class="jxr_linenumber" name="L46" href="#L46">46</a>          +362880.0d,                   <em class="jxr_comment">// 9</em>
-<a class="jxr_linenumber" name="L47" href="#L47">47</a>          +3628800.0d,                  <em class="jxr_comment">// 10</em>
-<a class="jxr_linenumber" name="L48" href="#L48">48</a>          +39916800.0d,                 <em class="jxr_comment">// 11</em>
-<a class="jxr_linenumber" name="L49" href="#L49">49</a>          +479001600.0d,                <em class="jxr_comment">// 12</em>
-<a class="jxr_linenumber" name="L50" href="#L50">50</a>          +6227020800.0d,               <em class="jxr_comment">// 13</em>
-<a class="jxr_linenumber" name="L51" href="#L51">51</a>          +87178291200.0d,              <em class="jxr_comment">// 14</em>
-<a class="jxr_linenumber" name="L52" href="#L52">52</a>          +1307674368000.0d,            <em class="jxr_comment">// 15</em>
-<a class="jxr_linenumber" name="L53" href="#L53">53</a>          +20922789888000.0d,           <em class="jxr_comment">// 16</em>
-<a class="jxr_linenumber" name="L54" href="#L54">54</a>          +355687428096000.0d,          <em class="jxr_comment">// 17</em>
-<a class="jxr_linenumber" name="L55" href="#L55">55</a>          +6402373705728000.0d,         <em class="jxr_comment">// 18</em>
-<a class="jxr_linenumber" name="L56" href="#L56">56</a>          +121645100408832000.0d,       <em class="jxr_comment">// 19</em>
-<a class="jxr_linenumber" name="L57" href="#L57">57</a>          };
-<a class="jxr_linenumber" name="L58" href="#L58">58</a>  
-<a class="jxr_linenumber" name="L59" href="#L59">59</a>      <em class="jxr_javadoccomment">/** Coefficients for slowLog. */</em>
-<a class="jxr_linenumber" name="L60" href="#L60">60</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> LN_SPLIT_COEF[][] = {
-<a class="jxr_linenumber" name="L61" href="#L61">61</a>          {2.0, 0.0},
-<a class="jxr_linenumber" name="L62" href="#L62">62</a>          {0.6666666269302368, 3.9736429850260626E-8},
-<a class="jxr_linenumber" name="L63" href="#L63">63</a>          {0.3999999761581421, 2.3841857910019882E-8},
-<a class="jxr_linenumber" name="L64" href="#L64">64</a>          {0.2857142686843872, 1.7029898543501842E-8},
-<a class="jxr_linenumber" name="L65" href="#L65">65</a>          {0.2222222089767456, 1.3245471311735498E-8},
-<a class="jxr_linenumber" name="L66" href="#L66">66</a>          {0.1818181574344635, 2.4384203044354907E-8},
-<a class="jxr_linenumber" name="L67" href="#L67">67</a>          {0.1538461446762085, 9.140260083262505E-9},
-<a class="jxr_linenumber" name="L68" href="#L68">68</a>          {0.13333332538604736, 9.220590270857665E-9},
-<a class="jxr_linenumber" name="L69" href="#L69">69</a>          {0.11764700710773468, 1.2393345855018391E-8},
-<a class="jxr_linenumber" name="L70" href="#L70">70</a>          {0.10526403784751892, 8.251545029714408E-9},
-<a class="jxr_linenumber" name="L71" href="#L71">71</a>          {0.0952233225107193, 1.2675934823758863E-8},
-<a class="jxr_linenumber" name="L72" href="#L72">72</a>          {0.08713622391223907, 1.1430250008909141E-8},
-<a class="jxr_linenumber" name="L73" href="#L73">73</a>          {0.07842259109020233, 2.404307984052299E-9},
-<a class="jxr_linenumber" name="L74" href="#L74">74</a>          {0.08371849358081818, 1.176342548272881E-8},
-<a class="jxr_linenumber" name="L75" href="#L75">75</a>          {0.030589580535888672, 1.2958646899018938E-9},
-<a class="jxr_linenumber" name="L76" href="#L76">76</a>          {0.14982303977012634, 1.225743062930824E-8},
-<a class="jxr_linenumber" name="L77" href="#L77">77</a>      };
-<a class="jxr_linenumber" name="L78" href="#L78">78</a>  
-<a class="jxr_linenumber" name="L79" href="#L79">79</a>      <em class="jxr_javadoccomment">/** Table start declaration. */</em>
-<a class="jxr_linenumber" name="L80" href="#L80">80</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> String TABLE_START_DECL = <span class="jxr_string">"    {"</span>;
-<a class="jxr_linenumber" name="L81" href="#L81">81</a>  
-<a class="jxr_linenumber" name="L82" href="#L82">82</a>      <em class="jxr_javadoccomment">/** Table end declaration. */</em>
-<a class="jxr_linenumber" name="L83" href="#L83">83</a>      <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> String TABLE_END_DECL   = <span class="jxr_string">"    };"</span>;
-<a class="jxr_linenumber" name="L84" href="#L84">84</a>  
-<a class="jxr_linenumber" name="L85" href="#L85">85</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L86" href="#L86">86</a>  <em class="jxr_javadoccomment">     * Private Constructor.</em>
-<a class="jxr_linenumber" name="L87" href="#L87">87</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L88" href="#L88">88</a>      <strong class="jxr_keyword">private</strong> <a href="../../../../../org/apache/commons/math3/util/FastMathCalc.html">FastMathCalc</a>() {
-<a class="jxr_linenumber" name="L89" href="#L89">89</a>      }
-<a class="jxr_linenumber" name="L90" href="#L90">90</a>  
-<a class="jxr_linenumber" name="L91" href="#L91">91</a>      <em class="jxr_javadoccomment">/** Build the sine and cosine tables.</em>
-<a class="jxr_linenumber" name="L92" href="#L92">92</a>  <em class="jxr_javadoccomment">     * @param SINE_TABLE_A table of the most significant part of the sines</em>
-<a class="jxr_linenumber" name="L93" href="#L93">93</a>  <em class="jxr_javadoccomment">     * @param SINE_TABLE_B table of the least significant part of the sines</em>
-<a class="jxr_linenumber" name="L94" href="#L94">94</a>  <em class="jxr_javadoccomment">     * @param COSINE_TABLE_A table of the most significant part of the cosines</em>
-<a class="jxr_linenumber" name="L95" href="#L95">95</a>  <em class="jxr_javadoccomment">     * @param COSINE_TABLE_B table of the most significant part of the cosines</em>
-<a class="jxr_linenumber" name="L96" href="#L96">96</a>  <em class="jxr_javadoccomment">     * @param SINE_TABLE_LEN length of the tables</em>
-<a class="jxr_linenumber" name="L97" href="#L97">97</a>  <em class="jxr_javadoccomment">     * @param TANGENT_TABLE_A table of the most significant part of the tangents</em>
-<a class="jxr_linenumber" name="L98" href="#L98">98</a>  <em class="jxr_javadoccomment">     * @param TANGENT_TABLE_B table of the most significant part of the tangents</em>
-<a class="jxr_linenumber" name="L99" href="#L99">99</a>  <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L100" href="#L100">100</a>     @SuppressWarnings(<span class="jxr_string">"unused"</span>)
-<a class="jxr_linenumber" name="L101" href="#L101">101</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> buildSinCosTables(<strong class="jxr_keyword">double</strong>[] SINE_TABLE_A, <strong class="jxr_keyword">double</strong>[] SINE_TABLE_B,
-<a class="jxr_linenumber" name="L102" href="#L102">102</a>                                           <strong class="jxr_keyword">double</strong>[] COSINE_TABLE_A, <strong class="jxr_keyword">double</strong>[] COSINE_TABLE_B,
-<a class="jxr_linenumber" name="L103" href="#L103">103</a>                                           <strong class="jxr_keyword">int</strong> SINE_TABLE_LEN, <strong class="jxr_keyword">double</strong>[] TANGENT_TABLE_A, <strong class="jxr_keyword">double</strong>[] TANGENT_TABLE_B) {
-<a class="jxr_linenumber" name="L104" href="#L104">104</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L105" href="#L105">105</a> 
-<a class="jxr_linenumber" name="L106" href="#L106">106</a>         <em class="jxr_comment">/* Use taylor series for 0 &lt;= x &lt;= 6/8 */</em>
-<a class="jxr_linenumber" name="L107" href="#L107">107</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; 7; i++) {
-<a class="jxr_linenumber" name="L108" href="#L108">108</a>             <strong class="jxr_keyword">double</strong> x = i / 8.0;
-<a class="jxr_linenumber" name="L109" href="#L109">109</a> 
-<a class="jxr_linenumber" name="L110" href="#L110">110</a>             slowSin(x, result);
-<a class="jxr_linenumber" name="L111" href="#L111">111</a>             SINE_TABLE_A[i] = result[0];
-<a class="jxr_linenumber" name="L112" href="#L112">112</a>             SINE_TABLE_B[i] = result[1];
-<a class="jxr_linenumber" name="L113" href="#L113">113</a> 
-<a class="jxr_linenumber" name="L114" href="#L114">114</a>             slowCos(x, result);
-<a class="jxr_linenumber" name="L115" href="#L115">115</a>             COSINE_TABLE_A[i] = result[0];
-<a class="jxr_linenumber" name="L116" href="#L116">116</a>             COSINE_TABLE_B[i] = result[1];
-<a class="jxr_linenumber" name="L117" href="#L117">117</a>         }
-<a class="jxr_linenumber" name="L118" href="#L118">118</a> 
-<a class="jxr_linenumber" name="L119" href="#L119">119</a>         <em class="jxr_comment">/* Use angle addition formula to complete table to 13/8, just beyond pi/2 */</em>
-<a class="jxr_linenumber" name="L120" href="#L120">120</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 7; i &lt; SINE_TABLE_LEN; i++) {
-<a class="jxr_linenumber" name="L121" href="#L121">121</a>             <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L122" href="#L122">122</a>             <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L123" href="#L123">123</a>             <strong class="jxr_keyword">double</strong> as[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L124" href="#L124">124</a>             <strong class="jxr_keyword">double</strong> bs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L125" href="#L125">125</a>             <strong class="jxr_keyword">double</strong> temps[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L126" href="#L126">126</a> 
-<a class="jxr_linenumber" name="L127" href="#L127">127</a>             <strong class="jxr_keyword">if</strong> ( (i &amp; 1) == 0) {
-<a class="jxr_linenumber" name="L128" href="#L128">128</a>                 <em class="jxr_comment">// Even, use double angle</em>
-<a class="jxr_linenumber" name="L129" href="#L129">129</a>                 xs[0] = SINE_TABLE_A[i/2];
-<a class="jxr_linenumber" name="L130" href="#L130">130</a>                 xs[1] = SINE_TABLE_B[i/2];
-<a class="jxr_linenumber" name="L131" href="#L131">131</a>                 ys[0] = COSINE_TABLE_A[i/2];
-<a class="jxr_linenumber" name="L132" href="#L132">132</a>                 ys[1] = COSINE_TABLE_B[i/2];
-<a class="jxr_linenumber" name="L133" href="#L133">133</a> 
-<a class="jxr_linenumber" name="L134" href="#L134">134</a>                 <em class="jxr_comment">/* compute sine */</em>
-<a class="jxr_linenumber" name="L135" href="#L135">135</a>                 splitMult(xs, ys, result);
-<a class="jxr_linenumber" name="L136" href="#L136">136</a>                 SINE_TABLE_A[i] = result[0] * 2.0;
-<a class="jxr_linenumber" name="L137" href="#L137">137</a>                 SINE_TABLE_B[i] = result[1] * 2.0;
-<a class="jxr_linenumber" name="L138" href="#L138">138</a> 
-<a class="jxr_linenumber" name="L139" href="#L139">139</a>                 <em class="jxr_comment">/* Compute cosine */</em>
-<a class="jxr_linenumber" name="L140" href="#L140">140</a>                 splitMult(ys, ys, as);
-<a class="jxr_linenumber" name="L141" href="#L141">141</a>                 splitMult(xs, xs, temps);
-<a class="jxr_linenumber" name="L142" href="#L142">142</a>                 temps[0] = -temps[0];
-<a class="jxr_linenumber" name="L143" href="#L143">143</a>                 temps[1] = -temps[1];
-<a class="jxr_linenumber" name="L144" href="#L144">144</a>                 splitAdd(as, temps, result);
-<a class="jxr_linenumber" name="L145" href="#L145">145</a>                 COSINE_TABLE_A[i] = result[0];
-<a class="jxr_linenumber" name="L146" href="#L146">146</a>                 COSINE_TABLE_B[i] = result[1];
-<a class="jxr_linenumber" name="L147" href="#L147">147</a>             } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L148" href="#L148">148</a>                 xs[0] = SINE_TABLE_A[i/2];
-<a class="jxr_linenumber" name="L149" href="#L149">149</a>                 xs[1] = SINE_TABLE_B[i/2];
-<a class="jxr_linenumber" name="L150" href="#L150">150</a>                 ys[0] = COSINE_TABLE_A[i/2];
-<a class="jxr_linenumber" name="L151" href="#L151">151</a>                 ys[1] = COSINE_TABLE_B[i/2];
-<a class="jxr_linenumber" name="L152" href="#L152">152</a>                 as[0] = SINE_TABLE_A[i/2+1];
-<a class="jxr_linenumber" name="L153" href="#L153">153</a>                 as[1] = SINE_TABLE_B[i/2+1];
-<a class="jxr_linenumber" name="L154" href="#L154">154</a>                 bs[0] = COSINE_TABLE_A[i/2+1];
-<a class="jxr_linenumber" name="L155" href="#L155">155</a>                 bs[1] = COSINE_TABLE_B[i/2+1];
-<a class="jxr_linenumber" name="L156" href="#L156">156</a> 
-<a class="jxr_linenumber" name="L157" href="#L157">157</a>                 <em class="jxr_comment">/* compute sine */</em>
-<a class="jxr_linenumber" name="L158" href="#L158">158</a>                 splitMult(xs, bs, temps);
-<a class="jxr_linenumber" name="L159" href="#L159">159</a>                 splitMult(ys, as, result);
-<a class="jxr_linenumber" name="L160" href="#L160">160</a>                 splitAdd(result, temps, result);
-<a class="jxr_linenumber" name="L161" href="#L161">161</a>                 SINE_TABLE_A[i] = result[0];
-<a class="jxr_linenumber" name="L162" href="#L162">162</a>                 SINE_TABLE_B[i] = result[1];
-<a class="jxr_linenumber" name="L163" href="#L163">163</a> 
-<a class="jxr_linenumber" name="L164" href="#L164">164</a>                 <em class="jxr_comment">/* Compute cosine */</em>
-<a class="jxr_linenumber" name="L165" href="#L165">165</a>                 splitMult(ys, bs, result);
-<a class="jxr_linenumber" name="L166" href="#L166">166</a>                 splitMult(xs, as, temps);
-<a class="jxr_linenumber" name="L167" href="#L167">167</a>                 temps[0] = -temps[0];
-<a class="jxr_linenumber" name="L168" href="#L168">168</a>                 temps[1] = -temps[1];
-<a class="jxr_linenumber" name="L169" href="#L169">169</a>                 splitAdd(result, temps, result);
-<a class="jxr_linenumber" name="L170" href="#L170">170</a>                 COSINE_TABLE_A[i] = result[0];
-<a class="jxr_linenumber" name="L171" href="#L171">171</a>                 COSINE_TABLE_B[i] = result[1];
-<a class="jxr_linenumber" name="L172" href="#L172">172</a>             }
-<a class="jxr_linenumber" name="L173" href="#L173">173</a>         }
-<a class="jxr_linenumber" name="L174" href="#L174">174</a> 
-<a class="jxr_linenumber" name="L175" href="#L175">175</a>         <em class="jxr_comment">/* Compute tangent = sine/cosine */</em>
-<a class="jxr_linenumber" name="L176" href="#L176">176</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; SINE_TABLE_LEN; i++) {
-<a class="jxr_linenumber" name="L177" href="#L177">177</a>             <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L178" href="#L178">178</a>             <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L179" href="#L179">179</a>             <strong class="jxr_keyword">double</strong> as[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L180" href="#L180">180</a> 
-<a class="jxr_linenumber" name="L181" href="#L181">181</a>             as[0] = COSINE_TABLE_A[i];
-<a class="jxr_linenumber" name="L182" href="#L182">182</a>             as[1] = COSINE_TABLE_B[i];
-<a class="jxr_linenumber" name="L183" href="#L183">183</a> 
-<a class="jxr_linenumber" name="L184" href="#L184">184</a>             splitReciprocal(as, ys);
-<a class="jxr_linenumber" name="L185" href="#L185">185</a> 
-<a class="jxr_linenumber" name="L186" href="#L186">186</a>             xs[0] = SINE_TABLE_A[i];
-<a class="jxr_linenumber" name="L187" href="#L187">187</a>             xs[1] = SINE_TABLE_B[i];
-<a class="jxr_linenumber" name="L188" href="#L188">188</a> 
-<a class="jxr_linenumber" name="L189" href="#L189">189</a>             splitMult(xs, ys, as);
-<a class="jxr_linenumber" name="L190" href="#L190">190</a> 
-<a class="jxr_linenumber" name="L191" href="#L191">191</a>             TANGENT_TABLE_A[i] = as[0];
-<a class="jxr_linenumber" name="L192" href="#L192">192</a>             TANGENT_TABLE_B[i] = as[1];
-<a class="jxr_linenumber" name="L193" href="#L193">193</a>         }
-<a class="jxr_linenumber" name="L194" href="#L194">194</a> 
-<a class="jxr_linenumber" name="L195" href="#L195">195</a>     }
-<a class="jxr_linenumber" name="L196" href="#L196">196</a> 
-<a class="jxr_linenumber" name="L197" href="#L197">197</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L198" href="#L198">198</a> <em class="jxr_javadoccomment">     *  For x between 0 and pi/4 compute cosine using Talor series</em>
-<a class="jxr_linenumber" name="L199" href="#L199">199</a> <em class="jxr_javadoccomment">     *  cos(x) = 1 - x^2/2! + x^4/4! ...</em>
-<a class="jxr_linenumber" name="L200" href="#L200">200</a> <em class="jxr_javadoccomment">     * @param x number from which cosine is requested</em>
-<a class="jxr_linenumber" name="L201" href="#L201">201</a> <em class="jxr_javadoccomment">     * @param result placeholder where to put the result in extended precision</em>
-<a class="jxr_linenumber" name="L202" href="#L202">202</a> <em class="jxr_javadoccomment">     * (may be null)</em>
-<a class="jxr_linenumber" name="L203" href="#L203">203</a> <em class="jxr_javadoccomment">     * @return cos(x)</em>
-<a class="jxr_linenumber" name="L204" href="#L204">204</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L205" href="#L205">205</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">double</strong> slowCos(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[]) {
-<a class="jxr_linenumber" name="L206" href="#L206">206</a> 
-<a class="jxr_linenumber" name="L207" href="#L207">207</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L208" href="#L208">208</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L209" href="#L209">209</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> facts[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L210" href="#L210">210</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> as[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L211" href="#L211">211</a>         split(x, xs);
-<a class="jxr_linenumber" name="L212" href="#L212">212</a>         ys[0] = ys[1] = 0.0;
-<a class="jxr_linenumber" name="L213" href="#L213">213</a> 
-<a class="jxr_linenumber" name="L214" href="#L214">214</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = FACT.length-1; i &gt;= 0; i--) {
-<a class="jxr_linenumber" name="L215" href="#L215">215</a>             splitMult(xs, ys, as);
-<a class="jxr_linenumber" name="L216" href="#L216">216</a>             ys[0] = as[0]; ys[1] = as[1];
-<a class="jxr_linenumber" name="L217" href="#L217">217</a> 
-<a class="jxr_linenumber" name="L218" href="#L218">218</a>             <strong class="jxr_keyword">if</strong> ( (i &amp; 1) != 0) { <em class="jxr_comment">// skip odd entries</em>
-<a class="jxr_linenumber" name="L219" href="#L219">219</a>                 <strong class="jxr_keyword">continue</strong>;
-<a class="jxr_linenumber" name="L220" href="#L220">220</a>             }
-<a class="jxr_linenumber" name="L221" href="#L221">221</a> 
-<a class="jxr_linenumber" name="L222" href="#L222">222</a>             split(FACT[i], as);
-<a class="jxr_linenumber" name="L223" href="#L223">223</a>             splitReciprocal(as, facts);
-<a class="jxr_linenumber" name="L224" href="#L224">224</a> 
-<a class="jxr_linenumber" name="L225" href="#L225">225</a>             <strong class="jxr_keyword">if</strong> ( (i &amp; 2) != 0 ) { <em class="jxr_comment">// alternate terms are negative</em>
-<a class="jxr_linenumber" name="L226" href="#L226">226</a>                 facts[0] = -facts[0];
-<a class="jxr_linenumber" name="L227" href="#L227">227</a>                 facts[1] = -facts[1];
-<a class="jxr_linenumber" name="L228" href="#L228">228</a>             }
-<a class="jxr_linenumber" name="L229" href="#L229">229</a> 
-<a class="jxr_linenumber" name="L230" href="#L230">230</a>             splitAdd(ys, facts, as);
-<a class="jxr_linenumber" name="L231" href="#L231">231</a>             ys[0] = as[0]; ys[1] = as[1];
-<a class="jxr_linenumber" name="L232" href="#L232">232</a>         }
-<a class="jxr_linenumber" name="L233" href="#L233">233</a> 
-<a class="jxr_linenumber" name="L234" href="#L234">234</a>         <strong class="jxr_keyword">if</strong> (result != <strong class="jxr_keyword">null</strong>) {
-<a class="jxr_linenumber" name="L235" href="#L235">235</a>             result[0] = ys[0];
-<a class="jxr_linenumber" name="L236" href="#L236">236</a>             result[1] = ys[1];
-<a class="jxr_linenumber" name="L237" href="#L237">237</a>         }
-<a class="jxr_linenumber" name="L238" href="#L238">238</a> 
-<a class="jxr_linenumber" name="L239" href="#L239">239</a>         <strong class="jxr_keyword">return</strong> ys[0] + ys[1];
-<a class="jxr_linenumber" name="L240" href="#L240">240</a>     }
-<a class="jxr_linenumber" name="L241" href="#L241">241</a> 
-<a class="jxr_linenumber" name="L242" href="#L242">242</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L243" href="#L243">243</a> <em class="jxr_javadoccomment">     * For x between 0 and pi/4 compute sine using Taylor expansion:</em>
-<a class="jxr_linenumber" name="L244" href="#L244">244</a> <em class="jxr_javadoccomment">     * sin(x) = x - x^3/3! + x^5/5! - x^7/7! ...</em>
-<a class="jxr_linenumber" name="L245" href="#L245">245</a> <em class="jxr_javadoccomment">     * @param x number from which sine is requested</em>
-<a class="jxr_linenumber" name="L246" href="#L246">246</a> <em class="jxr_javadoccomment">     * @param result placeholder where to put the result in extended precision</em>
-<a class="jxr_linenumber" name="L247" href="#L247">247</a> <em class="jxr_javadoccomment">     * (may be null)</em>
-<a class="jxr_linenumber" name="L248" href="#L248">248</a> <em class="jxr_javadoccomment">     * @return sin(x)</em>
-<a class="jxr_linenumber" name="L249" href="#L249">249</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L250" href="#L250">250</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">double</strong> slowSin(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[]) {
-<a class="jxr_linenumber" name="L251" href="#L251">251</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L252" href="#L252">252</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L253" href="#L253">253</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> facts[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L254" href="#L254">254</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> as[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L255" href="#L255">255</a>         split(x, xs);
-<a class="jxr_linenumber" name="L256" href="#L256">256</a>         ys[0] = ys[1] = 0.0;
-<a class="jxr_linenumber" name="L257" href="#L257">257</a> 
-<a class="jxr_linenumber" name="L258" href="#L258">258</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = FACT.length-1; i &gt;= 0; i--) {
-<a class="jxr_linenumber" name="L259" href="#L259">259</a>             splitMult(xs, ys, as);
-<a class="jxr_linenumber" name="L260" href="#L260">260</a>             ys[0] = as[0]; ys[1] = as[1];
-<a class="jxr_linenumber" name="L261" href="#L261">261</a> 
-<a class="jxr_linenumber" name="L262" href="#L262">262</a>             <strong class="jxr_keyword">if</strong> ( (i &amp; 1) == 0) { <em class="jxr_comment">// Ignore even numbers</em>
-<a class="jxr_linenumber" name="L263" href="#L263">263</a>                 <strong class="jxr_keyword">continue</strong>;
-<a class="jxr_linenumber" name="L264" href="#L264">264</a>             }
-<a class="jxr_linenumber" name="L265" href="#L265">265</a> 
-<a class="jxr_linenumber" name="L266" href="#L266">266</a>             split(FACT[i], as);
-<a class="jxr_linenumber" name="L267" href="#L267">267</a>             splitReciprocal(as, facts);
-<a class="jxr_linenumber" name="L268" href="#L268">268</a> 
-<a class="jxr_linenumber" name="L269" href="#L269">269</a>             <strong class="jxr_keyword">if</strong> ( (i &amp; 2) != 0 ) { <em class="jxr_comment">// alternate terms are negative</em>
-<a class="jxr_linenumber" name="L270" href="#L270">270</a>                 facts[0] = -facts[0];
-<a class="jxr_linenumber" name="L271" href="#L271">271</a>                 facts[1] = -facts[1];
-<a class="jxr_linenumber" name="L272" href="#L272">272</a>             }
-<a class="jxr_linenumber" name="L273" href="#L273">273</a> 
-<a class="jxr_linenumber" name="L274" href="#L274">274</a>             splitAdd(ys, facts, as);
-<a class="jxr_linenumber" name="L275" href="#L275">275</a>             ys[0] = as[0]; ys[1] = as[1];
-<a class="jxr_linenumber" name="L276" href="#L276">276</a>         }
-<a class="jxr_linenumber" name="L277" href="#L277">277</a> 
-<a class="jxr_linenumber" name="L278" href="#L278">278</a>         <strong class="jxr_keyword">if</strong> (result != <strong class="jxr_keyword">null</strong>) {
-<a class="jxr_linenumber" name="L279" href="#L279">279</a>             result[0] = ys[0];
-<a class="jxr_linenumber" name="L280" href="#L280">280</a>             result[1] = ys[1];
-<a class="jxr_linenumber" name="L281" href="#L281">281</a>         }
-<a class="jxr_linenumber" name="L282" href="#L282">282</a> 
-<a class="jxr_linenumber" name="L283" href="#L283">283</a>         <strong class="jxr_keyword">return</strong> ys[0] + ys[1];
-<a class="jxr_linenumber" name="L284" href="#L284">284</a>     }
-<a class="jxr_linenumber" name="L285" href="#L285">285</a> 
-<a class="jxr_linenumber" name="L286" href="#L286">286</a> 
-<a class="jxr_linenumber" name="L287" href="#L287">287</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L288" href="#L288">288</a> <em class="jxr_javadoccomment">     *  For x between 0 and 1, returns exp(x), uses extended precision</em>
-<a class="jxr_linenumber" name="L289" href="#L289">289</a> <em class="jxr_javadoccomment">     *  @param x argument of exponential</em>
-<a class="jxr_linenumber" name="L290" href="#L290">290</a> <em class="jxr_javadoccomment">     *  @param result placeholder where to place exp(x) split in two terms</em>
-<a class="jxr_linenumber" name="L291" href="#L291">291</a> <em class="jxr_javadoccomment">     *  for extra precision (i.e. exp(x) = result[0] + result[1]</em>
-<a class="jxr_linenumber" name="L292" href="#L292">292</a> <em class="jxr_javadoccomment">     *  @return exp(x)</em>
-<a class="jxr_linenumber" name="L293" href="#L293">293</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L294" href="#L294">294</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">double</strong> slowexp(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[]) {
-<a class="jxr_linenumber" name="L295" href="#L295">295</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L296" href="#L296">296</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L297" href="#L297">297</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> facts[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L298" href="#L298">298</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> as[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L299" href="#L299">299</a>         split(x, xs);
-<a class="jxr_linenumber" name="L300" href="#L300">300</a>         ys[0] = ys[1] = 0.0;
-<a class="jxr_linenumber" name="L301" href="#L301">301</a> 
-<a class="jxr_linenumber" name="L302" href="#L302">302</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = FACT.length-1; i &gt;= 0; i--) {
-<a class="jxr_linenumber" name="L303" href="#L303">303</a>             splitMult(xs, ys, as);
-<a class="jxr_linenumber" name="L304" href="#L304">304</a>             ys[0] = as[0];
-<a class="jxr_linenumber" name="L305" href="#L305">305</a>             ys[1] = as[1];
-<a class="jxr_linenumber" name="L306" href="#L306">306</a> 
-<a class="jxr_linenumber" name="L307" href="#L307">307</a>             split(FACT[i], as);
-<a class="jxr_linenumber" name="L308" href="#L308">308</a>             splitReciprocal(as, facts);
-<a class="jxr_linenumber" name="L309" href="#L309">309</a> 
-<a class="jxr_linenumber" name="L310" href="#L310">310</a>             splitAdd(ys, facts, as);
-<a class="jxr_linenumber" name="L311" href="#L311">311</a>             ys[0] = as[0];
-<a class="jxr_linenumber" name="L312" href="#L312">312</a>             ys[1] = as[1];
-<a class="jxr_linenumber" name="L313" href="#L313">313</a>         }
-<a class="jxr_linenumber" name="L314" href="#L314">314</a> 
-<a class="jxr_linenumber" name="L315" href="#L315">315</a>         <strong class="jxr_keyword">if</strong> (result != <strong class="jxr_keyword">null</strong>) {
-<a class="jxr_linenumber" name="L316" href="#L316">316</a>             result[0] = ys[0];
-<a class="jxr_linenumber" name="L317" href="#L317">317</a>             result[1] = ys[1];
-<a class="jxr_linenumber" name="L318" href="#L318">318</a>         }
-<a class="jxr_linenumber" name="L319" href="#L319">319</a> 
-<a class="jxr_linenumber" name="L320" href="#L320">320</a>         <strong class="jxr_keyword">return</strong> ys[0] + ys[1];
-<a class="jxr_linenumber" name="L321" href="#L321">321</a>     }
-<a class="jxr_linenumber" name="L322" href="#L322">322</a> 
-<a class="jxr_linenumber" name="L323" href="#L323">323</a>     <em class="jxr_javadoccomment">/** Compute split[0], split[1] such that their sum is equal to d,</em>
-<a class="jxr_linenumber" name="L324" href="#L324">324</a> <em class="jxr_javadoccomment">     * and split[0] has its 30 least significant bits as zero.</em>
-<a class="jxr_linenumber" name="L325" href="#L325">325</a> <em class="jxr_javadoccomment">     * @param d number to split</em>
-<a class="jxr_linenumber" name="L326" href="#L326">326</a> <em class="jxr_javadoccomment">     * @param split placeholder where to place the result</em>
-<a class="jxr_linenumber" name="L327" href="#L327">327</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L328" href="#L328">328</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> split(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> d, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> split[]) {
-<a class="jxr_linenumber" name="L329" href="#L329">329</a>         <strong class="jxr_keyword">if</strong> (d &lt; 8e298 &amp;&amp; d &gt; -8e298) {
-<a class="jxr_linenumber" name="L330" href="#L330">330</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a = d * HEX_40000000;
-<a class="jxr_linenumber" name="L331" href="#L331">331</a>             split[0] = (d + a) - a;
-<a class="jxr_linenumber" name="L332" href="#L332">332</a>             split[1] = d - split[0];
-<a class="jxr_linenumber" name="L333" href="#L333">333</a>         } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L334" href="#L334">334</a>             <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a = d * 9.31322574615478515625E-10;
-<a class="jxr_linenumber" name="L335" href="#L335">335</a>             split[0] = (d + a - d) * HEX_40000000;
-<a class="jxr_linenumber" name="L336" href="#L336">336</a>             split[1] = d - split[0];
-<a class="jxr_linenumber" name="L337" href="#L337">337</a>         }
-<a class="jxr_linenumber" name="L338" href="#L338">338</a>     }
-<a class="jxr_linenumber" name="L339" href="#L339">339</a> 
-<a class="jxr_linenumber" name="L340" href="#L340">340</a>     <em class="jxr_javadoccomment">/** Recompute a split.</em>
-<a class="jxr_linenumber" name="L341" href="#L341">341</a> <em class="jxr_javadoccomment">     * @param a input/out array containing the split, changed</em>
-<a class="jxr_linenumber" name="L342" href="#L342">342</a> <em class="jxr_javadoccomment">     * on output</em>
-<a class="jxr_linenumber" name="L343" href="#L343">343</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L344" href="#L344">344</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> resplit(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a[]) {
-<a class="jxr_linenumber" name="L345" href="#L345">345</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> c = a[0] + a[1];
-<a class="jxr_linenumber" name="L346" href="#L346">346</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> d = -(c - a[0] - a[1]);
-<a class="jxr_linenumber" name="L347" href="#L347">347</a> 
-<a class="jxr_linenumber" name="L348" href="#L348">348</a>         <strong class="jxr_keyword">if</strong> (c &lt; 8e298 &amp;&amp; c &gt; -8e298) { <em class="jxr_comment">// MAGIC NUMBER</em>
-<a class="jxr_linenumber" name="L349" href="#L349">349</a>             <strong class="jxr_keyword">double</strong> z = c * HEX_40000000;
-<a class="jxr_linenumber" name="L350" href="#L350">350</a>             a[0] = (c + z) - z;
-<a class="jxr_linenumber" name="L351" href="#L351">351</a>             a[1] = c - a[0] + d;
-<a class="jxr_linenumber" name="L352" href="#L352">352</a>         } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L353" href="#L353">353</a>             <strong class="jxr_keyword">double</strong> z = c * 9.31322574615478515625E-10;
-<a class="jxr_linenumber" name="L354" href="#L354">354</a>             a[0] = (c + z - c) * HEX_40000000;
-<a class="jxr_linenumber" name="L355" href="#L355">355</a>             a[1] = c - a[0] + d;
-<a class="jxr_linenumber" name="L356" href="#L356">356</a>         }
-<a class="jxr_linenumber" name="L357" href="#L357">357</a>     }
-<a class="jxr_linenumber" name="L358" href="#L358">358</a> 
-<a class="jxr_linenumber" name="L359" href="#L359">359</a>     <em class="jxr_javadoccomment">/** Multiply two numbers in split form.</em>
-<a class="jxr_linenumber" name="L360" href="#L360">360</a> <em class="jxr_javadoccomment">     * @param a first term of multiplication</em>
-<a class="jxr_linenumber" name="L361" href="#L361">361</a> <em class="jxr_javadoccomment">     * @param b second term of multiplication</em>
-<a class="jxr_linenumber" name="L362" href="#L362">362</a> <em class="jxr_javadoccomment">     * @param ans placeholder where to put the result</em>
-<a class="jxr_linenumber" name="L363" href="#L363">363</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L364" href="#L364">364</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> splitMult(<strong class="jxr_keyword">double</strong> a[], <strong class="jxr_keyword">double</strong> b[], <strong class="jxr_keyword">double</strong> ans[]) {
-<a class="jxr_linenumber" name="L365" href="#L365">365</a>         ans[0] = a[0] * b[0];
-<a class="jxr_linenumber" name="L366" href="#L366">366</a>         ans[1] = a[0] * b[1] + a[1] * b[0] + a[1] * b[1];
-<a class="jxr_linenumber" name="L367" href="#L367">367</a> 
-<a class="jxr_linenumber" name="L368" href="#L368">368</a>         <em class="jxr_comment">/* Resplit */</em>
-<a class="jxr_linenumber" name="L369" href="#L369">369</a>         resplit(ans);
-<a class="jxr_linenumber" name="L370" href="#L370">370</a>     }
-<a class="jxr_linenumber" name="L371" href="#L371">371</a> 
-<a class="jxr_linenumber" name="L372" href="#L372">372</a>     <em class="jxr_javadoccomment">/** Add two numbers in split form.</em>
-<a class="jxr_linenumber" name="L373" href="#L373">373</a> <em class="jxr_javadoccomment">     * @param a first term of addition</em>
-<a class="jxr_linenumber" name="L374" href="#L374">374</a> <em class="jxr_javadoccomment">     * @param b second term of addition</em>
-<a class="jxr_linenumber" name="L375" href="#L375">375</a> <em class="jxr_javadoccomment">     * @param ans placeholder where to put the result</em>
-<a class="jxr_linenumber" name="L376" href="#L376">376</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L377" href="#L377">377</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> splitAdd(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a[], <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> b[], <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> ans[]) {
-<a class="jxr_linenumber" name="L378" href="#L378">378</a>         ans[0] = a[0] + b[0];
-<a class="jxr_linenumber" name="L379" href="#L379">379</a>         ans[1] = a[1] + b[1];
-<a class="jxr_linenumber" name="L380" href="#L380">380</a> 
-<a class="jxr_linenumber" name="L381" href="#L381">381</a>         resplit(ans);
-<a class="jxr_linenumber" name="L382" href="#L382">382</a>     }
-<a class="jxr_linenumber" name="L383" href="#L383">383</a> 
-<a class="jxr_linenumber" name="L384" href="#L384">384</a>     <em class="jxr_javadoccomment">/** Compute the reciprocal of in.  Use the following algorithm.</em>
-<a class="jxr_linenumber" name="L385" href="#L385">385</a> <em class="jxr_javadoccomment">     *  in = c + d.</em>
-<a class="jxr_linenumber" name="L386" href="#L386">386</a> <em class="jxr_javadoccomment">     *  want to find x + y such that x+y = 1/(c+d) and x is much</em>
-<a class="jxr_linenumber" name="L387" href="#L387">387</a> <em class="jxr_javadoccomment">     *  larger than y and x has several zero bits on the right.</em>
-<a class="jxr_linenumber" name="L388" href="#L388">388</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L389" href="#L389">389</a> <em class="jxr_javadoccomment">     *  Set b = 1/(2^22),  a = 1 - b.  Thus (a+b) = 1.</em>
-<a class="jxr_linenumber" name="L390" href="#L390">390</a> <em class="jxr_javadoccomment">     *  Use following identity to compute (a+b)/(c+d)</em>
-<a class="jxr_linenumber" name="L391" href="#L391">391</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L392" href="#L392">392</a> <em class="jxr_javadoccomment">     *  (a+b)/(c+d)  =   a/c   +    (bc - ad) / (c^2 + cd)</em>
-<a class="jxr_linenumber" name="L393" href="#L393">393</a> <em class="jxr_javadoccomment">     *  set x = a/c  and y = (bc - ad) / (c^2 + cd)</em>
-<a class="jxr_linenumber" name="L394" href="#L394">394</a> <em class="jxr_javadoccomment">     *  This will be close to the right answer, but there will be</em>
-<a class="jxr_linenumber" name="L395" href="#L395">395</a> <em class="jxr_javadoccomment">     *  some rounding in the calculation of X.  So by carefully</em>
-<a class="jxr_linenumber" name="L396" href="#L396">396</a> <em class="jxr_javadoccomment">     *  computing 1 - (c+d)(x+y) we can compute an error and</em>
-<a class="jxr_linenumber" name="L397" href="#L397">397</a> <em class="jxr_javadoccomment">     *  add that back in.   This is done carefully so that terms</em>
-<a class="jxr_linenumber" name="L398" href="#L398">398</a> <em class="jxr_javadoccomment">     *  of similar size are subtracted first.</em>
-<a class="jxr_linenumber" name="L399" href="#L399">399</a> <em class="jxr_javadoccomment">     *  @param in initial number, in split form</em>
-<a class="jxr_linenumber" name="L400" href="#L400">400</a> <em class="jxr_javadoccomment">     *  @param result placeholder where to put the result</em>
-<a class="jxr_linenumber" name="L401" href="#L401">401</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L402" href="#L402">402</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> splitReciprocal(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> in[], <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[]) {
-<a class="jxr_linenumber" name="L403" href="#L403">403</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> b = 1.0/4194304.0;
-<a class="jxr_linenumber" name="L404" href="#L404">404</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a = 1.0 - b;
-<a class="jxr_linenumber" name="L405" href="#L405">405</a> 
-<a class="jxr_linenumber" name="L406" href="#L406">406</a>         <strong class="jxr_keyword">if</strong> (in[0] == 0.0) {
-<a class="jxr_linenumber" name="L407" href="#L407">407</a>             in[0] = in[1];
-<a class="jxr_linenumber" name="L408" href="#L408">408</a>             in[1] = 0.0;
-<a class="jxr_linenumber" name="L409" href="#L409">409</a>         }
-<a class="jxr_linenumber" name="L410" href="#L410">410</a> 
-<a class="jxr_linenumber" name="L411" href="#L411">411</a>         result[0] = a / in[0];
-<a class="jxr_linenumber" name="L412" href="#L412">412</a>         result[1] = (b*in[0]-a*in[1]) / (in[0]*in[0] + in[0]*in[1]);
-<a class="jxr_linenumber" name="L413" href="#L413">413</a> 
-<a class="jxr_linenumber" name="L414" href="#L414">414</a>         <strong class="jxr_keyword">if</strong> (result[1] != result[1]) { <em class="jxr_comment">// can happen if result[1] is NAN</em>
-<a class="jxr_linenumber" name="L415" href="#L415">415</a>             result[1] = 0.0;
-<a class="jxr_linenumber" name="L416" href="#L416">416</a>         }
-<a class="jxr_linenumber" name="L417" href="#L417">417</a> 
-<a class="jxr_linenumber" name="L418" href="#L418">418</a>         <em class="jxr_comment">/* Resplit */</em>
-<a class="jxr_linenumber" name="L419" href="#L419">419</a>         resplit(result);
-<a class="jxr_linenumber" name="L420" href="#L420">420</a> 
-<a class="jxr_linenumber" name="L421" href="#L421">421</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; 2; i++) {
-<a class="jxr_linenumber" name="L422" href="#L422">422</a>             <em class="jxr_comment">/* this may be overkill, probably once is enough */</em>
-<a class="jxr_linenumber" name="L423" href="#L423">423</a>             <strong class="jxr_keyword">double</strong> err = 1.0 - result[0] * in[0] - result[0] * in[1] -
-<a class="jxr_linenumber" name="L424" href="#L424">424</a>             result[1] * in[0] - result[1] * in[1];
-<a class="jxr_linenumber" name="L425" href="#L425">425</a>             <em class="jxr_comment">/*err = 1.0 - err; */</em>
-<a class="jxr_linenumber" name="L426" href="#L426">426</a>             err *= result[0] + result[1];
-<a class="jxr_linenumber" name="L427" href="#L427">427</a>             <em class="jxr_comment">/*printf("err = %16e\n", err); */</em>
-<a class="jxr_linenumber" name="L428" href="#L428">428</a>             result[1] += err;
-<a class="jxr_linenumber" name="L429" href="#L429">429</a>         }
-<a class="jxr_linenumber" name="L430" href="#L430">430</a>     }
-<a class="jxr_linenumber" name="L431" href="#L431">431</a> 
-<a class="jxr_linenumber" name="L432" href="#L432">432</a>     <em class="jxr_javadoccomment">/** Compute (a[0] + a[1]) * (b[0] + b[1]) in extended precision.</em>
-<a class="jxr_linenumber" name="L433" href="#L433">433</a> <em class="jxr_javadoccomment">     * @param a first term of the multiplication</em>
-<a class="jxr_linenumber" name="L434" href="#L434">434</a> <em class="jxr_javadoccomment">     * @param b second term of the multiplication</em>
-<a class="jxr_linenumber" name="L435" href="#L435">435</a> <em class="jxr_javadoccomment">     * @param result placeholder where to put the result</em>
-<a class="jxr_linenumber" name="L436" href="#L436">436</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L437" href="#L437">437</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> quadMult(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a[], <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> b[], <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[]) {
-<a class="jxr_linenumber" name="L438" href="#L438">438</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L439" href="#L439">439</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L440" href="#L440">440</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> zs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L441" href="#L441">441</a> 
-<a class="jxr_linenumber" name="L442" href="#L442">442</a>         <em class="jxr_comment">/* a[0] * b[0] */</em>
-<a class="jxr_linenumber" name="L443" href="#L443">443</a>         split(a[0], xs);
-<a class="jxr_linenumber" name="L444" href="#L444">444</a>         split(b[0], ys);
-<a class="jxr_linenumber" name="L445" href="#L445">445</a>         splitMult(xs, ys, zs);
-<a class="jxr_linenumber" name="L446" href="#L446">446</a> 
-<a class="jxr_linenumber" name="L447" href="#L447">447</a>         result[0] = zs[0];
-<a class="jxr_linenumber" name="L448" href="#L448">448</a>         result[1] = zs[1];
-<a class="jxr_linenumber" name="L449" href="#L449">449</a> 
-<a class="jxr_linenumber" name="L450" href="#L450">450</a>         <em class="jxr_comment">/* a[0] * b[1] */</em>
-<a class="jxr_linenumber" name="L451" href="#L451">451</a>         split(b[1], ys);
-<a class="jxr_linenumber" name="L452" href="#L452">452</a>         splitMult(xs, ys, zs);
-<a class="jxr_linenumber" name="L453" href="#L453">453</a> 
-<a class="jxr_linenumber" name="L454" href="#L454">454</a>         <strong class="jxr_keyword">double</strong> tmp = result[0] + zs[0];
-<a class="jxr_linenumber" name="L455" href="#L455">455</a>         result[1] -= tmp - result[0] - zs[0];
-<a class="jxr_linenumber" name="L456" href="#L456">456</a>         result[0] = tmp;
-<a class="jxr_linenumber" name="L457" href="#L457">457</a>         tmp = result[0] + zs[1];
-<a class="jxr_linenumber" name="L458" href="#L458">458</a>         result[1] -= tmp - result[0] - zs[1];
-<a class="jxr_linenumber" name="L459" href="#L459">459</a>         result[0] = tmp;
-<a class="jxr_linenumber" name="L460" href="#L460">460</a> 
-<a class="jxr_linenumber" name="L461" href="#L461">461</a>         <em class="jxr_comment">/* a[1] * b[0] */</em>
-<a class="jxr_linenumber" name="L462" href="#L462">462</a>         split(a[1], xs);
-<a class="jxr_linenumber" name="L463" href="#L463">463</a>         split(b[0], ys);
-<a class="jxr_linenumber" name="L464" href="#L464">464</a>         splitMult(xs, ys, zs);
-<a class="jxr_linenumber" name="L465" href="#L465">465</a> 
-<a class="jxr_linenumber" name="L466" href="#L466">466</a>         tmp = result[0] + zs[0];
-<a class="jxr_linenumber" name="L467" href="#L467">467</a>         result[1] -= tmp - result[0] - zs[0];
-<a class="jxr_linenumber" name="L468" href="#L468">468</a>         result[0] = tmp;
-<a class="jxr_linenumber" name="L469" href="#L469">469</a>         tmp = result[0] + zs[1];
-<a class="jxr_linenumber" name="L470" href="#L470">470</a>         result[1] -= tmp - result[0] - zs[1];
-<a class="jxr_linenumber" name="L471" href="#L471">471</a>         result[0] = tmp;
-<a class="jxr_linenumber" name="L472" href="#L472">472</a> 
-<a class="jxr_linenumber" name="L473" href="#L473">473</a>         <em class="jxr_comment">/* a[1] * b[0] */</em>
-<a class="jxr_linenumber" name="L474" href="#L474">474</a>         split(a[1], xs);
-<a class="jxr_linenumber" name="L475" href="#L475">475</a>         split(b[1], ys);
-<a class="jxr_linenumber" name="L476" href="#L476">476</a>         splitMult(xs, ys, zs);
-<a class="jxr_linenumber" name="L477" href="#L477">477</a> 
-<a class="jxr_linenumber" name="L478" href="#L478">478</a>         tmp = result[0] + zs[0];
-<a class="jxr_linenumber" name="L479" href="#L479">479</a>         result[1] -= tmp - result[0] - zs[0];
-<a class="jxr_linenumber" name="L480" href="#L480">480</a>         result[0] = tmp;
-<a class="jxr_linenumber" name="L481" href="#L481">481</a>         tmp = result[0] + zs[1];
-<a class="jxr_linenumber" name="L482" href="#L482">482</a>         result[1] -= tmp - result[0] - zs[1];
-<a class="jxr_linenumber" name="L483" href="#L483">483</a>         result[0] = tmp;
-<a class="jxr_linenumber" name="L484" href="#L484">484</a>     }
-<a class="jxr_linenumber" name="L485" href="#L485">485</a> 
-<a class="jxr_linenumber" name="L486" href="#L486">486</a>     <em class="jxr_javadoccomment">/** Compute exp(p) for a integer p in extended precision.</em>
-<a class="jxr_linenumber" name="L487" href="#L487">487</a> <em class="jxr_javadoccomment">     * @param p integer whose exponential is requested</em>
-<a class="jxr_linenumber" name="L488" href="#L488">488</a> <em class="jxr_javadoccomment">     * @param result placeholder where to put the result in extended precision</em>
-<a class="jxr_linenumber" name="L489" href="#L489">489</a> <em class="jxr_javadoccomment">     * @return exp(p) in standard precision (equal to result[0] + result[1])</em>
-<a class="jxr_linenumber" name="L490" href="#L490">490</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L491" href="#L491">491</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">double</strong> expint(<strong class="jxr_keyword">int</strong> p, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> result[]) {
-<a class="jxr_linenumber" name="L492" href="#L492">492</a>         <em class="jxr_comment">//double x = M_E;</em>
-<a class="jxr_linenumber" name="L493" href="#L493">493</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> xs[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L494" href="#L494">494</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> as[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L495" href="#L495">495</a>         <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> ys[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L496" href="#L496">496</a>         <em class="jxr_comment">//split(x, xs);</em>
-<a class="jxr_linenumber" name="L497" href="#L497">497</a>         <em class="jxr_comment">//xs[1] = (double)(2.7182818284590452353602874713526625L - xs[0]);</em>
-<a class="jxr_linenumber" name="L498" href="#L498">498</a>         <em class="jxr_comment">//xs[0] = 2.71827697753906250000;</em>
-<a class="jxr_linenumber" name="L499" href="#L499">499</a>         <em class="jxr_comment">//xs[1] = 4.85091998273542816811e-06;</em>
-<a class="jxr_linenumber" name="L500" href="#L500">500</a>         <em class="jxr_comment">//xs[0] = Double.longBitsToDouble(0x4005bf0800000000L);</em>
-<a class="jxr_linenumber" name="L501" href="#L501">501</a>         <em class="jxr_comment">//xs[1] = Double.longBitsToDouble(0x3ed458a2bb4a9b00L);</em>
-<a class="jxr_linenumber" name="L502" href="#L502">502</a> 
-<a class="jxr_linenumber" name="L503" href="#L503">503</a>         <em class="jxr_comment">/* E */</em>
-<a class="jxr_linenumber" name="L504" href="#L504">504</a>         xs[0] = 2.718281828459045;
-<a class="jxr_linenumber" name="L505" href="#L505">505</a>         xs[1] = 1.4456468917292502E-16;
-<a class="jxr_linenumber" name="L506" href="#L506">506</a> 
-<a class="jxr_linenumber" name="L507" href="#L507">507</a>         split(1.0, ys);
-<a class="jxr_linenumber" name="L508" href="#L508">508</a> 
-<a class="jxr_linenumber" name="L509" href="#L509">509</a>         <strong class="jxr_keyword">while</strong> (p &gt; 0) {
-<a class="jxr_linenumber" name="L510" href="#L510">510</a>             <strong class="jxr_keyword">if</strong> ((p &amp; 1) != 0) {
-<a class="jxr_linenumber" name="L511" href="#L511">511</a>                 quadMult(ys, xs, as);
-<a class="jxr_linenumber" name="L512" href="#L512">512</a>                 ys[0] = as[0]; ys[1] = as[1];
-<a class="jxr_linenumber" name="L513" href="#L513">513</a>             }
-<a class="jxr_linenumber" name="L514" href="#L514">514</a> 
-<a class="jxr_linenumber" name="L515" href="#L515">515</a>             quadMult(xs, xs, as);
-<a class="jxr_linenumber" name="L516" href="#L516">516</a>             xs[0] = as[0]; xs[1] = as[1];
-<a class="jxr_linenumber" name="L517" href="#L517">517</a> 
-<a class="jxr_linenumber" name="L518" href="#L518">518</a>             p &gt;&gt;= 1;
-<a class="jxr_linenumber" name="L519" href="#L519">519</a>         }
-<a class="jxr_linenumber" name="L520" href="#L520">520</a> 
-<a class="jxr_linenumber" name="L521" href="#L521">521</a>         <strong class="jxr_keyword">if</strong> (result != <strong class="jxr_keyword">null</strong>) {
-<a class="jxr_linenumber" name="L522" href="#L522">522</a>             result[0] = ys[0];
-<a class="jxr_linenumber" name="L523" href="#L523">523</a>             result[1] = ys[1];
-<a class="jxr_linenumber" name="L524" href="#L524">524</a> 
-<a class="jxr_linenumber" name="L525" href="#L525">525</a>             resplit(result);
-<a class="jxr_linenumber" name="L526" href="#L526">526</a>         }
-<a class="jxr_linenumber" name="L527" href="#L527">527</a> 
-<a class="jxr_linenumber" name="L528" href="#L528">528</a>         <strong class="jxr_keyword">return</strong> ys[0] + ys[1];
-<a class="jxr_linenumber" name="L529" href="#L529">529</a>     }
-<a class="jxr_linenumber" name="L530" href="#L530">530</a>     <em class="jxr_javadoccomment">/** xi in the range of [1, 2].</em>
-<a class="jxr_linenumber" name="L531" href="#L531">531</a> <em class="jxr_javadoccomment">     *                                3        5        7</em>
-<a class="jxr_linenumber" name="L532" href="#L532">532</a> <em class="jxr_javadoccomment">     *      x+1           /          x        x        x          \</em>
-<a class="jxr_linenumber" name="L533" href="#L533">533</a> <em class="jxr_javadoccomment">     *  ln ----- =   2 *  |  x  +   ----  +  ----  +  ---- + ...  |</em>
-<a class="jxr_linenumber" name="L534" href="#L534">534</a> <em class="jxr_javadoccomment">     *      1-x           \          3        5        7          /</em>
-<a class="jxr_linenumber" name="L535" href="#L535">535</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L536" href="#L536">536</a> <em class="jxr_javadoccomment">     * So, compute a Remez approximation of the following function</em>
-<a class="jxr_linenumber" name="L537" href="#L537">537</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L538" href="#L538">538</a> <em class="jxr_javadoccomment">     *  ln ((sqrt(x)+1)/(1-sqrt(x)))  /  x</em>
-<a class="jxr_linenumber" name="L539" href="#L539">539</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L540" href="#L540">540</a> <em class="jxr_javadoccomment">     * This will be an even function with only positive coefficents.</em>
-<a class="jxr_linenumber" name="L541" href="#L541">541</a> <em class="jxr_javadoccomment">     * x is in the range [0 - 1/3].</em>
-<a class="jxr_linenumber" name="L542" href="#L542">542</a> <em class="jxr_javadoccomment">     *</em>
-<a class="jxr_linenumber" name="L543" href="#L543">543</a> <em class="jxr_javadoccomment">     * Transform xi for input to the above function by setting</em>
-<a class="jxr_linenumber" name="L544" href="#L544">544</a> <em class="jxr_javadoccomment">     * x = (xi-1)/(xi+1).   Input to the polynomial is x^2, then</em>
-<a class="jxr_linenumber" name="L545" href="#L545">545</a> <em class="jxr_javadoccomment">     * the result is multiplied by x.</em>
-<a class="jxr_linenumber" name="L546" href="#L546">546</a> <em class="jxr_javadoccomment">     * @param xi number from which log is requested</em>
-<a class="jxr_linenumber" name="L547" href="#L547">547</a> <em class="jxr_javadoccomment">     * @return log(xi)</em>
-<a class="jxr_linenumber" name="L548" href="#L548">548</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L549" href="#L549">549</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">double</strong>[] slowLog(<strong class="jxr_keyword">double</strong> xi) {
-<a class="jxr_linenumber" name="L550" href="#L550">550</a>         <strong class="jxr_keyword">double</strong> x[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L551" href="#L551">551</a>         <strong class="jxr_keyword">double</strong> x2[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L552" href="#L552">552</a>         <strong class="jxr_keyword">double</strong> y[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L553" href="#L553">553</a>         <strong class="jxr_keyword">double</strong> a[] = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[2];
-<a class="jxr_linenumber" name="L554" href="#L554">554</a> 
-<a class="jxr_linenumber" name="L555" href="#L555">555</a>         split(xi, x);
-<a class="jxr_linenumber" name="L556" href="#L556">556</a> 
-<a class="jxr_linenumber" name="L557" href="#L557">557</a>         <em class="jxr_comment">/* Set X = (x-1)/(x+1) */</em>
-<a class="jxr_linenumber" name="L558" href="#L558">558</a>         x[0] += 1.0;
-<a class="jxr_linenumber" name="L559" href="#L559">559</a>         resplit(x);
-<a class="jxr_linenumber" name="L560" href="#L560">560</a>         splitReciprocal(x, a);
-<a class="jxr_linenumber" name="L561" href="#L561">561</a>         x[0] -= 2.0;
-<a class="jxr_linenumber" name="L562" href="#L562">562</a>         resplit(x);
-<a class="jxr_linenumber" name="L563" href="#L563">563</a>         splitMult(x, a, y);
-<a class="jxr_linenumber" name="L564" href="#L564">564</a>         x[0] = y[0];
-<a class="jxr_linenumber" name="L565" href="#L565">565</a>         x[1] = y[1];
-<a class="jxr_linenumber" name="L566" href="#L566">566</a> 
-<a class="jxr_linenumber" name="L567" href="#L567">567</a>         <em class="jxr_comment">/* Square X -&gt; X2*/</em>
-<a class="jxr_linenumber" name="L568" href="#L568">568</a>         splitMult(x, x, x2);
-<a class="jxr_linenumber" name="L569" href="#L569">569</a> 
-<a class="jxr_linenumber" name="L570" href="#L570">570</a> 
-<a class="jxr_linenumber" name="L571" href="#L571">571</a>         <em class="jxr_comment">//x[0] -= 1.0;</em>
-<a class="jxr_linenumber" name="L572" href="#L572">572</a>         <em class="jxr_comment">//resplit(x);</em>
-<a class="jxr_linenumber" name="L573" href="#L573">573</a> 
-<a class="jxr_linenumber" name="L574" href="#L574">574</a>         y[0] = LN_SPLIT_COEF[LN_SPLIT_COEF.length-1][0];
-<a class="jxr_linenumber" name="L575" href="#L575">575</a>         y[1] = LN_SPLIT_COEF[LN_SPLIT_COEF.length-1][1];
-<a class="jxr_linenumber" name="L576" href="#L576">576</a> 
-<a class="jxr_linenumber" name="L577" href="#L577">577</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = LN_SPLIT_COEF.length-2; i &gt;= 0; i--) {
-<a class="jxr_linenumber" name="L578" href="#L578">578</a>             splitMult(y, x2, a);
-<a class="jxr_linenumber" name="L579" href="#L579">579</a>             y[0] = a[0];
-<a class="jxr_linenumber" name="L580" href="#L580">580</a>             y[1] = a[1];
-<a class="jxr_linenumber" name="L581" href="#L581">581</a>             splitAdd(y, LN_SPLIT_COEF[i], a);
-<a class="jxr_linenumber" name="L582" href="#L582">582</a>             y[0] = a[0];
-<a class="jxr_linenumber" name="L583" href="#L583">583</a>             y[1] = a[1];
-<a class="jxr_linenumber" name="L584" href="#L584">584</a>         }
-<a class="jxr_linenumber" name="L585" href="#L585">585</a> 
-<a class="jxr_linenumber" name="L586" href="#L586">586</a>         splitMult(y, x, a);
-<a class="jxr_linenumber" name="L587" href="#L587">587</a>         y[0] = a[0];
-<a class="jxr_linenumber" name="L588" href="#L588">588</a>         y[1] = a[1];
-<a class="jxr_linenumber" name="L589" href="#L589">589</a> 
-<a class="jxr_linenumber" name="L590" href="#L590">590</a>         <strong class="jxr_keyword">return</strong> y;
-<a class="jxr_linenumber" name="L591" href="#L591">591</a>     }
-<a class="jxr_linenumber" name="L592" href="#L592">592</a> 
-<a class="jxr_linenumber" name="L593" href="#L593">593</a> 
-<a class="jxr_linenumber" name="L594" href="#L594">594</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L595" href="#L595">595</a> <em class="jxr_javadoccomment">     * Print an array.</em>
-<a class="jxr_linenumber" name="L596" href="#L596">596</a> <em class="jxr_javadoccomment">     * @param out text output stream where output should be printed</em>
-<a class="jxr_linenumber" name="L597" href="#L597">597</a> <em class="jxr_javadoccomment">     * @param name array name</em>
-<a class="jxr_linenumber" name="L598" href="#L598">598</a> <em class="jxr_javadoccomment">     * @param expectedLen expected length of the array</em>
-<a class="jxr_linenumber" name="L599" href="#L599">599</a> <em class="jxr_javadoccomment">     * @param array2d array data</em>
-<a class="jxr_linenumber" name="L600" href="#L600">600</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L601" href="#L601">601</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> printarray(PrintStream out, String name, <strong class="jxr_keyword">int</strong> expectedLen, <strong class="jxr_keyword">double</strong>[][] array2d) {
-<a class="jxr_linenumber" name="L602" href="#L602">602</a>         out.println(name);
-<a class="jxr_linenumber" name="L603" href="#L603">603</a>         checkLen(expectedLen, array2d.length);
-<a class="jxr_linenumber" name="L604" href="#L604">604</a>         out.println(TABLE_START_DECL + <span class="jxr_string">" "</span>);
-<a class="jxr_linenumber" name="L605" href="#L605">605</a>         <strong class="jxr_keyword">int</strong> i = 0;
-<a class="jxr_linenumber" name="L606" href="#L606">606</a>         <strong class="jxr_keyword">for</strong>(<strong class="jxr_keyword">double</strong>[] array : array2d) { <em class="jxr_comment">// "double array[]" causes PMD parsing error</em>
-<a class="jxr_linenumber" name="L607" href="#L607">607</a>             out.print(<span class="jxr_string">"        {"</span>);
-<a class="jxr_linenumber" name="L608" href="#L608">608</a>             <strong class="jxr_keyword">for</strong>(<strong class="jxr_keyword">double</strong> d : array) { <em class="jxr_comment">// assume inner array has very few entries</em>
-<a class="jxr_linenumber" name="L609" href="#L609">609</a>                 out.printf(<span class="jxr_string">"%-25.25s"</span>, format(d)); <em class="jxr_comment">// multiple entries per line</em>
-<a class="jxr_linenumber" name="L610" href="#L610">610</a>             }
-<a class="jxr_linenumber" name="L611" href="#L611">611</a>             out.println(<span class="jxr_string">"}, // "</span> + i++);
-<a class="jxr_linenumber" name="L612" href="#L612">612</a>         }
-<a class="jxr_linenumber" name="L613" href="#L613">613</a>         out.println(TABLE_END_DECL);
-<a class="jxr_linenumber" name="L614" href="#L614">614</a>     }
-<a class="jxr_linenumber" name="L615" href="#L615">615</a> 
-<a class="jxr_linenumber" name="L616" href="#L616">616</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L617" href="#L617">617</a> <em class="jxr_javadoccomment">     * Print an array.</em>
-<a class="jxr_linenumber" name="L618" href="#L618">618</a> <em class="jxr_javadoccomment">     * @param out text output stream where output should be printed</em>
-<a class="jxr_linenumber" name="L619" href="#L619">619</a> <em class="jxr_javadoccomment">     * @param name array name</em>
-<a class="jxr_linenumber" name="L620" href="#L620">620</a> <em class="jxr_javadoccomment">     * @param expectedLen expected length of the array</em>
-<a class="jxr_linenumber" name="L621" href="#L621">621</a> <em class="jxr_javadoccomment">     * @param array array data</em>
-<a class="jxr_linenumber" name="L622" href="#L622">622</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L623" href="#L623">623</a>     <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> printarray(PrintStream out, String name, <strong class="jxr_keyword">int</strong> expectedLen, <strong class="jxr_keyword">double</strong>[] array) {
-<a class="jxr_linenumber" name="L624" href="#L624">624</a>         out.println(name + <span class="jxr_string">"="</span>);
-<a class="jxr_linenumber" name="L625" href="#L625">625</a>         checkLen(expectedLen, array.length);
-<a class="jxr_linenumber" name="L626" href="#L626">626</a>         out.println(TABLE_START_DECL);
-<a class="jxr_linenumber" name="L627" href="#L627">627</a>         <strong class="jxr_keyword">for</strong>(<strong class="jxr_keyword">double</strong> d : array){
-<a class="jxr_linenumber" name="L628" href="#L628">628</a>             out.printf(<span class="jxr_string">"        %s%n"</span>, format(d)); <em class="jxr_comment">// one entry per line</em>
-<a class="jxr_linenumber" name="L629" href="#L629">629</a>         }
-<a class="jxr_linenumber" name="L630" href="#L630">630</a>         out.println(TABLE_END_DECL);
-<a class="jxr_linenumber" name="L631" href="#L631">631</a>     }
-<a class="jxr_linenumber" name="L632" href="#L632">632</a> 
-<a class="jxr_linenumber" name="L633" href="#L633">633</a>     <em class="jxr_javadoccomment">/** Format a double.</em>
-<a class="jxr_linenumber" name="L634" href="#L634">634</a> <em class="jxr_javadoccomment">     * @param d double number to format</em>
-<a class="jxr_linenumber" name="L635" href="#L635">635</a> <em class="jxr_javadoccomment">     * @return formatted number</em>
-<a class="jxr_linenumber" name="L636" href="#L636">636</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L637" href="#L637">637</a>     <strong class="jxr_keyword">static</strong> String format(<strong class="jxr_keyword">double</strong> d) {
-<a class="jxr_linenumber" name="L638" href="#L638">638</a>         <strong class="jxr_keyword">if</strong> (d != d) {
-<a class="jxr_linenumber" name="L639" href="#L639">639</a>             <strong class="jxr_keyword">return</strong> <span class="jxr_string">"Double.NaN,"</span>;
-<a class="jxr_linenumber" name="L640" href="#L640">640</a>         } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L641" href="#L641">641</a>             <strong class="jxr_keyword">return</strong> ((d &gt;= 0) ? <span class="jxr_string">"+"</span> : <span class="jxr_string">""</span>) + Double.toString(d) + <span class="jxr_string">"d,"</span>;
-<a class="jxr_linenumber" name="L642" href="#L642">642</a>         }
-<a class="jxr_linenumber" name="L643" href="#L643">643</a>     }
-<a class="jxr_linenumber" name="L644" href="#L644">644</a> 
-<a class="jxr_linenumber" name="L645" href="#L645">645</a>     <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L646" href="#L646">646</a> <em class="jxr_javadoccomment">     * Check two lengths are equal.</em>
-<a class="jxr_linenumber" name="L647" href="#L647">647</a> <em class="jxr_javadoccomment">     * @param expectedLen expected length</em>
-<a class="jxr_linenumber" name="L648" href="#L648">648</a> <em class="jxr_javadoccomment">     * @param actual actual length</em>
-<a class="jxr_linenumber" name="L649" href="#L649">649</a> <em class="jxr_javadoccomment">     * @exception DimensionMismatchException if the two lengths are not equal</em>
-<a class="jxr_linenumber" name="L650" href="#L650">650</a> <em class="jxr_javadoccomment">     */</em>
-<a class="jxr_linenumber" name="L651" href="#L651">651</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">void</strong> checkLen(<strong class="jxr_keyword">int</strong> expectedLen, <strong class="jxr_keyword">int</strong> actual)
-<a class="jxr_linenumber" name="L652" href="#L652">652</a>         <strong class="jxr_keyword">throws</strong> <a href="../../../../../org/apache/commons/math3/exception/DimensionMismatchException.html">DimensionMismatchException</a> {
-<a class="jxr_linenumber" name="L653" href="#L653">653</a>         <strong class="jxr_keyword">if</strong> (expectedLen != actual) {
-<a class="jxr_linenumber" name="L654" href="#L654">654</a>             <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> <a href="../../../../../org/apache/commons/math3/exception/DimensionMismatchException.html">DimensionMismatchException</a>(actual, expectedLen);
-<a class="jxr_linenumber" name="L655" href="#L655">655</a>         }
-<a class="jxr_linenumber" name="L656" href="#L656">656</a>     }
-<a class="jxr_linenumber" name="L657" href="#L657">657</a> 
-<a class="jxr_linenumber" name="L658" href="#L658">658</a> }
-</pre>
-<hr/>
-<div id="footer">Copyright &#169; 2003&#x2013;2016 <a href="http://www.apache.org/">The Apache Software Foundation</a>. All rights reserved.</div>
-</body>
-</html>

http://git-wip-us.apache.org/repos/asf/commons-complex/blob/b3576eeb/site-content/.svn/pristine/01/015eab5accf046f00e97e246bfa9f061afd64447.svn-base
----------------------------------------------------------------------
diff --git a/site-content/.svn/pristine/01/015eab5accf046f00e97e246bfa9f061afd64447.svn-base b/site-content/.svn/pristine/01/015eab5accf046f00e97e246bfa9f061afd64447.svn-base
deleted file mode 100644
index 63ff742..0000000
--- a/site-content/.svn/pristine/01/015eab5accf046f00e97e246bfa9f061afd64447.svn-base
+++ /dev/null
@@ -1,38 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml" lang="en"><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8"/><link rel="stylesheet" href="../.resources/report.css" type="text/css"/><link rel="shortcut icon" href="../.resources/report.gif" type="image/gif"/><title>SingularMatrixException.java</title><link rel="stylesheet" href="../.resources/prettify.css" type="text/css"/><script type="text/javascript" src="../.resources/prettify.js"></script></head><body onload="window['PR_TAB_WIDTH']=4;prettyPrint()"><div class="breadcrumb" id="breadcrumb"><span class="info"><a href="../.sessions.html" class="el_session">Sessions</a></span><a href="../index.html" class="el_report">Apache Commons Math</a> &gt; <a href="index.source.html" class="el_package">org.apache.commons.math3.linear</a> &gt; <span class="el_source">SingularMatrixExce
 ption.java</span></div><h1>SingularMatrixException.java</h1><pre class="source lang-java linenums">/*
- * Licensed to the Apache Software Foundation (ASF) under one or more
- * contributor license agreements.  See the NOTICE file distributed with
- * this work for additional information regarding copyright ownership.
- * The ASF licenses this file to You under the Apache License, Version 2.0
- * (the &quot;License&quot;); you may not use this file except in compliance with
- * the License.  You may obtain a copy of the License at
- *
- *      http://www.apache.org/licenses/LICENSE-2.0
- *
- * Unless required by applicable law or agreed to in writing, software
- * distributed under the License is distributed on an &quot;AS IS&quot; BASIS,
- * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
- * See the License for the specific language governing permissions and
- * limitations under the License.
- */
-package org.apache.commons.math3.linear;
-
-import org.apache.commons.math3.exception.MathIllegalArgumentException;
-import org.apache.commons.math3.exception.util.LocalizedFormats;
-
-/**
- * Exception to be thrown when a non-singular matrix is expected.
- *
- * @since 3.0
- */
-public class SingularMatrixException extends MathIllegalArgumentException {
-    /** Serializable version Id. */
-    private static final long serialVersionUID = -4206514844735401070L;
-
-    /**
-     * Construct an exception.
-     */
-    public SingularMatrixException() {
-<span class="fc" id="L35">        super(LocalizedFormats.SINGULAR_MATRIX);</span>
-<span class="fc" id="L36">    }</span>
-}
-</pre><div class="footer"><span class="right">Created with <a href="http://www.eclemma.org/jacoco">JaCoCo</a> 0.7.5.201505241946</span></div></body></html>
\ No newline at end of file


Mime
View raw message