Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ChiSquaredDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ChiSquaredDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ChiSquaredDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,167 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>ChiSquaredDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/ChiSquaredDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> +<a class="jxr_linenumber" name="L18" href="#L18">18</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L19" href="#L19">19</a> +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <strong class="jxr_keyword">import</strong> java.util.stream.Stream; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.Arguments; +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.MethodSource; +<a class="jxr_linenumber" name="L24" href="#L24">24</a> +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <em class="jxr_javadoccomment"> * Test cases for {@link ChiSquaredDistribution}.</em> +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em> +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <strong class="jxr_keyword">class</strong> <a name="ChiSquaredDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/ChiSquaredDistributionTest.html#ChiSquaredDistributionTest">ChiSquaredDistributionTest</a> <strong class="jxr_keyword">extends</strong> <a name="BaseContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/BaseContinuousDistributionTest.html#BaseContinuousDistributionTest">BaseContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L30" href="#L30">30</a> @Override +<a class="jxr_linenumber" name="L31" href="#L31">31</a> ContinuousDistribution makeDistribution(Object... parameters) { +<a class="jxr_linenumber" name="L32" href="#L32">32</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> df = (Double) parameters[0]; +<a class="jxr_linenumber" name="L33" href="#L33">33</a> <strong class="jxr_keyword">return</strong> ChiSquaredDistribution.of(df); +<a class="jxr_linenumber" name="L34" href="#L34">34</a> } +<a class="jxr_linenumber" name="L35" href="#L35">35</a> +<a class="jxr_linenumber" name="L36" href="#L36">36</a> @Override +<a class="jxr_linenumber" name="L37" href="#L37">37</a> Object[][] makeInvalidParameters() { +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] { +<a class="jxr_linenumber" name="L39" href="#L39">39</a> {0.0}, +<a class="jxr_linenumber" name="L40" href="#L40">40</a> {-0.1} +<a class="jxr_linenumber" name="L41" href="#L41">41</a> }; +<a class="jxr_linenumber" name="L42" href="#L42">42</a> } +<a class="jxr_linenumber" name="L43" href="#L43">43</a> +<a class="jxr_linenumber" name="L44" href="#L44">44</a> @Override +<a class="jxr_linenumber" name="L45" href="#L45">45</a> String[] getParameterNames() { +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"DegreesOfFreedom"</span>}; +<a class="jxr_linenumber" name="L47" href="#L47">47</a> } +<a class="jxr_linenumber" name="L48" href="#L48">48</a> +<a class="jxr_linenumber" name="L49" href="#L49">49</a> @Override +<a class="jxr_linenumber" name="L50" href="#L50">50</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() { +<a class="jxr_linenumber" name="L51" href="#L51">51</a> <strong class="jxr_keyword">return</strong> 5e-15; +<a class="jxr_linenumber" name="L52" href="#L52">52</a> } +<a class="jxr_linenumber" name="L53" href="#L53">53</a> +<a class="jxr_linenumber" name="L54" href="#L54">54</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em> +<a class="jxr_linenumber" name="L55" href="#L55">55</a> +<a class="jxr_linenumber" name="L56" href="#L56">56</a> @ParameterizedTest +<a class="jxr_linenumber" name="L57" href="#L57">57</a> @MethodSource +<a class="jxr_linenumber" name="L58" href="#L58">58</a> <strong class="jxr_keyword">void</strong> testAdditionalDensity(<strong class="jxr_keyword">double</strong> df, <strong class="jxr_keyword">double</strong>[] points, <strong class="jxr_keyword">double</strong>[] values, <a name="DoubleTolerance" href="../../../../../org/apache/commons/statistics/distribution/DoubleTolerance.html#DoubleTolerance">DoubleTolerance</a> tol) { +<a class="jxr_linenumber" name="L59" href="#L59">59</a> testDensity(ChiSquaredDistribution.of(df), points, values, tol); +<a class="jxr_linenumber" name="L60" href="#L60">60</a> } +<a class="jxr_linenumber" name="L61" href="#L61">61</a> +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <strong class="jxr_keyword">static</strong> Stream<Arguments> testAdditionalDensity() { +<a class="jxr_linenumber" name="L63" href="#L63">63</a> <em class="jxr_comment">// Values have many digits above the decimal point so use relative tolerance</em> +<a class="jxr_linenumber" name="L64" href="#L64">64</a> <strong class="jxr_keyword">final</strong> <a name="DoubleTolerance" href="../../../../../org/apache/commons/statistics/distribution/DoubleTolerance.html#DoubleTolerance">DoubleTolerance</a> tol = DoubleTolerances.relative(5e-14); +<a class="jxr_linenumber" name="L65" href="#L65">65</a> +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <em class="jxr_comment">// R 2.5:</em> +<a class="jxr_linenumber" name="L67" href="#L67">67</a> <em class="jxr_comment">// x <- c(-0.1, 1e-6, 0.5, 1, 2, 5)</em> +<a class="jxr_linenumber" name="L68" href="#L68">68</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] x = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[]{-0.1, 1e-6, 0.5, 1, 2, 5}; +<a class="jxr_linenumber" name="L69" href="#L69">69</a> <strong class="jxr_keyword">return</strong> Stream.of( +<a class="jxr_linenumber" name="L70" href="#L70">70</a> <em class="jxr_comment">// print(dchisq(x, df=1), digits=17)</em> +<a class="jxr_linenumber" name="L71" href="#L71">71</a> Arguments.of(1, x, <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L72" href="#L72">72</a> 0, 398.942080930342626743, 0.439391289467722435, 0.241970724519143365, +<a class="jxr_linenumber" name="L73" href="#L73">73</a> 0.103776874355148693, 0.014644982561926489}, tol), +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <em class="jxr_comment">// print(dchisq(x, df=0.1), digits=17)</em> +<a class="jxr_linenumber" name="L75" href="#L75">75</a> Arguments.of(0.1, x, <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L76" href="#L76">76</a> 0, 2.4864539972849805e+04, 7.4642387316120481e-02, +<a class="jxr_linenumber" name="L77" href="#L77">77</a> 3.0090777182393683e-02, 9.4472991589506262e-03, 8.8271993957607896e-04}, tol), +<a class="jxr_linenumber" name="L78" href="#L78">78</a> <em class="jxr_comment">// print(dchisq(x, df=2), digits=17)</em> +<a class="jxr_linenumber" name="L79" href="#L79">79</a> Arguments.of(2, x, <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L80" href="#L80">80</a> 0, 0.49999975000006253, 0.38940039153570244, +<a class="jxr_linenumber" name="L81" href="#L81">81</a> 0.30326532985631671, 0.18393972058572117, 0.04104249931194940}, tol), +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <em class="jxr_comment">// print(dchisq(x, df=10), digits=17)</em> +<a class="jxr_linenumber" name="L83" href="#L83">83</a> Arguments.of(10, x, <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L84" href="#L84">84</a> 0, 1.3020826822918329e-27, 6.3378969976514082e-05, +<a class="jxr_linenumber" name="L85" href="#L85">85</a> 7.8975346316749191e-04, 7.6641550244050524e-03, 6.6800942890542614e-02}, tol), +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <em class="jxr_comment">// print(dchisq(x, df=100), digits=17)</em> +<a class="jxr_linenumber" name="L87" href="#L87">87</a> Arguments.of(100, x, <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L88" href="#L88">88</a> 0, 0.0000000000000000e+00, 2.0200026568141969e-93, +<a class="jxr_linenumber" name="L89" href="#L89">89</a> 8.8562141121618944e-79, 3.0239224849774644e-64, 2.1290671364111626e-45}, tol), +<a class="jxr_linenumber" name="L90" href="#L90">90</a> +<a class="jxr_linenumber" name="L91" href="#L91">91</a> <em class="jxr_comment">// Progressively larger degrees of freedom (df) with values around the mean (df).</em> +<a class="jxr_linenumber" name="L92" href="#L92">92</a> <em class="jxr_comment">// Note that the CDF tends towards a step function (0 to 1) around the mean, the</em> +<a class="jxr_linenumber" name="L93" href="#L93">93</a> <em class="jxr_comment">// density is tiny and the computation has large cancellation leading to inaccuracy.</em> +<a class="jxr_linenumber" name="L94" href="#L94">94</a> <em class="jxr_comment">// Note that R's dchisq and matlab's chi2pdf computations are close and have a</em> +<a class="jxr_linenumber" name="L95" href="#L95">95</a> <em class="jxr_comment">// similar error to the scipy result as the current java code.</em> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> +<a class="jxr_linenumber" name="L97" href="#L97">97</a> <em class="jxr_comment">// scipy.stats 1.9.1</em> +<a class="jxr_linenumber" name="L98" href="#L98">98</a> <em class="jxr_comment">// chi2.pdf([250, 500, 1000, 1500, 2000, 2500], 1000)</em> +<a class="jxr_linenumber" name="L99" href="#L99">99</a> Arguments.of(1000, +<a class="jxr_linenumber" name="L100" href="#L100">100</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {250, 500, 1000, 1500, 2000, 2500}, +<a class="jxr_linenumber" name="L101" href="#L101">101</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L102" href="#L102">102</a> 2.4144472784936886e-140, 2.0416219443308211e-044, +<a class="jxr_linenumber" name="L103" href="#L103">103</a> 8.9191339347531283e-003, 1.7629519620803219e-023, +<a class="jxr_linenumber" name="L104" href="#L104">104</a> 1.0400388688836408e-069, 6.3317766286480137e-130 +<a class="jxr_linenumber" name="L105" href="#L105">105</a> }, DoubleTolerances.relative(5e-13)), +<a class="jxr_linenumber" name="L106" href="#L106">106</a> <em class="jxr_comment">// chi2.pdf([7000, 8000, 9000, 10000, 11000, 12000, 14000], 10000)</em> +<a class="jxr_linenumber" name="L107" href="#L107">107</a> Arguments.of(10000, +<a class="jxr_linenumber" name="L108" href="#L108">108</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {7000, 8000, 9000, 10000, 11000, 12000, 14000}, +<a class="jxr_linenumber" name="L109" href="#L109">109</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L110" href="#L110">110</a> 3.4451863344803051e-126, 1.9575583029092260e-053, +<a class="jxr_linenumber" name="L111" href="#L111">111</a> 7.1768433769351565e-015, 2.8209009023369056e-003, +<a class="jxr_linenumber" name="L112" href="#L112">112</a> 1.6794999068429000e-013, 9.6151246453108889e-042, +<a class="jxr_linenumber" name="L113" href="#L113">113</a> 2.2671193242727422e-141 +<a class="jxr_linenumber" name="L114" href="#L114">114</a> }, DoubleTolerances.relative(5e-12)), +<a class="jxr_linenumber" name="L115" href="#L115">115</a> <em class="jxr_comment">// chi2.pdf([90000, 93000, 97000, 100000, 103000, 106000, 110000], 100000)</em> +<a class="jxr_linenumber" name="L116" href="#L116">116</a> Arguments.of(100000, +<a class="jxr_linenumber" name="L117" href="#L117">117</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {90000, 93000, 97000, 100000, 103000, 106000, 110000}, +<a class="jxr_linenumber" name="L118" href="#L118">118</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L119" href="#L119">119</a> 3.9267505859463047e-120, 1.4455254881569810e-059, +<a class="jxr_linenumber" name="L120" href="#L120">120</a> 9.8188082087821045e-014, 8.9206057128873189e-004, +<a class="jxr_linenumber" name="L121" href="#L121">121</a> 2.2754616642908043e-013, 2.1622505071910235e-041, +<a class="jxr_linenumber" name="L122" href="#L122">122</a> 1.1771994592645341e-105 +<a class="jxr_linenumber" name="L123" href="#L123">123</a> }, DoubleTolerances.relative(2e-10)), +<a class="jxr_linenumber" name="L124" href="#L124">124</a> <em class="jxr_comment">// chi2.pdf([970000, 980000, 990000, 1000000, 1010000, 1020000, 1030000], 1000000)</em> +<a class="jxr_linenumber" name="L125" href="#L125">125</a> Arguments.of(1000000, +<a class="jxr_linenumber" name="L126" href="#L126">126</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {970000, 980000, 990000, 1000000, 1010000, 1020000, 1030000}, +<a class="jxr_linenumber" name="L127" href="#L127">127</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L128" href="#L128">128</a> 5.5973818150499247e-104, 2.7658855125007757e-048, +<a class="jxr_linenumber" name="L129" href="#L129">129</a> 3.3455543700508753e-015, 2.8209474455352760e-004, +<a class="jxr_linenumber" name="L130" href="#L130">130</a> 4.5767314194128692e-015, 3.8269849102866545e-047, +<a class="jxr_linenumber" name="L131" href="#L131">131</a> 4.2922183625107689e-100 +<a class="jxr_linenumber" name="L132" href="#L132">132</a> }, DoubleTolerances.relative(1e-9)), +<a class="jxr_linenumber" name="L133" href="#L133">133</a> <em class="jxr_comment">// chi2.pdf([9889000, 9890000, 9900000, 10000000, 10100000, 10150000], 10000000)</em> +<a class="jxr_linenumber" name="L134" href="#L134">134</a> Arguments.of(10000000, +<a class="jxr_linenumber" name="L135" href="#L135">135</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {9889000, 9890000, 9900000, 10000000, 10100000, 10150000}, +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L137" href="#L137">137</a> 1.5256172522921843e-139, 4.0708897874504053e-137, +<a class="jxr_linenumber" name="L138" href="#L138">138</a> 4.4858296398343578e-114, 8.9206205026501391e-005, +<a class="jxr_linenumber" name="L139" href="#L139">139</a> 1.2327974895005824e-112, 1.1722829145179983e-246 +<a class="jxr_linenumber" name="L140" href="#L140">140</a> }, DoubleTolerances.relative(2e-8)), +<a class="jxr_linenumber" name="L141" href="#L141">141</a> <em class="jxr_comment">// chi2.pdf([1e9-1e7, 1e9-1e6-1e5, 1e9-1e6-1e3, 1e9-1e6, 1e9, 1e9+1e6, 1e9+1e6+1e3, 1e9+1e6+1e5, 1e9+1e7], 1e9)</em> +<a class="jxr_linenumber" name="L142" href="#L142">142</a> Arguments.of(1e9, +<a class="jxr_linenumber" name="L143" href="#L143">143</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {1e9 - 1e7, 1e9 - 1e6 - 1e5, 1e9 - 1e6 - 1e3, 1e9 - 1e6, 1e9, +<a class="jxr_linenumber" name="L144" href="#L144">144</a> 1e9 + 1e6, 1e9 + 1e6 + 1e3, 1e9 + 1e6 + 1e5, 1e9 + 1e7}, +<a class="jxr_linenumber" name="L145" href="#L145">145</a> <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { +<a class="jxr_linenumber" name="L146" href="#L146">146</a> 0.0000000000000000e+000, 3.0225312085986285e-137, +<a class="jxr_linenumber" name="L147" href="#L147">147</a> 1.2226420014038719e-114, 2.0173074570506826e-114, +<a class="jxr_linenumber" name="L148" href="#L148">148</a> 8.9206087939511128e-006, 2.8097537346819163e-114, +<a class="jxr_linenumber" name="L149" href="#L149">149</a> 1.7046276559389835e-114, 4.7000061695431404e-137, +<a class="jxr_linenumber" name="L150" href="#L150">150</a> 0.0000000000000000e+000 +<a class="jxr_linenumber" name="L151" href="#L151">151</a> }, DoubleTolerances.relative(3e-6)) +<a class="jxr_linenumber" name="L152" href="#L152">152</a> ); +<a class="jxr_linenumber" name="L153" href="#L153">153</a> } +<a class="jxr_linenumber" name="L154" href="#L154">154</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html>
Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ContinuousDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ContinuousDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/ContinuousDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,166 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>ContinuousDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/ContinuousDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <strong class="jxr_keyword">import</strong> java.util.concurrent.ThreadLocalRandom; +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <strong class="jxr_keyword">import</strong> org.apache.commons.rng.UniformRandomProvider; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest; +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.ValueSource; +<a class="jxr_linenumber" name="L25" href="#L25">25</a> +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <em class="jxr_javadoccomment"> * Test default implementations in the {@link ContinuousDistribution} interface.</em> +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <strong class="jxr_keyword">class</strong> <a name="ContinuousDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/ContinuousDistributionTest.html#ContinuousDistributionTest">ContinuousDistributionTest</a> { +<a class="jxr_linenumber" name="L30" href="#L30">30</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <em class="jxr_javadoccomment"> * Test the default interface methods.</em> +<a class="jxr_linenumber" name="L32" href="#L32">32</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L33" href="#L33">33</a> @Test +<a class="jxr_linenumber" name="L34" href="#L34">34</a> <strong class="jxr_keyword">void</strong> testDefaultMethods() { +<a class="jxr_linenumber" name="L35" href="#L35">35</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> high = 0.54; +<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> low = 0.313; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">final</strong> ContinuousDistribution dist = <strong class="jxr_keyword">new</strong> ContinuousDistribution() { +<a class="jxr_linenumber" name="L39" href="#L39">39</a> @Override +<a class="jxr_linenumber" name="L40" href="#L40">40</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> inverseCumulativeProbability(<strong class="jxr_keyword">double</strong> p) { +<a class="jxr_linenumber" name="L41" href="#L41">41</a> <em class="jxr_comment">// For the default inverseSurvivalProbability(double) method</em> +<a class="jxr_linenumber" name="L42" href="#L42">42</a> <strong class="jxr_keyword">return</strong> 10 * p; +<a class="jxr_linenumber" name="L43" href="#L43">43</a> } +<a class="jxr_linenumber" name="L44" href="#L44">44</a> @Override +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getVariance() { +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <strong class="jxr_keyword">return</strong> 0; +<a class="jxr_linenumber" name="L47" href="#L47">47</a> } +<a class="jxr_linenumber" name="L48" href="#L48">48</a> @Override +<a class="jxr_linenumber" name="L49" href="#L49">49</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getSupportUpperBound() { +<a class="jxr_linenumber" name="L50" href="#L50">50</a> <strong class="jxr_keyword">return</strong> 0; +<a class="jxr_linenumber" name="L51" href="#L51">51</a> } +<a class="jxr_linenumber" name="L52" href="#L52">52</a> @Override +<a class="jxr_linenumber" name="L53" href="#L53">53</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getSupportLowerBound() { +<a class="jxr_linenumber" name="L54" href="#L54">54</a> <strong class="jxr_keyword">return</strong> 0; +<a class="jxr_linenumber" name="L55" href="#L55">55</a> } +<a class="jxr_linenumber" name="L56" href="#L56">56</a> @Override +<a class="jxr_linenumber" name="L57" href="#L57">57</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getMean() { +<a class="jxr_linenumber" name="L58" href="#L58">58</a> <strong class="jxr_keyword">return</strong> 0; +<a class="jxr_linenumber" name="L59" href="#L59">59</a> } +<a class="jxr_linenumber" name="L60" href="#L60">60</a> @Override +<a class="jxr_linenumber" name="L61" href="#L61">61</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> density(<strong class="jxr_keyword">double</strong> x) { +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <em class="jxr_comment">// Return input value for testing</em> +<a class="jxr_linenumber" name="L63" href="#L63">63</a> <strong class="jxr_keyword">return</strong> x; +<a class="jxr_linenumber" name="L64" href="#L64">64</a> } +<a class="jxr_linenumber" name="L65" href="#L65">65</a> @Override +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> cumulativeProbability(<strong class="jxr_keyword">double</strong> x) { +<a class="jxr_linenumber" name="L67" href="#L67">67</a> <em class="jxr_comment">// For the default probability(double, double) method</em> +<a class="jxr_linenumber" name="L68" href="#L68">68</a> <strong class="jxr_keyword">return</strong> x > 1 ? high : low; +<a class="jxr_linenumber" name="L69" href="#L69">69</a> } +<a class="jxr_linenumber" name="L70" href="#L70">70</a> @Override +<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">public</strong> Sampler createSampler(UniformRandomProvider rng) { +<a class="jxr_linenumber" name="L72" href="#L72">72</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">null</strong>; +<a class="jxr_linenumber" name="L73" href="#L73">73</a> } +<a class="jxr_linenumber" name="L74" href="#L74">74</a> }; +<a class="jxr_linenumber" name="L75" href="#L75">75</a> +<a class="jxr_linenumber" name="L76" href="#L76">76</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x : <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {Double.NaN, Double.POSITIVE_INFINITY, +<a class="jxr_linenumber" name="L77" href="#L77">77</a> Double.NEGATIVE_INFINITY, 0, 1, 0.123}) { +<a class="jxr_linenumber" name="L78" href="#L78">78</a> <em class="jxr_comment">// Return the log of the density</em> +<a class="jxr_linenumber" name="L79" href="#L79">79</a> Assertions.assertEquals(Math.log(x), dist.logDensity(x)); +<a class="jxr_linenumber" name="L80" href="#L80">80</a> } +<a class="jxr_linenumber" name="L81" href="#L81">81</a> +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <em class="jxr_comment">// Should throw for bad range</em> +<a class="jxr_linenumber" name="L83" href="#L83">83</a> Assertions.assertThrows(DistributionException.<strong class="jxr_keyword">class</strong>, () -> dist.probability(0.5, 0.4)); +<a class="jxr_linenumber" name="L84" href="#L84">84</a> Assertions.assertEquals(high - low, dist.probability(0.5, 1.5)); +<a class="jxr_linenumber" name="L85" href="#L85">85</a> Assertions.assertEquals(high - low, dist.probability(0.5, 1.5)); +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p : <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {0.2, 0.5, 0.7}) { +<a class="jxr_linenumber" name="L87" href="#L87">87</a> Assertions.assertEquals(dist.inverseCumulativeProbability(1 - p), +<a class="jxr_linenumber" name="L88" href="#L88">88</a> dist.inverseSurvivalProbability(p)); +<a class="jxr_linenumber" name="L89" href="#L89">89</a> } +<a class="jxr_linenumber" name="L90" href="#L90">90</a> } +<a class="jxr_linenumber" name="L91" href="#L91">91</a> +<a class="jxr_linenumber" name="L92" href="#L92">92</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L93" href="#L93">93</a> <em class="jxr_javadoccomment"> * Test the {@link ContinuousDistribution.Sampler} default stream methods.</em> +<a class="jxr_linenumber" name="L94" href="#L94">94</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L95" href="#L95">95</a> <em class="jxr_javadoccomment"> * @param streamSize Number of values to generate.</em> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L97" href="#L97">97</a> @ParameterizedTest +<a class="jxr_linenumber" name="L98" href="#L98">98</a> @ValueSource(longs = {0, 1, 13}) +<a class="jxr_linenumber" name="L99" href="#L99">99</a> <strong class="jxr_keyword">void</strong> testSamplerStreamMethods(<strong class="jxr_keyword">long</strong> streamSize) { +<a class="jxr_linenumber" name="L100" href="#L100">100</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> seed = ThreadLocalRandom.current().nextDouble(); +<a class="jxr_linenumber" name="L101" href="#L101">101</a> <strong class="jxr_keyword">final</strong> ContinuousDistribution.Sampler s1 = createIncrementSampler(seed); +<a class="jxr_linenumber" name="L102" href="#L102">102</a> <strong class="jxr_keyword">final</strong> ContinuousDistribution.Sampler s2 = createIncrementSampler(seed); +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <strong class="jxr_keyword">final</strong> ContinuousDistribution.Sampler s3 = createIncrementSampler(seed); +<a class="jxr_linenumber" name="L104" href="#L104">104</a> <em class="jxr_comment">// Get the reference output from the sample() method</em> +<a class="jxr_linenumber" name="L105" href="#L105">105</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] x = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[(<strong class="jxr_keyword">int</strong>) streamSize]; +<a class="jxr_linenumber" name="L106" href="#L106">106</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < x.length; i++) { +<a class="jxr_linenumber" name="L107" href="#L107">107</a> x[i] = s1.sample(); +<a class="jxr_linenumber" name="L108" href="#L108">108</a> } +<a class="jxr_linenumber" name="L109" href="#L109">109</a> <em class="jxr_comment">// Test default stream methods</em> +<a class="jxr_linenumber" name="L110" href="#L110">110</a> Assertions.assertArrayEquals(x, s2.samples().limit(streamSize).toArray(), <span class="jxr_string">"samples()"</span>); +<a class="jxr_linenumber" name="L111" href="#L111">111</a> Assertions.assertArrayEquals(x, s3.samples(streamSize).toArray(), <span class="jxr_string">"samples(long)"</span>); +<a class="jxr_linenumber" name="L112" href="#L112">112</a> } +<a class="jxr_linenumber" name="L113" href="#L113">113</a> +<a class="jxr_linenumber" name="L114" href="#L114">114</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L115" href="#L115">115</a> <em class="jxr_javadoccomment"> * Test the {@link ContinuousDistribution.Sampler} default stream method with a bad stream size.</em> +<a class="jxr_linenumber" name="L116" href="#L116">116</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L117" href="#L117">117</a> <em class="jxr_javadoccomment"> * @param streamSize Number of values to generate.</em> +<a class="jxr_linenumber" name="L118" href="#L118">118</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L119" href="#L119">119</a> @ParameterizedTest +<a class="jxr_linenumber" name="L120" href="#L120">120</a> @ValueSource(longs = {-1, -6576237846822L}) +<a class="jxr_linenumber" name="L121" href="#L121">121</a> <strong class="jxr_keyword">void</strong> testSamplerStreamMethodsThrow(<strong class="jxr_keyword">long</strong> streamSize) { +<a class="jxr_linenumber" name="L122" href="#L122">122</a> <strong class="jxr_keyword">final</strong> ContinuousDistribution.Sampler s = createIncrementSampler(42); +<a class="jxr_linenumber" name="L123" href="#L123">123</a> Assertions.assertThrows(IllegalArgumentException.<strong class="jxr_keyword">class</strong>, () -> s.samples(streamSize)); +<a class="jxr_linenumber" name="L124" href="#L124">124</a> } +<a class="jxr_linenumber" name="L125" href="#L125">125</a> +<a class="jxr_linenumber" name="L126" href="#L126">126</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L127" href="#L127">127</a> <em class="jxr_javadoccomment"> * Test the {@link ContinuousDistribution.Sampler} default stream methods are not parallel.</em> +<a class="jxr_linenumber" name="L128" href="#L128">128</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L129" href="#L129">129</a> @Test +<a class="jxr_linenumber" name="L130" href="#L130">130</a> <strong class="jxr_keyword">void</strong> testSamplerStreamMethodsNotParallel() { +<a class="jxr_linenumber" name="L131" href="#L131">131</a> <strong class="jxr_keyword">final</strong> ContinuousDistribution.Sampler s = createIncrementSampler(42); +<a class="jxr_linenumber" name="L132" href="#L132">132</a> Assertions.assertFalse(s.samples().isParallel(), <span class="jxr_string">"samples() should not be parallel"</span>); +<a class="jxr_linenumber" name="L133" href="#L133">133</a> Assertions.assertFalse(s.samples(11).isParallel(), <span class="jxr_string">"samples(long) should not be parallel"</span>); +<a class="jxr_linenumber" name="L134" href="#L134">134</a> } +<a class="jxr_linenumber" name="L135" href="#L135">135</a> +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_javadoccomment"> * Creates the sampler with a given seed value.</em> +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment"> * Each successive output sample will increment this value by 1.</em> +<a class="jxr_linenumber" name="L139" href="#L139">139</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L140" href="#L140">140</a> <em class="jxr_javadoccomment"> * @param seed Seed value.</em> +<a class="jxr_linenumber" name="L141" href="#L141">141</a> <em class="jxr_javadoccomment"> * @return the sampler</em> +<a class="jxr_linenumber" name="L142" href="#L142">142</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L143" href="#L143">143</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> ContinuousDistribution.Sampler createIncrementSampler(<strong class="jxr_keyword">double</strong> seed) { +<a class="jxr_linenumber" name="L144" href="#L144">144</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> ContinuousDistribution.Sampler() { +<a class="jxr_linenumber" name="L145" href="#L145">145</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> x = seed; +<a class="jxr_linenumber" name="L146" href="#L146">146</a> +<a class="jxr_linenumber" name="L147" href="#L147">147</a> @Override +<a class="jxr_linenumber" name="L148" href="#L148">148</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> sample() { +<a class="jxr_linenumber" name="L149" href="#L149">149</a> <strong class="jxr_keyword">return</strong> x += 1; +<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> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html> Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/DiscreteDistributionTest.html ============================================================================== --- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/DiscreteDistributionTest.html (added) +++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/DiscreteDistributionTest.html Thu Dec 1 16:47:12 2022 @@ -0,0 +1,264 @@ +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd";> +<html xmlns="http://www.w3.org/1999/xhtml"; xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>DiscreteDistributionTest xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/DiscreteDistributionTest.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0"; target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.statistics.distribution; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <strong class="jxr_keyword">import</strong> java.util.concurrent.ThreadLocalRandom; +<a class="jxr_linenumber" name="L20" href="#L20">20</a> <strong class="jxr_keyword">import</strong> org.apache.commons.rng.UniformRandomProvider; +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions; +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test; +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest; +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.ValueSource; +<a class="jxr_linenumber" name="L25" href="#L25">25</a> +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <em class="jxr_javadoccomment"> * Test default implementations in the {@link DiscreteDistribution} interface.</em> +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <strong class="jxr_keyword">class</strong> <a name="DiscreteDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/DiscreteDistributionTest.html#DiscreteDistributionTest">DiscreteDistributionTest</a> { +<a class="jxr_linenumber" name="L30" href="#L30">30</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <em class="jxr_javadoccomment"> * Test the default interface methods.</em> +<a class="jxr_linenumber" name="L32" href="#L32">32</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L33" href="#L33">33</a> @Test +<a class="jxr_linenumber" name="L34" href="#L34">34</a> <strong class="jxr_keyword">void</strong> testDefaultMethods() { +<a class="jxr_linenumber" name="L35" href="#L35">35</a> <em class="jxr_comment">// Implement methods called by the defaults</em> +<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution dist = <strong class="jxr_keyword">new</strong> InvalidDiscreteDistribution() { +<a class="jxr_linenumber" name="L37" href="#L37">37</a> @Override +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> probability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L39" href="#L39">39</a> <strong class="jxr_keyword">return</strong> x; +<a class="jxr_linenumber" name="L40" href="#L40">40</a> } +<a class="jxr_linenumber" name="L41" href="#L41">41</a> +<a class="jxr_linenumber" name="L42" href="#L42">42</a> @Override +<a class="jxr_linenumber" name="L43" href="#L43">43</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> cumulativeProbability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L44" href="#L44">44</a> <em class="jxr_comment">// Return some different values to allow the survival probability to be</em> +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <em class="jxr_comment">// tested</em> +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <strong class="jxr_keyword">if</strong> (x < 0) { +<a class="jxr_linenumber" name="L47" href="#L47">47</a> <strong class="jxr_keyword">return</strong> x < -5 ? 0.25 : 0.5; +<a class="jxr_linenumber" name="L48" href="#L48">48</a> } +<a class="jxr_linenumber" name="L49" href="#L49">49</a> <strong class="jxr_keyword">return</strong> x > 5 ? 1.0 : 0.75; +<a class="jxr_linenumber" name="L50" href="#L50">50</a> } +<a class="jxr_linenumber" name="L51" href="#L51">51</a> +<a class="jxr_linenumber" name="L52" href="#L52">52</a> @Override +<a class="jxr_linenumber" name="L53" href="#L53">53</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> inverseCumulativeProbability(<strong class="jxr_keyword">double</strong> p) { +<a class="jxr_linenumber" name="L54" href="#L54">54</a> <em class="jxr_comment">// For the default inverseSurvivalProbability(double) method</em> +<a class="jxr_linenumber" name="L55" href="#L55">55</a> <strong class="jxr_keyword">return</strong> (<strong class="jxr_keyword">int</strong>) (10 * p); +<a class="jxr_linenumber" name="L56" href="#L56">56</a> } +<a class="jxr_linenumber" name="L57" href="#L57">57</a> }; +<a class="jxr_linenumber" name="L58" href="#L58">58</a> +<a class="jxr_linenumber" name="L59" href="#L59">59</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> x : <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] {Integer.MIN_VALUE, -1, 0, 1, 2, Integer.MAX_VALUE}) { +<a class="jxr_linenumber" name="L60" href="#L60">60</a> <em class="jxr_comment">// Return the log of the density</em> +<a class="jxr_linenumber" name="L61" href="#L61">61</a> Assertions.assertEquals(Math.log(x), dist.logProbability(x)); +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <em class="jxr_comment">// Must return 1 - CDF(x)</em> +<a class="jxr_linenumber" name="L63" href="#L63">63</a> Assertions.assertEquals(1.0 - dist.cumulativeProbability(x), dist.survivalProbability(x)); +<a class="jxr_linenumber" name="L64" href="#L64">64</a> } +<a class="jxr_linenumber" name="L65" href="#L65">65</a> +<a class="jxr_linenumber" name="L66" href="#L66">66</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p : <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {0.2, 0.5, 0.7}) { +<a class="jxr_linenumber" name="L67" href="#L67">67</a> Assertions.assertEquals(dist.inverseCumulativeProbability(1 - p), +<a class="jxr_linenumber" name="L68" href="#L68">68</a> dist.inverseSurvivalProbability(p)); +<a class="jxr_linenumber" name="L69" href="#L69">69</a> } +<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">/**</em> +<a class="jxr_linenumber" name="L73" href="#L73">73</a> <em class="jxr_javadoccomment"> * Test the default implementation of probability in a range.</em> +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L75" href="#L75">75</a> @Test +<a class="jxr_linenumber" name="L76" href="#L76">76</a> <strong class="jxr_keyword">void</strong> testDefaultProbabilityRange() { +<a class="jxr_linenumber" name="L77" href="#L77">77</a> <em class="jxr_comment">// Return a marker probability. This should be unique for each input.</em> +<a class="jxr_linenumber" name="L78" href="#L78">78</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution dist = <strong class="jxr_keyword">new</strong> InvalidDiscreteDistribution() { +<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> <strong class="jxr_keyword">double</strong> probability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L81" href="#L81">81</a> <em class="jxr_comment">// Value >= 1</em> +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <strong class="jxr_keyword">return</strong> x + (1L << 31) + 1; +<a class="jxr_linenumber" name="L83" href="#L83">83</a> } +<a class="jxr_linenumber" name="L84" href="#L84">84</a> +<a class="jxr_linenumber" name="L85" href="#L85">85</a> @Override +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> cumulativeProbability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L87" href="#L87">87</a> <em class="jxr_comment">// Value in [-1, 1)</em> +<a class="jxr_linenumber" name="L88" href="#L88">88</a> <strong class="jxr_keyword">return</strong> x * 0x1.0p-31; +<a class="jxr_linenumber" name="L89" href="#L89">89</a> } +<a class="jxr_linenumber" name="L90" href="#L90">90</a> +<a class="jxr_linenumber" name="L91" href="#L91">91</a> @Override +<a class="jxr_linenumber" name="L92" href="#L92">92</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getSupportLowerBound() { +<a class="jxr_linenumber" name="L93" href="#L93">93</a> <strong class="jxr_keyword">return</strong> Integer.MIN_VALUE; +<a class="jxr_linenumber" name="L94" href="#L94">94</a> } +<a class="jxr_linenumber" name="L95" href="#L95">95</a> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> @Override +<a class="jxr_linenumber" name="L97" href="#L97">97</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getSupportUpperBound() { +<a class="jxr_linenumber" name="L98" href="#L98">98</a> <strong class="jxr_keyword">return</strong> Integer.MAX_VALUE; +<a class="jxr_linenumber" name="L99" href="#L99">99</a> } +<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> <em class="jxr_comment">// Test default implementation</em> +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] values = { +<a class="jxr_linenumber" name="L104" href="#L104">104</a> Integer.MIN_VALUE, Integer.MIN_VALUE + 1, +<a class="jxr_linenumber" name="L105" href="#L105">105</a> -3, -2, -1, 0, 1, 2, 3, +<a class="jxr_linenumber" name="L106" href="#L106">106</a> Integer.MAX_VALUE - 1, Integer.MAX_VALUE +<a class="jxr_linenumber" name="L107" href="#L107">107</a> }; +<a class="jxr_linenumber" name="L108" href="#L108">108</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> x0 : values) { +<a class="jxr_linenumber" name="L109" href="#L109">109</a> <em class="jxr_comment">// probability(x, x) == 0.0</em> +<a class="jxr_linenumber" name="L110" href="#L110">110</a> Assertions.assertEquals(0.0, dist.probability(x0, x0)); +<a class="jxr_linenumber" name="L111" href="#L111">111</a> <strong class="jxr_keyword">if</strong> (x0 < dist.getSupportUpperBound()) { +<a class="jxr_linenumber" name="L112" href="#L112">112</a> <em class="jxr_comment">// probability(x, x + 1) == probability(x + 1)</em> +<a class="jxr_linenumber" name="L113" href="#L113">113</a> Assertions.assertEquals(dist.probability(x0 + 1), dist.probability(x0, x0 + 1)); +<a class="jxr_linenumber" name="L114" href="#L114">114</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> x1 : values) { +<a class="jxr_linenumber" name="L115" href="#L115">115</a> <strong class="jxr_keyword">if</strong> (x1 > x0 + 1) { +<a class="jxr_linenumber" name="L116" href="#L116">116</a> <em class="jxr_comment">// probability(x0, x1) == cdf(x1) - cdf(x0)</em> +<a class="jxr_linenumber" name="L117" href="#L117">117</a> Assertions.assertEquals(dist.cumulativeProbability(x1) - dist.cumulativeProbability(x0), +<a class="jxr_linenumber" name="L118" href="#L118">118</a> dist.probability(x0, x1)); +<a class="jxr_linenumber" name="L119" href="#L119">119</a> } <strong class="jxr_keyword">else</strong> <strong class="jxr_keyword">if</strong> (x1 < x0) { +<a class="jxr_linenumber" name="L120" href="#L120">120</a> Assertions.assertThrows(IllegalArgumentException.<strong class="jxr_keyword">class</strong>, () -> dist.probability(x0, x1)); +<a class="jxr_linenumber" name="L121" href="#L121">121</a> } +<a class="jxr_linenumber" name="L122" href="#L122">122</a> } +<a class="jxr_linenumber" name="L123" href="#L123">123</a> } +<a class="jxr_linenumber" name="L124" href="#L124">124</a> } +<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> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L128" href="#L128">128</a> <em class="jxr_javadoccomment"> * Test the default implementation of probability in a range calls the probability function</em> +<a class="jxr_linenumber" name="L129" href="#L129">129</a> <em class="jxr_javadoccomment"> * when x+1 or x-1 would overflow.</em> +<a class="jxr_linenumber" name="L130" href="#L130">130</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L131" href="#L131">131</a> @Test +<a class="jxr_linenumber" name="L132" href="#L132">132</a> <strong class="jxr_keyword">void</strong> testDefaultProbabilityRangeOverflow() { +<a class="jxr_linenumber" name="L133" href="#L133">133</a> <em class="jxr_comment">// Return only the probability.</em> +<a class="jxr_linenumber" name="L134" href="#L134">134</a> <em class="jxr_comment">// The CDF should not be called.</em> +<a class="jxr_linenumber" name="L135" href="#L135">135</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution dist = <strong class="jxr_keyword">new</strong> InvalidDiscreteDistribution() { +<a class="jxr_linenumber" name="L136" href="#L136">136</a> @Override +<a class="jxr_linenumber" name="L137" href="#L137">137</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> probability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <strong class="jxr_keyword">return</strong> x; +<a class="jxr_linenumber" name="L139" href="#L139">139</a> } +<a class="jxr_linenumber" name="L140" href="#L140">140</a> }; +<a class="jxr_linenumber" name="L141" href="#L141">141</a> +<a class="jxr_linenumber" name="L142" href="#L142">142</a> <em class="jxr_comment">// Extreme x at the integer limits</em> +<a class="jxr_linenumber" name="L143" href="#L143">143</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> min = Integer.MIN_VALUE; +<a class="jxr_linenumber" name="L144" href="#L144">144</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> max = Integer.MAX_VALUE; +<a class="jxr_linenumber" name="L145" href="#L145">145</a> Assertions.assertEquals(0.0, dist.probability(min, min)); +<a class="jxr_linenumber" name="L146" href="#L146">146</a> Assertions.assertEquals(0.0, dist.probability(max, max)); +<a class="jxr_linenumber" name="L147" href="#L147">147</a> Assertions.assertEquals(min + 1, dist.probability(min, min + 1)); +<a class="jxr_linenumber" name="L148" href="#L148">148</a> Assertions.assertEquals(max, dist.probability(max - 1, max)); +<a class="jxr_linenumber" name="L149" href="#L149">149</a> } +<a class="jxr_linenumber" name="L150" href="#L150">150</a> +<a class="jxr_linenumber" name="L151" href="#L151">151</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L152" href="#L152">152</a> <em class="jxr_javadoccomment"> * Test the {@link DiscreteDistribution.Sampler} default stream methods.</em> +<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"> * @param streamSize Number of values to generate.</em> +<a class="jxr_linenumber" name="L155" href="#L155">155</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L156" href="#L156">156</a> @ParameterizedTest +<a class="jxr_linenumber" name="L157" href="#L157">157</a> @ValueSource(longs = {0, 1, 13}) +<a class="jxr_linenumber" name="L158" href="#L158">158</a> <strong class="jxr_keyword">void</strong> testSamplerStreamMethods(<strong class="jxr_keyword">long</strong> streamSize) { +<a class="jxr_linenumber" name="L159" href="#L159">159</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> seed = ThreadLocalRandom.current().nextInt(); +<a class="jxr_linenumber" name="L160" href="#L160">160</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler s1 = createIncrementSampler(seed); +<a class="jxr_linenumber" name="L161" href="#L161">161</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler s2 = createIncrementSampler(seed); +<a class="jxr_linenumber" name="L162" href="#L162">162</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler s3 = createIncrementSampler(seed); +<a class="jxr_linenumber" name="L163" href="#L163">163</a> <em class="jxr_comment">// Get the reference output from the sample() method</em> +<a class="jxr_linenumber" name="L164" href="#L164">164</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] x = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[(<strong class="jxr_keyword">int</strong>) streamSize]; +<a class="jxr_linenumber" name="L165" href="#L165">165</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i < x.length; i++) { +<a class="jxr_linenumber" name="L166" href="#L166">166</a> x[i] = s1.sample(); +<a class="jxr_linenumber" name="L167" href="#L167">167</a> } +<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_comment">// Test default stream methods</em> +<a class="jxr_linenumber" name="L169" href="#L169">169</a> Assertions.assertArrayEquals(x, s2.samples().limit(streamSize).toArray(), <span class="jxr_string">"samples()"</span>); +<a class="jxr_linenumber" name="L170" href="#L170">170</a> Assertions.assertArrayEquals(x, s3.samples(streamSize).toArray(), <span class="jxr_string">"samples(long)"</span>); +<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"> * Test the {@link DiscreteDistribution.Sampler} default stream method with a bad stream size.</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 streamSize Number of values to generate.</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> @ParameterizedTest +<a class="jxr_linenumber" name="L179" href="#L179">179</a> @ValueSource(longs = {-1, -6576237846822L}) +<a class="jxr_linenumber" name="L180" href="#L180">180</a> <strong class="jxr_keyword">void</strong> testSamplerStreamMethodsThrow(<strong class="jxr_keyword">long</strong> streamSize) { +<a class="jxr_linenumber" name="L181" href="#L181">181</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler s = createIncrementSampler(42); +<a class="jxr_linenumber" name="L182" href="#L182">182</a> Assertions.assertThrows(IllegalArgumentException.<strong class="jxr_keyword">class</strong>, () -> s.samples(streamSize)); +<a class="jxr_linenumber" name="L183" href="#L183">183</a> } +<a class="jxr_linenumber" name="L184" href="#L184">184</a> +<a class="jxr_linenumber" name="L185" href="#L185">185</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment"> * Test the {@link DiscreteDistribution.Sampler} default stream methods are not parallel.</em> +<a class="jxr_linenumber" name="L187" href="#L187">187</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L188" href="#L188">188</a> @Test +<a class="jxr_linenumber" name="L189" href="#L189">189</a> <strong class="jxr_keyword">void</strong> testSamplerStreamMethodsNotParallel() { +<a class="jxr_linenumber" name="L190" href="#L190">190</a> <strong class="jxr_keyword">final</strong> DiscreteDistribution.Sampler s = createIncrementSampler(42); +<a class="jxr_linenumber" name="L191" href="#L191">191</a> Assertions.assertFalse(s.samples().isParallel(), <span class="jxr_string">"samples() should not be parallel"</span>); +<a class="jxr_linenumber" name="L192" href="#L192">192</a> Assertions.assertFalse(s.samples(11).isParallel(), <span class="jxr_string">"samples(long) should not be parallel"</span>); +<a class="jxr_linenumber" name="L193" href="#L193">193</a> } +<a class="jxr_linenumber" name="L194" href="#L194">194</a> +<a class="jxr_linenumber" name="L195" href="#L195">195</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L196" href="#L196">196</a> <em class="jxr_javadoccomment"> * Creates the sampler with a given seed value.</em> +<a class="jxr_linenumber" name="L197" href="#L197">197</a> <em class="jxr_javadoccomment"> * Each successive output sample will increment this value by 1.</em> +<a class="jxr_linenumber" name="L198" href="#L198">198</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L199" href="#L199">199</a> <em class="jxr_javadoccomment"> * @param seed Seed value.</em> +<a class="jxr_linenumber" name="L200" href="#L200">200</a> <em class="jxr_javadoccomment"> * @return the sampler</em> +<a class="jxr_linenumber" name="L201" href="#L201">201</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L202" href="#L202">202</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> DiscreteDistribution.Sampler createIncrementSampler(<strong class="jxr_keyword">int</strong> seed) { +<a class="jxr_linenumber" name="L203" href="#L203">203</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> DiscreteDistribution.Sampler() { +<a class="jxr_linenumber" name="L204" href="#L204">204</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> x = seed; +<a class="jxr_linenumber" name="L205" href="#L205">205</a> +<a class="jxr_linenumber" name="L206" href="#L206">206</a> @Override +<a class="jxr_linenumber" name="L207" href="#L207">207</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> sample() { +<a class="jxr_linenumber" name="L208" href="#L208">208</a> <strong class="jxr_keyword">return</strong> x += 1; +<a class="jxr_linenumber" name="L209" href="#L209">209</a> } +<a class="jxr_linenumber" name="L210" href="#L210">210</a> }; +<a class="jxr_linenumber" name="L211" href="#L211">211</a> } +<a class="jxr_linenumber" name="L212" href="#L212">212</a> +<a class="jxr_linenumber" name="L213" href="#L213">213</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L214" href="#L214">214</a> <em class="jxr_javadoccomment"> * Invalid implementation of DiscreteDistribution that raise an exception for all methods.</em> +<a class="jxr_linenumber" name="L215" href="#L215">215</a> <em class="jxr_javadoccomment"> * Ensures the methods are not called.</em> +<a class="jxr_linenumber" name="L216" href="#L216">216</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L217" href="#L217">217</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">abstract</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">class</strong> InvalidDiscreteDistribution <strong class="jxr_keyword">implements</strong> DiscreteDistribution { +<a class="jxr_linenumber" name="L218" href="#L218">218</a> @Override +<a class="jxr_linenumber" name="L219" href="#L219">219</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> probability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L220" href="#L220">220</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L221" href="#L221">221</a> } +<a class="jxr_linenumber" name="L222" href="#L222">222</a> @Override +<a class="jxr_linenumber" name="L223" href="#L223">223</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> cumulativeProbability(<strong class="jxr_keyword">int</strong> x) { +<a class="jxr_linenumber" name="L224" href="#L224">224</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L225" href="#L225">225</a> } +<a class="jxr_linenumber" name="L226" href="#L226">226</a> @Override +<a class="jxr_linenumber" name="L227" href="#L227">227</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> inverseCumulativeProbability(<strong class="jxr_keyword">double</strong> p) { +<a class="jxr_linenumber" name="L228" href="#L228">228</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L229" href="#L229">229</a> } +<a class="jxr_linenumber" name="L230" href="#L230">230</a> @Override +<a class="jxr_linenumber" name="L231" href="#L231">231</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getMean() { +<a class="jxr_linenumber" name="L232" href="#L232">232</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L233" href="#L233">233</a> } +<a class="jxr_linenumber" name="L234" href="#L234">234</a> @Override +<a class="jxr_linenumber" name="L235" href="#L235">235</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">double</strong> getVariance() { +<a class="jxr_linenumber" name="L236" href="#L236">236</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L237" href="#L237">237</a> } +<a class="jxr_linenumber" name="L238" href="#L238">238</a> @Override +<a class="jxr_linenumber" name="L239" href="#L239">239</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getSupportLowerBound() { +<a class="jxr_linenumber" name="L240" href="#L240">240</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L241" href="#L241">241</a> } +<a class="jxr_linenumber" name="L242" href="#L242">242</a> @Override +<a class="jxr_linenumber" name="L243" href="#L243">243</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getSupportUpperBound() { +<a class="jxr_linenumber" name="L244" href="#L244">244</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L245" href="#L245">245</a> } +<a class="jxr_linenumber" name="L246" href="#L246">246</a> @Override +<a class="jxr_linenumber" name="L247" href="#L247">247</a> <strong class="jxr_keyword">public</strong> Sampler createSampler(UniformRandomProvider rng) { +<a class="jxr_linenumber" name="L248" href="#L248">248</a> <strong class="jxr_keyword">throw</strong> <strong class="jxr_keyword">new</strong> AssertionError(); +<a class="jxr_linenumber" name="L249" href="#L249">249</a> } +<a class="jxr_linenumber" name="L250" href="#L250">250</a> } +<a class="jxr_linenumber" name="L251" href="#L251">251</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2018–2022 <a href="https://www.apache.org/";>The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html>
