1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17 package org.apache.commons.rng.sampling.distribution;
18
19 import org.apache.commons.rng.UniformRandomProvider;
20 import org.apache.commons.rng.sampling.RandomAssert;
21 import org.junit.jupiter.api.Assertions;
22 import org.junit.jupiter.api.Test;
23
24
25
26
27 class AhrensDieterMarsagliaTsangGammaSamplerTest {
28
29
30
31 @Test
32 void testConstructorThrowsWithZeroAlpha() {
33 final UniformRandomProvider rng = RandomAssert.seededRNG();
34 final double alpha = 0;
35 final double theta = 1;
36 Assertions.assertThrows(IllegalArgumentException.class,
37 () -> AhrensDieterMarsagliaTsangGammaSampler.of(rng, alpha, theta));
38 }
39
40
41
42
43 @Test
44 void testConstructorThrowsWithZeroTheta() {
45 final UniformRandomProvider rng = RandomAssert.seededRNG();
46 final double alpha = 1;
47 final double theta = 0;
48 Assertions.assertThrows(IllegalArgumentException.class,
49 () -> AhrensDieterMarsagliaTsangGammaSampler.of(rng, alpha, theta));
50 }
51
52
53
54
55 @Test
56 void testSharedStateSamplerWithAlphaBelowOne() {
57 testSharedStateSampler(0.5, 3.456);
58 }
59
60
61
62
63 @Test
64 void testSharedStateSamplerWithAlphaAboveOne() {
65 testSharedStateSampler(3.5, 3.456);
66 }
67
68
69
70
71
72
73
74 private static void testSharedStateSampler(double alpha, double theta) {
75 final UniformRandomProvider rng1 = RandomAssert.seededRNG();
76 final UniformRandomProvider rng2 = RandomAssert.seededRNG();
77
78 final AhrensDieterMarsagliaTsangGammaSampler sampler1 =
79 new AhrensDieterMarsagliaTsangGammaSampler(rng1, alpha, theta);
80 final SharedStateContinuousSampler sampler2 = sampler1.withUniformRandomProvider(rng2);
81 RandomAssert.assertProduceSameSequence(sampler1, sampler2);
82 }
83
84
85
86
87
88 @Test
89 void testToString() {
90 final UniformRandomProvider rng = RandomAssert.seededRNG();
91 Assertions.assertTrue(new AhrensDieterMarsagliaTsangGammaSampler(rng, 1.0, 2.0).toString()
92 .toLowerCase().contains("gamma"));
93 }
94 }