commons-notifications mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From er...@apache.org
Subject svn commit: r1002658 [35/35] - in /websites/production/commons/content/proper/commons-rng: ./ commons-rng-client-api/ commons-rng-client-api/apidocs/ commons-rng-client-api/apidocs/org/apache/commons/rng/ commons-rng-client-api/apidocs/org/apache/commo...
Date Mon, 12 Dec 2016 16:27:20 GMT
Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/PoissonSampler.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/PoissonSampler.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/PoissonSampler.html
Mon Dec 12 16:27:09 2016
@@ -74,119 +74,121 @@
 <a class="jxr_linenumber" name="L66" href="#L66">66</a>  
 <a class="jxr_linenumber" name="L67" href="#L67">67</a>          gaussian = <strong
class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/rng/sampling/distribution/BoxMullerGaussianSampler.html">BoxMullerGaussianSampler</a>(rng,
0, 1);
 <a class="jxr_linenumber" name="L68" href="#L68">68</a>          exponential
= <strong class="jxr_keyword">new</strong> <a href="../../../../../../org/apache/commons/rng/sampling/distribution/AhrensDieterExponentialSampler.html">AhrensDieterExponentialSampler</a>(rng,
1);
-<a class="jxr_linenumber" name="L69" href="#L69">69</a>          factorialLog
= InternalUtils.FactorialLog.create();
-<a class="jxr_linenumber" name="L70" href="#L70">70</a>      }
-<a class="jxr_linenumber" name="L71" href="#L71">71</a>  
-<a class="jxr_linenumber" name="L72" href="#L72">72</a>      <em class="jxr_javadoccomment">/**
{@inheritDoc} */</em>
-<a class="jxr_linenumber" name="L73" href="#L73">73</a>      @Override
-<a class="jxr_linenumber" name="L74" href="#L74">74</a>      <strong class="jxr_keyword">public</strong>
<strong class="jxr_keyword">int</strong> sample() {
-<a class="jxr_linenumber" name="L75" href="#L75">75</a>          <strong class="jxr_keyword">return</strong>
(<strong class="jxr_keyword">int</strong>) Math.min(nextPoisson(mean), Integer.MAX_VALUE);
-<a class="jxr_linenumber" name="L76" href="#L76">76</a>      }
-<a class="jxr_linenumber" name="L77" href="#L77">77</a>  
-<a class="jxr_linenumber" name="L78" href="#L78">78</a>      <em class="jxr_javadoccomment">/**
{@inheritDoc} */</em>
-<a class="jxr_linenumber" name="L79" href="#L79">79</a>      @Override
-<a class="jxr_linenumber" name="L80" href="#L80">80</a>      <strong class="jxr_keyword">public</strong>
String toString() {
-<a class="jxr_linenumber" name="L81" href="#L81">81</a>          <strong class="jxr_keyword">return</strong>
<span class="jxr_string">"Poisson deviate ["</span> + <strong class="jxr_keyword">super</strong>.toString()
+ <span class="jxr_string">"]"</span>;
-<a class="jxr_linenumber" name="L82" href="#L82">82</a>      }
-<a class="jxr_linenumber" name="L83" href="#L83">83</a>  
-<a class="jxr_linenumber" name="L84" href="#L84">84</a>      <em class="jxr_javadoccomment">/**</em>
-<a class="jxr_linenumber" name="L85" href="#L85">85</a>  <em class="jxr_javadoccomment">
    * @param meanPoisson Mean.</em>
-<a class="jxr_linenumber" name="L86" href="#L86">86</a>  <em class="jxr_javadoccomment">
    * @return the next sample.</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>
<strong class="jxr_keyword">long</strong> nextPoisson(<strong class="jxr_keyword">double</strong>
meanPoisson) {
-<a class="jxr_linenumber" name="L89" href="#L89">89</a>          <strong class="jxr_keyword">if</strong>
(meanPoisson &lt; PIVOT) {
-<a class="jxr_linenumber" name="L90" href="#L90">90</a>              <strong
class="jxr_keyword">double</strong> p = Math.exp(-meanPoisson);
-<a class="jxr_linenumber" name="L91" href="#L91">91</a>              <strong
class="jxr_keyword">long</strong> n = 0;
-<a class="jxr_linenumber" name="L92" href="#L92">92</a>              <strong
class="jxr_keyword">double</strong> r = 1;
-<a class="jxr_linenumber" name="L93" href="#L93">93</a>  
-<a class="jxr_linenumber" name="L94" href="#L94">94</a>              <strong
class="jxr_keyword">while</strong> (n &lt; 1000 * meanPoisson) {
-<a class="jxr_linenumber" name="L95" href="#L95">95</a>                  r *=
nextDouble();
-<a class="jxr_linenumber" name="L96" href="#L96">96</a>                  <strong
class="jxr_keyword">if</strong> (r &gt;= p) {
-<a class="jxr_linenumber" name="L97" href="#L97">97</a>                     
n++;
-<a class="jxr_linenumber" name="L98" href="#L98">98</a>                  } <strong
class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L99" href="#L99">99</a>                     
<strong class="jxr_keyword">break</strong>;
-<a class="jxr_linenumber" name="L100" href="#L100">100</a>                 }
-<a class="jxr_linenumber" name="L101" href="#L101">101</a>             }
-<a class="jxr_linenumber" name="L102" href="#L102">102</a>             <strong
class="jxr_keyword">return</strong> n;
-<a class="jxr_linenumber" name="L103" href="#L103">103</a>         } <strong
class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L104" href="#L104">104</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
lambda = Math.floor(meanPoisson);
-<a class="jxr_linenumber" name="L105" href="#L105">105</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
lambdaFractional = meanPoisson - lambda;
-<a class="jxr_linenumber" name="L106" href="#L106">106</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
logLambda = Math.log(lambda);
-<a class="jxr_linenumber" name="L107" href="#L107">107</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
logLambdaFactorial = factorialLog((<strong class="jxr_keyword">int</strong>) lambda);
-<a class="jxr_linenumber" name="L108" href="#L108">108</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong>
y2 = lambdaFractional &lt; Double.MIN_VALUE ? 0 : nextPoisson(lambdaFractional);
-<a class="jxr_linenumber" name="L109" href="#L109">109</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
delta = Math.sqrt(lambda * Math.log(32 * lambda / Math.PI + 1));
-<a class="jxr_linenumber" name="L110" href="#L110">110</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
halfDelta = delta / 2;
-<a class="jxr_linenumber" name="L111" href="#L111">111</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
twolpd = 2 * lambda + delta;
-<a class="jxr_linenumber" name="L112" href="#L112">112</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
a1 = Math.sqrt(Math.PI * twolpd) * Math.exp(1 / (8 * lambda));
-<a class="jxr_linenumber" name="L113" href="#L113">113</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
a2 = (twolpd / delta) * Math.exp(-delta * (1 + delta) / twolpd);
-<a class="jxr_linenumber" name="L114" href="#L114">114</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
aSum = a1 + a2 + 1;
-<a class="jxr_linenumber" name="L115" href="#L115">115</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
p1 = a1 / aSum;
-<a class="jxr_linenumber" name="L116" href="#L116">116</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
p2 = a2 / aSum;
-<a class="jxr_linenumber" name="L117" href="#L117">117</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
c1 = 1 / (8 * lambda);
-<a class="jxr_linenumber" name="L118" href="#L118">118</a> 
-<a class="jxr_linenumber" name="L119" href="#L119">119</a>             <strong
class="jxr_keyword">double</strong> x = 0;
-<a class="jxr_linenumber" name="L120" href="#L120">120</a>             <strong
class="jxr_keyword">double</strong> y = 0;
-<a class="jxr_linenumber" name="L121" href="#L121">121</a>             <strong
class="jxr_keyword">double</strong> v = 0;
-<a class="jxr_linenumber" name="L122" href="#L122">122</a>             <strong
class="jxr_keyword">int</strong> a = 0;
-<a class="jxr_linenumber" name="L123" href="#L123">123</a>             <strong
class="jxr_keyword">double</strong> t = 0;
-<a class="jxr_linenumber" name="L124" href="#L124">124</a>             <strong
class="jxr_keyword">double</strong> qr = 0;
-<a class="jxr_linenumber" name="L125" href="#L125">125</a>             <strong
class="jxr_keyword">double</strong> qa = 0;
-<a class="jxr_linenumber" name="L126" href="#L126">126</a>             <strong
class="jxr_keyword">while</strong> (<strong class="jxr_keyword">true</strong>)
{
-<a class="jxr_linenumber" name="L127" href="#L127">127</a>                 <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
u = nextDouble();
-<a class="jxr_linenumber" name="L128" href="#L128">128</a>                 <strong
class="jxr_keyword">if</strong> (u &lt;= p1) {
-<a class="jxr_linenumber" name="L129" href="#L129">129</a>                  
  <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
n = gaussian.sample();
-<a class="jxr_linenumber" name="L130" href="#L130">130</a>                  
  x = n * Math.sqrt(lambda + halfDelta) - 0.5d;
-<a class="jxr_linenumber" name="L131" href="#L131">131</a>                  
  <strong class="jxr_keyword">if</strong> (x &gt; delta || x &lt; -lambda)
{
-<a class="jxr_linenumber" name="L132" href="#L132">132</a>                  
      <strong class="jxr_keyword">continue</strong>;
-<a class="jxr_linenumber" name="L133" href="#L133">133</a>                  
  }
-<a class="jxr_linenumber" name="L134" href="#L134">134</a>                  
  y = x &lt; 0 ? Math.floor(x) : Math.ceil(x);
-<a class="jxr_linenumber" name="L135" href="#L135">135</a>                  
  <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
e = exponential.sample();
-<a class="jxr_linenumber" name="L136" href="#L136">136</a>                  
  v = -e - (n * n / 2) + c1;
-<a class="jxr_linenumber" name="L137" href="#L137">137</a>                 }
<strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L138" href="#L138">138</a>                  
  <strong class="jxr_keyword">if</strong> (u &gt; p1 + p2) {
-<a class="jxr_linenumber" name="L139" href="#L139">139</a>                  
      y = lambda;
-<a class="jxr_linenumber" name="L140" href="#L140">140</a>                  
      <strong class="jxr_keyword">break</strong>;
-<a class="jxr_linenumber" name="L141" href="#L141">141</a>                  
  } <strong class="jxr_keyword">else</strong> {
-<a class="jxr_linenumber" name="L142" href="#L142">142</a>                  
      x = delta + (twolpd / delta) * exponential.sample();
-<a class="jxr_linenumber" name="L143" href="#L143">143</a>                  
      y = Math.ceil(x);
-<a class="jxr_linenumber" name="L144" href="#L144">144</a>                  
      v = -exponential.sample() - delta * (x + 1) / twolpd;
-<a class="jxr_linenumber" name="L145" href="#L145">145</a>                  
  }
-<a class="jxr_linenumber" name="L146" href="#L146">146</a>                 }
-<a class="jxr_linenumber" name="L147" href="#L147">147</a>                 a
= x &lt; 0 ? 1 : 0;
-<a class="jxr_linenumber" name="L148" href="#L148">148</a>                 t
= y * (y + 1) / (2 * lambda);
-<a class="jxr_linenumber" name="L149" href="#L149">149</a>                 <strong
class="jxr_keyword">if</strong> (v &lt; -t &amp;&amp; a == 0) {
-<a class="jxr_linenumber" name="L150" href="#L150">150</a>                  
  y = lambda + y;
-<a class="jxr_linenumber" name="L151" href="#L151">151</a>                  
  <strong class="jxr_keyword">break</strong>;
-<a class="jxr_linenumber" name="L152" href="#L152">152</a>                 }
-<a class="jxr_linenumber" name="L153" href="#L153">153</a>                 qr
= t * ((2 * y + 1) / (6 * lambda) - 1);
-<a class="jxr_linenumber" name="L154" href="#L154">154</a>                 qa
= qr - (t * t) / (3 * (lambda + a * (y + 1)));
-<a class="jxr_linenumber" name="L155" href="#L155">155</a>                 <strong
class="jxr_keyword">if</strong> (v &lt; qa) {
-<a class="jxr_linenumber" name="L156" href="#L156">156</a>                  
  y = lambda + y;
-<a class="jxr_linenumber" name="L157" href="#L157">157</a>                  
  <strong class="jxr_keyword">break</strong>;
-<a class="jxr_linenumber" name="L158" href="#L158">158</a>                 }
-<a class="jxr_linenumber" name="L159" href="#L159">159</a>                 <strong
class="jxr_keyword">if</strong> (v &gt; qr) {
-<a class="jxr_linenumber" name="L160" href="#L160">160</a>                  
  <strong class="jxr_keyword">continue</strong>;
-<a class="jxr_linenumber" name="L161" href="#L161">161</a>                 }
-<a class="jxr_linenumber" name="L162" href="#L162">162</a>                 <strong
class="jxr_keyword">if</strong> (v &lt; y * logLambda - factorialLog((<strong
class="jxr_keyword">int</strong>) (y + lambda)) + logLambdaFactorial) {
-<a class="jxr_linenumber" name="L163" href="#L163">163</a>                  
  y = lambda + y;
-<a class="jxr_linenumber" name="L164" href="#L164">164</a>                  
  <strong class="jxr_keyword">break</strong>;
-<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>             <strong
class="jxr_keyword">return</strong> y2 + (<strong class="jxr_keyword">long</strong>)
y;
-<a class="jxr_linenumber" name="L168" href="#L168">168</a>         }
-<a class="jxr_linenumber" name="L169" href="#L169">169</a>     }
-<a class="jxr_linenumber" name="L170" href="#L170">170</a> 
-<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">
    * Compute the natural logarithm of the factorial of {@code n}.</em>
-<a class="jxr_linenumber" name="L173" href="#L173">173</a> <em class="jxr_javadoccomment">
    *</em>
-<a class="jxr_linenumber" name="L174" href="#L174">174</a> <em class="jxr_javadoccomment">
    * @param n Argument.</em>
-<a class="jxr_linenumber" name="L175" href="#L175">175</a> <em class="jxr_javadoccomment">
    * @return {@code log(n!)}</em>
-<a class="jxr_linenumber" name="L176" href="#L176">176</a> <em class="jxr_javadoccomment">
    * @throws IllegalArgumentException if {@code n &lt; 0}.</em>
-<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment">
    */</em>
