@@ -151,26 +151,26 @@ subroutine test_expon_rvs_csp
151
151
complex (sp) :: res(10 ), scale
152
152
integer :: i, n, k = 5
153
153
integer :: seed, get
154
- complex (sp) :: ans(10 ) = [(0.640164494514465332031250000000000000_sp , &
155
- 0.268204987049102783203125000000000000_sp ), &
156
- (0.141667589545249938964843750000000000_sp , &
157
- 2.60843825340270996093750000000000000_sp ), &
158
- (0.116705909371376037597656250000000000_sp , &
159
- 1.04120576381683349609375000000000000_sp ), &
160
- (2.06742310523986816406250000000000000_sp , &
161
- 0.723668336868286132812500000000000000_sp ), &
162
- (0.852514624595642089843750000000000000_sp , &
163
- 0.789043843746185302734375000000000000_sp ), &
164
- (1.09098446369171142578125000000000000_sp , &
165
- 1.48569476604461669921875000000000000_sp ), &
166
- (4.29633092880249023437500000000000000_sp , &
167
- 0.338216394186019897460937500000000000_sp ), &
168
- (0.340462744235992431640625000000000000_sp , &
169
- 0.172319442033767700195312500000000000_sp ), &
170
- (6.932352483272552490234375000000000000E -0002_sp , &
171
- 6.742518395185470581054687500000000000E -0002_sp ), &
172
- (1.03231632709503173828125000000000000_sp , &
173
- 0.421413004398345947265625000000000000_sp )]
154
+ complex (sp) :: ans(10 ) = [(0.640164505354053137153869101894088070_sp , &
155
+ 0.268204995991944639133208738712710328_sp ), &
156
+ (0.141667583643866407117073435983911608_sp , &
157
+ 2.60843818343998017361684560455614716_sp ), &
158
+ (0.116705911955157926040982374615850854_sp , &
159
+ 1.04120572522889689448533090398996145_sp ), &
160
+ (2.06742299986664055388985161698656149_sp , &
161
+ 0.723668357086202707062483341360348315_sp ), &
162
+ (0.852514651480446872255924972705542983_sp , &
163
+ 0.789043817146060844081034701957833041_sp ), &
164
+ (1.09098446025458177821576555288629603_sp , &
165
+ 1.48569478096399141264782883808948111_sp ), &
166
+ (4.29633077423243321391055360436439499_sp , &
167
+ 0.338216396454583145825267820328008412_sp ), &
168
+ (0.340462754402863337910289942556119029_sp , &
169
+ 0.172319442815022222381671213042864120_sp ), &
170
+ (6.932352666201882229746189523211795805E -0002_sp , &
171
+ 6.742518436285274002761624956292507704E -0002_sp ), &
172
+ (1.03231628501970258415809666985296648_sp , &
173
+ 0.421413014732743429480166241773986277_sp )]
174
174
175
175
print * , " Test exponential_distribution_rvs_csp"
176
176
seed = 593742186
@@ -188,26 +188,26 @@ subroutine test_expon_rvs_cdp
188
188
complex (dp) :: res(10 ), scale
189
189
integer :: i, n, k = 5
190
190
integer :: seed, get
191
- complex (dp) :: ans(10 ) = [(0.640164494514465332031250000000000000_dp , &
192
- 0.268204987049102783203125000000000000_dp ), &
193
- (0.141667589545249938964843750000000000_dp , &
194
- 2.60843825340270996093750000000000000_dp ), &
195
- (0.116705909371376037597656250000000000_dp , &
196
- 1.04120576381683349609375000000000000_dp ), &
197
- (2.06742310523986816406250000000000000_dp , &
198
- 0.723668336868286132812500000000000000_dp ), &
199
- (0.852514624595642089843750000000000000_dp , &
200
- 0.789043843746185302734375000000000000_dp ), &
201
- (1.09098446369171142578125000000000000_dp , &
202
- 1.48569476604461669921875000000000000_dp ), &
203
- (4.29633092880249023437500000000000000_dp , &
204
- 0.338216394186019897460937500000000000_dp ), &
205
- (0.340462744235992431640625000000000000_dp , &
206
- 0.172319442033767700195312500000000000_dp ), &
207
- (6.932352483272552490234375000000000000E -0002_dp , &
208
- 6.742518395185470581054687500000000000E -0002_dp ), &
209
- (1.03231632709503173828125000000000000_dp , &
210
- 0.421413004398345947265625000000000000_dp )]
191
+ complex (dp) :: ans(10 ) = [(0.640164505354053137153869101894088070_dp , &
192
+ 0.268204995991944639133208738712710328_dp ), &
193
+ (0.141667583643866407117073435983911608_dp , &
194
+ 2.60843818343998017361684560455614716_dp ), &
195
+ (0.116705911955157926040982374615850854_dp , &
196
+ 1.04120572522889689448533090398996145_dp ), &
197
+ (2.06742299986664055388985161698656149_dp , &
198
+ 0.723668357086202707062483341360348315_dp ), &
199
+ (0.852514651480446872255924972705542983_dp , &
200
+ 0.789043817146060844081034701957833041_dp ), &
201
+ (1.09098446025458177821576555288629603_dp , &
202
+ 1.48569478096399141264782883808948111_dp ), &
203
+ (4.29633077423243321391055360436439499_dp , &
204
+ 0.338216396454583145825267820328008412_dp ), &
205
+ (0.340462754402863337910289942556119029_dp , &
206
+ 0.172319442815022222381671213042864120_dp ), &
207
+ (6.932352666201882229746189523211795805E -0002_dp , &
208
+ 6.742518436285274002761624956292507704E -0002_dp ), &
209
+ (1.03231628501970258415809666985296648_dp , &
210
+ 0.421413014732743429480166241773986277_dp )]
211
211
212
212
print * , " Test exponential_distribution_rvs_cdp"
213
213
seed = 593742186
@@ -225,26 +225,26 @@ subroutine test_expon_rvs_cqp
225
225
complex (qp) :: res(10 ), scale
226
226
integer :: i, n, k = 5
227
227
integer :: seed, get
228
- complex (qp) :: ans(10 ) = [(0.640164494514465332031250000000000000_qp , &
229
- 0.268204987049102783203125000000000000_qp ), &
230
- (0.141667589545249938964843750000000000_qp , &
231
- 2.60843825340270996093750000000000000_qp ), &
232
- (0.116705909371376037597656250000000000_qp , &
233
- 1.04120576381683349609375000000000000_qp ), &
234
- (2.06742310523986816406250000000000000_qp , &
235
- 0.723668336868286132812500000000000000_qp ), &
236
- (0.852514624595642089843750000000000000_qp , &
237
- 0.789043843746185302734375000000000000_qp ), &
238
- (1.09098446369171142578125000000000000_qp , &
239
- 1.48569476604461669921875000000000000_qp ), &
240
- (4.29633092880249023437500000000000000_qp , &
241
- 0.338216394186019897460937500000000000_qp ), &
242
- (0.340462744235992431640625000000000000_qp , &
243
- 0.172319442033767700195312500000000000_qp ), &
244
- (6.932352483272552490234375000000000000E -0002_qp , &
245
- 6.742518395185470581054687500000000000E -0002_qp ), &
246
- (1.03231632709503173828125000000000000_qp , &
247
- 0.421413004398345947265625000000000000_qp )]
228
+ complex (qp) :: ans(10 ) = [(0.640164505354053137153869101894088070_qp , &
229
+ 0.268204995991944639133208738712710328_qp ), &
230
+ (0.141667583643866407117073435983911608_qp , &
231
+ 2.60843818343998017361684560455614716_qp ), &
232
+ (0.116705911955157926040982374615850854_qp , &
233
+ 1.04120572522889689448533090398996145_qp ), &
234
+ (2.06742299986664055388985161698656149_qp , &
235
+ 0.723668357086202707062483341360348315_qp ), &
236
+ (0.852514651480446872255924972705542983_qp , &
237
+ 0.789043817146060844081034701957833041_qp ), &
238
+ (1.09098446025458177821576555288629603_qp , &
239
+ 1.48569478096399141264782883808948111_qp ), &
240
+ (4.29633077423243321391055360436439499_qp , &
241
+ 0.338216396454583145825267820328008412_qp ), &
242
+ (0.340462754402863337910289942556119029_qp , &
243
+ 0.172319442815022222381671213042864120_qp ), &
244
+ (6.932352666201882229746189523211795805E -0002_qp , &
245
+ 6.742518436285274002761624956292507704E -0002_qp ), &
246
+ (1.03231628501970258415809666985296648_qp , &
247
+ 0.421413014732743429480166241773986277_qp )]
248
248
249
249
print * , " Test exponential_distribution_rvs_cqp"
250
250
seed = 593742186
@@ -334,11 +334,11 @@ subroutine test_expon_pdf_csp
334
334
integer :: i, n
335
335
integer :: seed, get
336
336
real :: res(3 ,5 )
337
- real :: ans(15 ) = [0.386471182 , 0.386471182 ,0.386471182 , 2.79592816E -03 ,&
338
- 4.01333207E-03 , 0.317740440 , 0.385551631 , &
339
- 5.02163824E -03 , 0.372386932 , 6.09764457E-03 , &
340
- 0.273956627 , 0.407586545 , 1.59074634E -03 ,0.136133015 , &
341
- 0.399842113 ]
337
+ real :: ans(15 ) = [0.386471182 , 0.386471182 ,0.386471182 ,2.79592793E -03 , &
338
+ 4.01333207E-03 , 0.317740440 , 0.385551631 , &
339
+ 5.02163777E -03 , 0.372386932 , 6.09764457E-03 , &
340
+ 0.273956627 , 0.407586545 , 1.59074657E -03 , &
341
+ 0.136133000 , 0.399842113 ]
342
342
343
343
print * , " Test exponential_distribution_pdf_csp"
344
344
seed = 123987654
@@ -357,11 +357,11 @@ subroutine test_expon_pdf_cdp
357
357
integer :: i, n
358
358
integer :: seed, get
359
359
real :: res(3 ,5 )
360
- real :: ans(15 ) = [0.386471182 , 0.386471182 ,0.386471182 , 2.79592816E -03 ,&
361
- 4.01333207E-03 , 0.317740440 , 0.385551631 , &
362
- 5.02163824E -03 , 0.372386932 , 6.09764457E-03 , &
363
- 0.273956627 , 0.407586545 , 1.59074634E -03 ,0.136133015 , &
364
- 0.399842113 ]
360
+ real :: ans(15 ) = [0.386471182 , 0.386471182 ,0.386471182 ,2.79592793E -03 , &
361
+ 4.01333207E-03 , 0.317740440 , 0.385551631 , &
362
+ 5.02163777E -03 , 0.372386932 , 6.09764457E-03 , &
363
+ 0.273956627 , 0.407586545 , 1.59074657E -03 , &
364
+ 0.136133000 , 0.399842113 ]
365
365
366
366
print * , " Test exponential_distribution_pdf_cdp"
367
367
seed = 123987654
@@ -380,11 +380,11 @@ subroutine test_expon_pdf_cqp
380
380
integer :: i, n
381
381
integer :: seed, get
382
382
real :: res(3 ,5 )
383
- real :: ans(15 ) = [0.386471182 , 0.386471182 ,0.386471182 , 2.79592816E -03 ,&
384
- 4.01333207E-03 , 0.317740440 , 0.385551631 , &
385
- 5.02163824E -03 , 0.372386932 , 6.09764457E-03 , &
386
- 0.273956627 , 0.407586545 , 1.59074634E -03 ,0.136133015 , &
387
- 0.399842113 ]
383
+ real :: ans(15 ) = [0.386471182 , 0.386471182 ,0.386471182 ,2.79592793E -03 , &
384
+ 4.01333207E-03 , 0.317740440 , 0.385551631 , &
385
+ 5.02163777E -03 , 0.372386932 , 6.09764457E-03 , &
386
+ 0.273956627 , 0.407586545 , 1.59074657E -03 , &
387
+ 0.136133000 , 0.399842113 ]
388
388
389
389
print * , " Test exponential_distribution_pdf_cqp"
390
390
seed = 123987654
@@ -473,10 +473,10 @@ subroutine test_expon_cdf_csp
473
473
integer :: i, n
474
474
integer :: seed, get
475
475
real :: res(3 ,5 )
476
- real :: ans(15 ) = [0.176930442 ,0.176930442 ,0.176930442 , 5.98644949E -02 , &
477
- 0.981560826 ,0.135309443 ,0.617795825 , 7.55468607E -02 , &
478
- 0.118341736 , 0.484595388 , 0.794088185 , 0.912919402 , &
479
- 0.914170802 , 0.370377690 , 0.793968201 ]
476
+ real :: ans(15 ) = [0.176930442 , 0.176930442 , 0.176930442 ,5.98644912E -02 ,&
477
+ 0.981560826 , 0.135309443 , 0.617795825 ,7.55468532E -02 ,&
478
+ 0.118341751 , 0.484595388 , 0.794088185 , 0.912919402 , &
479
+ 0.914170802 , 0.370377690 , 0.793968141 ]
480
480
481
481
print * , " Test exponential_distribution_cdf_csp"
482
482
seed = 621957438
@@ -496,10 +496,10 @@ subroutine test_expon_cdf_cdp
496
496
integer :: i, n
497
497
integer :: seed, get
498
498
real :: res(3 ,5 )
499
- real :: ans(15 ) = [0.176930442 ,0.176930442 ,0.176930442 , 5.98644949E -02 , &
500
- 0.981560826 ,0.135309443 ,0.617795825 , 7.55468607E -02 , &
501
- 0.118341736 , 0.484595388 , 0.794088185 , 0.912919402 , &
502
- 0.914170802 , 0.370377690 , 0.793968201 ]
499
+ real :: ans(15 ) = [0.176930442 , 0.176930442 , 0.176930442 ,5.98644912E -02 ,&
500
+ 0.981560826 , 0.135309443 , 0.617795825 ,7.55468532E -02 ,&
501
+ 0.118341751 , 0.484595388 , 0.794088185 , 0.912919402 , &
502
+ 0.914170802 , 0.370377690 , 0.793968141 ]
503
503
504
504
print * , " Test exponential_distribution_cdf_cdp"
505
505
seed = 621957438
@@ -519,10 +519,10 @@ subroutine test_expon_cdf_cqp
519
519
integer :: i, n
520
520
integer :: seed, get
521
521
real :: res(3 ,5 )
522
- real :: ans(15 ) = [0.176930442 ,0.176930442 ,0.176930442 , 5.98644949E -02 , &
523
- 0.981560826 ,0.135309443 ,0.617795825 , 7.55468607E -02 , &
524
- 0.118341736 , 0.484595388 , 0.794088185 , 0.912919402 , &
525
- 0.914170802 , 0.370377690 , 0.793968201 ]
522
+ real :: ans(15 ) = [0.176930442 , 0.176930442 , 0.176930442 ,5.98644912E -02 ,&
523
+ 0.981560826 , 0.135309443 , 0.617795825 ,7.55468532E -02 ,&
524
+ 0.118341751 , 0.484595388 , 0.794088185 , 0.912919402 , &
525
+ 0.914170802 , 0.370377690 , 0.793968141 ]
526
526
527
527
print * , " Test exponential_distribution_cdf_cqp"
528
528
seed = 621957438
0 commit comments