You are viewing a plain text version of this content. The canonical link for it is here.
Posted to notifications@commons.apache.org by lu...@apache.org on 2015/07/27 21:42:10 UTC
svn commit: r959803 [30/32] - in
/websites/production/commons/content/proper/commons-math: jacoco/
jacoco/org.apache.commons.math3.fraction/
jacoco/org.apache.commons.math3.genetics/
jacoco/org.apache.commons.math3.geometry.euclidean.threed/ jacoco/org...
Modified: websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/special/GammaTest.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/special/GammaTest.html (original)
+++ websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/special/GammaTest.html Mon Jul 27 19:42:09 2015
@@ -133,847 +133,861 @@
<a class="jxr_linenumber" name="L125" href="#L125">125</a> }
<a class="jxr_linenumber" name="L126" href="#L126">126</a>
<a class="jxr_linenumber" name="L127" href="#L127">127</a> @Test
-<a class="jxr_linenumber" name="L128" href="#L128">128</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testTrigamma() {
-<a class="jxr_linenumber" name="L129" href="#L129">129</a> <strong class="jxr_keyword">double</strong> eps = 1e-8;
-<a class="jxr_linenumber" name="L130" href="#L130">130</a> <em class="jxr_comment">// computed using webMathematica. For example, to compute trigamma($i) = Polygamma(1, $i), use</em>
-<a class="jxr_linenumber" name="L131" href="#L131">131</a> <em class="jxr_comment">//</em>
-<a class="jxr_linenumber" name="L132" href="#L132">132</a> <em class="jxr_comment">// http://functions.wolfram.com/webMathematica/Evaluated.jsp?name=PolyGamma2&plottype=0&vars={%221%22,%22$i%22}&digits=20</em>
-<a class="jxr_linenumber" name="L133" href="#L133">133</a> <strong class="jxr_keyword">double</strong>[] data = {
-<a class="jxr_linenumber" name="L134" href="#L134">134</a> 1e-4, 1.0000000164469368793e8,
-<a class="jxr_linenumber" name="L135" href="#L135">135</a> 1e-3, 1.0000016425331958690e6,
-<a class="jxr_linenumber" name="L136" href="#L136">136</a> 1e-2, 10001.621213528313220,
-<a class="jxr_linenumber" name="L137" href="#L137">137</a> 1e-1, 101.43329915079275882,
-<a class="jxr_linenumber" name="L138" href="#L138">138</a> 1, 1.6449340668482264365,
-<a class="jxr_linenumber" name="L139" href="#L139">139</a> 2, 0.64493406684822643647,
-<a class="jxr_linenumber" name="L140" href="#L140">140</a> 3, 0.39493406684822643647,
-<a class="jxr_linenumber" name="L141" href="#L141">141</a> 4, 0.28382295573711532536,
-<a class="jxr_linenumber" name="L142" href="#L142">142</a> 5, 0.22132295573711532536,
-<a class="jxr_linenumber" name="L143" href="#L143">143</a> 10, 0.10516633568168574612,
-<a class="jxr_linenumber" name="L144" href="#L144">144</a> 20, 0.051270822935203119832,
-<a class="jxr_linenumber" name="L145" href="#L145">145</a> 50, 0.020201333226697125806,
-<a class="jxr_linenumber" name="L146" href="#L146">146</a> 100, 0.010050166663333571395
-<a class="jxr_linenumber" name="L147" href="#L147">147</a> };
-<a class="jxr_linenumber" name="L148" href="#L148">148</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = data.length - 2; i >= 0; i -= 2) {
-<a class="jxr_linenumber" name="L149" href="#L149">149</a> Assert.assertEquals(String.format(<span class="jxr_string">"trigamma %.0f"</span>, data[i]), data[i + 1], Gamma.trigamma(data[i]), eps);
-<a class="jxr_linenumber" name="L150" href="#L150">150</a> }
-<a class="jxr_linenumber" name="L151" href="#L151">151</a> }
-<a class="jxr_linenumber" name="L152" href="#L152">152</a>
-<a class="jxr_linenumber" name="L153" href="#L153">153</a> <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L154" href="#L154">154</a> <em class="jxr_javadoccomment"> * Reference data for the {@link Gamma#logGamma(double)} function. This data</em>
-<a class="jxr_linenumber" name="L155" href="#L155">155</a> <em class="jxr_javadoccomment"> * was generated with the following <a</em>
-<a class="jxr_linenumber" name="L156" href="#L156">156</a> <em class="jxr_javadoccomment"> * href="<a href="http://maxima.sourceforge.net/" target="alexandria_uri">http://maxima.sourceforge.net/</a>">Maxima</a> script.</em>
-<a class="jxr_linenumber" name="L157" href="#L157">157</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L158" href="#L158">158</a> <em class="jxr_javadoccomment"> * <pre></em>
-<a class="jxr_linenumber" name="L159" href="#L159">159</a> <em class="jxr_javadoccomment"> * kill(all);</em>
-<a class="jxr_linenumber" name="L160" href="#L160">160</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L161" href="#L161">161</a> <em class="jxr_javadoccomment"> * fpprec : 64;</em>
-<a class="jxr_linenumber" name="L162" href="#L162">162</a> <em class="jxr_javadoccomment"> * gamln(x) := log(gamma(x));</em>
-<a class="jxr_linenumber" name="L163" href="#L163">163</a> <em class="jxr_javadoccomment"> * x : append(makelist(bfloat(i / 8), i, 1, 80),</em>
-<a class="jxr_linenumber" name="L164" href="#L164">164</a> <em class="jxr_javadoccomment"> * [0.8b0, 1b2, 1b3, 1b4, 1b5, 1b6, 1b7, 1b8, 1b9, 1b10]);</em>
-<a class="jxr_linenumber" name="L165" href="#L165">165</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L166" href="#L166">166</a> <em class="jxr_javadoccomment"> * for i : 1 while i <= length(x) do</em>
-<a class="jxr_linenumber" name="L167" href="#L167">167</a> <em class="jxr_javadoccomment"> * print("{", float(x[i]), ",", float(gamln(x[i])), "},");</em>
-<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_javadoccomment"> * </pre></em>
-<a class="jxr_linenumber" name="L169" href="#L169">169</a> <em class="jxr_javadoccomment"> */</em>
-<a class="jxr_linenumber" name="L170" href="#L170">170</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>[][] LOG_GAMMA_REF = {
-<a class="jxr_linenumber" name="L171" href="#L171">171</a> { 0.125 , 2.019418357553796 },
-<a class="jxr_linenumber" name="L172" href="#L172">172</a> { 0.25 , 1.288022524698077 },
-<a class="jxr_linenumber" name="L173" href="#L173">173</a> { 0.375 , .8630739822706475 },
-<a class="jxr_linenumber" name="L174" href="#L174">174</a> { 0.5 , .5723649429247001 },
-<a class="jxr_linenumber" name="L175" href="#L175">175</a> { 0.625 , .3608294954889402 },
-<a class="jxr_linenumber" name="L176" href="#L176">176</a> { 0.75 , .2032809514312954 },
-<a class="jxr_linenumber" name="L177" href="#L177">177</a> { 0.875 , .08585870722533433 },
-<a class="jxr_linenumber" name="L178" href="#L178">178</a> { 0.890625 , .07353860936979656 },
-<a class="jxr_linenumber" name="L179" href="#L179">179</a> { 0.90625 , .06169536624059108 },
-<a class="jxr_linenumber" name="L180" href="#L180">180</a> { 0.921875 , .05031670080005688 },
-<a class="jxr_linenumber" name="L181" href="#L181">181</a> { 0.9375 , 0.0393909017345823 },
-<a class="jxr_linenumber" name="L182" href="#L182">182</a> { 0.953125 , .02890678734595923 },
-<a class="jxr_linenumber" name="L183" href="#L183">183</a> { 0.96875 , .01885367233441289 },
-<a class="jxr_linenumber" name="L184" href="#L184">184</a> { 0.984375 , .009221337197578781 },
-<a class="jxr_linenumber" name="L185" href="#L185">185</a> { 1.0 , 0.0 },
-<a class="jxr_linenumber" name="L186" href="#L186">186</a> { 1.015625 , - 0.00881970970573307 },
-<a class="jxr_linenumber" name="L187" href="#L187">187</a> { 1.03125 , - .01724677500176807 },
-<a class="jxr_linenumber" name="L188" href="#L188">188</a> { 1.046875 , - .02528981394675729 },
-<a class="jxr_linenumber" name="L189" href="#L189">189</a> { 1.0625 , - .03295710029357782 },
-<a class="jxr_linenumber" name="L190" href="#L190">190</a> { 1.078125 , - .04025658272400143 },
-<a class="jxr_linenumber" name="L191" href="#L191">191</a> { 1.09375 , - .04719590272716985 },
-<a class="jxr_linenumber" name="L192" href="#L192">192</a> { 1.109375 , - .05378241123619192 },
-<a class="jxr_linenumber" name="L193" href="#L193">193</a> { 1.125 , - .06002318412603958 },
-<a class="jxr_linenumber" name="L194" href="#L194">194</a> { 1.25 , - .09827183642181316 },
-<a class="jxr_linenumber" name="L195" href="#L195">195</a> { 1.375 , - .1177552707410788 },
-<a class="jxr_linenumber" name="L196" href="#L196">196</a> { 1.5 , - .1207822376352452 },
-<a class="jxr_linenumber" name="L197" href="#L197">197</a> { 1.625 , - .1091741337567954 },
-<a class="jxr_linenumber" name="L198" href="#L198">198</a> { 1.75 , - .08440112102048555 },
-<a class="jxr_linenumber" name="L199" href="#L199">199</a> { 1.875 , - 0.0476726853991883 },
-<a class="jxr_linenumber" name="L200" href="#L200">200</a> { 1.890625 , - .04229320615532515 },
-<a class="jxr_linenumber" name="L201" href="#L201">201</a> { 1.90625 , - .03674470657266143 },
-<a class="jxr_linenumber" name="L202" href="#L202">202</a> { 1.921875 , - .03102893865389552 },
-<a class="jxr_linenumber" name="L203" href="#L203">203</a> { 1.9375 , - .02514761940298887 },
-<a class="jxr_linenumber" name="L204" href="#L204">204</a> { 1.953125 , - .01910243184040138 },
-<a class="jxr_linenumber" name="L205" href="#L205">205</a> { 1.96875 , - .01289502598016741 },
-<a class="jxr_linenumber" name="L206" href="#L206">206</a> { 1.984375 , - .006527019770560387 },
-<a class="jxr_linenumber" name="L207" href="#L207">207</a> { 2.0 , 0.0 },
-<a class="jxr_linenumber" name="L208" href="#L208">208</a> { 2.015625 , .006684476830232185 },
-<a class="jxr_linenumber" name="L209" href="#L209">209</a> { 2.03125 , .01352488366498562 },
-<a class="jxr_linenumber" name="L210" href="#L210">210</a> { 2.046875 , .02051972208453692 },
-<a class="jxr_linenumber" name="L211" href="#L211">211</a> { 2.0625 , .02766752152285702 },
-<a class="jxr_linenumber" name="L212" href="#L212">212</a> { 2.078125 , 0.0349668385135861 },
-<a class="jxr_linenumber" name="L213" href="#L213">213</a> { 2.09375 , .04241625596251728 },
-<a class="jxr_linenumber" name="L214" href="#L214">214</a> { 2.109375 , .05001438244545164 },
-<a class="jxr_linenumber" name="L215" href="#L215">215</a> { 2.125 , .05775985153034387 },
-<a class="jxr_linenumber" name="L216" href="#L216">216</a> { 2.25 , .1248717148923966 },
-<a class="jxr_linenumber" name="L217" href="#L217">217</a> { 2.375 , .2006984603774558 },
-<a class="jxr_linenumber" name="L218" href="#L218">218</a> { 2.5 , .2846828704729192 },
-<a class="jxr_linenumber" name="L219" href="#L219">219</a> { 2.625 , .3763336820249054 },
-<a class="jxr_linenumber" name="L220" href="#L220">220</a> { 2.75 , .4752146669149371 },
-<a class="jxr_linenumber" name="L221" href="#L221">221</a> { 2.875 , .5809359740231859 },
-<a class="jxr_linenumber" name="L222" href="#L222">222</a> { 2.890625 , .5946142560817441 },
-<a class="jxr_linenumber" name="L223" href="#L223">223</a> { 2.90625 , .6083932548009232 },
-<a class="jxr_linenumber" name="L224" href="#L224">224</a> { 2.921875 , .6222723333588501 },
-<a class="jxr_linenumber" name="L225" href="#L225">225</a> { 2.9375 , .6362508628423761 },
-<a class="jxr_linenumber" name="L226" href="#L226">226</a> { 2.953125 , .6503282221022278 },
-<a class="jxr_linenumber" name="L227" href="#L227">227</a> { 2.96875 , .6645037976116387 },
-<a class="jxr_linenumber" name="L228" href="#L228">228</a> { 2.984375 , 0.678776983328359 },
-<a class="jxr_linenumber" name="L229" href="#L229">229</a> { 3.0 , .6931471805599453 },
-<a class="jxr_linenumber" name="L230" href="#L230">230</a> { 3.015625 , .7076137978322324 },
-<a class="jxr_linenumber" name="L231" href="#L231">231</a> { 3.03125 , .7221762507608962 },
-<a class="jxr_linenumber" name="L232" href="#L232">232</a> { 3.046875 , .7368339619260166 },
-<a class="jxr_linenumber" name="L233" href="#L233">233</a> { 3.0625 , 0.751586360749556 },
-<a class="jxr_linenumber" name="L234" href="#L234">234</a> { 3.078125 , .7664328833756681 },
-<a class="jxr_linenumber" name="L235" href="#L235">235</a> { 3.09375 , .7813729725537568 },
-<a class="jxr_linenumber" name="L236" href="#L236">236</a> { 3.109375 , .7964060775242092 },
-<a class="jxr_linenumber" name="L237" href="#L237">237</a> { 3.125 , 0.811531653906724 },
-<a class="jxr_linenumber" name="L238" href="#L238">238</a> { 3.25 , .9358019311087253 },
-<a class="jxr_linenumber" name="L239" href="#L239">239</a> { 3.375 , 1.06569589786406 },
-<a class="jxr_linenumber" name="L240" href="#L240">240</a> { 3.5 , 1.200973602347074 },
-<a class="jxr_linenumber" name="L241" href="#L241">241</a> { 3.625 , 1.341414578068493 },
-<a class="jxr_linenumber" name="L242" href="#L242">242</a> { 3.75 , 1.486815578593417 },
-<a class="jxr_linenumber" name="L243" href="#L243">243</a> { 3.875 , 1.6369886482725 },
-<a class="jxr_linenumber" name="L244" href="#L244">244</a> { 4.0 , 1.791759469228055 },
-<a class="jxr_linenumber" name="L245" href="#L245">245</a> { 4.125 , 1.950965937095089 },
-<a class="jxr_linenumber" name="L246" href="#L246">246</a> { 4.25 , 2.114456927450371 },
-<a class="jxr_linenumber" name="L247" href="#L247">247</a> { 4.375 , 2.282091222188554 },
-<a class="jxr_linenumber" name="L248" href="#L248">248</a> { 4.5 , 2.453736570842442 },
-<a class="jxr_linenumber" name="L249" href="#L249">249</a> { 4.625 , 2.62926886637513 },
-<a class="jxr_linenumber" name="L250" href="#L250">250</a> { 4.75 , 2.808571418575736 },
-<a class="jxr_linenumber" name="L251" href="#L251">251</a> { 4.875 , 2.99153431107781 },
-<a class="jxr_linenumber" name="L252" href="#L252">252</a> { 5.0 , 3.178053830347946 },
-<a class="jxr_linenumber" name="L253" href="#L253">253</a> { 5.125 , 3.368031956881733 },
-<a class="jxr_linenumber" name="L254" href="#L254">254</a> { 5.25 , 3.561375910386697 },
-<a class="jxr_linenumber" name="L255" href="#L255">255</a> { 5.375 , 3.757997741998131 },
-<a class="jxr_linenumber" name="L256" href="#L256">256</a> { 5.5 , 3.957813967618717 },
-<a class="jxr_linenumber" name="L257" href="#L257">257</a> { 5.625 , 4.160745237339519 },
-<a class="jxr_linenumber" name="L258" href="#L258">258</a> { 5.75 , 4.366716036622286 },
-<a class="jxr_linenumber" name="L259" href="#L259">259</a> { 5.875 , 4.57565441552762 },
-<a class="jxr_linenumber" name="L260" href="#L260">260</a> { 6.0 , 4.787491742782046 },
-<a class="jxr_linenumber" name="L261" href="#L261">261</a> { 6.125 , 5.002162481906205 },
-<a class="jxr_linenumber" name="L262" href="#L262">262</a> { 6.25 , 5.219603986990229 },
-<a class="jxr_linenumber" name="L263" href="#L263">263</a> { 6.375 , 5.439756316011858 },
-<a class="jxr_linenumber" name="L264" href="#L264">264</a> { 6.5 , 5.662562059857142 },
-<a class="jxr_linenumber" name="L265" href="#L265">265</a> { 6.625 , 5.887966185430003 },
-<a class="jxr_linenumber" name="L266" href="#L266">266</a> { 6.75 , 6.115915891431546 },
-<a class="jxr_linenumber" name="L267" href="#L267">267</a> { 6.875 , 6.346360475557843 },
-<a class="jxr_linenumber" name="L268" href="#L268">268</a> { 7.0 , 6.579251212010101 },
-<a class="jxr_linenumber" name="L269" href="#L269">269</a> { 7.125 , 6.814541238336996 },
-<a class="jxr_linenumber" name="L270" href="#L270">270</a> { 7.25 , 7.05218545073854 },
-<a class="jxr_linenumber" name="L271" href="#L271">271</a> { 7.375 , 7.292140407056348 },
-<a class="jxr_linenumber" name="L272" href="#L272">272</a> { 7.5 , 7.534364236758733 },
-<a class="jxr_linenumber" name="L273" href="#L273">273</a> { 7.625 , 7.778816557302289 },
-<a class="jxr_linenumber" name="L274" href="#L274">274</a> { 7.75 , 8.025458396315983 },
-<a class="jxr_linenumber" name="L275" href="#L275">275</a> { 7.875 , 8.274252119110479 },
-<a class="jxr_linenumber" name="L276" href="#L276">276</a> { 8.0 , 8.525161361065415 },
-<a class="jxr_linenumber" name="L277" href="#L277">277</a> { 8.125 , 8.77815096449171 },
-<a class="jxr_linenumber" name="L278" href="#L278">278</a> { 8.25 , 9.033186919605123 },
-<a class="jxr_linenumber" name="L279" href="#L279">279</a> { 8.375 , 9.290236309282232 },
-<a class="jxr_linenumber" name="L280" href="#L280">280</a> { 8.5 , 9.549267257300997 },
-<a class="jxr_linenumber" name="L281" href="#L281">281</a> { 8.625 , 9.810248879795765 },
-<a class="jxr_linenumber" name="L282" href="#L282">282</a> { 8.75 , 10.07315123968124 },
-<a class="jxr_linenumber" name="L283" href="#L283">283</a> { 8.875 , 10.33794530382217 },
-<a class="jxr_linenumber" name="L284" href="#L284">284</a> { 9.0 , 10.60460290274525 },
-<a class="jxr_linenumber" name="L285" href="#L285">285</a> { 9.125 , 10.87309669270751 },
-<a class="jxr_linenumber" name="L286" href="#L286">286</a> { 9.25 , 11.14340011995171 },
-<a class="jxr_linenumber" name="L287" href="#L287">287</a> { 9.375 , 11.41548738699336 },
-<a class="jxr_linenumber" name="L288" href="#L288">288</a> { 9.5 , 11.68933342079727 },
-<a class="jxr_linenumber" name="L289" href="#L289">289</a> { 9.625 , 11.96491384271319 },
-<a class="jxr_linenumber" name="L290" href="#L290">290</a> { 9.75 , 12.24220494005076 },
-<a class="jxr_linenumber" name="L291" href="#L291">291</a> { 9.875 , 12.52118363918365 },
-<a class="jxr_linenumber" name="L292" href="#L292">292</a> { 10.0 , 12.80182748008147 },
-<a class="jxr_linenumber" name="L293" href="#L293">293</a> { 0.8 , .1520596783998376 },
-<a class="jxr_linenumber" name="L294" href="#L294">294</a> { 100.0 , 359.1342053695754 },
-<a class="jxr_linenumber" name="L295" href="#L295">295</a> { 1000.0 , 5905.220423209181 },
-<a class="jxr_linenumber" name="L296" href="#L296">296</a> { 10000.0 , 82099.71749644238 },
-<a class="jxr_linenumber" name="L297" href="#L297">297</a> { 100000.0 , 1051287.708973657 },
-<a class="jxr_linenumber" name="L298" href="#L298">298</a> { 1000000.0 , 1.2815504569147612e+7 },
-<a class="jxr_linenumber" name="L299" href="#L299">299</a> { 10000000.0 , 1.511809493694739e+8 },
-<a class="jxr_linenumber" name="L300" href="#L300">300</a> { 1.e+8 , 1.7420680661038346e+9 },
-<a class="jxr_linenumber" name="L301" href="#L301">301</a> { 1.e+9 , 1.972326582750371e+10 },
-<a class="jxr_linenumber" name="L302" href="#L302">302</a> { 1.e+10 , 2.202585092888106e+11 },
-<a class="jxr_linenumber" name="L303" href="#L303">303</a> };
-<a class="jxr_linenumber" name="L304" href="#L304">304</a>
-<a class="jxr_linenumber" name="L305" href="#L305">305</a> @Test
-<a class="jxr_linenumber" name="L306" href="#L306">306</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGamma() {
-<a class="jxr_linenumber" name="L307" href="#L307">307</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> ulps = 3;
-<a class="jxr_linenumber" name="L308" href="#L308">308</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < LOG_GAMMA_REF.length; i++) {
-<a class="jxr_linenumber" name="L309" href="#L309">309</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] data = LOG_GAMMA_REF[i];
-<a class="jxr_linenumber" name="L310" href="#L310">310</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = data[0];
-<a class="jxr_linenumber" name="L311" href="#L311">311</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expected = data[1];
-<a class="jxr_linenumber" name="L312" href="#L312">312</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actual = Gamma.logGamma(x);
-<a class="jxr_linenumber" name="L313" href="#L313">313</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol;
-<a class="jxr_linenumber" name="L314" href="#L314">314</a> <strong class="jxr_keyword">if</strong> (expected == 0.0) {
-<a class="jxr_linenumber" name="L315" href="#L315">315</a> tol = 1E-15;
-<a class="jxr_linenumber" name="L316" href="#L316">316</a> } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L317" href="#L317">317</a> tol = ulps * FastMath.ulp(expected);
-<a class="jxr_linenumber" name="L318" href="#L318">318</a> }
-<a class="jxr_linenumber" name="L319" href="#L319">319</a> Assert.assertEquals(Double.toString(x), expected, actual, tol);
-<a class="jxr_linenumber" name="L320" href="#L320">320</a> }
-<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> @Test
-<a class="jxr_linenumber" name="L324" href="#L324">324</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGammaPrecondition1() {
-<a class="jxr_linenumber" name="L325" href="#L325">325</a> Assert.assertTrue(Double.isNaN(Gamma.logGamma(0.0)));
-<a class="jxr_linenumber" name="L326" href="#L326">326</a> }
-<a class="jxr_linenumber" name="L327" href="#L327">327</a>
-<a class="jxr_linenumber" name="L328" href="#L328">328</a> @Test
-<a class="jxr_linenumber" name="L329" href="#L329">329</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGammaPrecondition2() {
-<a class="jxr_linenumber" name="L330" href="#L330">330</a> Assert.assertTrue(Double.isNaN(Gamma.logGamma(-1.0)));
-<a class="jxr_linenumber" name="L331" href="#L331">331</a> }
-<a class="jxr_linenumber" name="L332" href="#L332">332</a>
-<a class="jxr_linenumber" name="L333" href="#L333">333</a> <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L334" href="#L334">334</a> <em class="jxr_javadoccomment"> * <p></em>
-<a class="jxr_linenumber" name="L335" href="#L335">335</a> <em class="jxr_javadoccomment"> * Reference values for the {@link Gamma#invGamma1pm1(double)} method.</em>
-<a class="jxr_linenumber" name="L336" href="#L336">336</a> <em class="jxr_javadoccomment"> * These values were generated with the following <a</em>
-<a class="jxr_linenumber" name="L337" href="#L337">337</a> <em class="jxr_javadoccomment"> * href="<a href="http://maxima.sourceforge.net/" target="alexandria_uri">http://maxima.sourceforge.net/</a>">Maxima</a> script</em>
-<a class="jxr_linenumber" name="L338" href="#L338">338</a> <em class="jxr_javadoccomment"> * </p></em>
-<a class="jxr_linenumber" name="L339" href="#L339">339</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L340" href="#L340">340</a> <em class="jxr_javadoccomment"> * <pre></em>
-<a class="jxr_linenumber" name="L341" href="#L341">341</a> <em class="jxr_javadoccomment"> * kill(all);</em>
-<a class="jxr_linenumber" name="L342" href="#L342">342</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L343" href="#L343">343</a> <em class="jxr_javadoccomment"> * fpprec : 64;</em>
-<a class="jxr_linenumber" name="L344" href="#L344">344</a> <em class="jxr_javadoccomment"> * gam1(x) := 1 / gamma(1 + x) - 1;</em>
-<a class="jxr_linenumber" name="L345" href="#L345">345</a> <em class="jxr_javadoccomment"> * x : makelist(bfloat(i / 8), i, -4, 12);</em>
-<a class="jxr_linenumber" name="L346" href="#L346">346</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L347" href="#L347">347</a> <em class="jxr_javadoccomment"> * for i : 1 while i <= length(x) do print("{",</em>
-<a class="jxr_linenumber" name="L348" href="#L348">348</a> <em class="jxr_javadoccomment"> * float(x[i]),</em>
-<a class="jxr_linenumber" name="L349" href="#L349">349</a> <em class="jxr_javadoccomment"> * ",",</em>
-<a class="jxr_linenumber" name="L350" href="#L350">350</a> <em class="jxr_javadoccomment"> * float(gam1(x[i])),</em>
-<a class="jxr_linenumber" name="L351" href="#L351">351</a> <em class="jxr_javadoccomment"> * "},");</em>
-<a class="jxr_linenumber" name="L352" href="#L352">352</a> <em class="jxr_javadoccomment"> * </pre></em>
-<a class="jxr_linenumber" name="L353" href="#L353">353</a> <em class="jxr_javadoccomment"> */</em>
-<a class="jxr_linenumber" name="L354" href="#L354">354</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>[][] INV_GAMMA1P_M1_REF = {
-<a class="jxr_linenumber" name="L355" href="#L355">355</a> { -0.5 , -.4358104164522437 },
-<a class="jxr_linenumber" name="L356" href="#L356">356</a> { -0.375 , -.3029021533379859 },
-<a class="jxr_linenumber" name="L357" href="#L357">357</a> { -0.25 , -0.183951060901737 },
-<a class="jxr_linenumber" name="L358" href="#L358">358</a> { -0.125 , -.08227611018520711 },
-<a class="jxr_linenumber" name="L359" href="#L359">359</a> { 0.0 , 0.0 },
-<a class="jxr_linenumber" name="L360" href="#L360">360</a> { 0.125 , .06186116458306091 },
-<a class="jxr_linenumber" name="L361" href="#L361">361</a> { 0.25 , .1032626513208373 },
-<a class="jxr_linenumber" name="L362" href="#L362">362</a> { 0.375 , .1249687649039041 },
-<a class="jxr_linenumber" name="L363" href="#L363">363</a> { 0.5 , .1283791670955126 },
-<a class="jxr_linenumber" name="L364" href="#L364">364</a> { 0.625 , .1153565546592225 },
-<a class="jxr_linenumber" name="L365" href="#L365">365</a> { 0.75 , 0.0880652521310173 },
-<a class="jxr_linenumber" name="L366" href="#L366">366</a> { 0.875 , .04882730264547758 },
-<a class="jxr_linenumber" name="L367" href="#L367">367</a> { 1.0 , 0.0 },
-<a class="jxr_linenumber" name="L368" href="#L368">368</a> { 1.125 , -.05612340925950141 },
-<a class="jxr_linenumber" name="L369" href="#L369">369</a> { 1.25 , -.1173898789433302 },
-<a class="jxr_linenumber" name="L370" href="#L370">370</a> { 1.375 , -.1818408982517061 },
-<a class="jxr_linenumber" name="L371" href="#L371">371</a> { 1.5 , -0.247747221936325 },
-<a class="jxr_linenumber" name="L372" href="#L372">372</a> };
-<a class="jxr_linenumber" name="L373" href="#L373">373</a>
-<a class="jxr_linenumber" name="L374" href="#L374">374</a> @Test
-<a class="jxr_linenumber" name="L375" href="#L375">375</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInvGamma1pm1() {
-<a class="jxr_linenumber" name="L376" href="#L376">376</a>
-<a class="jxr_linenumber" name="L377" href="#L377">377</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> ulps = 3;
-<a class="jxr_linenumber" name="L378" href="#L378">378</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < INV_GAMMA1P_M1_REF.length; i++) {
-<a class="jxr_linenumber" name="L379" href="#L379">379</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] ref = INV_GAMMA1P_M1_REF[i];
-<a class="jxr_linenumber" name="L380" href="#L380">380</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = ref[0];
-<a class="jxr_linenumber" name="L381" href="#L381">381</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expected = ref[1];
-<a class="jxr_linenumber" name="L382" href="#L382">382</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actual = Gamma.invGamma1pm1(x);
-<a class="jxr_linenumber" name="L383" href="#L383">383</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol = ulps * FastMath.ulp(expected);
-<a class="jxr_linenumber" name="L384" href="#L384">384</a> Assert.assertEquals(Double.toString(x), expected, actual, tol);
-<a class="jxr_linenumber" name="L385" href="#L385">385</a> }
-<a class="jxr_linenumber" name="L386" href="#L386">386</a> }
+<a class="jxr_linenumber" name="L128" href="#L128">128</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDigammaNonRealArgs() {
+<a class="jxr_linenumber" name="L129" href="#L129">129</a> Assert.assertTrue(Double.isNaN(Gamma.digamma(Double.NaN)));
+<a class="jxr_linenumber" name="L130" href="#L130">130</a> Assert.assertTrue(Double.isInfinite(Gamma.digamma(Double.POSITIVE_INFINITY)));
+<a class="jxr_linenumber" name="L131" href="#L131">131</a> Assert.assertTrue(Double.isInfinite(Gamma.digamma(Double.NEGATIVE_INFINITY)));
+<a class="jxr_linenumber" name="L132" href="#L132">132</a> }
+<a class="jxr_linenumber" name="L133" href="#L133">133</a>
+<a class="jxr_linenumber" name="L134" href="#L134">134</a> @Test
+<a class="jxr_linenumber" name="L135" href="#L135">135</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testTrigamma() {
+<a class="jxr_linenumber" name="L136" href="#L136">136</a> <strong class="jxr_keyword">double</strong> eps = 1e-8;
+<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_comment">// computed using webMathematica. For example, to compute trigamma($i) = Polygamma(1, $i), use</em>
+<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_comment">//</em>
+<a class="jxr_linenumber" name="L139" href="#L139">139</a> <em class="jxr_comment">// http://functions.wolfram.com/webMathematica/Evaluated.jsp?name=PolyGamma2&plottype=0&vars={%221%22,%22$i%22}&digits=20</em>
+<a class="jxr_linenumber" name="L140" href="#L140">140</a> <strong class="jxr_keyword">double</strong>[] data = {
+<a class="jxr_linenumber" name="L141" href="#L141">141</a> 1e-4, 1.0000000164469368793e8,
+<a class="jxr_linenumber" name="L142" href="#L142">142</a> 1e-3, 1.0000016425331958690e6,
+<a class="jxr_linenumber" name="L143" href="#L143">143</a> 1e-2, 10001.621213528313220,
+<a class="jxr_linenumber" name="L144" href="#L144">144</a> 1e-1, 101.43329915079275882,
+<a class="jxr_linenumber" name="L145" href="#L145">145</a> 1, 1.6449340668482264365,
+<a class="jxr_linenumber" name="L146" href="#L146">146</a> 2, 0.64493406684822643647,
+<a class="jxr_linenumber" name="L147" href="#L147">147</a> 3, 0.39493406684822643647,
+<a class="jxr_linenumber" name="L148" href="#L148">148</a> 4, 0.28382295573711532536,
+<a class="jxr_linenumber" name="L149" href="#L149">149</a> 5, 0.22132295573711532536,
+<a class="jxr_linenumber" name="L150" href="#L150">150</a> 10, 0.10516633568168574612,
+<a class="jxr_linenumber" name="L151" href="#L151">151</a> 20, 0.051270822935203119832,
+<a class="jxr_linenumber" name="L152" href="#L152">152</a> 50, 0.020201333226697125806,
+<a class="jxr_linenumber" name="L153" href="#L153">153</a> 100, 0.010050166663333571395
+<a class="jxr_linenumber" name="L154" href="#L154">154</a> };
+<a class="jxr_linenumber" name="L155" href="#L155">155</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = data.length - 2; i >= 0; i -= 2) {
+<a class="jxr_linenumber" name="L156" href="#L156">156</a> Assert.assertEquals(String.format(<span class="jxr_string">"trigamma %.0f"</span>, data[i]), data[i + 1], Gamma.trigamma(data[i]), eps);
+<a class="jxr_linenumber" name="L157" href="#L157">157</a> }
+<a class="jxr_linenumber" name="L158" href="#L158">158</a> }
+<a class="jxr_linenumber" name="L159" href="#L159">159</a>
+<a class="jxr_linenumber" name="L160" href="#L160">160</a> @Test
+<a class="jxr_linenumber" name="L161" href="#L161">161</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testTrigammaNonRealArgs() {
+<a class="jxr_linenumber" name="L162" href="#L162">162</a> Assert.assertTrue(Double.isNaN(Gamma.trigamma(Double.NaN)));
+<a class="jxr_linenumber" name="L163" href="#L163">163</a> Assert.assertTrue(Double.isInfinite(Gamma.trigamma(Double.POSITIVE_INFINITY)));
+<a class="jxr_linenumber" name="L164" href="#L164">164</a> Assert.assertTrue(Double.isInfinite(Gamma.trigamma(Double.NEGATIVE_INFINITY)));
+<a class="jxr_linenumber" name="L165" href="#L165">165</a> }
+<a class="jxr_linenumber" name="L166" href="#L166">166</a>
+<a class="jxr_linenumber" name="L167" href="#L167">167</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_javadoccomment"> * Reference data for the {@link Gamma#logGamma(double)} function. This data</em>
+<a class="jxr_linenumber" name="L169" href="#L169">169</a> <em class="jxr_javadoccomment"> * was generated with the following <a</em>
+<a class="jxr_linenumber" name="L170" href="#L170">170</a> <em class="jxr_javadoccomment"> * href="<a href="http://maxima.sourceforge.net/" target="alexandria_uri">http://maxima.sourceforge.net/</a>">Maxima</a> script.</em>
+<a class="jxr_linenumber" name="L171" href="#L171">171</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L172" href="#L172">172</a> <em class="jxr_javadoccomment"> * <pre></em>
+<a class="jxr_linenumber" name="L173" href="#L173">173</a> <em class="jxr_javadoccomment"> * kill(all);</em>
+<a class="jxr_linenumber" name="L174" href="#L174">174</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L175" href="#L175">175</a> <em class="jxr_javadoccomment"> * fpprec : 64;</em>
+<a class="jxr_linenumber" name="L176" href="#L176">176</a> <em class="jxr_javadoccomment"> * gamln(x) := log(gamma(x));</em>
+<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment"> * x : append(makelist(bfloat(i / 8), i, 1, 80),</em>
+<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_javadoccomment"> * [0.8b0, 1b2, 1b3, 1b4, 1b5, 1b6, 1b7, 1b8, 1b9, 1b10]);</em>
+<a class="jxr_linenumber" name="L179" href="#L179">179</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L180" href="#L180">180</a> <em class="jxr_javadoccomment"> * for i : 1 while i <= length(x) do</em>
+<a class="jxr_linenumber" name="L181" href="#L181">181</a> <em class="jxr_javadoccomment"> * print("{", float(x[i]), ",", float(gamln(x[i])), "},");</em>
+<a class="jxr_linenumber" name="L182" href="#L182">182</a> <em class="jxr_javadoccomment"> * </pre></em>
+<a class="jxr_linenumber" name="L183" href="#L183">183</a> <em class="jxr_javadoccomment"> */</em>
+<a class="jxr_linenumber" name="L184" href="#L184">184</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>[][] LOG_GAMMA_REF = {
+<a class="jxr_linenumber" name="L185" href="#L185">185</a> { 0.125 , 2.019418357553796 },
+<a class="jxr_linenumber" name="L186" href="#L186">186</a> { 0.25 , 1.288022524698077 },
+<a class="jxr_linenumber" name="L187" href="#L187">187</a> { 0.375 , .8630739822706475 },
+<a class="jxr_linenumber" name="L188" href="#L188">188</a> { 0.5 , .5723649429247001 },
+<a class="jxr_linenumber" name="L189" href="#L189">189</a> { 0.625 , .3608294954889402 },
+<a class="jxr_linenumber" name="L190" href="#L190">190</a> { 0.75 , .2032809514312954 },
+<a class="jxr_linenumber" name="L191" href="#L191">191</a> { 0.875 , .08585870722533433 },
+<a class="jxr_linenumber" name="L192" href="#L192">192</a> { 0.890625 , .07353860936979656 },
+<a class="jxr_linenumber" name="L193" href="#L193">193</a> { 0.90625 , .06169536624059108 },
+<a class="jxr_linenumber" name="L194" href="#L194">194</a> { 0.921875 , .05031670080005688 },
+<a class="jxr_linenumber" name="L195" href="#L195">195</a> { 0.9375 , 0.0393909017345823 },
+<a class="jxr_linenumber" name="L196" href="#L196">196</a> { 0.953125 , .02890678734595923 },
+<a class="jxr_linenumber" name="L197" href="#L197">197</a> { 0.96875 , .01885367233441289 },
+<a class="jxr_linenumber" name="L198" href="#L198">198</a> { 0.984375 , .009221337197578781 },
+<a class="jxr_linenumber" name="L199" href="#L199">199</a> { 1.0 , 0.0 },
+<a class="jxr_linenumber" name="L200" href="#L200">200</a> { 1.015625 , - 0.00881970970573307 },
+<a class="jxr_linenumber" name="L201" href="#L201">201</a> { 1.03125 , - .01724677500176807 },
+<a class="jxr_linenumber" name="L202" href="#L202">202</a> { 1.046875 , - .02528981394675729 },
+<a class="jxr_linenumber" name="L203" href="#L203">203</a> { 1.0625 , - .03295710029357782 },
+<a class="jxr_linenumber" name="L204" href="#L204">204</a> { 1.078125 , - .04025658272400143 },
+<a class="jxr_linenumber" name="L205" href="#L205">205</a> { 1.09375 , - .04719590272716985 },
+<a class="jxr_linenumber" name="L206" href="#L206">206</a> { 1.109375 , - .05378241123619192 },
+<a class="jxr_linenumber" name="L207" href="#L207">207</a> { 1.125 , - .06002318412603958 },
+<a class="jxr_linenumber" name="L208" href="#L208">208</a> { 1.25 , - .09827183642181316 },
+<a class="jxr_linenumber" name="L209" href="#L209">209</a> { 1.375 , - .1177552707410788 },
+<a class="jxr_linenumber" name="L210" href="#L210">210</a> { 1.5 , - .1207822376352452 },
+<a class="jxr_linenumber" name="L211" href="#L211">211</a> { 1.625 , - .1091741337567954 },
+<a class="jxr_linenumber" name="L212" href="#L212">212</a> { 1.75 , - .08440112102048555 },
+<a class="jxr_linenumber" name="L213" href="#L213">213</a> { 1.875 , - 0.0476726853991883 },
+<a class="jxr_linenumber" name="L214" href="#L214">214</a> { 1.890625 , - .04229320615532515 },
+<a class="jxr_linenumber" name="L215" href="#L215">215</a> { 1.90625 , - .03674470657266143 },
+<a class="jxr_linenumber" name="L216" href="#L216">216</a> { 1.921875 , - .03102893865389552 },
+<a class="jxr_linenumber" name="L217" href="#L217">217</a> { 1.9375 , - .02514761940298887 },
+<a class="jxr_linenumber" name="L218" href="#L218">218</a> { 1.953125 , - .01910243184040138 },
+<a class="jxr_linenumber" name="L219" href="#L219">219</a> { 1.96875 , - .01289502598016741 },
+<a class="jxr_linenumber" name="L220" href="#L220">220</a> { 1.984375 , - .006527019770560387 },
+<a class="jxr_linenumber" name="L221" href="#L221">221</a> { 2.0 , 0.0 },
+<a class="jxr_linenumber" name="L222" href="#L222">222</a> { 2.015625 , .006684476830232185 },
+<a class="jxr_linenumber" name="L223" href="#L223">223</a> { 2.03125 , .01352488366498562 },
+<a class="jxr_linenumber" name="L224" href="#L224">224</a> { 2.046875 , .02051972208453692 },
+<a class="jxr_linenumber" name="L225" href="#L225">225</a> { 2.0625 , .02766752152285702 },
+<a class="jxr_linenumber" name="L226" href="#L226">226</a> { 2.078125 , 0.0349668385135861 },
+<a class="jxr_linenumber" name="L227" href="#L227">227</a> { 2.09375 , .04241625596251728 },
+<a class="jxr_linenumber" name="L228" href="#L228">228</a> { 2.109375 , .05001438244545164 },
+<a class="jxr_linenumber" name="L229" href="#L229">229</a> { 2.125 , .05775985153034387 },
+<a class="jxr_linenumber" name="L230" href="#L230">230</a> { 2.25 , .1248717148923966 },
+<a class="jxr_linenumber" name="L231" href="#L231">231</a> { 2.375 , .2006984603774558 },
+<a class="jxr_linenumber" name="L232" href="#L232">232</a> { 2.5 , .2846828704729192 },
+<a class="jxr_linenumber" name="L233" href="#L233">233</a> { 2.625 , .3763336820249054 },
+<a class="jxr_linenumber" name="L234" href="#L234">234</a> { 2.75 , .4752146669149371 },
+<a class="jxr_linenumber" name="L235" href="#L235">235</a> { 2.875 , .5809359740231859 },
+<a class="jxr_linenumber" name="L236" href="#L236">236</a> { 2.890625 , .5946142560817441 },
+<a class="jxr_linenumber" name="L237" href="#L237">237</a> { 2.90625 , .6083932548009232 },
+<a class="jxr_linenumber" name="L238" href="#L238">238</a> { 2.921875 , .6222723333588501 },
+<a class="jxr_linenumber" name="L239" href="#L239">239</a> { 2.9375 , .6362508628423761 },
+<a class="jxr_linenumber" name="L240" href="#L240">240</a> { 2.953125 , .6503282221022278 },
+<a class="jxr_linenumber" name="L241" href="#L241">241</a> { 2.96875 , .6645037976116387 },
+<a class="jxr_linenumber" name="L242" href="#L242">242</a> { 2.984375 , 0.678776983328359 },
+<a class="jxr_linenumber" name="L243" href="#L243">243</a> { 3.0 , .6931471805599453 },
+<a class="jxr_linenumber" name="L244" href="#L244">244</a> { 3.015625 , .7076137978322324 },
+<a class="jxr_linenumber" name="L245" href="#L245">245</a> { 3.03125 , .7221762507608962 },
+<a class="jxr_linenumber" name="L246" href="#L246">246</a> { 3.046875 , .7368339619260166 },
+<a class="jxr_linenumber" name="L247" href="#L247">247</a> { 3.0625 , 0.751586360749556 },
+<a class="jxr_linenumber" name="L248" href="#L248">248</a> { 3.078125 , .7664328833756681 },
+<a class="jxr_linenumber" name="L249" href="#L249">249</a> { 3.09375 , .7813729725537568 },
+<a class="jxr_linenumber" name="L250" href="#L250">250</a> { 3.109375 , .7964060775242092 },
+<a class="jxr_linenumber" name="L251" href="#L251">251</a> { 3.125 , 0.811531653906724 },
+<a class="jxr_linenumber" name="L252" href="#L252">252</a> { 3.25 , .9358019311087253 },
+<a class="jxr_linenumber" name="L253" href="#L253">253</a> { 3.375 , 1.06569589786406 },
+<a class="jxr_linenumber" name="L254" href="#L254">254</a> { 3.5 , 1.200973602347074 },
+<a class="jxr_linenumber" name="L255" href="#L255">255</a> { 3.625 , 1.341414578068493 },
+<a class="jxr_linenumber" name="L256" href="#L256">256</a> { 3.75 , 1.486815578593417 },
+<a class="jxr_linenumber" name="L257" href="#L257">257</a> { 3.875 , 1.6369886482725 },
+<a class="jxr_linenumber" name="L258" href="#L258">258</a> { 4.0 , 1.791759469228055 },
+<a class="jxr_linenumber" name="L259" href="#L259">259</a> { 4.125 , 1.950965937095089 },
+<a class="jxr_linenumber" name="L260" href="#L260">260</a> { 4.25 , 2.114456927450371 },
+<a class="jxr_linenumber" name="L261" href="#L261">261</a> { 4.375 , 2.282091222188554 },
+<a class="jxr_linenumber" name="L262" href="#L262">262</a> { 4.5 , 2.453736570842442 },
+<a class="jxr_linenumber" name="L263" href="#L263">263</a> { 4.625 , 2.62926886637513 },
+<a class="jxr_linenumber" name="L264" href="#L264">264</a> { 4.75 , 2.808571418575736 },
+<a class="jxr_linenumber" name="L265" href="#L265">265</a> { 4.875 , 2.99153431107781 },
+<a class="jxr_linenumber" name="L266" href="#L266">266</a> { 5.0 , 3.178053830347946 },
+<a class="jxr_linenumber" name="L267" href="#L267">267</a> { 5.125 , 3.368031956881733 },
+<a class="jxr_linenumber" name="L268" href="#L268">268</a> { 5.25 , 3.561375910386697 },
+<a class="jxr_linenumber" name="L269" href="#L269">269</a> { 5.375 , 3.757997741998131 },
+<a class="jxr_linenumber" name="L270" href="#L270">270</a> { 5.5 , 3.957813967618717 },
+<a class="jxr_linenumber" name="L271" href="#L271">271</a> { 5.625 , 4.160745237339519 },
+<a class="jxr_linenumber" name="L272" href="#L272">272</a> { 5.75 , 4.366716036622286 },
+<a class="jxr_linenumber" name="L273" href="#L273">273</a> { 5.875 , 4.57565441552762 },
+<a class="jxr_linenumber" name="L274" href="#L274">274</a> { 6.0 , 4.787491742782046 },
+<a class="jxr_linenumber" name="L275" href="#L275">275</a> { 6.125 , 5.002162481906205 },
+<a class="jxr_linenumber" name="L276" href="#L276">276</a> { 6.25 , 5.219603986990229 },
+<a class="jxr_linenumber" name="L277" href="#L277">277</a> { 6.375 , 5.439756316011858 },
+<a class="jxr_linenumber" name="L278" href="#L278">278</a> { 6.5 , 5.662562059857142 },
+<a class="jxr_linenumber" name="L279" href="#L279">279</a> { 6.625 , 5.887966185430003 },
+<a class="jxr_linenumber" name="L280" href="#L280">280</a> { 6.75 , 6.115915891431546 },
+<a class="jxr_linenumber" name="L281" href="#L281">281</a> { 6.875 , 6.346360475557843 },
+<a class="jxr_linenumber" name="L282" href="#L282">282</a> { 7.0 , 6.579251212010101 },
+<a class="jxr_linenumber" name="L283" href="#L283">283</a> { 7.125 , 6.814541238336996 },
+<a class="jxr_linenumber" name="L284" href="#L284">284</a> { 7.25 , 7.05218545073854 },
+<a class="jxr_linenumber" name="L285" href="#L285">285</a> { 7.375 , 7.292140407056348 },
+<a class="jxr_linenumber" name="L286" href="#L286">286</a> { 7.5 , 7.534364236758733 },
+<a class="jxr_linenumber" name="L287" href="#L287">287</a> { 7.625 , 7.778816557302289 },
+<a class="jxr_linenumber" name="L288" href="#L288">288</a> { 7.75 , 8.025458396315983 },
+<a class="jxr_linenumber" name="L289" href="#L289">289</a> { 7.875 , 8.274252119110479 },
+<a class="jxr_linenumber" name="L290" href="#L290">290</a> { 8.0 , 8.525161361065415 },
+<a class="jxr_linenumber" name="L291" href="#L291">291</a> { 8.125 , 8.77815096449171 },
+<a class="jxr_linenumber" name="L292" href="#L292">292</a> { 8.25 , 9.033186919605123 },
+<a class="jxr_linenumber" name="L293" href="#L293">293</a> { 8.375 , 9.290236309282232 },
+<a class="jxr_linenumber" name="L294" href="#L294">294</a> { 8.5 , 9.549267257300997 },
+<a class="jxr_linenumber" name="L295" href="#L295">295</a> { 8.625 , 9.810248879795765 },
+<a class="jxr_linenumber" name="L296" href="#L296">296</a> { 8.75 , 10.07315123968124 },
+<a class="jxr_linenumber" name="L297" href="#L297">297</a> { 8.875 , 10.33794530382217 },
+<a class="jxr_linenumber" name="L298" href="#L298">298</a> { 9.0 , 10.60460290274525 },
+<a class="jxr_linenumber" name="L299" href="#L299">299</a> { 9.125 , 10.87309669270751 },
+<a class="jxr_linenumber" name="L300" href="#L300">300</a> { 9.25 , 11.14340011995171 },
+<a class="jxr_linenumber" name="L301" href="#L301">301</a> { 9.375 , 11.41548738699336 },
+<a class="jxr_linenumber" name="L302" href="#L302">302</a> { 9.5 , 11.68933342079727 },
+<a class="jxr_linenumber" name="L303" href="#L303">303</a> { 9.625 , 11.96491384271319 },
+<a class="jxr_linenumber" name="L304" href="#L304">304</a> { 9.75 , 12.24220494005076 },
+<a class="jxr_linenumber" name="L305" href="#L305">305</a> { 9.875 , 12.52118363918365 },
+<a class="jxr_linenumber" name="L306" href="#L306">306</a> { 10.0 , 12.80182748008147 },
+<a class="jxr_linenumber" name="L307" href="#L307">307</a> { 0.8 , .1520596783998376 },
+<a class="jxr_linenumber" name="L308" href="#L308">308</a> { 100.0 , 359.1342053695754 },
+<a class="jxr_linenumber" name="L309" href="#L309">309</a> { 1000.0 , 5905.220423209181 },
+<a class="jxr_linenumber" name="L310" href="#L310">310</a> { 10000.0 , 82099.71749644238 },
+<a class="jxr_linenumber" name="L311" href="#L311">311</a> { 100000.0 , 1051287.708973657 },
+<a class="jxr_linenumber" name="L312" href="#L312">312</a> { 1000000.0 , 1.2815504569147612e+7 },
+<a class="jxr_linenumber" name="L313" href="#L313">313</a> { 10000000.0 , 1.511809493694739e+8 },
+<a class="jxr_linenumber" name="L314" href="#L314">314</a> { 1.e+8 , 1.7420680661038346e+9 },
+<a class="jxr_linenumber" name="L315" href="#L315">315</a> { 1.e+9 , 1.972326582750371e+10 },
+<a class="jxr_linenumber" name="L316" href="#L316">316</a> { 1.e+10 , 2.202585092888106e+11 },
+<a class="jxr_linenumber" name="L317" href="#L317">317</a> };
+<a class="jxr_linenumber" name="L318" href="#L318">318</a>
+<a class="jxr_linenumber" name="L319" href="#L319">319</a> @Test
+<a class="jxr_linenumber" name="L320" href="#L320">320</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGamma() {
+<a class="jxr_linenumber" name="L321" href="#L321">321</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> ulps = 3;
+<a class="jxr_linenumber" name="L322" href="#L322">322</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < LOG_GAMMA_REF.length; i++) {
+<a class="jxr_linenumber" name="L323" href="#L323">323</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] data = LOG_GAMMA_REF[i];
+<a class="jxr_linenumber" name="L324" href="#L324">324</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = data[0];
+<a class="jxr_linenumber" name="L325" href="#L325">325</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expected = data[1];
+<a class="jxr_linenumber" name="L326" href="#L326">326</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actual = Gamma.logGamma(x);
+<a class="jxr_linenumber" name="L327" href="#L327">327</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol;
+<a class="jxr_linenumber" name="L328" href="#L328">328</a> <strong class="jxr_keyword">if</strong> (expected == 0.0) {
+<a class="jxr_linenumber" name="L329" href="#L329">329</a> tol = 1E-15;
+<a class="jxr_linenumber" name="L330" href="#L330">330</a> } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L331" href="#L331">331</a> tol = ulps * FastMath.ulp(expected);
+<a class="jxr_linenumber" name="L332" href="#L332">332</a> }
+<a class="jxr_linenumber" name="L333" href="#L333">333</a> Assert.assertEquals(Double.toString(x), expected, actual, tol);
+<a class="jxr_linenumber" name="L334" href="#L334">334</a> }
+<a class="jxr_linenumber" name="L335" href="#L335">335</a> }
+<a class="jxr_linenumber" name="L336" href="#L336">336</a>
+<a class="jxr_linenumber" name="L337" href="#L337">337</a> @Test
+<a class="jxr_linenumber" name="L338" href="#L338">338</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGammaPrecondition1() {
+<a class="jxr_linenumber" name="L339" href="#L339">339</a> Assert.assertTrue(Double.isNaN(Gamma.logGamma(0.0)));
+<a class="jxr_linenumber" name="L340" href="#L340">340</a> }
+<a class="jxr_linenumber" name="L341" href="#L341">341</a>
+<a class="jxr_linenumber" name="L342" href="#L342">342</a> @Test
+<a class="jxr_linenumber" name="L343" href="#L343">343</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGammaPrecondition2() {
+<a class="jxr_linenumber" name="L344" href="#L344">344</a> Assert.assertTrue(Double.isNaN(Gamma.logGamma(-1.0)));
+<a class="jxr_linenumber" name="L345" href="#L345">345</a> }
+<a class="jxr_linenumber" name="L346" href="#L346">346</a>
+<a class="jxr_linenumber" name="L347" href="#L347">347</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L348" href="#L348">348</a> <em class="jxr_javadoccomment"> * <p></em>
+<a class="jxr_linenumber" name="L349" href="#L349">349</a> <em class="jxr_javadoccomment"> * Reference values for the {@link Gamma#invGamma1pm1(double)} method.</em>
+<a class="jxr_linenumber" name="L350" href="#L350">350</a> <em class="jxr_javadoccomment"> * These values were generated with the following <a</em>
+<a class="jxr_linenumber" name="L351" href="#L351">351</a> <em class="jxr_javadoccomment"> * href="<a href="http://maxima.sourceforge.net/" target="alexandria_uri">http://maxima.sourceforge.net/</a>">Maxima</a> script</em>
+<a class="jxr_linenumber" name="L352" href="#L352">352</a> <em class="jxr_javadoccomment"> * </p></em>
+<a class="jxr_linenumber" name="L353" href="#L353">353</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L354" href="#L354">354</a> <em class="jxr_javadoccomment"> * <pre></em>
+<a class="jxr_linenumber" name="L355" href="#L355">355</a> <em class="jxr_javadoccomment"> * kill(all);</em>
+<a class="jxr_linenumber" name="L356" href="#L356">356</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L357" href="#L357">357</a> <em class="jxr_javadoccomment"> * fpprec : 64;</em>
+<a class="jxr_linenumber" name="L358" href="#L358">358</a> <em class="jxr_javadoccomment"> * gam1(x) := 1 / gamma(1 + x) - 1;</em>
+<a class="jxr_linenumber" name="L359" href="#L359">359</a> <em class="jxr_javadoccomment"> * x : makelist(bfloat(i / 8), i, -4, 12);</em>
+<a class="jxr_linenumber" name="L360" href="#L360">360</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L361" href="#L361">361</a> <em class="jxr_javadoccomment"> * for i : 1 while i <= length(x) do print("{",</em>
+<a class="jxr_linenumber" name="L362" href="#L362">362</a> <em class="jxr_javadoccomment"> * float(x[i]),</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> <em class="jxr_javadoccomment"> * float(gam1(x[i])),</em>
+<a class="jxr_linenumber" name="L365" href="#L365">365</a> <em class="jxr_javadoccomment"> * "},");</em>
+<a class="jxr_linenumber" name="L366" href="#L366">366</a> <em class="jxr_javadoccomment"> * </pre></em>
+<a class="jxr_linenumber" name="L367" href="#L367">367</a> <em class="jxr_javadoccomment"> */</em>
+<a class="jxr_linenumber" name="L368" href="#L368">368</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>[][] INV_GAMMA1P_M1_REF = {
+<a class="jxr_linenumber" name="L369" href="#L369">369</a> { -0.5 , -.4358104164522437 },
+<a class="jxr_linenumber" name="L370" href="#L370">370</a> { -0.375 , -.3029021533379859 },
+<a class="jxr_linenumber" name="L371" href="#L371">371</a> { -0.25 , -0.183951060901737 },
+<a class="jxr_linenumber" name="L372" href="#L372">372</a> { -0.125 , -.08227611018520711 },
+<a class="jxr_linenumber" name="L373" href="#L373">373</a> { 0.0 , 0.0 },
+<a class="jxr_linenumber" name="L374" href="#L374">374</a> { 0.125 , .06186116458306091 },
+<a class="jxr_linenumber" name="L375" href="#L375">375</a> { 0.25 , .1032626513208373 },
+<a class="jxr_linenumber" name="L376" href="#L376">376</a> { 0.375 , .1249687649039041 },
+<a class="jxr_linenumber" name="L377" href="#L377">377</a> { 0.5 , .1283791670955126 },
+<a class="jxr_linenumber" name="L378" href="#L378">378</a> { 0.625 , .1153565546592225 },
+<a class="jxr_linenumber" name="L379" href="#L379">379</a> { 0.75 , 0.0880652521310173 },
+<a class="jxr_linenumber" name="L380" href="#L380">380</a> { 0.875 , .04882730264547758 },
+<a class="jxr_linenumber" name="L381" href="#L381">381</a> { 1.0 , 0.0 },
+<a class="jxr_linenumber" name="L382" href="#L382">382</a> { 1.125 , -.05612340925950141 },
+<a class="jxr_linenumber" name="L383" href="#L383">383</a> { 1.25 , -.1173898789433302 },
+<a class="jxr_linenumber" name="L384" href="#L384">384</a> { 1.375 , -.1818408982517061 },
+<a class="jxr_linenumber" name="L385" href="#L385">385</a> { 1.5 , -0.247747221936325 },
+<a class="jxr_linenumber" name="L386" href="#L386">386</a> };
<a class="jxr_linenumber" name="L387" href="#L387">387</a>
-<a class="jxr_linenumber" name="L388" href="#L388">388</a> @Test(expected = NumberIsTooSmallException.<strong class="jxr_keyword">class</strong>)
-<a class="jxr_linenumber" name="L389" href="#L389">389</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInvGamma1pm1Precondition1() {
+<a class="jxr_linenumber" name="L388" href="#L388">388</a> @Test
+<a class="jxr_linenumber" name="L389" href="#L389">389</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInvGamma1pm1() {
<a class="jxr_linenumber" name="L390" href="#L390">390</a>
-<a class="jxr_linenumber" name="L391" href="#L391">391</a> Gamma.invGamma1pm1(-0.51);
-<a class="jxr_linenumber" name="L392" href="#L392">392</a> }
-<a class="jxr_linenumber" name="L393" href="#L393">393</a>
-<a class="jxr_linenumber" name="L394" href="#L394">394</a> @Test(expected = NumberIsTooLargeException.<strong class="jxr_keyword">class</strong>)
-<a class="jxr_linenumber" name="L395" href="#L395">395</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInvGamma1pm1Precondition2() {
-<a class="jxr_linenumber" name="L396" href="#L396">396</a>
-<a class="jxr_linenumber" name="L397" href="#L397">397</a> Gamma.invGamma1pm1(1.51);
-<a class="jxr_linenumber" name="L398" href="#L398">398</a> }
-<a class="jxr_linenumber" name="L399" href="#L399">399</a>
-<a class="jxr_linenumber" name="L400" href="#L400">400</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>[][] LOG_GAMMA1P_REF = {
-<a class="jxr_linenumber" name="L401" href="#L401">401</a> { - 0.5 , .5723649429247001 },
-<a class="jxr_linenumber" name="L402" href="#L402">402</a> { - 0.375 , .3608294954889402 },
-<a class="jxr_linenumber" name="L403" href="#L403">403</a> { - 0.25 , .2032809514312954 },
-<a class="jxr_linenumber" name="L404" href="#L404">404</a> { - 0.125 , .08585870722533433 },
-<a class="jxr_linenumber" name="L405" href="#L405">405</a> { 0.0 , 0.0 },
-<a class="jxr_linenumber" name="L406" href="#L406">406</a> { 0.125 , - .06002318412603958 },
-<a class="jxr_linenumber" name="L407" href="#L407">407</a> { 0.25 , - .09827183642181316 },
-<a class="jxr_linenumber" name="L408" href="#L408">408</a> { 0.375 , - .1177552707410788 },
-<a class="jxr_linenumber" name="L409" href="#L409">409</a> { 0.5 , - .1207822376352452 },
-<a class="jxr_linenumber" name="L410" href="#L410">410</a> { 0.625 , - .1091741337567954 },
-<a class="jxr_linenumber" name="L411" href="#L411">411</a> { 0.75 , - .08440112102048555 },
-<a class="jxr_linenumber" name="L412" href="#L412">412</a> { 0.875 , - 0.0476726853991883 },
-<a class="jxr_linenumber" name="L413" href="#L413">413</a> { 1.0 , 0.0 },
-<a class="jxr_linenumber" name="L414" href="#L414">414</a> { 1.125 , .05775985153034387 },
-<a class="jxr_linenumber" name="L415" href="#L415">415</a> { 1.25 , .1248717148923966 },
-<a class="jxr_linenumber" name="L416" href="#L416">416</a> { 1.375 , .2006984603774558 },
-<a class="jxr_linenumber" name="L417" href="#L417">417</a> { 1.5 , .2846828704729192 },
-<a class="jxr_linenumber" name="L418" href="#L418">418</a> };
-<a class="jxr_linenumber" name="L419" href="#L419">419</a>
-<a class="jxr_linenumber" name="L420" href="#L420">420</a> @Test
-<a class="jxr_linenumber" name="L421" href="#L421">421</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGamma1p() {
-<a class="jxr_linenumber" name="L422" href="#L422">422</a>
-<a class="jxr_linenumber" name="L423" href="#L423">423</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> ulps = 3;
-<a class="jxr_linenumber" name="L424" href="#L424">424</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < LOG_GAMMA1P_REF.length; i++) {
-<a class="jxr_linenumber" name="L425" href="#L425">425</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] ref = LOG_GAMMA1P_REF[i];
-<a class="jxr_linenumber" name="L426" href="#L426">426</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = ref[0];
-<a class="jxr_linenumber" name="L427" href="#L427">427</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expected = ref[1];
-<a class="jxr_linenumber" name="L428" href="#L428">428</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actual = Gamma.logGamma1p(x);
-<a class="jxr_linenumber" name="L429" href="#L429">429</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol = ulps * FastMath.ulp(expected);
-<a class="jxr_linenumber" name="L430" href="#L430">430</a> Assert.assertEquals(Double.toString(x), expected, actual, tol);
-<a class="jxr_linenumber" name="L431" href="#L431">431</a> }
-<a class="jxr_linenumber" name="L432" href="#L432">432</a> }
+<a class="jxr_linenumber" name="L391" href="#L391">391</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> ulps = 3;
+<a class="jxr_linenumber" name="L392" href="#L392">392</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < INV_GAMMA1P_M1_REF.length; i++) {
+<a class="jxr_linenumber" name="L393" href="#L393">393</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] ref = INV_GAMMA1P_M1_REF[i];
+<a class="jxr_linenumber" name="L394" href="#L394">394</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x = ref[0];
+<a class="jxr_linenumber" name="L395" href="#L395">395</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> expected = ref[1];
+<a class="jxr_linenumber" name="L396" href="#L396">396</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> actual = Gamma.invGamma1pm1(x);
+<a class="jxr_linenumber" name="L397" href="#L397">397</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> tol = ulps * FastMath.ulp(expected);
+<a class="jxr_linenumber" name="L398" href="#L398">398</a> Assert.assertEquals(Double.toString(x), expected, actual, tol);
+<a class="jxr_linenumber" name="L399" href="#L399">399</a> }
+<a class="jxr_linenumber" name="L400" href="#L400">400</a> }
+<a class="jxr_linenumber" name="L401" href="#L401">401</a>
+<a class="jxr_linenumber" name="L402" href="#L402">402</a> @Test(expected = NumberIsTooSmallException.<strong class="jxr_keyword">class</strong>)
+<a class="jxr_linenumber" name="L403" href="#L403">403</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInvGamma1pm1Precondition1() {
+<a class="jxr_linenumber" name="L404" href="#L404">404</a>
+<a class="jxr_linenumber" name="L405" href="#L405">405</a> Gamma.invGamma1pm1(-0.51);
+<a class="jxr_linenumber" name="L406" href="#L406">406</a> }
+<a class="jxr_linenumber" name="L407" href="#L407">407</a>
+<a class="jxr_linenumber" name="L408" href="#L408">408</a> @Test(expected = NumberIsTooLargeException.<strong class="jxr_keyword">class</strong>)
+<a class="jxr_linenumber" name="L409" href="#L409">409</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInvGamma1pm1Precondition2() {
+<a class="jxr_linenumber" name="L410" href="#L410">410</a>
+<a class="jxr_linenumber" name="L411" href="#L411">411</a> Gamma.invGamma1pm1(1.51);
+<a class="jxr_linenumber" name="L412" href="#L412">412</a> }
+<a class="jxr_linenumber" name="L413" href="#L413">413</a>
+<a class="jxr_linenumber" name="L414" href="#L414">414</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>[][] LOG_GAMMA1P_REF = {
+<a class="jxr_linenumber" name="L415" href="#L415">415</a> { - 0.5 , .5723649429247001 },
+<a class="jxr_linenumber" name="L416" href="#L416">416</a> { - 0.375 , .3608294954889402 },
+<a class="jxr_linenumber" name="L417" href="#L417">417</a> { - 0.25 , .2032809514312954 },
+<a class="jxr_linenumber" name="L418" href="#L418">418</a> { - 0.125 , .08585870722533433 },
+<a class="jxr_linenumber" name="L419" href="#L419">419</a> { 0.0 , 0.0 },
+<a class="jxr_linenumber" name="L420" href="#L420">420</a> { 0.125 , - .06002318412603958 },
+<a class="jxr_linenumber" name="L421" href="#L421">421</a> { 0.25 , - .09827183642181316 },
+<a class="jxr_linenumber" name="L422" href="#L422">422</a> { 0.375 , - .1177552707410788 },
+<a class="jxr_linenumber" name="L423" href="#L423">423</a> { 0.5 , - .1207822376352452 },
+<a class="jxr_linenumber" name="L424" href="#L424">424</a> { 0.625 , - .1091741337567954 },
+<a class="jxr_linenumber" name="L425" href="#L425">425</a> { 0.75 , - .08440112102048555 },
+<a class="jxr_linenumber" name="L426" href="#L426">426</a> { 0.875 , - 0.0476726853991883 },
+<a class="jxr_linenumber" name="L427" href="#L427">427</a> { 1.0 , 0.0 },
+<a class="jxr_linenumber" name="L428" href="#L428">428</a> { 1.125 , .05775985153034387 },
+<a class="jxr_linenumber" name="L429" href="#L429">429</a> { 1.25 , .1248717148923966 },
+<a class="jxr_linenumber" name="L430" href="#L430">430</a> { 1.375 , .2006984603774558 },
+<a class="jxr_linenumber" name="L431" href="#L431">431</a> { 1.5 , .2846828704729192 },
+<a class="jxr_linenumber" name="L432" href="#L432">432</a> };
<a class="jxr_linenumber" name="L433" href="#L433">433</a>
-<a class="jxr_linenumber" name="L434" href="#L434">434</a> @Test(expected = NumberIsTooSmallException.<strong class="jxr_keyword">class</strong>)
-<a class="jxr_linenumber" name="L435" href="#L435">435</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGamma1pPrecondition1() {
+<a class="jxr_linenumber" name="L434" href="#L434">434</a> @Test
+<a class="jxr_linenumber" name="L435" href="#L435">435</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGamma1p() {
<a class="jxr_linenumber" name="L436" href="#L436">436</a>
-<a class="jxr_linenumber" name="L437" href="#L437">437</a> Gamma.logGamma1p(-0.51);
-<a class="jxr_linenumber" name="L438" href="#L438">438</a> }
-<a class="jxr_linenumber" name="L439" href="#L439">439</a>
-<a class="jxr_linenumber" name="L440" href="#L440">440</a> @Test(expected = NumberIsTooLargeException.<strong class="jxr_keyword">class</strong>)
-<a class="jxr_linenumber" name="L441" href="#L441">441</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testLogGamma1pPrecondition2() {
-<a class="jxr_linenumber" name="L442" href="#L442">442</a>
-<a class="jxr_linenumber" name="L443" href="#L443">443</a> Gamma.logGamma1p(1.51);
-<a class="jxr_linenumber" name="L444" href="#L444">444</a> }
-<a class="jxr_linenumber" name="L445" href="#L445">445</a>
-<a class="jxr_linenumber" name="L446" href="#L446">446</a> <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L447" href="#L447">447</a> <em class="jxr_javadoccomment"> * Reference data for the {@link Gamma#gamma(double)} function. This</em>
-<a class="jxr_linenumber" name="L448" href="#L448">448</a> <em class="jxr_javadoccomment"> * data was generated with the following <a</em>
-<a class="jxr_linenumber" name="L449" href="#L449">449</a> <em class="jxr_javadoccomment"> * href="<a href="http://maxima.sourceforge.net/" target="alexandria_uri">http://maxima.sourceforge.net/</a>">Maxima</a> script.</em>
-<a class="jxr_linenumber" name="L450" href="#L450">450</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L451" href="#L451">451</a> <em class="jxr_javadoccomment"> * <pre></em>
-<a class="jxr_linenumber" name="L452" href="#L452">452</a> <em class="jxr_javadoccomment"> * kill(all);</em>
-<a class="jxr_linenumber" name="L453" href="#L453">453</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L454" href="#L454">454</a> <em class="jxr_javadoccomment"> * fpprec : 64;</em>
-<a class="jxr_linenumber" name="L455" href="#L455">455</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L456" href="#L456">456</a> <em class="jxr_javadoccomment"> * EPSILON : 10**(-fpprec + 1);</em>
-<a class="jxr_linenumber" name="L457" href="#L457">457</a> <em class="jxr_javadoccomment"> * isInteger(x) := abs(x - floor(x)) <= EPSILON * abs(x);</em>
-<a class="jxr_linenumber" name="L458" href="#L458">458</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L459" href="#L459">459</a> <em class="jxr_javadoccomment"> * x : makelist(bfloat(i / 8), i, -160, 160);</em>
-<a class="jxr_linenumber" name="L460" href="#L460">460</a> <em class="jxr_javadoccomment"> * x : append(x, makelist(bfloat(i / 2), i, 41, 200));</em>
-<a class="jxr_linenumber" name="L461" href="#L461">461</a> <em class="jxr_javadoccomment"> *</em>
-<a class="jxr_linenumber" name="L462" href="#L462">462</a> <em class="jxr_javadoccomment"> * for i : 1 while i <= length(x) do if not(isInteger(x[i])) then</em>
-<a class="jxr_linenumber" name="L463" href="#L463">463</a> <em class="jxr_javadoccomment"> * print("{", float(x[i]), ",", float(gamma(x[i])), "},");</em>
-<a class="jxr_linenumber" name="L464" href="#L464">464</a> <em class="jxr_javadoccomment"> * </pre></em>
-<a class="jxr_linenumber" name="L465" href="#L465">465</a> <em class="jxr_javadoccomment"> */</em>
-<a class="jxr_linenumber" name="L466" href="#L466">466</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>[][] GAMMA_REF = {
-<a class="jxr_linenumber" name="L467" href="#L467">467</a> { - 19.875 , 4.920331854832504e-18 },
-<a class="jxr_linenumber" name="L468" href="#L468">468</a> { - 19.75 , 3.879938752480031e-18 },
-<a class="jxr_linenumber" name="L469" href="#L469">469</a> { - 19.625 , 4.323498423815027e-18 },
-<a class="jxr_linenumber" name="L470" href="#L470">470</a> { - 19.5 , 5.811045977502237e-18 },
-<a class="jxr_linenumber" name="L471" href="#L471">471</a> { - 19.375 , 9.14330910942125e-18 },
-<a class="jxr_linenumber" name="L472" href="#L472">472</a> { - 19.25 , 1.735229114436739e-17 },
-<a class="jxr_linenumber" name="L473" href="#L473">473</a> { - 19.125 , 4.653521565668223e-17 },
-<a class="jxr_linenumber" name="L474" href="#L474">474</a> { - 18.875 , - 9.779159561479603e-17 },
-<a class="jxr_linenumber" name="L475" href="#L475">475</a> { - 18.75 , - 7.662879036148061e-17 },
-<a class="jxr_linenumber" name="L476" href="#L476">476</a> { - 18.625 , - 8.484865656736991e-17 },
-<a class="jxr_linenumber" name="L477" href="#L477">477</a> { - 18.5 , - 1.133153965612936e-16 },
-<a class="jxr_linenumber" name="L478" href="#L478">478</a> { - 18.375 , - 1.771516139950367e-16 },
-<a class="jxr_linenumber" name="L479" href="#L479">479</a> { - 18.25 , - 3.340316045290721e-16 },
-<a class="jxr_linenumber" name="L480" href="#L480">480</a> { - 18.125 , - 8.899859994340475e-16 },
-<a class="jxr_linenumber" name="L481" href="#L481">481</a> { - 17.875 , 1.845816367229275e-15 },
-<a class="jxr_linenumber" name="L482" href="#L482">482</a> { - 17.75 , 1.436789819277761e-15 },
-<a class="jxr_linenumber" name="L483" href="#L483">483</a> { - 17.625 , 1.580306228567265e-15 },
-<a class="jxr_linenumber" name="L484" href="#L484">484</a> { - 17.5 , 2.096334836383932e-15 },
-<a class="jxr_linenumber" name="L485" href="#L485">485</a> { - 17.375 , 3.255160907158799e-15 },
-<a class="jxr_linenumber" name="L486" href="#L486">486</a> { - 17.25 , 6.096076782655566e-15 },
-<a class="jxr_linenumber" name="L487" href="#L487">487</a> { - 17.125 , 1.613099623974211e-14 },
-<a class="jxr_linenumber" name="L488" href="#L488">488</a> { - 16.875 , - 3.29939675642233e-14 },
-<a class="jxr_linenumber" name="L489" href="#L489">489</a> { - 16.75 , - 2.550301929218027e-14 },
-<a class="jxr_linenumber" name="L490" href="#L490">490</a> { - 16.625 , - 2.785289727849803e-14 },
-<a class="jxr_linenumber" name="L491" href="#L491">491</a> { - 16.5 , - 3.66858596367188e-14 },
-<a class="jxr_linenumber" name="L492" href="#L492">492</a> { - 16.375 , - 5.655842076188414e-14 },
-<a class="jxr_linenumber" name="L493" href="#L493">493</a> { - 16.25 , - 1.051573245008085e-13 },
-<a class="jxr_linenumber" name="L494" href="#L494">494</a> { - 16.125 , - 2.762433106055837e-13 },
-<a class="jxr_linenumber" name="L495" href="#L495">495</a> { - 15.875 , 5.567732026462681e-13 },
-<a class="jxr_linenumber" name="L496" href="#L496">496</a> { - 15.75 , 4.271755731440195e-13 },
-<a class="jxr_linenumber" name="L497" href="#L497">497</a> { - 15.625 , 4.630544172550298e-13 },
-<a class="jxr_linenumber" name="L498" href="#L498">498</a> { - 15.5 , 6.053166840058604e-13 },
-<a class="jxr_linenumber" name="L499" href="#L499">499</a> { - 15.375 , 9.261441399758529e-13 },
-<a class="jxr_linenumber" name="L500" href="#L500">500</a> { - 15.25 , 1.708806523138138e-12 },
-<a class="jxr_linenumber" name="L501" href="#L501">501</a> { - 15.125 , 4.454423383515037e-12 },
-<a class="jxr_linenumber" name="L502" href="#L502">502</a> { - 14.875 , - 8.838774592009505e-12 },
-<a class="jxr_linenumber" name="L503" href="#L503">503</a> { - 14.75 , - 6.728015277018307e-12 },
-<a class="jxr_linenumber" name="L504" href="#L504">504</a> { - 14.625 , - 7.235225269609841e-12 },
-<a class="jxr_linenumber" name="L505" href="#L505">505</a> { - 14.5 , - 9.382408602090835e-12 },
-<a class="jxr_linenumber" name="L506" href="#L506">506</a> { - 14.375 , - 1.423946615212874e-11 },
-<a class="jxr_linenumber" name="L507" href="#L507">507</a> { - 14.25 , - 2.605929947785661e-11 },
-<a class="jxr_linenumber" name="L508" href="#L508">508</a> { - 14.125 , - 6.737315367566492e-11 },
-<a class="jxr_linenumber" name="L509" href="#L509">509</a> { - 13.875 , 1.314767720561414e-10 },
-<a class="jxr_linenumber" name="L510" href="#L510">510</a> { - 13.75 , 9.923822533602004e-11 },
-<a class="jxr_linenumber" name="L511" href="#L511">511</a> { - 13.625 , 1.058151695680439e-10 },
-<a class="jxr_linenumber" name="L512" href="#L512">512</a> { - 13.5 , 1.360449247303171e-10 },
-<a class="jxr_linenumber" name="L513" href="#L513">513</a> { - 13.375 , 2.046923259368506e-10 },
-<a class="jxr_linenumber" name="L514" href="#L514">514</a> { - 13.25 , 3.713450175594567e-10 },
-<a class="jxr_linenumber" name="L515" href="#L515">515</a> { - 13.125 , 9.516457956687671e-10 },
-<a class="jxr_linenumber" name="L516" href="#L516">516</a> { - 12.875 , - 1.8242402122789617e-9 },
-<a class="jxr_linenumber" name="L517" href="#L517">517</a> { - 12.75 , - 1.3645255983702756e-9 },
-<a class="jxr_linenumber" name="L518" href="#L518">518</a> { - 12.625 , - 1.4417316853645984e-9 },
-<a class="jxr_linenumber" name="L519" href="#L519">519</a> { - 12.5 , - 1.836606483859281e-9 },
-<a class="jxr_linenumber" name="L520" href="#L520">520</a> { - 12.375 , - 2.7377598594053765e-9 },
-<a class="jxr_linenumber" name="L521" href="#L521">521</a> { - 12.25 , - 4.9203214826628017e-9 },
-<a class="jxr_linenumber" name="L522" href="#L522">522</a> { - 12.125 , - 1.2490351068152569e-8 },
-<a class="jxr_linenumber" name="L523" href="#L523">523</a> { - 11.875 , 2.3487092733091633e-8 },
-<a class="jxr_linenumber" name="L524" href="#L524">524</a> { - 11.75 , 1.7397701379221012e-8 },
-<a class="jxr_linenumber" name="L525" href="#L525">525</a> { - 11.625 , 1.8201862527728055e-8 },
-<a class="jxr_linenumber" name="L526" href="#L526">526</a> { - 11.5 , 2.295758104824101e-8 },
-<a class="jxr_linenumber" name="L527" href="#L527">527</a> { - 11.375 , 3.3879778260141535e-8 },
-<a class="jxr_linenumber" name="L528" href="#L528">528</a> { - 11.25 , 6.027393816261931e-8 },
-<a class="jxr_linenumber" name="L529" href="#L529">529</a> { - 11.125 , 1.5144550670134987e-7 },
-<a class="jxr_linenumber" name="L530" href="#L530">530</a> { - 10.875 , - 2.7890922620546316e-7 },
-<a class="jxr_linenumber" name="L531" href="#L531">531</a> { - 10.75 , - 2.044229912058469e-7 },
-<a class="jxr_linenumber" name="L532" href="#L532">532</a> { - 10.625 , - 2.1159665188483867e-7 },
-<a class="jxr_linenumber" name="L533" href="#L533">533</a> { - 10.5 , - 2.640121820547716e-7 },
-<a class="jxr_linenumber" name="L534" href="#L534">534</a> { - 10.375 , - 3.8538247770911e-7 },
-<a class="jxr_linenumber" name="L535" href="#L535">535</a> { - 10.25 , - 6.780818043294673e-7 },
-<a class="jxr_linenumber" name="L536" href="#L536">536</a> { - 10.125 , - 1.6848312620525174e-6 },
-<a class="jxr_linenumber" name="L537" href="#L537">537</a> { - 9.875 , 3.0331378349844124e-6 },
-<a class="jxr_linenumber" name="L538" href="#L538">538</a> { - 9.75 , 2.1975471554628537e-6 },
-<a class="jxr_linenumber" name="L539" href="#L539">539</a> { - 9.625 , 2.2482144262764103e-6 },
-<a class="jxr_linenumber" name="L540" href="#L540">540</a> { - 9.5 , 2.772127911575102e-6 },
-<a class="jxr_linenumber" name="L541" href="#L541">541</a> { - 9.375 , 3.998343206232017e-6 },
-<a class="jxr_linenumber" name="L542" href="#L542">542</a> { - 9.25 , 6.95033849437704e-6 },
-<a class="jxr_linenumber" name="L543" href="#L543">543</a> { - 9.125 , 1.7058916528281737e-5 },
-<a class="jxr_linenumber" name="L544" href="#L544">544</a> { - 8.875 , - 2.9952236120471065e-5 },
-<a class="jxr_linenumber" name="L545" href="#L545">545</a> { - 8.75 , - 2.1426084765762826e-5 },
-<a class="jxr_linenumber" name="L546" href="#L546">546</a> { - 8.625 , - 2.163906385291045e-5 },
-<a class="jxr_linenumber" name="L547" href="#L547">547</a> { - 8.5 , - 2.633521515996347e-5 },
-<a class="jxr_linenumber" name="L548" href="#L548">548</a> { - 8.375 , - 3.748446755842515e-5 },
-<a class="jxr_linenumber" name="L549" href="#L549">549</a> { - 8.25 , - 6.429063107298763e-5 },
-<a class="jxr_linenumber" name="L550" href="#L550">550</a> { - 8.125 , - 1.5566261332057085e-4 },
-<a class="jxr_linenumber" name="L551" href="#L551">551</a> { - 7.875 , 2.658260955691807e-4 },
-<a class="jxr_linenumber" name="L552" href="#L552">552</a> { - 7.75 , 1.874782417004247e-4 },
-<a class="jxr_linenumber" name="L553" href="#L553">553</a> { - 7.625 , 1.8663692573135265e-4 },
-<a class="jxr_linenumber" name="L554" href="#L554">554</a> { - 7.5 , 2.238493288596895e-4 },
-<a class="jxr_linenumber" name="L555" href="#L555">555</a> { - 7.375 , 3.1393241580181064e-4 },
-<a class="jxr_linenumber" name="L556" href="#L556">556</a> { - 7.25 , 5.303977063521479e-4 },
-<a class="jxr_linenumber" name="L557" href="#L557">557</a> { - 7.125 , .001264758733229638 },
-<a class="jxr_linenumber" name="L558" href="#L558">558</a> { - 6.875 , - .002093380502607298 },
-<a class="jxr_linenumber" name="L559" href="#L559">559</a> { - 6.75 , - .001452956373178292 },
-<a class="jxr_linenumber" name="L560" href="#L560">560</a> { - 6.625 , - .001423106558701564 },
-<a class="jxr_linenumber" name="L561" href="#L561">561</a> { - 6.5 , - .001678869966447671 },
-<a class="jxr_linenumber" name="L562" href="#L562">562</a> { - 6.375 , - .002315251566538353 },
-<a class="jxr_linenumber" name="L563" href="#L563">563</a> { - 6.25 , - .003845383371053072 },
-<a class="jxr_linenumber" name="L564" href="#L564">564</a> { - 6.125 , - .009011405974261174 },
-<a class="jxr_linenumber" name="L565" href="#L565">565</a> { - 5.875 , .01439199095542518 },
-<a class="jxr_linenumber" name="L566" href="#L566">566</a> { - 5.75 , .009807455518953468 },
-<a class="jxr_linenumber" name="L567" href="#L567">567</a> { - 5.625 , .009428080951397862 },
-<a class="jxr_linenumber" name="L568" href="#L568">568</a> { - 5.5 , .01091265478190986 },
-<a class="jxr_linenumber" name="L569" href="#L569">569</a> { - 5.375 , 0.014759728736682 },
-<a class="jxr_linenumber" name="L570" href="#L570">570</a> { - 5.25 , 0.0240336460690817 },
-<a class="jxr_linenumber" name="L571" href="#L571">571</a> { - 5.125 , .05519486159234969 },
-<a class="jxr_linenumber" name="L572" href="#L572">572</a> { - 4.875 , - 0.0845529468631229 },
-<a class="jxr_linenumber" name="L573" href="#L573">573</a> { - 4.75 , - .05639286923398244 },
-<a class="jxr_linenumber" name="L574" href="#L574">574</a> { - 4.625 , - .05303295535161297 },
-<a class="jxr_linenumber" name="L575" href="#L575">575</a> { - 4.5 , - .06001960130050425 },
-<a class="jxr_linenumber" name="L576" href="#L576">576</a> { - 4.375 , - .07933354195966577 },
-<a class="jxr_linenumber" name="L577" href="#L577">577</a> { - 4.25 , - .1261766418626789 },
-<a class="jxr_linenumber" name="L578" href="#L578">578</a> { - 4.125 , - .2828736656607921 },
-<a class="jxr_linenumber" name="L579" href="#L579">579</a> { - 3.875 , .4121956159577241 },
-<a class="jxr_linenumber" name="L580" href="#L580">580</a> { - 3.75 , .2678661288614166 },
-<a class="jxr_linenumber" name="L581" href="#L581">581</a> { - 3.625 , 0.24527741850121 },
-<a class="jxr_linenumber" name="L582" href="#L582">582</a> { - 3.5 , .2700882058522691 },
-<a class="jxr_linenumber" name="L583" href="#L583">583</a> { - 3.375 , .3470842460735378 },
-<a class="jxr_linenumber" name="L584" href="#L584">584</a> { - 3.25 , .5362507279163854 },
-<a class="jxr_linenumber" name="L585" href="#L585">585</a> { - 3.125 , 1.166853870850768 },
-<a class="jxr_linenumber" name="L586" href="#L586">586</a> { - 2.875 , - 1.597258011836181 },
-<a class="jxr_linenumber" name="L587" href="#L587">587</a> { - 2.75 , - 1.004497983230312 },
-<a class="jxr_linenumber" name="L588" href="#L588">588</a> { - 2.625 , - .8891306420668862 },
-<a class="jxr_linenumber" name="L589" href="#L589">589</a> { - 2.5 , - .9453087204829419 },
-<a class="jxr_linenumber" name="L590" href="#L590">590</a> { - 2.375 , - 1.17140933049819 },
-<a class="jxr_linenumber" name="L591" href="#L591">591</a> { - 2.25 , - 1.742814865728253 },
-<a class="jxr_linenumber" name="L592" href="#L592">592</a> { - 2.125 , - 3.646418346408649 },
-<a class="jxr_linenumber" name="L593" href="#L593">593</a> { - 1.875 , 4.59211678402902 },
-<a class="jxr_linenumber" name="L594" href="#L594">594</a> { - 1.75 , 2.762369453883359 },
-<a class="jxr_linenumber" name="L595" href="#L595">595</a> { - 1.625 , 2.333967935425576 },
-<a class="jxr_linenumber" name="L596" href="#L596">596</a> { - 1.5 , 2.363271801207355 },
-<a class="jxr_linenumber" name="L597" href="#L597">597</a> { - 1.375 , 2.782097159933201 },
-<a class="jxr_linenumber" name="L598" href="#L598">598</a> { - 1.25 , 3.921333447888569 },
-<a class="jxr_linenumber" name="L599" href="#L599">599</a> { - 1.125 , 7.748638986118379 },
-<a class="jxr_linenumber" name="L600" href="#L600">600</a> { - 0.875 , - 8.610218970054413 },
-<a class="jxr_linenumber" name="L601" href="#L601">601</a> { - 0.75 , - 4.834146544295877 },
-<a class="jxr_linenumber" name="L602" href="#L602">602</a> { - 0.625 , - 3.792697895066561 },
-<a class="jxr_linenumber" name="L603" href="#L603">603</a> { - 0.5 , - 3.544907701811032 },
-<a class="jxr_linenumber" name="L604" href="#L604">604</a> { - 0.375 , - 3.825383594908152 },
-<a class="jxr_linenumber" name="L605" href="#L605">605</a> { - 0.25 , - 4.901666809860711 },
-<a class="jxr_linenumber" name="L606" href="#L606">606</a> { - 0.125 , - 8.717218859383175 },
-<a class="jxr_linenumber" name="L607" href="#L607">607</a> { 0.125 , 7.533941598797612 },
-<a class="jxr_linenumber" name="L608" href="#L608">608</a> { 0.25 , 3.625609908221908 },
-<a class="jxr_linenumber" name="L609" href="#L609">609</a> { 0.375 , 2.370436184416601 },
-<a class="jxr_linenumber" name="L610" href="#L610">610</a> { 0.5 , 1.772453850905516 },
-<a class="jxr_linenumber" name="L611" href="#L611">611</a> { 0.625 , 1.434518848090557 },
-<a class="jxr_linenumber" name="L612" href="#L612">612</a> { 0.75 , 1.225416702465178 },
-<a class="jxr_linenumber" name="L613" href="#L613">613</a> { 0.875 , 1.089652357422897 },
-<a class="jxr_linenumber" name="L614" href="#L614">614</a> { 1.0 , 1.0 },
-<a class="jxr_linenumber" name="L615" href="#L615">615</a> { 1.125 , .9417426998497015 },
-<a class="jxr_linenumber" name="L616" href="#L616">616</a> { 1.25 , 0.906402477055477 },
-<a class="jxr_linenumber" name="L617" href="#L617">617</a> { 1.375 , .8889135691562253 },
-<a class="jxr_linenumber" name="L618" href="#L618">618</a> { 1.5 , 0.886226925452758 },
-<a class="jxr_linenumber" name="L619" href="#L619">619</a> { 1.625 , 0.896574280056598 },
-<a class="jxr_linenumber" name="L620" href="#L620">620</a> { 1.75 , .9190625268488832 },
-<a class="jxr_linenumber" name="L621" href="#L621">621</a> { 1.875 , .9534458127450348 },
-<a class="jxr_linenumber" name="L622" href="#L622">622</a> { 2.0 , 1.0 },
-<a class="jxr_linenumber" name="L623" href="#L623">623</a> { 2.125 , 1.059460537330914 },
-<a class="jxr_linenumber" name="L624" href="#L624">624</a> { 2.25 , 1.133003096319346 },
-<a class="jxr_linenumber" name="L625" href="#L625">625</a> { 2.375 , 1.22225615758981 },
-<a class="jxr_linenumber" name="L626" href="#L626">626</a> { 2.5 , 1.329340388179137 },
-<a class="jxr_linenumber" name="L627" href="#L627">627</a> { 2.625 , 1.456933205091972 },
-<a class="jxr_linenumber" name="L628" href="#L628">628</a> { 2.75 , 1.608359421985546 },
-<a class="jxr_linenumber" name="L629" href="#L629">629</a> { 2.875 , 1.78771089889694 },
-<a class="jxr_linenumber" name="L630" href="#L630">630</a> { 3.0 , 2.0 },
-<a class="jxr_linenumber" name="L631" href="#L631">631</a> { 3.125 , 2.251353641828193 },
[... 886 lines stripped ...]