-<a class="jxr_linenumber" name="L178" href="#L178">178</a>     <strong class="jxr_keyword">private</strong>
<strong class="jxr_keyword">double</strong> factorialLog(<strong class="jxr_keyword">int</strong>
n) {
-<a class="jxr_linenumber" name="L179" href="#L179">179</a>         <strong
class="jxr_keyword">return</strong> factorialLog.value(n);
-<a class="jxr_linenumber" name="L180" href="#L180">180</a>     }
-<a class="jxr_linenumber" name="L181" href="#L181">181</a> }
+<a class="jxr_linenumber" name="L69" href="#L69">69</a>          factorialLog
= mean &lt; PIVOT ?
+<a class="jxr_linenumber" name="L70" href="#L70">70</a>              <strong
class="jxr_keyword">null</strong> : <em class="jxr_comment">// Not used.</em>
+<a class="jxr_linenumber" name="L71" href="#L71">71</a>              InternalUtils.FactorialLog.create().withCache((<strong
class="jxr_keyword">int</strong>) Math.min(mean, 2 * PIVOT));
+<a class="jxr_linenumber" name="L72" href="#L72">72</a>      }
+<a class="jxr_linenumber" name="L73" href="#L73">73</a>  
+<a class="jxr_linenumber" name="L74" href="#L74">74</a>      <em class="jxr_javadoccomment">/**
{@inheritDoc} */</em>
+<a class="jxr_linenumber" name="L75" href="#L75">75</a>      @Override
+<a class="jxr_linenumber" name="L76" href="#L76">76</a>      <strong class="jxr_keyword">public</strong>
<strong class="jxr_keyword">int</strong> sample() {
+<a class="jxr_linenumber" name="L77" href="#L77">77</a>          <strong class="jxr_keyword">return</strong>
(<strong class="jxr_keyword">int</strong>) Math.min(nextPoisson(mean), Integer.MAX_VALUE);
+<a class="jxr_linenumber" name="L78" href="#L78">78</a>      }
+<a class="jxr_linenumber" name="L79" href="#L79">79</a>  
+<a class="jxr_linenumber" name="L80" href="#L80">80</a>      <em class="jxr_javadoccomment">/**
{@inheritDoc} */</em>
+<a class="jxr_linenumber" name="L81" href="#L81">81</a>      @Override
+<a class="jxr_linenumber" name="L82" href="#L82">82</a>      <strong class="jxr_keyword">public</strong>
String toString() {
+<a class="jxr_linenumber" name="L83" href="#L83">83</a>          <strong class="jxr_keyword">return</strong>
<span class="jxr_string">"Poisson deviate ["</span> + <strong class="jxr_keyword">super</strong>.toString()
+ <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>  
+<a class="jxr_linenumber" name="L86" href="#L86">86</a>      <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L87" href="#L87">87</a>  <em class="jxr_javadoccomment">
    * @param meanPoisson Mean.</em>
+<a class="jxr_linenumber" name="L88" href="#L88">88</a>  <em class="jxr_javadoccomment">
    * @return the next sample.</em>
+<a class="jxr_linenumber" name="L89" href="#L89">89</a>  <em class="jxr_javadoccomment">
    */</em>
+<a class="jxr_linenumber" name="L90" href="#L90">90</a>      <strong class="jxr_keyword">private</strong>
<strong class="jxr_keyword">long</strong> nextPoisson(<strong class="jxr_keyword">double</strong>
meanPoisson) {
+<a class="jxr_linenumber" name="L91" href="#L91">91</a>          <strong class="jxr_keyword">if</strong>
(meanPoisson &lt; PIVOT) {
+<a class="jxr_linenumber" name="L92" href="#L92">92</a>              <strong
class="jxr_keyword">double</strong> p = Math.exp(-meanPoisson);
+<a class="jxr_linenumber" name="L93" href="#L93">93</a>              <strong
class="jxr_keyword">long</strong> n = 0;
+<a class="jxr_linenumber" name="L94" href="#L94">94</a>              <strong
class="jxr_keyword">double</strong> r = 1;
+<a class="jxr_linenumber" name="L95" href="#L95">95</a>  
+<a class="jxr_linenumber" name="L96" href="#L96">96</a>              <strong
class="jxr_keyword">while</strong> (n &lt; 1000 * meanPoisson) {
+<a class="jxr_linenumber" name="L97" href="#L97">97</a>                  r *=
nextDouble();
+<a class="jxr_linenumber" name="L98" href="#L98">98</a>                  <strong
class="jxr_keyword">if</strong> (r &gt;= p) {
+<a class="jxr_linenumber" name="L99" href="#L99">99</a>                     
n++;
+<a class="jxr_linenumber" name="L100" href="#L100">100</a>                 }
<strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L101" href="#L101">101</a>                  
  <strong class="jxr_keyword">break</strong>;
+<a class="jxr_linenumber" name="L102" href="#L102">102</a>                 }
+<a class="jxr_linenumber" name="L103" href="#L103">103</a>             }
+<a class="jxr_linenumber" name="L104" href="#L104">104</a>             <strong
class="jxr_keyword">return</strong> n;
+<a class="jxr_linenumber" name="L105" href="#L105">105</a>         } <strong
class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L106" href="#L106">106</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
lambda = Math.floor(meanPoisson);
+<a class="jxr_linenumber" name="L107" href="#L107">107</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
lambdaFractional = meanPoisson - lambda;
+<a class="jxr_linenumber" name="L108" href="#L108">108</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
logLambda = Math.log(lambda);
+<a class="jxr_linenumber" name="L109" href="#L109">109</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
logLambdaFactorial = factorialLog((<strong class="jxr_keyword">int</strong>) lambda);
+<a class="jxr_linenumber" name="L110" href="#L110">110</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong>
y2 = lambdaFractional &lt; Double.MIN_VALUE ? 0 : nextPoisson(lambdaFractional);
+<a class="jxr_linenumber" name="L111" href="#L111">111</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
delta = Math.sqrt(lambda * Math.log(32 * lambda / Math.PI + 1));
+<a class="jxr_linenumber" name="L112" href="#L112">112</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
halfDelta = delta / 2;
+<a class="jxr_linenumber" name="L113" href="#L113">113</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
twolpd = 2 * lambda + delta;
+<a class="jxr_linenumber" name="L114" href="#L114">114</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
a1 = Math.sqrt(Math.PI * twolpd) * Math.exp(1 / (8 * lambda));
+<a class="jxr_linenumber" name="L115" href="#L115">115</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
a2 = (twolpd / delta) * Math.exp(-delta * (1 + delta) / twolpd);
+<a class="jxr_linenumber" name="L116" href="#L116">116</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
aSum = a1 + a2 + 1;
+<a class="jxr_linenumber" name="L117" href="#L117">117</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
p1 = a1 / aSum;
+<a class="jxr_linenumber" name="L118" href="#L118">118</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
p2 = a2 / aSum;
+<a class="jxr_linenumber" name="L119" href="#L119">119</a>             <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
c1 = 1 / (8 * lambda);
+<a class="jxr_linenumber" name="L120" href="#L120">120</a> 
+<a class="jxr_linenumber" name="L121" href="#L121">121</a>             <strong
class="jxr_keyword">double</strong> x = 0;
+<a class="jxr_linenumber" name="L122" href="#L122">122</a>             <strong
class="jxr_keyword">double</strong> y = 0;
+<a class="jxr_linenumber" name="L123" href="#L123">123</a>             <strong
class="jxr_keyword">double</strong> v = 0;
+<a class="jxr_linenumber" name="L124" href="#L124">124</a>             <strong
class="jxr_keyword">int</strong> a = 0;
+<a class="jxr_linenumber" name="L125" href="#L125">125</a>             <strong
class="jxr_keyword">double</strong> t = 0;
+<a class="jxr_linenumber" name="L126" href="#L126">126</a>             <strong
class="jxr_keyword">double</strong> qr = 0;
+<a class="jxr_linenumber" name="L127" href="#L127">127</a>             <strong
class="jxr_keyword">double</strong> qa = 0;
+<a class="jxr_linenumber" name="L128" href="#L128">128</a>             <strong
class="jxr_keyword">while</strong> (<strong class="jxr_keyword">true</strong>)
{
+<a class="jxr_linenumber" name="L129" href="#L129">129</a>                 <strong
class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
u = nextDouble();
+<a class="jxr_linenumber" name="L130" href="#L130">130</a>                 <strong
class="jxr_keyword">if</strong> (u &lt;= p1) {
+<a class="jxr_linenumber" name="L131" href="#L131">131</a>                  
  <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
n = gaussian.sample();
+<a class="jxr_linenumber" name="L132" href="#L132">132</a>                  
  x = n * Math.sqrt(lambda + halfDelta) - 0.5d;
+<a class="jxr_linenumber" name="L133" href="#L133">133</a>                  
  <strong class="jxr_keyword">if</strong> (x &gt; delta || x &lt; -lambda)
{
+<a class="jxr_linenumber" name="L134" href="#L134">134</a>                  
      <strong class="jxr_keyword">continue</strong>;
+<a class="jxr_linenumber" name="L135" href="#L135">135</a>                  
  }
+<a class="jxr_linenumber" name="L136" href="#L136">136</a>                  
  y = x &lt; 0 ? Math.floor(x) : Math.ceil(x);
+<a class="jxr_linenumber" name="L137" href="#L137">137</a>                  
  <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>
e = exponential.sample();
+<a class="jxr_linenumber" name="L138" href="#L138">138</a>                  
  v = -e - 0.5 * n * n + c1;
+<a class="jxr_linenumber" name="L139" href="#L139">139</a>                 }
<strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L140" href="#L140">140</a>                  
  <strong class="jxr_keyword">if</strong> (u &gt; p1 + p2) {
+<a class="jxr_linenumber" name="L141" href="#L141">141</a>                  
      y = lambda;
+<a class="jxr_linenumber" name="L142" href="#L142">142</a>                  
      <strong class="jxr_keyword">break</strong>;
+<a class="jxr_linenumber" name="L143" href="#L143">143</a>                  
  } <strong class="jxr_keyword">else</strong> {
+<a class="jxr_linenumber" name="L144" href="#L144">144</a>                  
      x = delta + (twolpd / delta) * exponential.sample();
+<a class="jxr_linenumber" name="L145" href="#L145">145</a>                  
      y = Math.ceil(x);
+<a class="jxr_linenumber" name="L146" href="#L146">146</a>                  
      v = -exponential.sample() - delta * (x + 1) / twolpd;
+<a class="jxr_linenumber" name="L147" href="#L147">147</a>                  
  }
+<a class="jxr_linenumber" name="L148" href="#L148">148</a>                 }
+<a class="jxr_linenumber" name="L149" href="#L149">149</a>                 a
= x &lt; 0 ? 1 : 0;
+<a class="jxr_linenumber" name="L150" href="#L150">150</a>                 t
= y * (y + 1) / (2 * lambda);
+<a class="jxr_linenumber" name="L151" href="#L151">151</a>                 <strong
class="jxr_keyword">if</strong> (v &lt; -t &amp;&amp; a == 0) {
+<a class="jxr_linenumber" name="L152" href="#L152">152</a>                  
  y = lambda + y;
+<a class="jxr_linenumber" name="L153" href="#L153">153</a>                  
  <strong class="jxr_keyword">break</strong>;
+<a class="jxr_linenumber" name="L154" href="#L154">154</a>                 }
+<a class="jxr_linenumber" name="L155" href="#L155">155</a>                 qr
= t * ((2 * y + 1) / (6 * lambda) - 1);
+<a class="jxr_linenumber" name="L156" href="#L156">156</a>                 qa
= qr - (t * t) / (3 * (lambda + a * (y + 1)));
+<a class="jxr_linenumber" name="L157" href="#L157">157</a>                 <strong
class="jxr_keyword">if</strong> (v &lt; qa) {
+<a class="jxr_linenumber" name="L158" href="#L158">158</a>                  
  y = lambda + y;
+<a class="jxr_linenumber" name="L159" href="#L159">159</a>                  
  <strong class="jxr_keyword">break</strong>;
+<a class="jxr_linenumber" name="L160" href="#L160">160</a>                 }
+<a class="jxr_linenumber" name="L161" href="#L161">161</a>                 <strong
class="jxr_keyword">if</strong> (v &gt; qr) {
+<a class="jxr_linenumber" name="L162" href="#L162">162</a>                  
  <strong class="jxr_keyword">continue</strong>;
+<a class="jxr_linenumber" name="L163" href="#L163">163</a>                 }
+<a class="jxr_linenumber" name="L164" href="#L164">164</a>                 <strong
class="jxr_keyword">if</strong> (v &lt; y * logLambda - factorialLog((<strong
class="jxr_keyword">int</strong>) (y + lambda)) + logLambdaFactorial) {
+<a class="jxr_linenumber" name="L165" href="#L165">165</a>                  
  y = lambda + y;
+<a class="jxr_linenumber" name="L166" href="#L166">166</a>                  
  <strong class="jxr_keyword">break</strong>;
+<a class="jxr_linenumber" name="L167" href="#L167">167</a>                 }
+<a class="jxr_linenumber" name="L168" href="#L168">168</a>             }
+<a class="jxr_linenumber" name="L169" href="#L169">169</a>             <strong
class="jxr_keyword">return</strong> y2 + (<strong class="jxr_keyword">long</strong>)
y;
+<a class="jxr_linenumber" name="L170" href="#L170">170</a>         }
+<a class="jxr_linenumber" name="L171" href="#L171">171</a>     }
+<a class="jxr_linenumber" name="L172" href="#L172">172</a> 
+<a class="jxr_linenumber" name="L173" href="#L173">173</a>     <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L174" href="#L174">174</a> <em class="jxr_javadoccomment">
    * Compute the natural logarithm of the factorial of {@code n}.</em>
