Skip to content

Commit a605a89

Browse files
Jim-215-FisherJim-215-Fisher
authored
chg. complex number with kinds
1 parent 374c4b9 commit a605a89

File tree

1 file changed

+84
-84
lines changed

1 file changed

+84
-84
lines changed

‎src/tests/stats/test_distribution_normal.f90

Lines changed: 84 additions & 84 deletions
Original file line numberDiff line numberDiff line change
@@ -158,26 +158,26 @@ subroutine test_nor_rvs_csp
158158
complex(sp) :: res(10), loc, scale
159159
integer :: i, n, k = 5
160160
integer :: seed, get
161-
complex(sp) :: ans(10) = [(2.12531018257141113281250000000000000_sp, &
162-
1.46507704257965087890625000000000000_sp), &
163-
(1.08284163475036621093750000000000000_sp, &
164-
0.277168631553649902343750000000000000_sp), &
165-
(1.41924941539764404296875000000000000_sp, &
166-
0.498445570468902587890625000000000000_sp), &
167-
(1.72639131546020507812500000000000000_sp, &
168-
0.715802907943725585937500000000000000_sp), &
169-
(1.98950588703155517578125000000000000_sp, &
170-
0.115721315145492553710937500000000000_sp), &
171-
(-1.16929018497467041015625000000000000_sp, &
172-
0.250744730234146118164062500000000000_sp), &
173-
(1.57160544395446777343750000000000000_sp, &
174-
0.638282597064971923828125000000000000_sp), &
175-
(-1.36106109619140625000000000000000000_sp, &
176-
0.166259199380874633789062500000000000_sp), &
177-
(1.13403093814849853515625000000000000_sp, &
178-
1.04232621192932128906250000000000000_sp), &
179-
(-1.68220531940460205078125000000000000_sp, &
180-
1.63361442089080810546875000000000000_sp)]
161+
complex(sp) :: ans(10) = [(2.12531029488530509574673033057479188_sp, &
162+
1.46507698734032082432676702410390135_sp), &
163+
(1.08284164094813181722365413861552952_sp, &
164+
0.277168639672963013076412153168348595_sp), &
165+
(1.41924946329521489696290359461272601_sp, &
166+
0.498445561155580918466512230224907398_sp), &
167+
(1.72639126368764062036120776610914618_sp, &
168+
0.715802936564464420410303091557580046_sp), &
169+
(1.98950590834134349860207180427096318_sp, &
170+
0.115721315405046931701349421928171068_sp), &
171+
(-1.16929014824793620075382705181255005_sp, &
172+
0.250744737486995217246033007540972903_sp), &
173+
(1.57160542831869509683428987045772374_sp, &
174+
0.638282596371312238581197107123443857_sp), &
175+
(-1.36106107654239116833139178197598085_sp, &
176+
0.166259201494369124318950525776017457_sp), &
177+
(1.13403096805387920698038328737311531_sp, &
178+
1.04232618148691447146347854868508875_sp), &
179+
(-1.68220535920475811053620418533682823_sp, &
180+
1.63361446685040256898702182297711261_sp)]
181181

182182
print *, "Test normal_distribution_rvs_csp"
183183
seed = 25836914
@@ -197,26 +197,26 @@ subroutine test_nor_rvs_cdp
197197
complex(dp) :: res(10), loc, scale
198198
integer :: i, n, k = 5
199199
integer :: seed, get
200-
complex(dp) :: ans(10) = [(2.12531018257141113281250000000000000_dp, &
201-
1.46507704257965087890625000000000000_dp), &
202-
(1.08284163475036621093750000000000000_dp, &
203-
0.277168631553649902343750000000000000_dp), &
204-
(1.41924941539764404296875000000000000_dp, &
205-
0.498445570468902587890625000000000000_dp), &
206-
(1.72639131546020507812500000000000000_dp, &
207-
0.715802907943725585937500000000000000_dp), &
208-
(1.98950588703155517578125000000000000_dp, &
209-
0.115721315145492553710937500000000000_dp), &
210-
(-1.16929018497467041015625000000000000_dp, &
211-
0.250744730234146118164062500000000000_dp), &
212-
(1.57160544395446777343750000000000000_dp, &
213-
0.638282597064971923828125000000000000_dp), &
214-
(-1.36106109619140625000000000000000000_dp, &
215-
0.166259199380874633789062500000000000_dp), &
216-
(1.13403093814849853515625000000000000_dp, &
217-
1.04232621192932128906250000000000000_dp), &
218-
(-1.68220531940460205078125000000000000_dp, &
219-
1.63361442089080810546875000000000000_dp)]
200+
complex(dp) :: ans(10) = [(2.12531029488530509574673033057479188_dp, &
201+
1.46507698734032082432676702410390135_dp), &
202+
(1.08284164094813181722365413861552952_dp, &
203+
0.277168639672963013076412153168348595_dp), &
204+
(1.41924946329521489696290359461272601_dp, &
205+
0.498445561155580918466512230224907398_dp), &
206+
(1.72639126368764062036120776610914618_dp, &
207+
0.715802936564464420410303091557580046_dp), &
208+
(1.98950590834134349860207180427096318_dp, &
209+
0.115721315405046931701349421928171068_dp), &
210+
(-1.16929014824793620075382705181255005_dp, &
211+
0.250744737486995217246033007540972903_dp), &
212+
(1.57160542831869509683428987045772374_dp, &
213+
0.638282596371312238581197107123443857_dp), &
214+
(-1.36106107654239116833139178197598085_dp, &
215+
0.166259201494369124318950525776017457_dp), &
216+
(1.13403096805387920698038328737311531_dp, &
217+
1.04232618148691447146347854868508875_dp), &
218+
(-1.68220535920475811053620418533682823_dp, &
219+
1.63361446685040256898702182297711261_dp)]
220220

