16 integer(IK),
parameter :: nsam
= 2
17 integer(IK),
allocatable :: iweight(:)
18 integer(IK) :: isam, ndim, lb(nsam), ub(nsam)
19 integer(IK) :: dim, itry, ntry
= 5
20 type(display_type) :: disp
24 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
25 call disp%show(
"!Compute the biased merged variance and mean of a univariate sample.")
26 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
30 real(TKG) :: mean(
0:nsam), meanMerged
31 real(TKG),
allocatable :: sample(:)
32 real(TKG) :: var(
0:nsam), varMerged
35 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
38 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
39 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
42 lb(isam)
= ub(isam
- 1)
+ 1
49 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
53 call disp%show(
"mean(0) = getMean(sample)")
55 call disp%show(
"mean(0) ! reference")
57 call disp%show(
"var(0) = getVar(sample)")
62 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)))")
63 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)))")
66 var(isam)
= getVar(sample(lb(isam):ub(isam)))
67 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
69 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(2), mean(2), var(1), mean(1), ub(1) / real(ub(2), TKG))")
70 call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
75 call disp%show(
"call setVarMeanMerged(var(2), mean(2), var(1), mean(1), ub(1) / real(ub(2), TKG))")
76 call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
79 call disp%show(
"mean(0) ! reference")
90 complex(TKG) :: mean(
0:nsam), meanMerged
91 complex(TKG),
allocatable :: sample(:)
92 real(TKG) :: var(
0:nsam), varMerged
95 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
98 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
99 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
102 lb(isam)
= ub(isam
- 1)
+ 1
109 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ub(nsam)), getUnifRand(-1., +1., ub(nsam)), TKG)")
113 call disp%show(
"mean(0) = getMean(sample)")
115 call disp%show(
"mean(0) ! reference")
117 call disp%show(
"var(0) = getVar(sample)")
119 call disp%show(
"var(0) ! reference")
122 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)))")
123 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)))")
126 var(isam)
= getVar(sample(lb(isam):ub(isam)))
127 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
129 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(2), mean(2), var(1), mean(1), ub(1) / real(ub(2), TKG))")
130 call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
135 call disp%show(
"call setVarMeanMerged(var(2), mean(2), var(1), mean(1), ub(1) / real(ub(2), TKG))")
136 call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
139 call disp%show(
"mean(0) ! reference")
143 call disp%show(
"var(0) ! reference")
150 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
151 call disp%show(
"!Compute the biased merged variance and mean of a frequency weighted univariate sample.")
152 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
156 real(TKG) :: mean(
0:nsam), meanMerged
157 real(TKG),
allocatable :: sample(:)
158 real(TKG) :: var(
0:nsam), varMerged
161 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
164 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
165 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
168 lb(isam)
= ub(isam
- 1)
+ 1
175 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
179 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, 1, IK))")
183 call disp%show(
"mean(0) = getMean(sample, iweight)")
184 mean(
0)
= getMean(sample, iweight)
185 call disp%show(
"mean(0) ! reference")
187 call disp%show(
"var(0) = getVar(sample, iweight)")
188 var(
0)
= getVar(sample, iweight)
189 call disp%show(
"var(0) ! reference")
192 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
193 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
196 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
197 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
199 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(2), mean(2), var(1), mean(1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
200 call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
205 call disp%show(
"call setVarMeanMerged(var(2), mean(2), var(1), mean(1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
206 call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
209 call disp%show(
"mean(0) ! reference")
213 call disp%show(
"var(0) ! reference")
220 complex(TKG) :: mean(
0:nsam), meanMerged
221 complex(TKG),
allocatable :: sample(:)
222 real(TKG) :: var(
0:nsam), varMerged
225 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
228 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
229 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
232 lb(isam)
= ub(isam
- 1)
+ 1
239 call disp%show(
"sample = cmplx(getUnifRand(0., 1., ub(nsam)), -getUnifRand(0., 1., ub(nsam)), TKG)")
243 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, 1, IK))")
247 call disp%show(
"mean(0) = getMean(sample, iweight)")
248 mean(
0)
= getMean(sample, iweight)
249 call disp%show(
"mean(0) ! reference")
251 call disp%show(
"var(0) = getVar(sample, iweight)")
252 var(
0)
= getVar(sample, iweight)
253 call disp%show(
"var(0) ! reference")
256 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
257 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
260 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
261 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
263 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(2), mean(2), var(1), mean(1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
264 call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
269 call disp%show(
"call setVarMeanMerged(var(2), mean(2), var(1), mean(1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
270 call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
273 call disp%show(
"mean(0) ! reference")
277 call disp%show(
"var(0) ! reference")
284 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
285 call disp%show(
"!Compute the biased merged variance and mean of a reliability weighted univariate sample.")
286 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
290 real(TKG) :: mean(
0:nsam), meanMerged
291 real(TKG),
allocatable :: sample(:)
292 real(TKG),
allocatable :: rweight(:)
293 real(TKG) :: var(
0:nsam), varMerged
296 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
299 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
300 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
303 lb(isam)
= ub(isam
- 1)
+ 1
310 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
314 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, 1, IK))")
318 call disp%show(
"mean(0) = getMean(sample, rweight)")
319 mean(
0)
= getMean(sample, rweight)
320 call disp%show(
"mean(0) ! reference")
322 call disp%show(
"var(0) = getVar(sample, rweight)")
323 var(
0)
= getVar(sample, rweight)
324 call disp%show(
"var(0) ! reference")
327 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
328 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
331 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
332 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
334 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(2), mean(2), var(1), mean(1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
335 call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
340 call disp%show(
"call setVarMeanMerged(var(2), mean(2), var(1), mean(1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
341 call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
344 call disp%show(
"mean(0) ! reference")
348 call disp%show(
"var(0) ! reference")
355 complex(TKG) :: mean(
0:nsam), meanMerged
356 complex(TKG),
allocatable :: sample(:)
357 real(TKG),
allocatable :: rweight(:)
358 real(TKG) :: var(
0:nsam), varMerged
361 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
364 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
365 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
368 lb(isam)
= ub(isam
- 1)
+ 1
375 call disp%show(
"sample = cmplx(getUnifRand(0., 1., ub(nsam)), -getUnifRand(0., 1., ub(nsam)), TKG)")
379 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, 1, IK))")
383 call disp%show(
"mean(0) = getMean(sample, rweight)")
384 mean(
0)
= getMean(sample, rweight)
385 call disp%show(
"mean(0) ! reference")
387 call disp%show(
"var(0) = getVar(sample, rweight)")
388 var(
0)
= getVar(sample, rweight)
389 call disp%show(
"var(0) ! reference")
392 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
393 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
396 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
397 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
399 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(2), mean(2), var(1), mean(1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
400 call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
405 call disp%show(
"call setVarMeanMerged(var(2), mean(2), var(1), mean(1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
406 call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
409 call disp%show(
"mean(0) ! reference")
413 call disp%show(
"var(0) ! reference")
420 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
421 call disp%show(
"!Compute the biased merged variance and mean of a multivariate sample.")
422 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
426 real(TKG),
allocatable :: sample(:,:)
427 real(TKG),
allocatable :: mean(:,:), meanMerged(:)
428 real(TKG),
allocatable :: var(:,:), varMerged(:)
431 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
434 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
435 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
438 lb(isam)
= ub(isam
- 1)
+ 1
445 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
447 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
448 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
449 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
450 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
451 call disp%show(
"call setResized(meanMerged, ndim)")
453 call disp%show(
"call setResized(varMerged, ndim)")
455 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
459 call disp%show(
"mean(:,0) = getVar(sample, dim)")
460 mean(:,
0)
= getVar(sample, dim)
461 call disp%show(
"mean(:,0) ! reference")
463 call disp%show(
"var(:,0) = getVar(sample, dim)")
464 var(:,
0)
= getVar(sample, dim)
465 call disp%show(
"var(:,0) ! reference")
468 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), dim)")
469 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
472 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
473 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
475 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(:,2), mean(:,2), var(:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG))")
476 call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
481 call disp%show(
"call setVarMeanMerged(var(:,2), mean(:,2), var(:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG))")
482 call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
485 call disp%show(
"mean(:,0) ! reference")
489 call disp%show(
"var(:,0) ! reference")
496 complex(TKG),
allocatable :: sample(:,:)
497 complex(TKG),
allocatable :: mean(:,:), meanMerged(:)
498 real(TKG),
allocatable :: var(:,:), varMerged(:)
501 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
504 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
505 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
508 lb(isam)
= ub(isam
- 1)
+ 1
515 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
517 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
518 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
519 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
520 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
521 call disp%show(
"call setResized(meanMerged, ndim)")
523 call disp%show(
"call setResized(varMerged, ndim)")
525 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ndim, ub(nsam)), -getUnifRand(-1., +1., ndim, ub(nsam)), TKG)")
529 call disp%show(
"mean(:,0) = getVar(sample, dim)")
530 mean(:,
0)
= getVar(sample, dim)
531 call disp%show(
"mean(:,0) ! reference")
533 call disp%show(
"var(:,0) = getVar(sample, dim)")
534 var(:,
0)
= getVar(sample, dim)
535 call disp%show(
"var(:,0) ! reference")
538 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), dim)")
539 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
542 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
543 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
545 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(:,2), mean(:,2), var(:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG))")
546 call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
551 call disp%show(
"call setVarMeanMerged(var(:,2), mean(:,2), var(:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG))")
552 call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
555 call disp%show(
"mean(:,0) ! reference")
559 call disp%show(
"var(:,0) ! reference")
566 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
567 call disp%show(
"!Compute the biased merged variance and mean of a frequency weighted multivariate sample.")
568 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
572 real(TKG),
allocatable :: sample(:,:)
573 real(TKG),
allocatable :: mean(:,:), meanMerged(:)
574 real(TKG),
allocatable :: var(:,:), varMerged(:)
577 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
580 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
581 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
584 lb(isam)
= ub(isam
- 1)
+ 1
591 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
593 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
594 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
595 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
596 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
597 call disp%show(
"call setResized(meanMerged, ndim)")
599 call disp%show(
"call setResized(varMerged, ndim)")
601 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
605 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
609 call disp%show(
"mean(:,0) = getVar(sample, 2_IK, iweight)")
610 mean(:,
0)
= getVar(sample,
2_IK, iweight)
611 call disp%show(
"mean(:,0) ! reference")
613 call disp%show(
"var(:,0) = getVar(sample, 2_IK, iweight)")
614 var(:,
0)
= getVar(sample,
2_IK, iweight)
615 call disp%show(
"var(:,0) ! reference")
618 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
619 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
622 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
623 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
625 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
626 call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
631 call disp%show(
"call setVarMeanMerged(var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
632 call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
635 call disp%show(
"mean(:,0) ! reference")
639 call disp%show(
"var(:,0) ! reference")
646 complex(TKG),
allocatable :: sample(:,:)
647 complex(TKG),
allocatable :: mean(:,:), meanMerged(:)
648 real(TKG),
allocatable :: var(:,:), varMerged(:)
651 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
654 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
655 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
658 lb(isam)
= ub(isam
- 1)
+ 1
665 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
667 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
668 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
669 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
670 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
671 call disp%show(
"call setResized(meanMerged, ndim)")
673 call disp%show(
"call setResized(varMerged, ndim)")
675 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ndim, ub(nsam)), getUnifRand(-1., +1., ndim, ub(nsam)), TKG)")
679 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
683 call disp%show(
"mean(:,0) = getVar(sample, 2_IK, iweight)")
684 mean(:,
0)
= getVar(sample,
2_IK, iweight)
685 call disp%show(
"mean(:,0) ! reference")
687 call disp%show(
"var(:,0) = getVar(sample, 2_IK, iweight)")
688 var(:,
0)
= getVar(sample,
2_IK, iweight)
689 call disp%show(
"var(:,0) ! reference")
692 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
693 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
696 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
697 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
699 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
700 call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
705 call disp%show(
"call setVarMeanMerged(var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
706 call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
709 call disp%show(
"mean(:,0) ! reference")
713 call disp%show(
"var(:,0) ! reference")
720 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
721 call disp%show(
"!Compute the biased merged variance and mean of a reliability weighted multivariate sample.")
722 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
726 real(TKG),
allocatable :: sample(:,:)
727 real(TKG),
allocatable :: mean(:,:), meanMerged(:)
728 real(TKG),
allocatable :: var(:,:), varMerged(:)
729 real(TKG),
allocatable :: rweight(:)
732 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
735 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
736 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
739 lb(isam)
= ub(isam
- 1)
+ 1
746 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
748 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
749 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
750 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
751 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
752 call disp%show(
"call setResized(meanMerged, ndim)")
754 call disp%show(
"call setResized(varMerged, ndim)")
756 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
760 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
764 call disp%show(
"mean(:,0) = getVar(sample, 2_IK, rweight)")
765 mean(:,
0)
= getVar(sample,
2_IK, rweight)
766 call disp%show(
"mean(:,0) ! reference")
768 call disp%show(
"var(:,0) = getVar(sample, 2_IK, rweight)")
769 var(:,
0)
= getVar(sample,
2_IK, rweight)
770 call disp%show(
"var(:,0) ! reference")
773 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
774 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
777 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
778 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
780 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
781 call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
786 call disp%show(
"call setVarMeanMerged(var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
787 call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
790 call disp%show(
"mean(:,0) ! reference")
794 call disp%show(
"var(:,0) ! reference")
801 complex(TKG),
allocatable :: sample(:,:)
802 complex(TKG),
allocatable :: mean(:,:), meanMerged(:)
803 real(TKG),
allocatable :: var(:,:), varMerged(:)
804 real(TKG),
allocatable :: rweight(:)
807 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
810 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
811 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
814 lb(isam)
= ub(isam
- 1)
+ 1
821 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
823 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
824 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
825 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
826 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
827 call disp%show(
"call setResized(meanMerged, ndim)")
829 call disp%show(
"call setResized(varMerged, ndim)")
831 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ndim, ub(nsam)), getUnifRand(-1., +1., ndim, ub(nsam)), TKG)")
835 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
839 call disp%show(
"mean(:,0) = getVar(sample, 2_IK, rweight)")
840 mean(:,
0)
= getVar(sample,
2_IK, rweight)
841 call disp%show(
"mean(:,0) ! reference")
843 call disp%show(
"var(:,0) = getVar(sample, 2_IK, rweight)")
844 var(:,
0)
= getVar(sample,
2_IK, rweight)
845 call disp%show(
"var(:,0) ! reference")
848 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
849 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
852 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
853 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
855 call disp%show(
"call setVarMeanMerged(varMerged, meanMerged, var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
856 call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
861 call disp%show(
"call setVarMeanMerged(var(:,2), mean(:,2), var(:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
862 call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
865 call disp%show(
"mean(:,0) ! reference")
869 call disp%show(
"var(:,0) ! reference")
Generate minimally-spaced character, integer, real sequences or sequences at fixed intervals of size ...
Resize (shrink or expand) an input allocatable array of rank 1..3 to arbitrary new lower and upper bo...
Allocate or resize (shrink or expand) an input allocatable scalar string or array of rank 1....
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
This is a generic method of the derived type display_type with pass attribute.
This is a generic method of the derived type display_type with pass attribute.
Generate and return the (weighted) mean of an input sample of nsam observations with ndim = 1 or 2 at...
Generate and return the variance of the input sample of type complex or real of shape (nsam) or (ndim...
This module contains procedures and generic interfaces for generating ranges of discrete character,...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for input/output (IO) or generic display operations on st...
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
This module contains classes and procedures for computing the first moment (i.e., the statistical mea...
Generate and return an object of type display_type.
9 lb(isam)
= ub(isam
- 1)
+ 1
18+0.119940042E-1,
+0.656136096,
+0.920551598,
+0.100176871,
+0.860603511,
+0.201040685,
+0.459808290,
+0.532367587,
+0.672253072,
+0.793510675E-2,
+0.413463652
26 var(isam)
= getVar(sample(lb(isam):ub(isam)))
27 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
29call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
34call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
47 lb(isam)
= ub(isam
- 1)
+ 1
56+0.883841753,
+0.354862809E-1,
+0.707520783,
+0.353002489,
+0.971128941E-1,
+0.679025769,
+0.887165725,
+0.346522450,
+0.944431603,
+0.259654522E-1,
+0.854691327
64 var(isam)
= getVar(sample(lb(isam):ub(isam)))
65 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
67call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
72call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
85 lb(isam)
= ub(isam
- 1)
+ 1
94+0.133238673,
+0.704561889,
+0.235955894,
+0.434056163,
+0.439929366,
+0.597409070,
+0.800524533,
+0.882233262,
+0.512249291,
+0.932722628
102 var(isam)
= getVar(sample(lb(isam):ub(isam)))
103 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
105call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
110call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
123 lb(isam)
= ub(isam
- 1)
+ 1
132+0.311637640,
+0.286906958,
+0.195507348,
+0.811418474,
+0.552980423,
+0.104915917,
+0.750227213,
+0.781923532E-1,
+0.117942452,
+0.849288762
140 var(isam)
= getVar(sample(lb(isam):ub(isam)))
141 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
143call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
148call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
161 lb(isam)
= ub(isam
- 1)
+ 1
170+0.769666314,
+0.289621472,
+0.782481730,
+0.229598761,
+0.572728693,
+0.294855356,
+0.365101993,
+0.518420577
178 var(isam)
= getVar(sample(lb(isam):ub(isam)))
179 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
181call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
186call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
199 lb(isam)
= ub(isam
- 1)
+ 1
208(
-0.385736227E-1,
-0.971312165), (
-0.164464593,
+0.457289696), (
-0.552359462,
-0.583754778), (
-0.244892836E-1,
-0.118603349), (
+0.934637308,
+0.229700208), (
-0.168417931,
+0.209864140), (
-0.785303235,
-0.622492790), (
-0.141442776,
-0.474508405), (
+0.145599246,
+0.854253769E-3)
211(
-0.883127078E-1,
-0.208107024)
216 var(isam)
= getVar(sample(lb(isam):ub(isam)))
217 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
219call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
221(
-0.883127004E-1,
-0.208107024)
224call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
226(
-0.883127004E-1,
-0.208107024)
228(
-0.883127078E-1,
-0.208107024)
237 lb(isam)
= ub(isam
- 1)
+ 1
246(
-0.645243645,
+0.117631197), (
-0.365102530,
-0.439955711), (
+0.737687469,
+0.750073314), (
-0.790288687,
-0.346800208), (
-0.858775258,
-0.964322329), (
+0.373319030,
+0.493835330), (
+0.363248467,
-0.768578768), (
-0.687694550,
-0.548529267), (
+0.472329974,
+0.716293812), (
-0.322200656,
+0.764347076), (
-0.346867085,
+0.717595339), (
+0.556713343,
-0.145169020)
249(
-0.126072839,
+0.288683977E-1)
254 var(isam)
= getVar(sample(lb(isam):ub(isam)))
255 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
257call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
259(
-0.126072854,
+0.288683958E-1)
262call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
264(
-0.126072854,
+0.288683958E-1)
266(
-0.126072839,
+0.288683977E-1)
275 lb(isam)
= ub(isam
- 1)
+ 1
284(
-0.463739991,
-0.740990281), (
-0.887680292,
+0.577617884E-1), (
-0.696246147,
-0.493150949), (
+0.914370537,
+0.327253819), (
-0.237545729,
-0.638078570), (
+0.699199677,
-0.951303482)
287(
-0.111940347,
-0.406417936)
292 var(isam)
= getVar(sample(lb(isam):ub(isam)))
293 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
295call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
297(
-0.111940339,
-0.406417936)
300call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
302(
-0.111940339,
-0.406417936)
304(
-0.111940347,
-0.406417936)
313 lb(isam)
= ub(isam
- 1)
+ 1
322(
+0.271996379,
-0.778391004), (
+0.362214565,
-0.723257422), (
-0.615412474,
-0.175877213), (
-0.438934922,
-0.869427085), (
-0.199115157,
-0.280784488), (
-0.546283126,
-0.614637136E-1), (
-0.222089887,
+0.420297980), (
-0.617734909,
+0.170706987), (
-0.783601522,
-0.256677985), (
-0.881325006E-1,
-0.971402764)
325(
-0.287709385,
-0.352627665)
330 var(isam)
= getVar(sample(lb(isam):ub(isam)))
331 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
333call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
335(
-0.287709385,
-0.352627695)
338call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
340(
-0.287709385,
-0.352627695)
342(
-0.287709385,
-0.352627665)
351 lb(isam)
= ub(isam
- 1)
+ 1
360(
+0.230885625,
-0.781304359), (
-0.107927680,
+0.313308120), (
-0.275286317,
+0.807674170), (
-0.492192030,
+0.894426107E-1), (
+0.900775075,
-0.742510319), (
+0.691864014,
+0.187800765), (
-0.853224277,
-0.652338982)
363(
+0.135563444E-1,
-0.111132570)
368 var(isam)
= getVar(sample(lb(isam):ub(isam)))
369 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
371call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
373(
+0.135563482E-1,
-0.111132562)
376call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1), ub(
1)
/ real(ub(
2), TKG))
378(
+0.135563482E-1,
-0.111132562)
380(
+0.135563444E-1,
-0.111132570)
394 lb(isam)
= ub(isam
- 1)
+ 1
403+0.877293646,
+0.872215331,
+0.125901878,
+0.301837325,
+0.908573806,
+0.414544940,
+0.434481144,
+0.148264587,
+0.611916482,
+0.848079145
406+3,
+1,
+5,
+3,
+7,
+9,
+4,
+3,
+6,
+7
407mean(
0)
= getMean(sample, iweight)
410var(
0)
= getVar(sample, iweight)
414 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
415 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
417call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
422call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
435 lb(isam)
= ub(isam
- 1)
+ 1
444+0.700306714,
+0.292092562,
+0.974595904,
+0.867057383,
+0.441709697,
+0.668324947,
+0.509613633,
+0.879361153
447+2,
+5,
+8,
+1,
+10,
+7,
+1,
+4
448mean(
0)
= getMean(sample, iweight)
451var(
0)
= getVar(sample, iweight)
455 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
456 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
458call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
463call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
476 lb(isam)
= ub(isam
- 1)
+ 1
485+0.958505034,
+0.345331252,
+0.319890797,
+0.160288930,
+0.762265980,
+0.535349429,
+0.990700126E-1,
+0.201173484,
+0.212709844,
+0.776189625,
+0.778449178,
+0.504525006,
+0.270502448,
+0.664611638
488+1,
+8,
+4,
+4,
+3,
+9,
+6,
+8,
+6,
+8,
+5,
+3,
+2,
+2
489mean(
0)
= getMean(sample, iweight)
492var(
0)
= getVar(sample, iweight)
496 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
497 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
499call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
504call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
517 lb(isam)
= ub(isam
- 1)
+ 1
526+0.812932909,
+0.840605795,
+0.753042042,
+0.138438940,
+0.214533925,
+0.625773966,
+0.594727159,
+0.524592102,
+0.151322603,
+0.410685182,
+0.898426771E-1,
+0.861285269
529+8,
+10,
+8,
+3,
+9,
+1,
+9,
+10,
+1,
+8,
+2,
+9
530mean(
0)
= getMean(sample, iweight)
533var(
0)
= getVar(sample, iweight)
537 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
538 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
540call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
545call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
558 lb(isam)
= ub(isam
- 1)
+ 1
567+0.450314939,
+0.303190827,
+0.667677045,
+0.492998481,
+0.766684115,
+0.107900679,
+0.386415839,
+0.668759227,
+0.317474127,
+0.284253597,
+0.542022586,
+0.207650661E-1,
+0.941341877
570+6,
+1,
+6,
+1,
+1,
+5,
+9,
+9,
+1,
+5,
+2,
+6,
+5
571mean(
0)
= getMean(sample, iweight)
574var(
0)
= getVar(sample, iweight)
578 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
579 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
581call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
586call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
599 lb(isam)
= ub(isam
- 1)
+ 1
608(
+0.733839154,
-0.736027062), (
+0.238857985,
-0.570147872), (
+0.161682308,
-0.531322896), (
+0.824257672,
-0.677505434), (
+0.780953765E-1,
-0.651413560), (
+0.426721990,
-0.748082161), (
+0.683678210,
-0.855036020), (
+0.214939833,
-0.527190447), (
+0.978096664,
-0.140259802)
611+5,
+1,
+8,
+6,
+9,
+4,
+1,
+3,
+6
612mean(
0)
= getMean(sample, iweight)
614(
+0.459392637,
-0.574398041)
615var(
0)
= getVar(sample, iweight)
619 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
620 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
622call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
624(
+0.459392577,
-0.574398041)
627call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
629(
+0.459392577,
-0.574398041)
631(
+0.459392637,
-0.574398041)
640 lb(isam)
= ub(isam
- 1)
+ 1
649(
+0.953109980,
-0.314150333), (
+0.780628860,
-0.188690245), (
+0.485206187,
-0.691465676), (
+0.236192822,
-0.561198115), (
+0.709286392,
-0.338161469), (
+0.389771342,
-0.542102456), (
+0.531400800,
-0.437221527E-1), (
+0.149240315,
-0.982882679), (
+0.512519538,
-0.503458142), (
+0.533620477,
-0.602279365), (
+0.660914183,
-0.493254721)
652+4,
+5,
+5,
+4,
+2,
+9,
+8,
+8,
+9,
+5,
+4
653mean(
0)
= getMean(sample, iweight)
655(
+0.498134375,
-0.495012641)
656var(
0)
= getVar(sample, iweight)
660 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
661 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
663call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
665(
+0.498134375,
-0.495012641)
668call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
670(
+0.498134375,
-0.495012641)
672(
+0.498134375,
-0.495012641)
681 lb(isam)
= ub(isam
- 1)
+ 1
690(
+0.608615577,
-0.418816626), (
+0.619734049,
-0.964758635), (
+0.437514067,
-0.212158978), (
+0.248523414,
-0.206597209), (
+0.731221914,
-0.797430575)
694mean(
0)
= getMean(sample, iweight)
696(
+0.534081519,
-0.469987184)
697var(
0)
= getVar(sample, iweight)
701 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
702 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
704call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
706(
+0.534081519,
-0.469987184)
709call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
711(
+0.534081519,
-0.469987184)
713(
+0.534081519,
-0.469987184)
722 lb(isam)
= ub(isam
- 1)
+ 1
731(
+0.288896084,
-0.653880894), (
+0.387549162,
-0.979286194), (
+0.399118066,
-0.748419583), (
+0.423218250,
-0.312043726), (
+0.623215854,
-0.630982041), (
+0.744935870E-1,
-0.158288479E-1), (
+0.171287358,
-0.475338459), (
+0.967443705,
-0.975493848), (
+0.253366828E-1,
-0.505101979), (
+0.576567054E-1,
-0.335952044), (
+0.344425440E-2,
-0.901057124E-1)
734+9,
+6,
+7,
+3,
+9,
+2,
+7,
+2,
+2,
+3,
+3
735mean(
0)
= getMean(sample, iweight)
737(
+0.341785759,
-0.588923454)
738var(
0)
= getVar(sample, iweight)
742 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
743 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
745call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
747(
+0.341785759,
-0.588923454)
750call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
752(
+0.341785759,
-0.588923454)
754(
+0.341785759,
-0.588923454)
763 lb(isam)
= ub(isam
- 1)
+ 1
772(
+0.657687306,
-0.670729876E-1), (
+0.958387852,
-0.722694397E-2), (
+0.988820493,
-0.263041496), (
+0.481504560,
-0.756367266), (
+0.410300493,
-0.245879471), (
+0.795342386,
-0.992211699), (
+0.929423094,
-0.182306767E-1), (
+0.234219849,
-0.333952785)
775+7,
+10,
+9,
+5,
+7,
+3,
+7,
+2
776mean(
0)
= getMean(sample, iweight)
778(
+0.754542530,
-0.243686050)
779var(
0)
= getVar(sample, iweight)
783 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
784 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
786call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
788(
+0.754542530,
-0.243686050)
791call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
793(
+0.754542530,
-0.243686050)
795(
+0.754542530,
-0.243686050)
809 lb(isam)
= ub(isam
- 1)
+ 1
818+0.674505651,
+0.884572268E-1,
+0.778975904,
+0.651336432,
+0.980571866,
+0.141536474,
+0.291891158,
+0.321805716,
+0.400037169
821+1.75373650,
+1.03356385,
+1.52788496,
+1.01779795,
+1.20727324,
+1.60602701,
+1.30024660,
+1.08889067,
+1.73037934
822mean(
0)
= getMean(sample, rweight)
825var(
0)
= getVar(sample, rweight)
829 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
830 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
832call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
837call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
850 lb(isam)
= ub(isam
- 1)
+ 1
859+0.614780426,
+0.470804274,
+0.192424715,
+0.478024364,
+0.483757198,
+0.961663127E-1,
+0.713184476E-1,
+0.392000079,
+0.621323168,
+0.438278556,
+0.204148531
862+1.01158857,
+1.01859975,
+1.73260868,
+1.48818684,
+1.35715485,
+1.78745222,
+1.86245263,
+1.75748146,
+1.96912062,
+1.42044401,
+1.42447019
863mean(
0)
= getMean(sample, rweight)
866var(
0)
= getVar(sample, rweight)
870 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
871 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
873call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
878call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
891 lb(isam)
= ub(isam
- 1)
+ 1
900+0.209337831,
+0.230150700,
+0.859026134,
+0.810046315,
+0.476253629,
+0.970555782,
+0.960021973,
+0.910234153,
+0.870920360,
+0.699376225,
+0.661750972,
+0.782559395
903+1.94758272,
+1.76901412,
+1.87601042,
+1.63177502,
+1.02643919,
+1.19849062,
+1.03706694,
+1.53485990,
+1.15975714,
+1.80401039,
+1.67372990,
+1.41304803
904mean(
0)
= getMean(sample, rweight)
907var(
0)
= getVar(sample, rweight)
911 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
912 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
914call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
919call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
932 lb(isam)
= ub(isam
- 1)
+ 1
941+0.331320763E-1,
+0.997135401,
+0.416052222,
+0.908958912E-2,
+0.999018729,
+0.264216542,
+0.818529725E-1,
+0.427850306,
+0.120869100,
+0.107904911,
+0.470962107
944+1.49814856,
+1.64824224,
+1.89569652,
+1.23003721,
+1.67444062,
+1.46046174,
+1.02513421,
+1.33302343,
+1.79053736,
+1.49233651,
+1.19293594
945mean(
0)
= getMean(sample, rweight)
948var(
0)
= getVar(sample, rweight)
952 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
953 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
955call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
960call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
973 lb(isam)
= ub(isam
- 1)
+ 1
982+0.479176342,
+0.651354790,
+0.979081631,
+0.643035412,
+0.646971166,
+0.728912175,
+0.733106077,
+0.726554751,
+0.807007492,
+0.740020156,
+0.853646338
985+1.01996756,
+1.20847034,
+1.86048603,
+1.59032822,
+1.16948724,
+1.67066908,
+1.58738041,
+1.76641357,
+1.98581839,
+1.71012497,
+1.34003353
986mean(
0)
= getMean(sample, rweight)
989var(
0)
= getVar(sample, rweight)
993 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
994 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
996call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1001call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1014 lb(isam)
= ub(isam
- 1)
+ 1
1023(
+0.689120829,
-0.234999061), (
+0.402820230,
-0.445851266), (
+0.928521395,
-0.256282091E-1), (
+0.458933353,
-0.315659881), (
+0.141881645,
-0.894450724), (
+0.667115092,
-0.426235199)
1026+1.25882101,
+1.01175332,
+1.05331218,
+1.30983901,
+1.80418277,
+1.32992387
1027mean(
0)
= getMean(sample, rweight)
1029(
+0.514606774,
-0.433581024)
1030var(
0)
= getVar(sample, rweight)
1034 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1035 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1037call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1039(
+0.514606833,
-0.433580995)
1042call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1044(
+0.514606833,
-0.433580995)
1046(
+0.514606774,
-0.433581024)
1055 lb(isam)
= ub(isam
- 1)
+ 1
1064(
+0.365207911,
-0.388750851), (
+0.475437462,
-0.850406229), (
+0.600203335,
-0.738627851), (
+0.238107026,
-0.970534921), (
+0.744269729,
-0.638918519), (
+0.544522703,
-0.445047200), (
+0.359141469,
-0.507187247)
1067+1.51620340,
+1.85358036,
+1.72813451,
+1.21773052,
+1.07183599,
+1.50133920,
+1.03417802
1068mean(
0)
= getMean(sample, rweight)
1070(
+0.478568673,
-0.655197263)
1071var(
0)
= getVar(sample, rweight)
1075 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1076 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1078call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1080(
+0.478568703,
-0.655197382)
1083call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1085(
+0.478568703,
-0.655197382)
1087(
+0.478568673,
-0.655197263)
1096 lb(isam)
= ub(isam
- 1)
+ 1
1105(
+0.739192128,
-0.842893600), (
+0.557888150E-1,
-0.790932655), (
+0.493421853,
-0.256256402), (
+0.133751631E-1,
-0.820576131)
1108+1.61477947,
+1.66125584,
+1.36738837,
+1.45803833
1109mean(
0)
= getMean(sample, rweight)
1111(
+0.324596375,
-0.691942751)
1112var(
0)
= getVar(sample, rweight)
1116 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1117 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1119call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1121(
+0.324596405,
-0.691942692)
1124call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1126(
+0.324596405,
-0.691942692)
1128(
+0.324596375,
-0.691942751)
1137 lb(isam)
= ub(isam
- 1)
+ 1
1146(
+0.778856456,
-0.643124104), (
+0.929328263,
-0.953751445), (
+0.306445360,
-0.764832199), (
+0.653658748,
-0.180498183), (
+0.167709410,
-0.995094776), (
+0.274486303,
-0.459129333), (
+0.697831631,
-0.748343527), (
+0.999758244,
-0.207869589), (
+0.638001621,
-0.797863543)
1149+1.46506858,
+1.14070487,
+1.34108961,
+1.08089876,
+1.27528250,
+1.24871540,
+1.94184399,
+1.76558363,
+1.08629298
1150mean(
0)
= getMean(sample, rweight)
1152(
+0.622786820,
-0.630207837)
1153var(
0)
= getVar(sample, rweight)
1157 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1158 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1160call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1162(
+0.622786880,
-0.630207837)
1165call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1167(
+0.622786880,
-0.630207837)
1169(
+0.622786820,
-0.630207837)
1178 lb(isam)
= ub(isam
- 1)
+ 1
1187(
+0.498375833,
-0.341944516), (
+0.785785675,
-0.838963985), (
+0.803011954,
-0.905782282), (
+0.308905840,
-0.354638219), (
+0.672143102,
-0.715029776), (
+0.299706280,
-0.650145113), (
+0.229217350,
-0.692987084), (
+0.968931794,
-0.803514123), (
+0.308839679E-1,
-0.159259737), (
+0.426950395,
-0.564006925)
1190+1.24266434,
+1.08515835,
+1.61370671,
+1.62555397,
+1.66845918,
+1.37063169,
+1.32197857,
+1.62536025,
+1.23835897,
+1.98431218
1191mean(
0)
= getMean(sample, rweight)
1193(
+0.512009501,
-0.608823538)
1194var(
0)
= getVar(sample, rweight)
1198 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1199 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
1201call setVarMeanMerged(varMerged, meanMerged, var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1203(
+0.512009561,
-0.608823657)
1206call setVarMeanMerged(var(
2), mean(
2), var(
1), mean(
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1208(
+0.512009561,
-0.608823657)
1210(
+0.512009501,
-0.608823538)
1224 lb(isam)
= ub(isam
- 1)
+ 1
1232call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1233call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1238-0.439300776,
-0.947754502,
+0.633023024,
-0.410256982,
+0.196280479E-1,
+0.104475260,
+0.199488044
1239+0.713191628,
-0.763127685,
+0.995138764,
-0.165531993,
+0.520690560,
-0.845351934,
-0.545616984
1240-0.311024070,
+0.587350607,
-0.697266936,
+0.609308243,
+0.846387744,
-0.507608891,
-0.164826751
1241mean(:,
0)
= getVar(sample, dim)
1243+0.230054438,
+0.484424531,
+0.325943947
1244var(:,
0)
= getVar(sample, dim)
1246+0.230054438,
+0.484424531,
+0.325943947
1248 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1249 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1251call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1253-0.120099694,
-0.129439235E-1,
+0.517599843E-1
1255+0.230054438,
+0.484424531,
+0.325943947
1256call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1258-0.120099694,
-0.129439235E-1,
+0.517599843E-1
1260+0.230054438,
+0.484424531,
+0.325943947
1262+0.230054438,
+0.484424531,
+0.325943947
1264+0.230054438,
+0.484424531,
+0.325943947
1269 lb(isam)
= ub(isam
- 1)
+ 1
1277call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1278call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1283-0.320945621,
-0.938246965,
+0.846944928,
+0.792569995,
+0.607313871,
-0.270295858,
+0.727086663,
+0.571539164,
-0.546037793,
+0.907956958,
-0.270005107,
+0.442611098,
+0.785134315,
-0.528756261
1284-0.502276540,
+0.382720113,
-0.453468084,
-0.497294307,
-0.457640767,
+0.142306089,
-0.484768867,
+0.934381008,
+0.372496843,
+0.547003627,
-0.256724358E-1,
+0.137520432,
-0.172556996,
+0.629006028
1285-0.250745058,
+0.875776172,
-0.220005393,
-0.611425161,
-0.606422901,
+0.772572875,
-0.392596841,
+0.325260401,
+0.205267668E-1,
-0.626538157,
-0.697822332,
-0.849208832E-1,
+0.591000795,
+0.219418526
1286-0.803479671,
-0.860755444,
+0.858319998,
-0.183054447,
-0.511176705,
+0.575032234,
+0.788634658,
-0.472276092,
+0.826205134,
+0.577724218,
+0.620198846,
+0.164198637,
+0.447374940,
-0.961143017
1287+0.495147109,
+0.682114124,
+0.967505336,
-0.212088943,
-0.533003807,
-0.331007957,
+0.313996315,
-0.430609703,
-0.972972035,
-0.365960240,
-0.638307333,
+0.685997009E-1,
-0.162395000,
+0.450876117
1288mean(:,
0)
= getVar(sample, dim)
1290+0.382186949,
+0.217904523,
+0.266257435,
+0.433174640,
+0.289978445
1291var(:,
0)
= getVar(sample, dim)
1293+0.382186949,
+0.217904523,
+0.266257435,
+0.433174640,
+0.289978445
1295 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1296 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1298call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1300+0.200490683,
+0.394111574E-1,
-0.489943698E-1,
+0.761288106E-1,
-0.477218777E-1
1302+0.382186979,
+0.217904538,
+0.266257405,
+0.433174670,
+0.289978445
1303call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1305+0.200490683,
+0.394111574E-1,
-0.489943698E-1,
+0.761288106E-1,
-0.477218777E-1
1307+0.382186949,
+0.217904523,
+0.266257435,
+0.433174640,
+0.289978445
1309+0.382186979,
+0.217904538,
+0.266257405,
+0.433174670,
+0.289978445
1311+0.382186949,
+0.217904523,
+0.266257435,
+0.433174640,
+0.289978445
1316 lb(isam)
= ub(isam
- 1)
+ 1
1324call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1325call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1330-0.158958673,
-0.233285904,
+0.215965033,
+0.425416350,
+0.384628773E-1,
+0.931401491,
+0.661632657,
-0.852002025,
-0.514157057,
-0.905305386,
+0.842246532,
+0.581661105,
-0.678930759,
+0.708365679
1331+0.612628102,
-0.732807994,
+0.537922740,
+0.945479155,
-0.975803494,
-0.592299104,
+0.775635004,
-0.655461788,
+0.231520295,
+0.361016989E-1,
+0.717669725E-1,
-0.547278643,
+0.535398364,
-0.925517559
1332-0.776354194,
+0.749330759,
-0.920078278,
+0.378347516,
-0.644398093,
-0.191198587E-1,
-0.209176540E-1,
+0.346883535,
-0.817781091,
-0.427559733,
+0.572685242,
-0.118050814,
+0.513249993,
-0.196097612
1333-0.583040833,
+0.264287472,
+0.673619270,
+0.223239183,
+0.348567963E-3,
+0.741120100,
+0.108615398,
-0.387412429,
+0.531768322,
-0.454299450E-1,
+0.911919117,
-0.243890882,
-0.744576693,
-0.913441896
1334+0.858712792,
+0.453865528E-2,
+0.789859772,
-0.710540056,
+0.407098532,
+0.545049787,
+0.747678757,
+0.262606144E-1,
+0.601000190,
+0.868038774,
+0.895237327,
-0.853961110,
+0.547324777,
-0.596507549
1335+0.960217714,
+0.204936266E-1,
-0.542162061,
+0.990568399,
-0.520812154,
+0.102156758,
-0.196271896,
+0.140751600,
+0.674492240,
-0.441921473,
+0.668405175,
-0.731903791,
-0.227602005,
-0.443277359
1336mean(:,
0)
= getVar(sample, dim)
1338+0.382407039,
+0.421383917,
+0.289260745,
+0.297758013,
+0.355087638,
+0.314566612
1339var(:,
0)
= getVar(sample, dim)
1341+0.382407039,
+0.421383917,
+0.289260745,
+0.297758013,
+0.355087638,
+0.314566612
1343 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1344 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1346call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1348+0.758937150E-1,
-0.487654470E-1,
-0.985614508E-1,
+0.383660495E-1,
+0.294985086,
+0.323667750E-1
1350+0.382407039,
+0.421383798,
+0.289260745,
+0.297758013,
+0.355087638,
+0.314566642
1351call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1353+0.758937150E-1,
-0.487654470E-1,
-0.985614508E-1,
+0.383660495E-1,
+0.294985086,
+0.323667750E-1
1355+0.382407039,
+0.421383917,
+0.289260745,
+0.297758013,
+0.355087638,
+0.314566612
1357+0.382407039,
+0.421383798,
+0.289260745,
+0.297758013,
+0.355087638,
+0.314566642
1359+0.382407039,
+0.421383917,
+0.289260745,
+0.297758013,
+0.355087638,
+0.314566612
1364 lb(isam)
= ub(isam
- 1)
+ 1
1372call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1373call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1378-0.150514126,
+0.711577535,
+0.216513395,
-0.497512102,
-0.351622343,
+0.579126358,
-0.851261616,
-0.146335959,
+0.136786580
1379+0.482598782,
+0.549315572,
+0.190539002,
+0.628557324,
-0.234736681,
+0.438307881,
+0.216528893,
+0.240388513,
-0.900430679E-1
1380mean(:,
0)
= getVar(sample, dim)
1382+0.225925073,
+0.749487057E-1
1383var(:,
0)
= getVar(sample, dim)
1385+0.225925073,
+0.749487057E-1
1387 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1388 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1390call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1392-0.392491333E-1,
+0.269050717
1394+0.225925088,
+0.749487132E-1
1395call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1397-0.392491333E-1,
+0.269050717
1399+0.225925073,
+0.749487057E-1
1401+0.225925088,
+0.749487132E-1
1403+0.225925073,
+0.749487057E-1
1408 lb(isam)
= ub(isam
- 1)
+ 1
1416call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1417call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1422-0.416686893,
+0.125025034,
+0.971437693E-1,
+0.872558594,
-0.932730079,
-0.610853791,
-0.327465177,
+0.807903886,
-0.467597485,
-0.484535336,
-0.489243031,
+0.680097699
1423-0.712541223,
-0.225304365E-1,
-0.437398672,
+0.991033673,
+0.837484241,
-0.442648530,
-0.636169076,
+0.866333365,
-0.572261095,
+0.764882565E-1,
+0.776366115,
+0.119745493
1424-0.567359447,
+0.674515128,
+0.717805386,
-0.231692195,
-0.158733845,
+0.956128836E-1,
+0.611454248,
-0.425406694E-1,
+0.107408643,
+0.727708817,
-0.991427898E-1,
+0.138234138
1425mean(:,
0)
= getVar(sample, dim)
1427+0.334076822,
+0.385443270,
+0.166780770
1428var(:,
0)
= getVar(sample, dim)
1430+0.334076822,
+0.385443270,
+0.166780770
1432 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1433 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1435call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1437-0.955319032E-1,
+0.703251809E-1,
+0.164439201
1439+0.334076852,
+0.385443270,
+0.166780785
1440call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1442-0.955319032E-1,
+0.703251809E-1,
+0.164439201
1444+0.334076822,
+0.385443270,
+0.166780770
1446+0.334076852,
+0.385443270,
+0.166780785
1448+0.334076822,
+0.385443270,
+0.166780770
1453 lb(isam)
= ub(isam
- 1)
+ 1
1461call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1462call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1467(
-0.827937245,
+0.573437333), (
+0.925697446,
+0.444543123), (
-0.545472503,
+0.460077405), (
+0.542983532,
+0.839096665), (
-0.698645949,
-0.978960872), (
-0.829718947,
-0.274715185), (
+0.344073534,
+0.619020820), (
-0.545555353,
-0.879950881), (
+0.922487855,
-0.102486610E-1), (
-0.782257557,
+0.229459167), (
-0.352678061,
+0.930501103), (
+0.945703983E-1,
-0.688433766)
1468(
-0.229400754,
-0.530421972), (
-0.187813878,
-0.410598278), (
-0.578216076,
+0.802461863), (
-0.357069135,
-0.391425610), (
-0.215777159E-1,
+0.666727185), (
+0.110720992,
+0.697784662), (
-0.652069211,
+0.738842010), (
+0.602635145,
+0.525061369), (
+0.671922803,
-0.732921004), (
+0.378312230,
+0.658055305), (
-0.484220505,
-0.505790710E-1), (
+0.102074623,
+0.438763976)
1469(
-0.753292322,
+0.666397333), (
+0.759135485,
-0.437500358), (
+0.327015281,
-0.366309047), (
-0.800070524,
-0.372627974E-1), (
-0.243941903,
+0.702132821), (
-0.195869327,
+0.763610005), (
+0.675318122,
+0.841464162), (
-0.962815523,
-0.837390423), (
-0.657793403,
+0.451066613), (
+0.109459758,
-0.885191917), (
-0.926375985,
+0.122865438), (
+0.521078110E-1,
-0.803279042)
1470(
+0.856962085,
-0.500666499), (
-0.539519429,
+0.220919132), (
+0.658990622,
+0.426123738), (
-0.418452263,
+0.682990313), (
-0.398637891,
-0.402785659), (
-0.780998707,
-0.519223452), (
-0.633221865,
+0.487448215), (
+0.831326365,
-0.270884156), (
+0.986258864,
+0.271868110), (
+0.192654490,
-0.669101238), (
+0.743857861,
-0.871444941), (
+0.926627398,
-0.228209496)
1471(
+0.605386257,
+0.778706074), (
-0.909575820,
+0.703624725), (
-0.902975798E-1,
+0.787807584), (
-0.952966690,
+0.398319840), (
+0.708492041,
-0.832051992), (
-0.817630649,
+0.484323502E-1), (
+0.566180348,
+0.942621827), (
-0.115485191E-1,
+0.641840696), (
-0.673574328,
-0.224295139), (
+0.938803673,
+0.584035873), (
+0.369288564,
-0.372304797), (
-0.424251437,
-0.685607910)
1472mean(:,
0)
= getVar(sample, dim)
1474(
+0.830049872,
+0.00000000), (
+0.486362100,
+0.00000000), (
+0.750634313,
+0.00000000), (
+0.691127717,
+0.00000000), (
+0.787409186,
+0.00000000)
1475var(:,
0)
= getVar(sample, dim)
1477+0.830049872,
+0.486362100,
+0.750634313,
+0.691127717,
+0.787409186
1479 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1480 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1482call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1484(
-0.146037757,
+0.105318844), (
-0.537251234E-1,
+0.200979233), (
-0.218093559,
+0.150502305E-1), (
+0.202153966,
-0.114413843), (
-0.576411784E-1,
+0.230927408)
1486+0.830049813,
+0.486362129,
+0.750634313,
+0.691127717,
+0.787409186
1487call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1489(
-0.146037757,
+0.105318844), (
-0.537251234E-1,
+0.200979233), (
-0.218093559,
+0.150502305E-1), (
+0.202153966,
-0.114413843), (
-0.576411784E-1,
+0.230927408)
1491(
+0.830049872,
+0.00000000), (
+0.486362100,
+0.00000000), (
+0.750634313,
+0.00000000), (
+0.691127717,
+0.00000000), (
+0.787409186,
+0.00000000)
1493+0.830049813,
+0.486362129,
+0.750634313,
+0.691127717,
+0.787409186
1495+0.830049872,
+0.486362100,
+0.750634313,
+0.691127717,
+0.787409186
1500 lb(isam)
= ub(isam
- 1)
+ 1
1508call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1509call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1514(
+0.763936281,
-0.335784316), (
-0.886662841,
-0.659366846), (
+0.536572099,
-0.479399920), (
-0.999518752,
+0.110117078), (
-0.633312106,
+0.477422714), (
+0.800188899,
+0.630500555), (
-0.910677910,
+0.722780228), (
+0.152210832,
-0.267416716)
1515(
+0.735235929,
+0.807211280), (
-0.207703352,
-0.672501922), (
-0.972353816,
+0.939628839), (
+0.375396967,
-0.284383893), (
+0.754493475,
-0.344394684), (
-0.274041772,
+0.834401011), (
+0.427717209,
-0.783475161), (
-0.821393967,
+0.604826927)
1516(
+0.711423874,
-0.476529598), (
-0.962325215,
+0.443145037), (
+0.319249392,
-0.112603664), (
+0.697380662,
-0.562021971), (
-0.123974681,
+0.795188785), (
-0.102982521E-1,
+0.909578085), (
-0.118997455,
+0.215515137), (
-0.345981956,
-0.839719772E-1)
1517mean(:,
0)
= getVar(sample, dim)
1519(
+0.797658920,
+0.00000000), (
+0.860319376,
+0.00000000), (
+0.534044802,
+0.00000000)
1520var(:,
0)
= getVar(sample, dim)
1522+0.797658920,
+0.860319376,
+0.534044802
1524 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1525 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1527call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1529(
-0.147157937,
+0.248565972E-1), (
+0.216883421E-2,
+0.137664050), (
+0.208095461E-1,
+0.141037494)
1531+0.797658980,
+0.860319436,
+0.534044802
1532call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1534(
-0.147157937,
+0.248565972E-1), (
+0.216883421E-2,
+0.137664050), (
+0.208095461E-1,
+0.141037494)
1536(
+0.797658920,
+0.00000000), (
+0.860319376,
+0.00000000), (
+0.534044802,
+0.00000000)
1538+0.797658980,
+0.860319436,
+0.534044802
1540+0.797658920,
+0.860319376,
+0.534044802
1545 lb(isam)
= ub(isam
- 1)
+ 1
1553call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1554call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1559(
-0.836043477,
-0.367998004), (
-0.479369164,
+0.980343938), (
+0.860828400,
-0.351111531), (
+0.929677367,
-0.589284778), (
+0.497790575E-1,
+0.316970348E-1), (
+0.678704977,
+0.514211655), (
+0.705893397,
-0.767410636), (
+0.400087714,
+0.483188987), (
-0.247502208,
+0.870597243)
1560(
-0.887015700,
-0.778934598), (
+0.147802830,
-0.194910169), (
+0.853671670,
-0.390640378), (
+0.582563281,
+0.289507151), (
+0.364023924,
-0.667358041), (
-0.245419741E-1,
-0.967603207), (
-0.599232554,
+0.612211585), (
+0.683482051,
+0.394176841), (
+0.821595907,
-0.396433115)
1561mean(:,
0)
= getVar(sample, dim)
1563(
+0.731698155,
+0.00000000), (
+0.616895556,
+0.00000000)
1564var(:,
0)
= getVar(sample, dim)
1566+0.731698155,
+0.616895556
1568 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1569 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1571call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1573(
+0.229117364,
+0.893593282E-1), (
+0.215816647,
-0.233331561)
1575+0.731698155,
+0.616895616
1576call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1578(
+0.229117364,
+0.893593282E-1), (
+0.215816647,
-0.233331561)
1580(
+0.731698155,
+0.00000000), (
+0.616895556,
+0.00000000)
1582+0.731698155,
+0.616895616
1584+0.731698155,
+0.616895556
1589 lb(isam)
= ub(isam
- 1)
+ 1
1597call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1598call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1603(
+0.760060430,
-0.140191555), (
-0.665151834,
+0.211441636), (
-0.727792501,
-0.572162271), (
-0.764912724,
-0.163124800E-1), (
-0.475780606,
+0.635030746), (
+0.545262575,
-0.845309258), (
-0.168730378,
+0.390247345), (
+0.591508627,
-0.761119723), (
+0.592610598,
-0.199554801), (
+0.147854209,
-0.260728717), (
+0.186126947,
-0.711760640)
1604(
+0.454974413,
+0.391404867), (
+0.594176054,
-0.802765131), (
-0.490337610,
+0.853772163), (
-0.951193810,
+0.395136595), (
+0.249110699,
+0.922598958), (
+0.637888789,
-0.498852849), (
-0.768915415,
+0.138136148E-1), (
-0.730981469,
+0.457081318), (
-0.630509853E-1,
+0.943181992), (
-0.674258828,
+0.958418846), (
+0.367269635,
-0.575836539)
1605mean(:,
0)
= getVar(sample, dim)
1607(
+0.530160487,
+0.00000000), (
+0.727527142,
+0.00000000)
1608var(:,
0)
= getVar(sample, dim)
1610+0.530160487,
+0.727527142
1612 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1613 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1615call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1617(
+0.191412866E-2,
-0.206401795), (
-0.125028968,
+0.277995795)
1619+0.530160487,
+0.727527082
1620call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1622(
+0.191412866E-2,
-0.206401795), (
-0.125028968,
+0.277995795)
1624(
+0.530160487,
+0.00000000), (
+0.727527142,
+0.00000000)
1626+0.530160487,
+0.727527082
1628+0.530160487,
+0.727527142
1633 lb(isam)
= ub(isam
- 1)
+ 1
1641call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1642call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1647(
-0.785773754,
-0.413294792), (
-0.858735442,
+0.948035359), (
-0.114555955,
-0.920747757), (
+0.371067405,
+0.178279400), (
+0.198675632,
-0.188819528), (
+0.416556716,
-0.977314234), (
-0.810105085,
-0.971097708), (
-0.300586104,
-0.569568872E-1), (
+0.851434469,
+0.126441717)
1648(
-0.859861732,
-0.498287916), (
-0.979881883,
-0.874804974), (
+0.447566867,
-0.669304371), (
+0.958600998,
+0.789242983E-1), (
-0.479042530,
+0.346481204), (
+0.685823798,
-0.626841903), (
+0.925171733,
-0.699363351), (
+0.122157335,
-0.405358434), (
-0.935447931,
-0.206273317)
1649mean(:,
0)
= getVar(sample, dim)
1651(
+0.710976303,
+0.00000000), (
+0.726188362,
+0.00000000)
1652var(:,
0)
= getVar(sample, dim)
1654+0.710976303,
+0.726188362
1656 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1657 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1659call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1661(
-0.114669114,
-0.252830505), (
-0.127681494E-1,
-0.394981027)
1663+0.710976362,
+0.726188302
1664call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1666(
-0.114669114,
-0.252830505), (
-0.127681494E-1,
-0.394981027)
1668(
+0.710976303,
+0.00000000), (
+0.726188362,
+0.00000000)
1670+0.710976362,
+0.726188302
1672+0.710976303,
+0.726188362
1682 lb(isam)
= ub(isam
- 1)
+ 1
1690call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1691call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1696-0.643719077,
-0.855207324,
+0.450945616,
+0.362075090,
+0.551812768,
-0.690978646,
-0.744119763,
-0.518817544,
+0.198843718,
-0.131254196E-1
1697+0.834832788,
+0.807408810,
-0.139973044,
-0.257270098,
+0.667196393,
+0.625988841,
-0.250999451,
+0.990952849,
-0.267226100,
+0.290871620
1698+0.691654921,
+0.641442657,
-0.802313805,
+0.266432524,
+0.531755209,
-0.415774941,
-0.703921318,
-0.232786417,
+0.812192678,
-0.336648345
1701+1,
+9,
+3,
+1,
+3,
+7,
+10,
+6,
+5,
+1
1702mean(:,
0)
= getVar(sample,
2_IK, iweight)
1704+0.243014738,
+0.261459768,
+0.366754174
1705var(:,
0)
= getVar(sample,
2_IK, iweight)
1707+0.243014738,
+0.261459768,
+0.366754174
1709 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1710 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1712call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1714-0.421306401,
+0.352136910,
-0.370140523E-1
1716+0.243014738,
+0.261459768,
+0.366754144
1717call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1719-0.421306401,
+0.352136910,
-0.370140523E-1
1721+0.243014738,
+0.261459768,
+0.366754174
1723+0.243014738,
+0.261459768,
+0.366754144
1725+0.243014738,
+0.261459768,
+0.366754174
1730 lb(isam)
= ub(isam
- 1)
+ 1
1738call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1739call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1744+0.927004576,
-0.853206158,
+0.827405572,
+0.298742056E-1,
+0.560175776,
+0.499327660,
+0.257999063
1745-0.682084918,
+0.831631422,
-0.733739734,
-0.617742062,
+0.636015892,
+0.125550628,
-0.543881178
1748+8,
+6,
+6,
+10,
+3,
+4,
+4
1749mean(:,
0)
= getVar(sample,
2_IK, iweight)
1751+0.338830203,
+0.354304641
1752var(:,
0)
= getVar(sample,
2_IK, iweight)
1754+0.338830203,
+0.354304641
1756 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1757 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1759call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1761+0.299263656,
-0.263707906
1763+0.338830203,
+0.354304641
1764call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1766+0.299263656,
-0.263707906
1768+0.338830203,
+0.354304641
1770+0.338830203,
+0.354304641
1772+0.338830203,
+0.354304641
1777 lb(isam)
= ub(isam
- 1)
+ 1
1785call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1786call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1791+0.936405182,
-0.521316290,
+0.864827275,
+0.745933056E-1,
+0.207241893,
-0.220037937,
-0.933571935
1792+0.457517266,
+0.780159712,
-0.906429291,
-0.200025678,
+0.142363310,
+0.566818476,
-0.980188251
1795+1,
+10,
+10,
+3,
+3,
+1,
+8
1796mean(:,
0)
= getVar(sample,
2_IK, iweight)
1798+0.501959622,
+0.578007996
1799var(:,
0)
= getVar(sample,
2_IK, iweight)
1801+0.501959622,
+0.578007996
1803 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1804 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1806call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1808-0.686553568E-1,
-0.229245931
1810+0.501959562,
+0.578008056
1811call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1813-0.686553568E-1,
-0.229245931
1815+0.501959622,
+0.578007996
1817+0.501959562,
+0.578008056
1819+0.501959622,
+0.578007996
1824 lb(isam)
= ub(isam
- 1)
+ 1
1832call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1833call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1838+0.905704856,
+0.308492184E-1,
-0.787090182,
-0.392057419,
+0.813247800,
+0.993899107E-1,
-0.794128418,
-0.892188549,
-0.657222867,
-0.591577291,
+0.912078500,
-0.286368251
1839-0.236485481,
-0.755500674,
-0.697118640,
+0.698418975,
+0.742473006,
-0.841741204,
+0.574502468,
-0.772358775,
-0.755270481,
-0.685257316,
+0.461509228,
+0.685600042E-1
1840+0.939602852E-1,
+0.785410881,
+0.382186413,
-0.408938646,
+0.382983685E-1,
-0.322721720,
+0.682053804,
-0.377744317,
+0.374070883,
-0.204478264,
-0.925498724,
+0.671899199
1843+10,
+9,
+1,
+7,
+6,
+8,
+3,
+9,
+7,
+10,
+6,
+1
1844mean(:,
0)
= getVar(sample,
2_IK, iweight)
1846+0.446371555,
+0.373269498,
+0.225802064
1847var(:,
0)
= getVar(sample,
2_IK, iweight)
1849+0.446371555,
+0.373269498,
+0.225802064
1851 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1852 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1854call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1856-0.553831905E-1,
-0.282872856,
-0.322725177E-1
1858+0.446371585,
+0.373269498,
+0.225802049
1859call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1861-0.553831905E-1,
-0.282872856,
-0.322725177E-1
1863+0.446371555,
+0.373269498,
+0.225802064
1865+0.446371585,
+0.373269498,
+0.225802049
1867+0.446371555,
+0.373269498,
+0.225802064
1872 lb(isam)
= ub(isam
- 1)
+ 1
1880call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1881call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1886-0.995423555,
-0.503279090,
-0.822163343,
-0.510252118,
+0.437843800E-2,
-0.813341379,
-0.670410872
1887-0.237360716,
+0.333773136,
+0.279461145E-1,
+0.336447239,
+0.130962133,
-0.482149363,
+0.598000169
1888+0.213051677,
+0.889516592,
+0.729994297,
-0.711555481,
+0.967219949,
-0.897472262,
+0.326425195
1891+2,
+8,
+9,
+7,
+5,
+10,
+2
1892mean(:,
0)
= getVar(sample,
2_IK, iweight)
1894+0.737216994E-1,
+0.112813115,
+0.627091587
1895var(:,
0)
= getVar(sample,
2_IK, iweight)
1897+0.737216994E-1,
+0.112813115,
+0.627091587
1899 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1900 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1902call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1904-0.614899039,
+0.425913110E-1,
+0.131291240
1906+0.737217069E-1,
+0.112813100,
+0.627091646
1907call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1909-0.614899039,
+0.425913110E-1,
+0.131291240
1911+0.737216994E-1,
+0.112813115,
+0.627091587
1913+0.737217069E-1,
+0.112813100,
+0.627091646
1915+0.737216994E-1,
+0.112813115,
+0.627091587
1920 lb(isam)
= ub(isam
- 1)
+ 1
1928call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1929call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1934(
-0.851281643,
-0.986930609), (
-0.964846611E-1,
+0.146651864), (
+0.518817663,
-0.443830609), (
+0.414522409,
+0.234564304), (
-0.373918414,
+0.259147525), (
-0.220223665E-1,
+0.808806539), (
+0.650298595,
-0.275458097), (
+0.749525309,
-0.521147013), (
+0.445556283,
+0.493030787), (
+0.272839546,
+0.752681136), (
-0.876741409E-1,
-0.746885180), (
+0.912574291,
-0.379040122)
1935(
+0.666175961,
+0.480147839), (
+0.610077024,
-0.457082272), (
-0.747153401,
-0.776517391E-1), (
+0.693219185,
-0.693808079), (
+0.703310966E-2,
-0.215394735), (
+0.521205425,
-0.777920485), (
-0.946692944,
-0.929207087), (
+0.738097310,
+0.282233119), (
+0.731930733,
-0.426679850E-1), (
+0.479791164,
+0.105232000E-1), (
+0.217943907,
+0.529036999), (
+0.178094745,
+0.217967391)
1938+2,
+7,
+6,
+8,
+4,
+10,
+3,
+9,
+6,
+1,
+2,
+2
1939mean(:,
0)
= getVar(sample,
2_IK, iweight)
1941(
+0.445248663,
+0.00000000), (
+0.479920685,
+0.00000000)
1942var(:,
0)
= getVar(sample,
2_IK, iweight)
1944+0.445248663,
+0.479920685
1946 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1947 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1949call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1951(
+0.260463893,
+0.555534363E-1), (
+0.356202722,
-0.264923990)
1953+0.445248634,
+0.479920685
1954call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1956(
+0.260463893,
+0.555534363E-1), (
+0.356202722,
-0.264923990)
1958(
+0.445248663,
+0.00000000), (
+0.479920685,
+0.00000000)
1960+0.445248634,
+0.479920685
1962+0.445248663,
+0.479920685
1967 lb(isam)
= ub(isam
- 1)
+ 1
1975call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1976call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1981(
-0.604471922,
-0.722172856), (
-0.586372614E-1,
-0.765535116), (
-0.992283702,
+0.807943702), (
-0.487717509,
+0.104156017), (
-0.254726410,
+0.939801216), (
-0.909863949,
+0.103176475), (
-0.543208003,
-0.837131381)
1982(
-0.587000966,
+0.841384172), (
+0.663629413,
-0.793049216), (
+0.719133019,
-0.723772764), (
+0.986236334,
+0.142232656), (
+0.924685240,
+0.738291502), (
+0.523684144,
+0.216694236), (
-0.194679022,
-0.102457047)
1985+9,
+3,
+8,
+2,
+8,
+2,
+5
1986mean(:,
0)
= getVar(sample,
2_IK, iweight)
1988(
+0.689130783,
+0.00000000), (
+0.822354198,
+0.00000000)
1989var(:,
0)
= getVar(sample,
2_IK, iweight)
1991+0.689130783,
+0.822354198
1993 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1994 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1996call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1998(
-0.570363402,
+0.382380113E-1), (
+0.321753293,
+0.149054825)
2000+0.689130843,
+0.822354257
2001call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2003(
-0.570363402,
+0.382380113E-1), (
+0.321753293,
+0.149054825)
2005(
+0.689130783,
+0.00000000), (
+0.822354198,
+0.00000000)
2007+0.689130843,
+0.822354257
2009+0.689130783,
+0.822354198
2014 lb(isam)
= ub(isam
- 1)
+ 1
2022call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2023call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2028(
-0.356173277,
-0.950295568), (
+0.430824757E-1,
+0.240799665), (
-0.874439478E-1,
-0.805547595), (
+0.775285125,
+0.601713896), (
+0.617304564,
-0.802394032), (
+0.148570895,
+0.196672559), (
-0.715219021,
-0.606914997), (
-0.609425902,
-0.936508060), (
+0.829522252,
-0.111786127E-1), (
-0.263325572,
-0.558312058), (
+0.844980478E-1,
-0.950670958), (
-0.373999357,
-0.936996341)
2029(
+0.642017007,
-0.639532804E-1), (
+0.226229787,
-0.699286222), (
-0.474665284,
+0.890850902), (
-0.124120235,
+0.727239609), (
+0.800080299,
-0.591337681E-2), (
-0.535969973,
-0.212570429), (
-0.552075267,
-0.281301737), (
+0.138498902,
-0.416224837), (
+0.197449327,
-0.216166854), (
+0.338417053,
-0.694341421), (
-0.940167069,
-0.481247663), (
+0.136414766E-1,
-0.129129410)
2030(
+0.382838368,
+0.366681337), (
+0.558368206,
+0.170911431), (
+0.474322319,
+0.131924510), (
+0.449027538,
+0.712810040), (
-0.588616729,
+0.254755616), (
-0.953671813,
+0.753389239), (
-0.768563390,
-0.538604498), (
-0.574646950,
-0.400382161), (
-0.804148674,
+0.734889984), (
+0.356234908,
+0.105269432), (
-0.239333272,
-0.844217062), (
-0.183472633,
-0.910837770)
2033+7,
+8,
+10,
+8,
+8,
+7,
+5,
+3,
+9,
+2,
+4,
+5
2034mean(:,
0)
= getVar(sample,
2_IK, iweight)
2036(
+0.560278833,
+0.00000000), (
+0.502727032,
+0.00000000), (
+0.619186282,
+0.00000000)
2037var(:,
0)
= getVar(sample,
2_IK, iweight)
2039+0.560278833,
+0.502727032,
+0.619186282
2041 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
2042 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
2044call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2046(
+0.120531492,
-0.375774711), (
-0.487191975E-2,
-0.185637102E-1), (
-0.129851326,
+0.174564198)
2048+0.560278773,
+0.502727032,
+0.619186163
2049call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2051(
+0.120531492,
-0.375774711), (
-0.487191975E-2,
-0.185637102E-1), (
-0.129851326,
+0.174564198)
2053(
+0.560278833,
+0.00000000), (
+0.502727032,
+0.00000000), (
+0.619186282,
+0.00000000)
2055+0.560278773,
+0.502727032,
+0.619186163
2057+0.560278833,
+0.502727032,
+0.619186282
2062 lb(isam)
= ub(isam
- 1)
+ 1
2070call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2071call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2076(
+0.789366961E-1,
-0.459667087), (
+0.555796146,
-0.999973059), (
+0.645811200,
-0.851449728), (
-0.376818180,
-0.988416433), (
-0.295202613,
-0.753653407), (
-0.652528405,
+0.149038315), (
-0.281677485,
-0.990754247), (
+0.296763778,
-0.147832155), (
+0.535935044,
+0.262844443), (
-0.715446830,
-0.838607669)
2079+8,
+5,
+7,
+2,
+8,
+8,
+1,
+8,
+2,
+2
2080mean(:,
0)
= getVar(sample,
2_IK, iweight)
2082(
+0.399099320,
+0.00000000)
2083var(:,
0)
= getVar(sample,
2_IK, iweight)
2087 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
2088 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
2090call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2092(
+0.260603391E-1,
-0.485804737)
2095call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2097(
+0.260603391E-1,
-0.485804737)
2099(
+0.399099320,
+0.00000000)
2108 lb(isam)
= ub(isam
- 1)
+ 1
2116call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2117call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2122(
+0.926758051E-1,
-0.655876875), (
+0.843337893,
+0.290679097), (
-0.752840757,
+0.496578574), (
+0.751332760,
-0.224754930), (
-0.154177785,
-0.299545169), (
+0.343067050,
+0.270943999), (
+0.474354029E-1,
+0.884594440), (
-0.530583382,
+0.971300602E-1)
2125+9,
+7,
+10,
+2,
+9,
+10,
+10,
+3
2126mean(:,
0)
= getVar(sample,
2_IK, iweight)
2128(
+0.481507272,
+0.00000000)
2129var(:,
0)
= getVar(sample,
2_IK, iweight)
2133 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
2134 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
2136call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2138(
+0.272896588E-1,
+0.163316786)
2141call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
2143(
+0.272896588E-1,
+0.163316786)
2145(
+0.481507272,
+0.00000000)
2159 lb(isam)
= ub(isam
- 1)
+ 1
2167call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2168call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2173-0.456859827,
+0.846340775,
-0.529965758,
-0.529519439,
-0.531724572,
-0.953049302,
+0.453287363,
+0.871036887
2174+0.166800976,
-0.741932273,
-0.113613844,
-0.748641491,
-0.368432879,
+0.924904704,
+0.704228640,
+0.286103725
2175-0.549233437,
-0.343745589,
+0.656392455,
+0.232125282,
+0.167057753,
+0.423360467,
+0.834948659,
+0.533569694
2178+1.91578388,
+1.15440679,
+1.63607156,
+1.16868019,
+1.48426163,
+1.07659340,
+1.26548696,
+1.85025275
2179mean(:,
0)
= getVar(sample,
2_IK, rweight)
2181+0.431653619,
+0.282213181,
+0.218195528
2182var(:,
0)
= getVar(sample,
2_IK, rweight)
2184+0.431653619,
+0.282213181,
+0.218195528
2186 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2187 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2189call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2191-0.877910405E-1,
+0.235217065E-1,
+0.228865653
2193+0.431653559,
+0.282213181,
+0.218195498
2194call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2196-0.877910405E-1,
+0.235217065E-1,
+0.228865653
2198+0.431653619,
+0.282213181,
+0.218195528
2200+0.431653559,
+0.282213181,
+0.218195498
2202+0.431653619,
+0.282213181,
+0.218195528
2207 lb(isam)
= ub(isam
- 1)
+ 1
2215call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2216call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2221-0.952730060,
+0.417450070,
+0.382472157,
-0.969472051,
+0.620872378,
-0.354057431,
+0.287482977,
+0.840806365,
-0.448304653,
-0.916323543,
-0.241936564
2222+0.355140805,
-0.192204118,
-0.904233813,
-0.706048250,
+0.524552584,
-0.972798228,
+0.333295703,
+0.913730741,
-0.999574184,
+0.770824075,
+0.752492905
2225+1.89920223,
+1.93254256,
+1.61339211,
+1.76530468,
+1.03323865,
+1.51040733,
+1.81978631,
+1.25476003,
+1.21311593,
+1.65833759,
+1.27826726
2226mean(:,
0)
= getVar(sample,
2_IK, rweight)
2228+0.402585149,
+0.493393302
2229var(:,
0)
= getVar(sample,
2_IK, rweight)
2231+0.402585149,
+0.493393302
2233 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2234 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2236call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2238-0.164020360,
-0.323324800E-1
2240+0.402585179,
+0.493393302
2241call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2243-0.164020360,
-0.323324800E-1
2245+0.402585149,
+0.493393302
2247+0.402585179,
+0.493393302
2249+0.402585149,
+0.493393302
2254 lb(isam)
= ub(isam
- 1)
+ 1
2262call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2263call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2268+0.678837419,
-0.104835033,
-0.141752720,
+0.305764675E-1,
-0.168594837,
+0.702507615,
+0.115168214
2271+1.00804949,
+1.87534475,
+1.66213238,
+1.35118973,
+1.35291338,
+1.44975519,
+1.79031503
2272mean(:,
0)
= getVar(sample,
2_IK, rweight)
2275var(:,
0)
= getVar(sample,
2_IK, rweight)
2279 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2280 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2282call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2287call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2300 lb(isam)
= ub(isam
- 1)
+ 1
2308call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2309call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2314-0.590085983E-1,
+0.665111303,
-0.860111475,
+0.600747705,
+0.477968574,
+0.938823223,
+0.506662369,
-0.664431095,
+0.306898355,
+0.769069195
2315+0.866360784,
-0.552640080,
-0.592116594,
-0.741476417,
-0.945961356,
-0.169301748,
-0.820216656,
-0.336416125,
+0.696911812,
-0.724162936
2316-0.841434360,
+0.328338504,
+0.995761633,
+0.221701026,
+0.866441846,
-0.921578050,
+0.230515003E-2,
-0.475273013,
-0.357847810,
+0.283270717
2317-0.622297764,
-0.670937300,
-0.403313398,
+0.171689391,
+0.372068524,
+0.741893411,
+0.860774159,
+0.515783548,
-0.331782818,
+0.646715164E-1
2320+1.31561375,
+1.56880832,
+1.33455181,
+1.83362663,
+1.84589934,
+1.50811875,
+1.30193865,
+1.60850406,
+1.88360548,
+1.17562199
2321mean(:,
0)
= getVar(sample,
2_IK, rweight)
2323+0.313359886,
+0.363156915,
+0.382107615,
+0.266435415
2324var(:,
0)
= getVar(sample,
2_IK, rweight)
2326+0.313359886,
+0.363156915,
+0.382107615,
+0.266435415
2328 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2329 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2331call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2333+0.279048830,
-0.326873511,
+0.162927955E-1,
+0.723423734E-1
2335+0.313359886,
+0.363156915,
+0.382107615,
+0.266435444
2336call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2338+0.279048830,
-0.326873511,
+0.162927955E-1,
+0.723423734E-1
2340+0.313359886,
+0.363156915,
+0.382107615,
+0.266435415
2342+0.313359886,
+0.363156915,
+0.382107615,
+0.266435444
2344+0.313359886,
+0.363156915,
+0.382107615,
+0.266435415
2349 lb(isam)
= ub(isam
- 1)
+ 1
2357call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2358call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2363-0.380760312,
-0.181761146,
-0.551533222,
-0.252892613,
-0.533695459,
+0.153751612,
-0.348669410,
+0.561230779,
+0.658128858,
+0.937129021,
-0.522557139,
-0.703782201,
-0.673147082
2364-0.850821376,
-0.344029188,
-0.978334785,
+0.912092090,
+0.477890134,
-0.179162741,
+0.784522176,
+0.158455610,
-0.765434980,
-0.778558731,
+0.305698276,
+0.176838994,
-0.454311371E-1
2365+0.159700632,
+0.382051468E-1,
+0.950282097,
+0.138051152,
-0.958918810,
-0.166980743,
-0.432868838,
-0.311475277,
+0.492194653,
-0.528613329,
+0.688581467E-1,
+0.440961599,
+0.911751628
2366-0.524776340,
-0.390508533,
+0.791839242,
+0.441127181,
-0.270178318,
-0.740340233,
-0.234954000,
-0.272068262,
+0.622387528,
+0.566700697E-1,
+0.550194502,
-0.596840262,
+0.568737864
2367-0.596731782,
-0.980297446,
-0.880720139,
+0.675351501,
+0.860759020E-1,
+0.284641027,
-0.702136040,
-0.713813305E-2,
+0.894187808,
-0.124615192,
+0.582032084,
+0.486657143,
-0.363450646
2368-0.361077070,
+0.867906809,
-0.863591671,
-0.945906401,
-0.593127728,
-0.277345061,
+0.572922230E-1,
-0.306186318,
-0.822239757,
+0.786304474E-2,
-0.126777172,
-0.142967105,
+0.477281809
2371+1.51350498,
+1.89345598,
+1.29765248,
+1.30025578,
+1.67977345,
+1.35092080,
+1.91176558,
+1.49859715,
+1.47617102,
+1.05615437,
+1.90248680,
+1.73263478,
+1.75387597
2372mean(:,
0)
= getVar(sample,
2_IK, rweight)
2374+0.243654490,
+0.343116134,
+0.279995084,
+0.254535228,
+0.364152372,
+0.262150049
2375var(:,
0)
= getVar(sample,
2_IK, rweight)
2377+0.243654490,
+0.343116134,
+0.279995084,
+0.254535228,
+0.364152372,
+0.262150049
2379 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2380 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2382call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2384-0.192081332,
-0.426470488E-1,
+0.617861599E-1,
-0.125053227E-1,
-0.661122724E-1,
-0.188110501
2386+0.243654519,
+0.343116164,
+0.279995114,
+0.254535198,
+0.364152431,
+0.262150049
2387call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2389-0.192081332,
-0.426470488E-1,
+0.617861599E-1,
-0.125053227E-1,
-0.661122724E-1,
-0.188110501
2391+0.243654490,
+0.343116134,
+0.279995084,
+0.254535228,
+0.364152372,
+0.262150049
2393+0.243654519,
+0.343116164,
+0.279995114,
+0.254535198,
+0.364152431,
+0.262150049
2395+0.243654490,
+0.343116134,
+0.279995084,
+0.254535228,
+0.364152372,
+0.262150049
2400 lb(isam)
= ub(isam
- 1)
+ 1
2408call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2409call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2414(
-0.559929132,
+0.107594609), (
-0.127242565,
+0.910096884), (
-0.265121579,
-0.456588626), (
+0.228326917,
+0.470838547), (
+0.895473957,
+0.338093758), (
-0.510994196,
-0.222503781), (
+0.996068716,
-0.968947411E-1), (
+0.387076139E-1,
+0.746616483), (
+0.962647915,
+0.732249379)
2415(
-0.224885702,
-0.549366951), (
-0.468231559,
+0.669173479), (
-0.333083391,
+0.280504227), (
+0.291358113,
+0.147912502E-1), (
+0.125867486,
+0.545435309), (
-0.259938598,
-0.162636042E-1), (
-0.648244262,
-0.566329598), (
-0.656138062,
+0.851830006), (
-0.681183696,
-0.794881105)
2418+1.95570564,
+1.48555732,
+1.48572564,
+1.94028366,
+1.03725314,
+1.65249991,
+1.59643817,
+1.66730642,
+1.70330381
2419mean(:,
0)
= getVar(sample,
2_IK, rweight)
2421(
+0.539648235,
+0.00000000), (
+0.416557014,
+0.00000000)
2422var(:,
0)
= getVar(sample,
2_IK, rweight)
2424+0.539648235,
+0.416557014
2426 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2427 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2429call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2431(
+0.147605628,
+0.283530354), (
-0.320369333,
+0.456161937E-2)
2433+0.539648175,
+0.416557014
2434call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2436(
+0.147605628,
+0.283530354), (
-0.320369333,
+0.456161937E-2)
2438(
+0.539648235,
+0.00000000), (
+0.416557014,
+0.00000000)
2440+0.539648175,
+0.416557014
2442+0.539648235,
+0.416557014
2447 lb(isam)
= ub(isam
- 1)
+ 1
2455call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2456call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2461(
-0.475404501,
-0.745295286), (
+0.654235840,
+0.180504441), (
-0.477253795,
+0.838284016), (
+0.623126030E-1,
-0.580994964), (
-0.336308122,
+0.918545365), (
-0.723591685,
+0.406844139), (
+0.261555791,
-0.496661901), (
-0.395248413,
-0.209601998), (
-0.286167622,
-0.860378027), (
-0.880562067,
+0.718285918), (
+0.230759501,
+0.133298635E-1)
2462(
-0.330305934,
+0.869467258E-1), (
+0.257432103,
-0.260700941), (
-0.948377490,
-0.154860735), (
-0.262456775,
+0.644871235), (
-0.481317997,
+0.409146786), (
-0.253222704,
-0.978476048), (
+0.807442904,
+0.433811545), (
-0.865841627,
+0.677737474), (
-0.306224704,
-0.165075779), (
-0.543186426,
-0.639907598), (
+0.608236074,
-0.456133366)
2465+1.50786710,
+1.90293026,
+1.22824669,
+1.23750782,
+1.61801374,
+1.95698977,
+1.52990723,
+1.46897864,
+1.81417561,
+1.13496017,
+1.33197188
2466mean(:,
0)
= getVar(sample,
2_IK, rweight)
2468(
+0.577200532,
+0.00000000), (
+0.538022697,
+0.00000000)
2469var(:,
0)
= getVar(sample,
2_IK, rweight)
2471+0.577200532,
+0.538022697
2473 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2474 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2476call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2478(
-0.199193686,
+0.102087110E-2), (
-0.189499885,
-0.588145107E-1)
2480+0.577200532,
+0.538022637
2481call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2483(
-0.199193686,
+0.102087110E-2), (
-0.189499885,
-0.588145107E-1)
2485(
+0.577200532,
+0.00000000), (
+0.538022697,
+0.00000000)
2487+0.577200532,
+0.538022637
2489+0.577200532,
+0.538022697
2494 lb(isam)
= ub(isam
- 1)
+ 1
2502call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2503call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2508(
-0.483351827,
-0.942086101), (
-0.153762460,
-0.816159725), (
-0.215927482,
+0.623725653), (
-0.635021091,
-0.687877059), (
-0.409031868,
-0.519144535), (
+0.138303041E-1,
+0.846430302), (
-0.792427897,
-0.760880470), (
+0.173882008,
+0.578499794), (
-0.126905680,
-0.144201875), (
+0.482061744,
+0.434282184)
2511+1.93162394,
+1.85349691,
+1.93623066,
+1.37309849,
+1.87251413,
+1.37600076,
+1.30573750,
+1.88095713,
+1.51429009,
+1.76100111
2512mean(:,
0)
= getVar(sample,
2_IK, rweight)
2514(
+0.559749186,
+0.00000000)
2515var(:,
0)
= getVar(sample,
2_IK, rweight)
2519 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2520 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2522call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2524(
-0.196755737,
-0.133040398)
2527call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2529(
-0.196755737,
-0.133040398)
2531(
+0.559749186,
+0.00000000)
2540 lb(isam)
= ub(isam
- 1)
+ 1
2548call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2549call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2554(
-0.694346547,
+0.443354607), (
+0.155163646,
-0.194434881), (
+0.678552389,
+0.993994474), (
+0.276175737E-1,
+0.479347944), (
-0.740348339,
+0.664105058), (
+0.147382617,
+0.890783191), (
-0.662088156,
-0.781434774), (
-0.323125124,
-0.418649197)
2557+1.81097853,
+1.31777596,
+1.65954888,
+1.55192065,
+1.80724859,
+1.72039008,
+1.18610168,
+1.86958385
2558mean(:,
0)
= getVar(sample,
2_IK, rweight)
2560(
+0.571343303,
+0.00000000)
2561var(:,
0)
= getVar(sample,
2_IK, rweight)
2565 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2566 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2568call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2570(
-0.182448030,
+0.306673706)
2573call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2575(
-0.182448030,
+0.306673706)
2577(
+0.571343303,
+0.00000000)
2586 lb(isam)
= ub(isam
- 1)
+ 1
2594call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2595call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
2600(
-0.522180438,
-0.978707910), (
+0.913498282,
-0.639264941), (
+0.510934591,
-0.460795403), (
+0.566257238,
-0.191065907), (
+0.802366972,
+0.688931942), (
+0.823998809,
+0.839356780), (
+0.753705502E-1,
+0.415634751)
2603+1.62331700,
+1.67917943,
+1.08532345,
+1.72158051,
+1.65279627,
+1.24655652,
+1.58064413
2604mean(:,
0)
= getVar(sample,
2_IK, rweight)
2606(
+0.656648755,
+0.00000000)
2607var(:,
0)
= getVar(sample,
2_IK, rweight)
2611 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2612 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2614call setVarMeanMerged(varMerged, meanMerged, var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2616(
+0.442715883,
-0.613161027E-1)
2619call setVarMeanMerged(var(:,
2), mean(:,
2), var(:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2621(
+0.442715883,
-0.613161027E-1)
2623(
+0.656648755,
+0.00000000)