+<a class="jxr_linenumber" name="L175" href="#L175">175</a> <em class="jxr_javadoccomment">
    *</em>
+<a class="jxr_linenumber" name="L176" href="#L176">176</a> <em class="jxr_javadoccomment">
    * @param n Argument.</em>
+<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment">
    * @return {@code log(n!)}</em>
+<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_javadoccomment">
    * @throws IllegalArgumentException if {@code n &lt; 0}.</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>     <strong class="jxr_keyword">private</strong>
<strong class="jxr_keyword">double</strong> factorialLog(<strong class="jxr_keyword">int</strong>
n) {
+<a class="jxr_linenumber" name="L181" href="#L181">181</a>         <strong
class="jxr_keyword">return</strong> factorialLog.value(n);
+<a class="jxr_linenumber" name="L182" href="#L182">182</a>     }
+<a class="jxr_linenumber" name="L183" href="#L183">183</a> }
 </pre>
 <hr/>
 <div id="footer">Copyright &#169; 2016 <a href="https://www.apache.org/">The
Apache Software Foundation</a>. All rights reserved.</div>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/RejectionInversionZipfSampler.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/RejectionInversionZipfSampler.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/RejectionInversionZipfSampler.html
Mon Dec 12 16:27:09 2016
@@ -77,7 +77,7 @@
 <a class="jxr_linenumber" name="L69" href="#L69">69</a>  <em class="jxr_javadoccomment">
    * &lt;blockquote&gt;</em>
 <a class="jxr_linenumber" name="L70" href="#L70">70</a>  <em class="jxr_javadoccomment">
    *   Wolfgang Hörmann and Gerhard Derflinger.</em>
 <a class="jxr_linenumber" name="L71" href="#L71">71</a>  <em class="jxr_javadoccomment">
    *   &lt;i&gt;"Rejection-inversion to generate variates from monotone discrete</em>