221221
print *, "Test normal_distribution_rvs_cdp"
222222
seed = 25836914
@@ -236,26 +236,26 @@ subroutine test_nor_rvs_cqp
236236
complex(qp) :: res(10), loc, scale
237237
integer :: i, n, k = 5
238238
integer :: seed, get
239-
complex(qp) :: ans(10) = [(2.12531018257141113281250000000000000_qp, &
240-
1.46507704257965087890625000000000000_qp), &
241-
(1.08284163475036621093750000000000000_qp, &
242-
0.277168631553649902343750000000000000_qp), &
243-
(1.41924941539764404296875000000000000_qp, &
244-
0.498445570468902587890625000000000000_qp), &
245-
(1.72639131546020507812500000000000000_qp, &
246-
0.715802907943725585937500000000000000_qp), &
247-
(1.98950588703155517578125000000000000_qp, &
248-
0.115721315145492553710937500000000000_qp), &
249-
(-1.16929018497467041015625000000000000_qp, &
250-
0.250744730234146118164062500000000000_qp), &
251-
(1.57160544395446777343750000000000000_qp, &
252-
0.638282597064971923828125000000000000_qp), &
253-
(-1.36106109619140625000000000000000000_qp, &
254-
0.166259199380874633789062500000000000_qp), &
255-
(1.13403093814849853515625000000000000_qp, &
256-
1.04232621192932128906250000000000000_qp), &
257-
(-1.68220531940460205078125000000000000_qp, &
258-
1.63361442089080810546875000000000000_qp)]
239+
complex(qp) :: ans(10) = [(2.12531029488530509574673033057479188_qp, &
240+
1.46507698734032082432676702410390135_qp), &
241+
(1.08284164094813181722365413861552952_qp, &
242+
0.277168639672963013076412153168348595_qp), &
243+
(1.41924946329521489696290359461272601_qp, &
244+
0.498445561155580918466512230224907398_qp), &
245+
(1.72639126368764062036120776610914618_qp, &
246+
0.715802936564464420410303091557580046_qp), &
247+
(1.98950590834134349860207180427096318_qp, &
248+
0.115721315405046931701349421928171068_qp), &
249+
(-1.16929014824793620075382705181255005_qp, &
250+
0.250744737486995217246033007540972903_qp), &
251+
(1.57160542831869509683428987045772374_qp, &
252+
0.638282596371312238581197107123443857_qp), &
253+
(-1.36106107654239116833139178197598085_qp, &
254+
0.166259201494369124318950525776017457_qp), &
255+
(1.13403096805387920698038328737311531_qp, &
256+
1.04232618148691447146347854868508875_qp), &
257+
(-1.68220535920475811053620418533682823_qp, &
258+
1.63361446685040256898702182297711261_qp)]
259259

260260
print *, "Test normal_distribution_rvs_cqp"
261261
seed = 25836914
@@ -347,10 +347,10 @@ subroutine test_nor_pdf_csp
347347
integer :: i, n, k = 5
348348
integer :: seed, get
349349
real :: res(3,5)
350-
real :: ans(15) = [0.129377306, 0.129377306, 0.129377306,4.05915640E-02,&
351-
0.209143385, 2.98881028E-02,0.128679410, 0.177484736,&
352-
3.82205285E-02, 7.09915683E-02, 4.56126593E-02, &
353-
6.57454208E-02,0.165161029,3.86104845E-02,0.196922958]
350+
real :: ans(15) = [0.129377320, 0.129377320,0.129377320,4.05915640E-02, &
351+
0.209143385,2.98881028E-02, 0.128679410, 0.177484736,&
352+
3.82205322E-02, 7.09915683E-02, 4.56126593E-02, &
353+
6.57454133E-02,0.165161043,3.86104807E-02,0.196922958]
354354

