@@ -788,21 +788,25 @@ def test_randomdist_noncentral_f(randomdist):
788
788
np .testing .assert_allclose (actual , desired , atol = 1e-7 , rtol = 1e-10 )
789
789
790
790
791
- def test_randomdist_normal (randomdist ):
791
+ def test_randomdist_normal_case1 (randomdist ):
792
792
rnd .seed (randomdist .seed , brng = randomdist .brng )
793
793
actual = rnd .normal (loc = .123456789 , scale = 2.0 , size = (3 , 2 ))
794
794
desired = np .array ([[4.405778774782659 , - 4.020116569348963 ],
795
795
[- 1.732103577371005 , 1.2282652034983546 ],
796
796
[0.21648943171034918 , 4.625591634211608 ]])
797
797
np .testing .assert_allclose (actual , desired , atol = 1e-7 , rtol = 1e-10 )
798
798
799
+
800
+ def test_randomdist_normal_case2 (randomdist ):
799
801
rnd .seed (randomdist .seed , brng = randomdist .brng )
800
802
actual = rnd .normal (loc = .123456789 , scale = 2.0 , size = (3 , 2 ), method = "BoxMuller" )
801
803
desired = np .array ([[0.16673479781277187 , - 3.4809986872165952 ],
802
804
[- 0.05193761082535492 , 3.249201213154922 ],
803
805
[- 0.11915582299214138 , 3.555636100927892 ]])
804
806
np .testing .assert_allclose (actual , desired , atol = 1e-8 , rtol = 1e-8 )
805
807
808
+
809
+ def test_randomdist_normal_case3 (randomdist ):
806
810
rnd .seed (randomdist .seed , brng = randomdist .brng )
807
811
actual = rnd .normal (loc = .123456789 , scale = 2.0 , size = (3 , 2 ), method = "BoxMuller2" )
808
812
desired = np .array ([[0.16673479781277187 , 0.48153966449249175 ],
@@ -893,14 +897,16 @@ def test_randomdist_standard_gamma(randomdist):
893
897
np .testing .assert_allclose (actual , desired , atol = 1e-7 , rtol = 1e-10 )
894
898
895
899
896
- def test_randomdist_standard_normal (randomdist ):
900
+ def test_randomdist_standard_normal_case1 (randomdist ):
897
901
rnd .seed (randomdist .seed , brng = randomdist .brng )
898
902
actual = rnd .standard_normal (size = (3 , 2 ))
899
903
desired = np .array ([[2.1411609928913298 , - 2.071786679174482 ],
900
904
[- 0.9277801831855025 , 0.5524042072491773 ],
901
905
[0.04651632135517459 , 2.2510674226058036 ]])
902
906
np .testing .assert_allclose (actual , desired , atol = 1e-7 , rtol = 1e-10 )
903
907
908
+
909
+ def test_randomdist_standard_normal_case2 (randomdist ):
904
910
rnd .seed (randomdist .seed , brng = randomdist .brng )
905
911
actual = rnd .standard_normal (size = (3 , 2 ), method = 'BoxMuller2' )
906
912
desired = np .array ([[0.021639004406385935 , 0.17904143774624587 ],
0 commit comments