-<a class="jxr_linenumber" name="L72" href="#L72">72</a>  <em class="jxr_javadoccomment">
    *    distributions",&lt;i&gt;&lt;br&gt;</em>
+<a class="jxr_linenumber" name="L72" href="#L72">72</a>  <em class="jxr_javadoccomment">
    *    distributions",&lt;/i&gt;&lt;br&gt;</em>
 <a class="jxr_linenumber" name="L73" href="#L73">73</a>  <em class="jxr_javadoccomment">
    *   &lt;strong&gt;ACM Transactions on Modeling and Computer Simulation&lt;/strong&gt;
(TOMACS) 6.3 (1996): 169-184.</em>
 <a class="jxr_linenumber" name="L74" href="#L74">74</a>  <em class="jxr_javadoccomment">
    * &lt;/blockquote&gt;</em>
 <a class="jxr_linenumber" name="L75" href="#L75">75</a>  <em class="jxr_javadoccomment">
    */</em>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/package-frame.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/package-frame.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.sampling.distribution</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.sampling.distribution</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/package-summary.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/distribution/package-summary.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.sampling.distribution</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.sampling.distribution</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/package-frame.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/package-frame.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.sampling</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.sampling</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/package-summary.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/sampling/package-summary.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.sampling</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.sampling</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/internal/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/internal/package-frame.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/internal/package-frame.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.simple.internal</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.simple.internal</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/internal/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/internal/package-summary.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/internal/package-summary.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.simple.internal</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.simple.internal</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/package-frame.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/package-frame.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.simple</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.simple</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/package-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/package-summary.html
