You are viewing a plain text version of this content. The canonical link for it is here.
Posted to commits@commons.apache.org by ah...@apache.org on 2022/12/01 16:47:13 UTC
svn commit: r58392 [28/50] - in /dev/commons/statistics/1.0-RC1: ./ binaries/ site/ site/css/ site/images/ site/images/logos/ site/img/ site/js/ site/release-notes/ site/style/ site/userguide/ site/xref-test/ site/xref-test/org/ site/xref-test/org/apac...
Added: dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/TrapezoidalDistributionTest.html
==============================================================================
--- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/TrapezoidalDistributionTest.html (added)
+++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/TrapezoidalDistributionTest.html Thu Dec 1 16:47:12 2022
@@ -0,0 +1,202 @@
+<!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>TrapezoidalDistributionTest xref</title>
+<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" />
+</head>
+<body>
+<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/TrapezoidalDistributionTest.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.math.BigDecimal;
+<a class="jxr_linenumber" name="L21" href="#L21">21</a> <strong class="jxr_keyword">import</strong> java.math.MathContext;
+<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> java.util.stream.Stream;
+<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Assertions;
+<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.api.Test;
+<a class="jxr_linenumber" name="L25" href="#L25">25</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest;
+<a class="jxr_linenumber" name="L26" href="#L26">26</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.Arguments;
+<a class="jxr_linenumber" name="L27" href="#L27">27</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.MethodSource;
+<a class="jxr_linenumber" name="L28" href="#L28">28</a>
+<a class="jxr_linenumber" name="L29" href="#L29">29</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L30" href="#L30">30</a> <em class="jxr_javadoccomment"> * Test cases for {@link TrapezoidalDistribution}.</em>
+<a class="jxr_linenumber" name="L31" href="#L31">31</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</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> <strong class="jxr_keyword">class</strong> <a name="TrapezoidalDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/TrapezoidalDistributionTest.html#TrapezoidalDistributionTest">TrapezoidalDistributionTest</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="L34" href="#L34">34</a> @Override
+<a class="jxr_linenumber" name="L35" href="#L35">35</a> ContinuousDistribution makeDistribution(Object... parameters) {
+<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a = (Double) parameters[0];
+<a class="jxr_linenumber" name="L37" href="#L37">37</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> b = (Double) parameters[1];
+<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> c = (Double) parameters[2];
+<a class="jxr_linenumber" name="L39" href="#L39">39</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> d = (Double) parameters[3];
+<a class="jxr_linenumber" name="L40" href="#L40">40</a> <strong class="jxr_keyword">return</strong> TrapezoidalDistribution.of(a, b, c, d);
+<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> @Override
+<a class="jxr_linenumber" name="L44" href="#L44">44</a> Object[][] makeInvalidParameters() {
+<a class="jxr_linenumber" name="L45" href="#L45">45</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] {
+<a class="jxr_linenumber" name="L46" href="#L46">46</a> {0.0, 0.0, 0.0, 0.0},
+<a class="jxr_linenumber" name="L47" href="#L47">47</a> <em class="jxr_comment">// 1.0, 2.0, 3.0, 4.0 is OK - move points to incorrect locations</em>
+<a class="jxr_linenumber" name="L48" href="#L48">48</a> {5.0, 2.0, 3.0, 4.0}, <em class="jxr_comment">// a > d</em>
+<a class="jxr_linenumber" name="L49" href="#L49">49</a> {1.0, 5.0, 3.0, 4.0}, <em class="jxr_comment">// b > d</em>
+<a class="jxr_linenumber" name="L50" href="#L50">50</a> {1.0, 2.0, 5.0, 4.0}, <em class="jxr_comment">// c > d</em>
+<a class="jxr_linenumber" name="L51" href="#L51">51</a> {3.5, 2.0, 3.0, 4.0}, <em class="jxr_comment">// a > c</em>
+<a class="jxr_linenumber" name="L52" href="#L52">52</a> {1.0, 3.5, 3.0, 4.0}, <em class="jxr_comment">// b > c</em>
+<a class="jxr_linenumber" name="L53" href="#L53">53</a> {2.5, 2.0, 3.0, 4.0}, <em class="jxr_comment">// a > b</em>
+<a class="jxr_linenumber" name="L54" href="#L54">54</a> {1.0, 2.0, 3.0, 0.0}, <em class="jxr_comment">// d < a</em>
+<a class="jxr_linenumber" name="L55" href="#L55">55</a> {1.0, 2.0, 0.0, 4.0}, <em class="jxr_comment">// c < a</em>
+<a class="jxr_linenumber" name="L56" href="#L56">56</a> {1.0, 0.0, 3.0, 4.0}, <em class="jxr_comment">// b < a</em>
+<a class="jxr_linenumber" name="L57" href="#L57">57</a> {1.0, 2.0, 3.0, 1.5}, <em class="jxr_comment">// d < b</em>
+<a class="jxr_linenumber" name="L58" href="#L58">58</a> {1.0, 2.0, 1.5, 4.0}, <em class="jxr_comment">// c < b</em>
+<a class="jxr_linenumber" name="L59" href="#L59">59</a> {1.0, 2.0, 3.0, 2.5}, <em class="jxr_comment">// d < c</em>
+<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>
+<a class="jxr_linenumber" name="L63" href="#L63">63</a> @Override
+<a class="jxr_linenumber" name="L64" href="#L64">64</a> String[] getParameterNames() {
+<a class="jxr_linenumber" name="L65" href="#L65">65</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"SupportLowerBound"</span>, <span class="jxr_string">"B"</span>, <span class="jxr_string">"C"</span>, <span class="jxr_string">"SupportUpperBound"</span>};
+<a class="jxr_linenumber" name="L66" href="#L66">66</a> }
+<a class="jxr_linenumber" name="L67" href="#L67">67</a>
+<a class="jxr_linenumber" name="L68" href="#L68">68</a> @Override
+<a class="jxr_linenumber" name="L69" href="#L69">69</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() {
+<a class="jxr_linenumber" name="L70" href="#L70">70</a> <em class="jxr_comment">// Tolerance is 4.440892098500626E-15.</em>
+<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">return</strong> 20 * RELATIVE_EPS;
+<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_comment">//-------------------- Additional test cases -------------------------------</em>
+<a class="jxr_linenumber" name="L75" href="#L75">75</a>
+<a class="jxr_linenumber" name="L76" href="#L76">76</a> @ParameterizedTest
+<a class="jxr_linenumber" name="L77" href="#L77">77</a> @MethodSource
+<a class="jxr_linenumber" name="L78" href="#L78">78</a> <strong class="jxr_keyword">void</strong> testAdditionalMoments(<strong class="jxr_keyword">double</strong> a, <strong class="jxr_keyword">double</strong> b, <strong class="jxr_keyword">double</strong> c, <strong class="jxr_keyword">double</strong> d, <strong class="jxr_keyword">double</strong> mean, <strong class="jxr_keyword">double</strong> variance, <strong class="jxr_keyword">int</strong> ulps) {
+<a class="jxr_linenumber" name="L79" href="#L79">79</a> <strong class="jxr_keyword">final</strong> TrapezoidalDistribution dist = TrapezoidalDistribution.of(a, b, c, d);
+<a class="jxr_linenumber" name="L80" href="#L80">80</a> testMoments(dist, mean, variance, DoubleTolerances.ulps(ulps));
+<a class="jxr_linenumber" name="L81" href="#L81">81</a> }
+<a class="jxr_linenumber" name="L82" href="#L82">82</a>
+<a class="jxr_linenumber" name="L83" href="#L83">83</a> <strong class="jxr_keyword">static</strong> Stream<Arguments> testAdditionalMoments() {
+<a class="jxr_linenumber" name="L84" href="#L84">84</a> <strong class="jxr_keyword">return</strong> Stream.of(
+<a class="jxr_linenumber" name="L85" href="#L85">85</a> <em class="jxr_comment">// Computed using scipy.stats.trapezoid</em>
+<a class="jxr_linenumber" name="L86" href="#L86">86</a> <em class="jxr_comment">// Up slope, then flat</em>
+<a class="jxr_linenumber" name="L87" href="#L87">87</a> Arguments.of(0, 0.1, 1, 1, 0.5245614035087719, 0.07562480763311791, 8),
+<a class="jxr_linenumber" name="L88" href="#L88">88</a> Arguments.of(0, 1e-3, 1, 1, 0.5002499583124894, 0.08325006249999839, 8),
+<a class="jxr_linenumber" name="L89" href="#L89">89</a> Arguments.of(0, 1e-6, 1, 1, 0.5000002499999582, 0.08333325000006259, 8),
+<a class="jxr_linenumber" name="L90" href="#L90">90</a> Arguments.of(0, 1e-9, 1, 1, 0.50000000025, 0.08333333324999997, 4),
+<a class="jxr_linenumber" name="L91" href="#L91">91</a> Arguments.of(0, 1e-12, 1, 1, 0.50000000000025, 0.08333333333324999, 4),
+<a class="jxr_linenumber" name="L92" href="#L92">92</a> Arguments.of(0, 1e-15, 1, 1, 0.5000000000000003, 0.0833333333333332, 4),
+<a class="jxr_linenumber" name="L93" href="#L93">93</a> Arguments.of(0, 0, 1, 1, 0.5, 0.08333333333333331, 1),
+<a class="jxr_linenumber" name="L94" href="#L94">94</a> <em class="jxr_comment">// Flat, then down slope</em>
+<a class="jxr_linenumber" name="L95" href="#L95">95</a> Arguments.of(0, 0, 0.9, 1, 0.47543859649122816, 0.07562480763311777, 4),
+<a class="jxr_linenumber" name="L96" href="#L96">96</a> Arguments.of(0, 0, 0.999, 1, 0.49975004168751025, 0.08325006249999842, 4),
+<a class="jxr_linenumber" name="L97" href="#L97">97</a> Arguments.of(0, 0, 0.999999, 1, 0.4999997500000417, 0.08333325000006248, 2),
+<a class="jxr_linenumber" name="L98" href="#L98">98</a> Arguments.of(0, 0, 0.999999999, 1, 0.49999999975000003, 0.08333333325000003, 1),
+<a class="jxr_linenumber" name="L99" href="#L99">99</a> Arguments.of(0, 0, 0.999999999999, 1, 0.49999999999975003, 0.08333333333324999, 1),
+<a class="jxr_linenumber" name="L100" href="#L100">100</a> Arguments.of(0, 0, 0.999999999999999, 1, 0.4999999999999998, 0.08333333333333326, 1),
+<a class="jxr_linenumber" name="L101" href="#L101">101</a> Arguments.of(0, 0, 1, 1, 0.5, 0.08333333333333331, 1),
+<a class="jxr_linenumber" name="L102" href="#L102">102</a> <em class="jxr_comment">// Computed</em>
+<a class="jxr_linenumber" name="L103" href="#L103">103</a> moments(-3, 1, 4, 12, 10),
+<a class="jxr_linenumber" name="L104" href="#L104">104</a> moments(1, 2, 3, 4, 4),
+<a class="jxr_linenumber" name="L105" href="#L105">105</a> moments(11, 15, 22, 30, 4),
+<a class="jxr_linenumber" name="L106" href="#L106">106</a> moments(0, 1, 9, 10, 4),
+<a class="jxr_linenumber" name="L107" href="#L107">107</a> moments(-12, -10, -2, -1, 8)
+<a class="jxr_linenumber" name="L108" href="#L108">108</a> );
+<a class="jxr_linenumber" name="L109" href="#L109">109</a> }
+<a class="jxr_linenumber" name="L110" href="#L110">110</a>
+<a class="jxr_linenumber" name="L111" href="#L111">111</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L112" href="#L112">112</a> <em class="jxr_javadoccomment"> * Compute the mean and variance and return the arguments: [a, b, c, d, mean, variance, ulps].</em>
+<a class="jxr_linenumber" name="L113" href="#L113">113</a> <em class="jxr_javadoccomment"> * The ulps are the expected units of least precision error for the test. This is typically</em>
+<a class="jxr_linenumber" name="L114" href="#L114">114</a> <em class="jxr_javadoccomment"> * limited by the variance as it depends on accuracy of the mean squared.</em>
+<a class="jxr_linenumber" name="L115" href="#L115">115</a> <em class="jxr_javadoccomment"> */</em>
+<a class="jxr_linenumber" name="L116" href="#L116">116</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> Arguments moments(<strong class="jxr_keyword">double</strong> a, <strong class="jxr_keyword">double</strong> b, <strong class="jxr_keyword">double</strong> c, <strong class="jxr_keyword">double</strong> d, <strong class="jxr_keyword">int</strong> ulps) {
+<a class="jxr_linenumber" name="L117" href="#L117">117</a> <strong class="jxr_keyword">final</strong> BigDecimal aa = <strong class="jxr_keyword">new</strong> BigDecimal(a);
+<a class="jxr_linenumber" name="L118" href="#L118">118</a> <strong class="jxr_keyword">final</strong> BigDecimal bb = <strong class="jxr_keyword">new</strong> BigDecimal(b);
+<a class="jxr_linenumber" name="L119" href="#L119">119</a> <strong class="jxr_keyword">final</strong> BigDecimal cc = <strong class="jxr_keyword">new</strong> BigDecimal(c);
+<a class="jxr_linenumber" name="L120" href="#L120">120</a> <strong class="jxr_keyword">final</strong> BigDecimal dd = <strong class="jxr_keyword">new</strong> BigDecimal(d);
+<a class="jxr_linenumber" name="L121" href="#L121">121</a> <strong class="jxr_keyword">final</strong> MathContext mc = MathContext.DECIMAL128;
+<a class="jxr_linenumber" name="L122" href="#L122">122</a> <em class="jxr_comment">// Mean</em>
+<a class="jxr_linenumber" name="L123" href="#L123">123</a> BigDecimal divisor = dd.add(cc).subtract(aa).subtract(bb).multiply(BigDecimal.valueOf(3));
+<a class="jxr_linenumber" name="L124" href="#L124">124</a> BigDecimal dc = dd.pow(3).subtract(cc.pow(3)).divide(dd.subtract(cc), mc);
+<a class="jxr_linenumber" name="L125" href="#L125">125</a> BigDecimal ba = bb.pow(3).subtract(aa.pow(3)).divide(bb.subtract(aa), mc);
+<a class="jxr_linenumber" name="L126" href="#L126">126</a> <strong class="jxr_keyword">final</strong> BigDecimal mu = dc.subtract(ba).divide(divisor, mc);
+<a class="jxr_linenumber" name="L127" href="#L127">127</a> <em class="jxr_comment">// Variance</em>
+<a class="jxr_linenumber" name="L128" href="#L128">128</a> divisor = divisor.multiply(BigDecimal.valueOf(2));
+<a class="jxr_linenumber" name="L129" href="#L129">129</a> dc = dd.pow(4).subtract(cc.pow(4)).divide(dd.subtract(cc), mc);
+<a class="jxr_linenumber" name="L130" href="#L130">130</a> ba = bb.pow(4).subtract(aa.pow(4)).divide(bb.subtract(aa), mc);
+<a class="jxr_linenumber" name="L131" href="#L131">131</a> <strong class="jxr_keyword">final</strong> BigDecimal <strong class="jxr_keyword">var</strong> = dc.subtract(ba).divide(divisor, mc).subtract(mu.pow(2));
+<a class="jxr_linenumber" name="L132" href="#L132">132</a> <strong class="jxr_keyword">return</strong> Arguments.of(a, b, c, d, mu.doubleValue(), <strong class="jxr_keyword">var</strong>.doubleValue(), ulps);
+<a class="jxr_linenumber" name="L133" href="#L133">133</a> }
+<a class="jxr_linenumber" name="L134" href="#L134">134</a>
+<a class="jxr_linenumber" name="L135" href="#L135">135</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L136" href="#L136">136</a> <em class="jxr_javadoccomment"> * Create a trapezoid with a very long upper tail to explicitly test the survival</em>
+<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_javadoccomment"> * probability is high precision.</em>
+<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment"> */</em>
+<a class="jxr_linenumber" name="L139" href="#L139">139</a> @Test
+<a class="jxr_linenumber" name="L140" href="#L140">140</a> <strong class="jxr_keyword">void</strong> testAdditionalSurvivalProbabilityHighPrecision() {
+<a class="jxr_linenumber" name="L141" href="#L141">141</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> a = 0;
+<a class="jxr_linenumber" name="L142" href="#L142">142</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> b = 0;
+<a class="jxr_linenumber" name="L143" href="#L143">143</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> c = 1;
+<a class="jxr_linenumber" name="L144" href="#L144">144</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> d = 1 << 14;
+<a class="jxr_linenumber" name="L145" href="#L145">145</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x1 = Math.nextDown(d);
+<a class="jxr_linenumber" name="L146" href="#L146">146</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> x2 = Math.nextDown(x1);
+<a class="jxr_linenumber" name="L147" href="#L147">147</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p1 = survivalProbability(a, b, c, d, x1);
+<a class="jxr_linenumber" name="L148" href="#L148">148</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> p2 = survivalProbability(a, b, c, d, x2);
+<a class="jxr_linenumber" name="L149" href="#L149">149</a> <strong class="jxr_keyword">final</strong> TrapezoidalDistribution dist = TrapezoidalDistribution.of(a, b, c, d);
+<a class="jxr_linenumber" name="L150" href="#L150">150</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] points = {x1, x2};
+<a class="jxr_linenumber" name="L151" href="#L151">151</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong>[] probabilities = {p1, p2};
+<a class="jxr_linenumber" name="L152" href="#L152">152</a>
+<a class="jxr_linenumber" name="L153" href="#L153">153</a> <em class="jxr_comment">// This fails if the sf(x) = 1 - cdf(x)</em>
+<a class="jxr_linenumber" name="L154" href="#L154">154</a> testSurvivalProbabilityHighPrecision(
+<a class="jxr_linenumber" name="L155" href="#L155">155</a> dist,
+<a class="jxr_linenumber" name="L156" href="#L156">156</a> points,
+<a class="jxr_linenumber" name="L157" href="#L157">157</a> probabilities,
+<a class="jxr_linenumber" name="L158" href="#L158">158</a> DoubleTolerances.relative(1e-15));
+<a class="jxr_linenumber" name="L159" href="#L159">159</a>
+<a class="jxr_linenumber" name="L160" href="#L160">160</a> <em class="jxr_comment">// This fails if the isf(p) = icdf(1 - p)</em>
+<a class="jxr_linenumber" name="L161" href="#L161">161</a> testInverseSurvivalProbability(
+<a class="jxr_linenumber" name="L162" href="#L162">162</a> dist,
+<a class="jxr_linenumber" name="L163" href="#L163">163</a> probabilities,
+<a class="jxr_linenumber" name="L164" href="#L164">164</a> points,
+<a class="jxr_linenumber" name="L165" href="#L165">165</a> DoubleTolerances.ulps(0));
+<a class="jxr_linenumber" name="L166" href="#L166">166</a> }
+<a class="jxr_linenumber" name="L167" href="#L167">167</a>
+<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L169" href="#L169">169</a> <em class="jxr_javadoccomment"> * Compute the trapezoid distribution survival probability for the value {@code x}</em>
+<a class="jxr_linenumber" name="L170" href="#L170">170</a> <em class="jxr_javadoccomment"> * in the region {@code [c, d]}.</em>
+<a class="jxr_linenumber" name="L171" href="#L171">171</a> <em class="jxr_javadoccomment"> *</em>
+<a class="jxr_linenumber" name="L172" href="#L172">172</a> <em class="jxr_javadoccomment"> * @param a Lower limit of the distribution (inclusive).</em>
+<a class="jxr_linenumber" name="L173" href="#L173">173</a> <em class="jxr_javadoccomment"> * @param b Start of the trapezoid constant density.</em>
+<a class="jxr_linenumber" name="L174" href="#L174">174</a> <em class="jxr_javadoccomment"> * @param c End of the trapezoid constant density.</em>
+<a class="jxr_linenumber" name="L175" href="#L175">175</a> <em class="jxr_javadoccomment"> * @param d Upper limit of the distribution (inclusive).</em>
+<a class="jxr_linenumber" name="L176" href="#L176">176</a> <em class="jxr_javadoccomment"> * @param x Value in [c, d].</em>
+<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment"> * @return the probability</em>
+<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_javadoccomment"> */</em>
+<a class="jxr_linenumber" name="L179" href="#L179">179</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">double</strong> survivalProbability(<strong class="jxr_keyword">double</strong> a, <strong class="jxr_keyword">double</strong> b, <strong class="jxr_keyword">double</strong> c, <strong class="jxr_keyword">double</strong> d, <strong class="jxr_keyword">double</strong> x) {
+<a class="jxr_linenumber" name="L180" href="#L180">180</a> Assertions.assertTrue(c <= x && x <= d, <span class="jxr_string">"Domain error"</span>);
+<a class="jxr_linenumber" name="L181" href="#L181">181</a> <strong class="jxr_keyword">final</strong> BigDecimal aa = <strong class="jxr_keyword">new</strong> BigDecimal(a);
+<a class="jxr_linenumber" name="L182" href="#L182">182</a> <strong class="jxr_keyword">final</strong> BigDecimal bb = <strong class="jxr_keyword">new</strong> BigDecimal(b);
+<a class="jxr_linenumber" name="L183" href="#L183">183</a> <strong class="jxr_keyword">final</strong> BigDecimal cc = <strong class="jxr_keyword">new</strong> BigDecimal(c);
+<a class="jxr_linenumber" name="L184" href="#L184">184</a> <strong class="jxr_keyword">final</strong> BigDecimal dd = <strong class="jxr_keyword">new</strong> BigDecimal(d);
+<a class="jxr_linenumber" name="L185" href="#L185">185</a> <strong class="jxr_keyword">final</strong> BigDecimal divisor = dd.add(cc).subtract(aa).subtract(bb).multiply(dd.subtract(cc));
+<a class="jxr_linenumber" name="L186" href="#L186">186</a> <strong class="jxr_keyword">return</strong> dd.subtract(<strong class="jxr_keyword">new</strong> BigDecimal(x)).pow(2)
+<a class="jxr_linenumber" name="L187" href="#L187">187</a> .divide(divisor, MathContext.DECIMAL128).doubleValue();
+<a class="jxr_linenumber" name="L188" href="#L188">188</a> }
+<a class="jxr_linenumber" name="L189" href="#L189">189</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/TriangularDistributionTest.html
==============================================================================
--- dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/TriangularDistributionTest.html (added)
+++ dev/commons/statistics/1.0-RC1/site/xref-test/org/apache/commons/statistics/distribution/TriangularDistributionTest.html Thu Dec 1 16:47:12 2022
@@ -0,0 +1,108 @@
+<!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>TriangularDistributionTest xref</title>
+<link type="text/css" rel="stylesheet" href="../../../../../stylesheet.css" />
+</head>
+<body>
+<div id="overview"><a href="../../../../../../testapidocs/org/apache/commons/statistics/distribution/TriangularDistributionTest.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.api.Assertions;
+<a class="jxr_linenumber" name="L22" href="#L22">22</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.ParameterizedTest;
+<a class="jxr_linenumber" name="L23" href="#L23">23</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.Arguments;
+<a class="jxr_linenumber" name="L24" href="#L24">24</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.CsvSource;
+<a class="jxr_linenumber" name="L25" href="#L25">25</a> <strong class="jxr_keyword">import</strong> org.junit.jupiter.params.provider.MethodSource;
+<a class="jxr_linenumber" name="L26" href="#L26">26</a>
+<a class="jxr_linenumber" name="L27" href="#L27">27</a> <em class="jxr_javadoccomment">/**</em>
+<a class="jxr_linenumber" name="L28" href="#L28">28</a> <em class="jxr_javadoccomment"> * Test cases for {@link TriangularDistribution}.</em>
+<a class="jxr_linenumber" name="L29" href="#L29">29</a> <em class="jxr_javadoccomment"> * Extends {@link BaseContinuousDistributionTest}. See javadoc of that class for details.</em>
+<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> <strong class="jxr_keyword">class</strong> <a name="TriangularDistributionTest" href="../../../../../org/apache/commons/statistics/distribution/TriangularDistributionTest.html#TriangularDistributionTest">TriangularDistributionTest</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="L32" href="#L32">32</a> @Override
+<a class="jxr_linenumber" name="L33" href="#L33">33</a> ContinuousDistribution makeDistribution(Object... parameters) {
+<a class="jxr_linenumber" name="L34" href="#L34">34</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> lower = (Double) parameters[0];
+<a class="jxr_linenumber" name="L35" href="#L35">35</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> mode = (Double) parameters[1];
+<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">double</strong> upper = (Double) parameters[2];
+<a class="jxr_linenumber" name="L37" href="#L37">37</a> <strong class="jxr_keyword">return</strong> TriangularDistribution.of(lower, mode, upper);
+<a class="jxr_linenumber" name="L38" href="#L38">38</a> }
+<a class="jxr_linenumber" name="L39" href="#L39">39</a>
+<a class="jxr_linenumber" name="L40" href="#L40">40</a>
+<a class="jxr_linenumber" name="L41" href="#L41">41</a> @Override
+<a class="jxr_linenumber" name="L42" href="#L42">42</a> Object[][] makeInvalidParameters() {
+<a class="jxr_linenumber" name="L43" href="#L43">43</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> Object[][] {
+<a class="jxr_linenumber" name="L44" href="#L44">44</a> {0.0, 0.0, 0.0},
+<a class="jxr_linenumber" name="L45" href="#L45">45</a> <em class="jxr_comment">// 1.0, 2.0, 3 is OK - move points to incorrect locations</em>
+<a class="jxr_linenumber" name="L46" href="#L46">46</a> {4.0, 2.0, 3.0},
+<a class="jxr_linenumber" name="L47" href="#L47">47</a> {3.0, 2.0, 3.0},
+<a class="jxr_linenumber" name="L48" href="#L48">48</a> {2.5, 2.0, 3.0},
+<a class="jxr_linenumber" name="L49" href="#L49">49</a> {1.0, 0.0, 3.0},
+<a class="jxr_linenumber" name="L50" href="#L50">50</a> {1.0, 4.0, 3.0},
+<a class="jxr_linenumber" name="L51" href="#L51">51</a> {1.0, 2.0, -1.0},
+<a class="jxr_linenumber" name="L52" href="#L52">52</a> {1.0, 2.0, 1.5},
+<a class="jxr_linenumber" name="L53" href="#L53">53</a> };
+<a class="jxr_linenumber" name="L54" href="#L54">54</a> }
+<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> String[] getParameterNames() {
+<a class="jxr_linenumber" name="L58" href="#L58">58</a> <strong class="jxr_keyword">return</strong> <strong class="jxr_keyword">new</strong> String[] {<span class="jxr_string">"SupportLowerBound"</span>, <span class="jxr_string">"Mode"</span>, <span class="jxr_string">"SupportUpperBound"</span>};
+<a class="jxr_linenumber" name="L59" href="#L59">59</a> }
+<a class="jxr_linenumber" name="L60" href="#L60">60</a>
+<a class="jxr_linenumber" name="L61" href="#L61">61</a> @Override
+<a class="jxr_linenumber" name="L62" href="#L62">62</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">double</strong> getRelativeTolerance() {
+<a class="jxr_linenumber" name="L63" href="#L63">63</a> <em class="jxr_comment">// Tolerance is 4.440892098500626E-15.</em>
+<a class="jxr_linenumber" name="L64" href="#L64">64</a> <strong class="jxr_keyword">return</strong> 20 * RELATIVE_EPS;
+<a class="jxr_linenumber" name="L65" href="#L65">65</a> }
+<a class="jxr_linenumber" name="L66" href="#L66">66</a>
+<a class="jxr_linenumber" name="L67" href="#L67">67</a> <em class="jxr_comment">//-------------------- Additional test cases -------------------------------</em>
+<a class="jxr_linenumber" name="L68" href="#L68">68</a>
+<a class="jxr_linenumber" name="L69" href="#L69">69</a> @ParameterizedTest
+<a class="jxr_linenumber" name="L70" href="#L70">70</a> @MethodSource
+<a class="jxr_linenumber" name="L71" href="#L71">71</a> <strong class="jxr_keyword">void</strong> testAdditionalMoments(<strong class="jxr_keyword">double</strong> a, <strong class="jxr_keyword">double</strong> b, <strong class="jxr_keyword">double</strong> c, <strong class="jxr_keyword">double</strong> mean, <strong class="jxr_keyword">double</strong> variance) {
+<a class="jxr_linenumber" name="L72" href="#L72">72</a> <strong class="jxr_keyword">final</strong> TriangularDistribution dist = TriangularDistribution.of(a, b, c);
+<a class="jxr_linenumber" name="L73" href="#L73">73</a> testMoments(dist, mean, variance, DoubleTolerances.equals());
+<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">static</strong> Stream<Arguments> testAdditionalMoments() {
+<a class="jxr_linenumber" name="L77" href="#L77">77</a> <strong class="jxr_keyword">return</strong> Stream.of(
+<a class="jxr_linenumber" name="L78" href="#L78">78</a> Arguments.of(0, 0.5, 1.0, 0.5, 1 / 24.0),
+<a class="jxr_linenumber" name="L79" href="#L79">79</a> Arguments.of(0, 1, 1, 2 / 3.0, 1 / 18.0),
+<a class="jxr_linenumber" name="L80" href="#L80">80</a> Arguments.of(-3, 2, 12, 3 + (2 / 3.0), 175 / 18.0)
+<a class="jxr_linenumber" name="L81" href="#L81">81</a> );
+<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> @ParameterizedTest
+<a class="jxr_linenumber" name="L85" href="#L85">85</a> @CsvSource({
+<a class="jxr_linenumber" name="L86" href="#L86">86</a> <span class="jxr_string">"1, 2, 3"</span>,
+<a class="jxr_linenumber" name="L87" href="#L87">87</a> <span class="jxr_string">"0.12, 3.45, 12.56"</span>,
+<a class="jxr_linenumber" name="L88" href="#L88">88</a> })
+<a class="jxr_linenumber" name="L89" href="#L89">89</a> <strong class="jxr_keyword">void</strong> testAdditionalParameterAccessors(<strong class="jxr_keyword">double</strong> lower, <strong class="jxr_keyword">double</strong> mode, <strong class="jxr_keyword">double</strong> upper) {
+<a class="jxr_linenumber" name="L90" href="#L90">90</a> <strong class="jxr_keyword">final</strong> TriangularDistribution dist = TriangularDistribution.of(lower, mode, upper);
+<a class="jxr_linenumber" name="L91" href="#L91">91</a> Assertions.assertEquals(lower, dist.getSupportLowerBound());
+<a class="jxr_linenumber" name="L92" href="#L92">92</a> Assertions.assertEquals(mode, dist.getMode());
+<a class="jxr_linenumber" name="L93" href="#L93">93</a> Assertions.assertEquals(upper, dist.getSupportUpperBound());
+<a class="jxr_linenumber" name="L94" href="#L94">94</a> }
+<a class="jxr_linenumber" name="L95" href="#L95">95</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>