@@ -24,46 +24,38 @@ TEST(MathFunctions, inv_Phi) {
2424TEST (MathFunctions, Equal) {
2525 using stan::math::inv_Phi;
2626 // test output generated with WolframAlpha
27- constexpr std::array p_tails{
28- 0.000000001 ,
29- 0.0000001 ,
30- 0.99999 ,
31- 0.9999999 ,
32- 0.999999999
33- };
34- constexpr std::array p{0.00001 , 0.001 , 0.05 , 0.15 , 0.25 ,
35- 0.35 , 0.45 , 0.55 , 0.65 , 0.75 , 0.85 , 0.95 ,
36- 0.999 };
27+ constexpr std::array p_tails{0.000000001 , 0.0000001 , 0.99999 , 0.9999999 ,
28+ 0.999999999 };
29+ constexpr std::array p{0.00001 , 0.001 , 0.05 , 0.15 , 0.25 , 0.35 , 0.45 ,
30+ 0.55 , 0.65 , 0.75 , 0.85 , 0.95 , 0.999 };
3731 constexpr std::array exact_tails{
38- -5.9978070150076868715623102049115374195951202210145432633059905935 ,
39- -5.199337582192816931587347266962336866509737160238716454318950531 ,
40- 4.2648907939228246284985246989063446293560532226954907262010508062 ,
41- 5.1993375821928169315873472669623368665097371602387164543189505310 ,
42- 5.9978070150076868715623102049115374195951202210145432633059905935
43- };
32+ -5.9978070150076868715623102049115374195951202210145432633059905935 ,
33+ -5.199337582192816931587347266962336866509737160238716454318950531 ,
34+ 4.2648907939228246284985246989063446293560532226954907262010508062 ,
35+ 5.1993375821928169315873472669623368665097371602387164543189505310 ,
36+ 5.9978070150076868715623102049115374195951202210145432633059905935 };
4437 constexpr std::array exact{
45- -4.264890793922824628498524698906344629356053222695490726201050806 ,
46- -3.090232306167813541540399830107379205491008491865808855697171108 ,
47- -1.644853626951472714863848907991632136083195744275322071769672094 ,
48- -1.036433389493789579713244074673503366134740595985917627904454866 ,
49- -0.674489750196081743202227014541307185386904415049861895662093788 ,
50- -0.385320466407567623810762390270473572721091620003945341485947546 ,
51- -0.125661346855074034210184388300799303397350646690021834224555394 ,
52- 0.1256613468550740342101843883007993033973506466900218342245553948 ,
53- 0.3853204664075676238107623902704735727210916200039453414859475461 ,
54- 0.6744897501960817432022270145413071853869044150498618956620937885 ,
55- 1.0364333894937895797132440746735033661347405959859176279044548663 ,
56- 1.6448536269514727148638489079916321360831957442753220717696720944 ,
57- 3.0902323061678135415403998301073792054910084918658088556971711085 };
58-
38+ -4.264890793922824628498524698906344629356053222695490726201050806 ,
39+ -3.090232306167813541540399830107379205491008491865808855697171108 ,
40+ -1.644853626951472714863848907991632136083195744275322071769672094 ,
41+ -1.036433389493789579713244074673503366134740595985917627904454866 ,
42+ -0.674489750196081743202227014541307185386904415049861895662093788 ,
43+ -0.385320466407567623810762390270473572721091620003945341485947546 ,
44+ -0.125661346855074034210184388300799303397350646690021834224555394 ,
45+ 0.1256613468550740342101843883007993033973506466900218342245553948 ,
46+ 0.3853204664075676238107623902704735727210916200039453414859475461 ,
47+ 0.6744897501960817432022270145413071853869044150498618956620937885 ,
48+ 1.0364333894937895797132440746735033661347405959859176279044548663 ,
49+ 1.6448536269514727148638489079916321360831957442753220717696720944 ,
50+ 3.0902323061678135415403998301073792054910084918658088556971711085 };
5951
6052 for (int i = 0 ; i < p_tails.size (); ++i) {
61- EXPECT_NEAR (exact_tails[i], inv_Phi (p_tails[i]), 4.6e-9 ) <<
62- " at p_tails(" << i << " ) = " << p_tails[i];
53+ EXPECT_NEAR (exact_tails[i], inv_Phi (p_tails[i]), 4.6e-9 )
54+ << " at p_tails(" << i << " ) = " << p_tails[i];
6355 }
6456 for (int i = 0 ; i < p.size (); ++i) {
65- EXPECT_NEAR (exact[i], inv_Phi (p[i]), 1.5e-15 ) <<
66- " at p(" << i << " ) = " << p[i];
57+ EXPECT_NEAR (exact[i], inv_Phi (p[i]), 1.5e-15 )
58+ << " at p(" << i << " ) = " << p[i];
6759 }
6860}
6961
0 commit comments