(original)
+++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/simple/package-summary.html
Mon Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference Package org.apache.commons.rng.simple</title>
+		<title>Apache Commons RNG 1.0 Reference Package org.apache.commons.rng.simple</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../stylesheet.css" title="style"
/>
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/overview-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/overview-frame.html (original)
+++ websites/production/commons/content/proper/commons-rng/xref/overview-frame.html Mon Dec
12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference</title>
+		<title>Apache Commons RNG 1.0 Reference</title>
 		<link rel="stylesheet" type="text/css" href="stylesheet.css" title="style" />
 	</head>
 	<body>

Modified: websites/production/commons/content/proper/commons-rng/xref/overview-summary.html
==============================================================================
--- websites/production/commons/content/proper/commons-rng/xref/overview-summary.html (original)
+++ websites/production/commons/content/proper/commons-rng/xref/overview-summary.html Mon
Dec 12 16:27:09 2016
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Apache Commons Parent RNG 1.0-SNAPSHOT Reference</title>
+		<title>Apache Commons RNG 1.0 Reference</title>
 		<link rel="stylesheet" type="text/css" href="stylesheet.css" title="style" />
 	</head>
 	<body>
@@ -24,7 +24,7 @@
 	        </ul>
       	</div>
 		
-		      	<h2>Apache Commons Parent RNG 1.0-SNAPSHOT Reference</h2>
+		      	<h2>Apache Commons RNG 1.0 Reference</h2>
 
 		<table class="summary">
         	<thead>



Mime
View raw message