355355
print *, "Test normal_distribution_pdf_csp"
356356
seed = 741852963
@@ -370,10 +370,10 @@ subroutine test_nor_pdf_cdp
370370
integer :: i, n, k = 5
371371
integer :: seed, get
372372
real :: res(3,5)
373-
real :: ans(15) = [0.129377306, 0.129377306, 0.129377306,4.05915640E-02,&
374-
0.209143385, 2.98881028E-02,0.128679410, 0.177484736,&
375-
3.82205285E-02, 7.09915683E-02, 4.56126593E-02, &
376-
6.57454208E-02,0.165161029,3.86104845E-02,0.196922958]
373+
real :: ans(15) = [0.129377320, 0.129377320,0.129377320,4.05915640E-02, &
374+
0.209143385,2.98881028E-02, 0.128679410, 0.177484736,&
375+
3.82205322E-02, 7.09915683E-02, 4.56126593E-02, &
376+
6.57454133E-02,0.165161043,3.86104807E-02,0.196922958]
377377

378378
print *, "Test normal_distribution_pdf_cdp"
379379
seed = 741852963
@@ -393,10 +393,10 @@ subroutine test_nor_pdf_cqp
393393
integer :: i, n, k = 5
394394
integer :: seed, get
395395
real :: res(3,5)
396-
real :: ans(15) = [0.129377306, 0.129377306, 0.129377306,4.05915640E-02,&
397-
0.209143385, 2.98881028E-02,0.128679410, 0.177484736,&
398-
3.82205285E-02, 7.09915683E-02, 4.56126593E-02, &
399-
6.57454208E-02,0.165161029,3.86104845E-02,0.196922958]
396+
real :: ans(15) = [0.129377320, 0.129377320,0.129377320,4.05915640E-02, &
397+
0.209143385,2.98881028E-02, 0.128679410, 0.177484736,&
398+
3.82205322E-02, 7.09915683E-02, 4.56126593E-02, &
399+
6.57454133E-02,0.165161043,3.86104807E-02,0.196922958]
400400

401401
print *, "Test normal_distribution_pdf_cqp"
402402
seed = 741852963
@@ -486,11 +486,11 @@ subroutine test_nor_cdf_csp
486486
integer :: i, n
487487
integer :: seed, get
488488
real :: res(3,5)
489-
real :: ans(15) = [1.07458141E-02, 1.07458141E-02, 1.07458141E-02, &
490-
6.86483160E-02, 7.95486569E-02, 2.40523368E-02, &
491-
3.35096717E-02, 0.315778911,0.446311295,0.102010213, &
489+
real :: ans(15) = [1.07458131E-02, 1.07458131E-02, 1.07458131E-02, &
490+
6.86483234E-02, 7.95486644E-02, 2.40523387E-02, &
491+
3.35096754E-02,0.315778911,0.446311295, 0.102010213, &
492492
7.66918957E-02, 0.564691007, 0.708769500, &
493-
6.40553832E-02, 5.39999120E-02]
493+
6.40553832E-02, 5.39999157E-02]
494494

495495
print *, "Test normal_distribution_cdf_csp"
496496
seed = 369147582
@@ -510,11 +510,11 @@ subroutine test_nor_cdf_cdp
510510
integer :: i, n
511511
integer :: seed, get
512512
real :: res(3,5)
513-
real :: ans(15) = [1.07458141E-02, 1.07458141E-02, 1.07458141E-02, &
514-
6.86483160E-02, 7.95486569E-02, 2.40523368E-02, &
515-
3.35096717E-02, 0.315778911,0.446311295,0.102010213, &
513+
real :: ans(15) = [1.07458131E-02, 1.07458131E-02, 1.07458131E-02, &
514+
6.86483234E-02, 7.95486644E-02, 2.40523387E-02, &
515+
3.35096754E-02,0.315778911,0.446311295, 0.102010213, &
516516
7.66918957E-02, 0.564691007, 0.708769500, &
517-
6.40553832E-02, 5.39999120E-02]
517+
6.40553832E-02, 5.39999157E-02]
518518

519519
print *, "Test normal_distribution_cdf_cdp"
520520
seed = 369147582
@@ -534,11 +534,11 @@ subroutine test_nor_cdf_cqp
534534
integer :: i, n
535535
integer :: seed, get
536536
real :: res(3,5)
537-
real :: ans(15) = [1.07458141E-02, 1.07458141E-02, 1.07458141E-02, &
538-
6.86483160E-02, 7.95486569E-02, 2.40523368E-02, &
539-
3.35096717E-02, 0.315778911,0.446311295,0.102010213, &
537+
real :: ans(15) = [1.07458131E-02, 1.07458131E-02, 1.07458131E-02, &
538+
6.86483234E-02, 7.95486644E-02, 2.40523387E-02, &
539+
3.35096754E-02,0.315778911,0.446311295, 0.102010213, &
540540
7.66918957E-02, 0.564691007, 0.708769500, &
541-
6.40553832E-02, 5.39999120E-02]
541+
6.40553832E-02, 5.39999157E-02]
542542

543543
print *, "Test normal_distribution_cdf_cqp"
544544
seed = 369147582

0 commit comments

Comments
 (0)