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 of a univariate sample.")
26 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
30 real(TKG) :: mean(
1:nsam)
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(
"var(0) = getVar(sample)")
58 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)))")
59 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)))")
62 var(isam)
= getVar(sample(lb(isam):ub(isam)))
63 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
65 call disp%show(
"call setVarMerged(varMerged, var(2), var(1), mean(1) - mean(2), ub(1) / real(ub(2), TKG))")
66 call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
69 call disp%show(
"call setVarMerged(var(2), var(1), mean(1) - mean(2), ub(1) / real(ub(2), TKG))")
70 call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
80 complex(TKG) :: mean(
1:nsam)
81 complex(TKG),
allocatable :: sample(:)
82 real(TKG) :: var(
0:nsam), varMerged
85 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
88 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
89 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
92 lb(isam)
= ub(isam
- 1)
+ 1
99 call disp%show(
"sample = cmplx(getUnifRand(0., 1., ub(nsam)), -getUnifRand(0., 1., ub(nsam)), TKG)")
103 call disp%show(
"var(0) = getVar(sample)")
105 call disp%show(
"var(0) ! reference")
108 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)))")
109 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)))")
112 var(isam)
= getVar(sample(lb(isam):ub(isam)))
113 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
115 call disp%show(
"call setVarMerged(varMerged, var(2), var(1), mean(1) - mean(2), ub(1) / real(ub(2), TKG))")
116 call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
119 call disp%show(
"call setVarMerged(var(2), var(1), mean(1) - mean(2), ub(1) / real(ub(2), TKG))")
120 call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
123 call disp%show(
"var(0) ! reference")
130 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
131 call disp%show(
"!Compute the biased merged variance of a frequency weighted univariate sample.")
132 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
136 real(TKG) :: mean(
1:nsam)
137 real(TKG),
allocatable :: sample(:)
138 real(TKG) :: var(
0:nsam), varMerged
141 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
144 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
145 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
148 lb(isam)
= ub(isam
- 1)
+ 1
155 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
159 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, 1, IK))")
163 call disp%show(
"var(0) = getVar(sample, iweight)")
164 var(
0)
= getVar(sample, iweight)
165 call disp%show(
"var(0) ! reference")
168 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
169 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
172 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
173 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
175 call disp%show(
"call setVarMerged(varMerged, var(2), var(1), mean(1) - mean(2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
176 call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
179 call disp%show(
"call setVarMerged(var(2), var(1), mean(1) - mean(2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
180 call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
183 call disp%show(
"var(0) ! reference")
190 complex(TKG) :: mean(
1:nsam)
191 complex(TKG),
allocatable :: sample(:)
192 real(TKG) :: var(
0:nsam), varMerged
195 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
198 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
199 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
202 lb(isam)
= ub(isam
- 1)
+ 1
209 call disp%show(
"sample = cmplx(getUnifRand(0., 1., ub(nsam)), -getUnifRand(0., 1., ub(nsam)), TKG)")
213 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, 1, IK))")
217 call disp%show(
"var(0) = getVar(sample, iweight)")
218 var(
0)
= getVar(sample, iweight)
219 call disp%show(
"var(0) ! reference")
222 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
223 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
226 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
227 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
229 call disp%show(
"call setVarMerged(varMerged, var(2), var(1), mean(1) - mean(2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
230 call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
233 call disp%show(
"call setVarMerged(var(2), var(1), mean(1) - mean(2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
234 call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
237 call disp%show(
"var(0) ! reference")
244 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
245 call disp%show(
"!Compute the biased merged variance of a reliability weighted univariate sample.")
246 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
250 real(TKG) :: mean(
1:nsam)
251 real(TKG),
allocatable :: sample(:)
252 real(TKG),
allocatable :: rweight(:)
253 real(TKG) :: var(
0:nsam), varMerged
256 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
259 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
260 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
263 lb(isam)
= ub(isam
- 1)
+ 1
270 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
274 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, 1, IK))")
278 call disp%show(
"var(0) = getVar(sample, rweight)")
279 var(
0)
= getVar(sample, rweight)
280 call disp%show(
"var(0) ! reference")
283 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
284 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
287 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
288 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
290 call disp%show(
"call setVarMerged(varMerged, var(2), var(1), mean(1) - mean(2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
291 call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
294 call disp%show(
"call setVarMerged(var(2), var(1), mean(1) - mean(2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
295 call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
298 call disp%show(
"var(0) ! reference")
305 complex(TKG) :: mean(
1:nsam)
306 complex(TKG),
allocatable :: sample(:)
307 real(TKG),
allocatable :: rweight(:)
308 real(TKG) :: var(
0:nsam), varMerged
311 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
314 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
315 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
318 lb(isam)
= ub(isam
- 1)
+ 1
325 call disp%show(
"sample = cmplx(getUnifRand(0., 1., ub(nsam)), -getUnifRand(0., 1., ub(nsam)), TKG)")
329 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, 1, IK))")
333 call disp%show(
"var(0) = getVar(sample, rweight)")
334 var(
0)
= getVar(sample, rweight)
335 call disp%show(
"var(0) ! reference")
338 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
339 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
342 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
343 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
345 call disp%show(
"call setVarMerged(varMerged, var(2), var(1), mean(1) - mean(2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
346 call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
349 call disp%show(
"call setVarMerged(var(2), var(1), mean(1) - mean(2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
350 call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
353 call disp%show(
"var(0) ! reference")
360 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
361 call disp%show(
"!Compute the biased merged variance of a multivariate sample.")
362 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
366 real(TKG),
allocatable :: mean(:,:)
367 real(TKG),
allocatable :: sample(:,:)
368 real(TKG),
allocatable :: var(:,:), varMerged(:)
371 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
374 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
375 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
378 lb(isam)
= ub(isam
- 1)
+ 1
385 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
387 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
388 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
389 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
390 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
391 call disp%show(
"call setResized(varMerged, ndim)")
393 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
397 call disp%show(
"var(:,0) = getVar(sample, dim)")
398 var(:,
0)
= getVar(sample, dim)
399 call disp%show(
"var(:,0) ! reference")
402 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), dim)")
403 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
406 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
407 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
409 call disp%show(
"call setVarMerged(varMerged, var(:,2), var(:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
410 call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
413 call disp%show(
"call setVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
414 call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
417 call disp%show(
"var(:,0) ! reference")
424 complex(TKG),
allocatable :: mean(:,:)
425 complex(TKG),
allocatable :: sample(:,:)
426 real(TKG),
allocatable :: var(:,:), varMerged(:)
429 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
432 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
433 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
436 lb(isam)
= ub(isam
- 1)
+ 1
443 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
445 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
446 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
447 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
448 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
449 call disp%show(
"call setResized(varMerged, ndim)")
451 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ndim, ub(nsam)), -getUnifRand(-1., +1., ndim, ub(nsam)), TKG)")
455 call disp%show(
"var(:,0) = getVar(sample, dim)")
456 var(:,
0)
= getVar(sample, dim)
457 call disp%show(
"var(:,0) ! reference")
460 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), dim)")
461 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
464 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
465 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
467 call disp%show(
"call setVarMerged(varMerged, var(:,2), var(:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
468 call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
471 call disp%show(
"call setVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
472 call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
475 call disp%show(
"var(:,0) ! reference")
482 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
483 call disp%show(
"!Compute the biased merged variance of a frequency weighted multivariate sample.")
484 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
488 real(TKG),
allocatable :: mean(:,:)
489 real(TKG),
allocatable :: sample(:,:)
490 real(TKG),
allocatable :: var(:,:), varMerged(:)
493 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
496 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
497 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
500 lb(isam)
= ub(isam
- 1)
+ 1
507 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
509 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
510 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
511 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
512 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
513 call disp%show(
"call setResized(varMerged, ndim)")
515 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
519 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
523 call disp%show(
"var(:,0) = getVar(sample, 2_IK, iweight)")
524 var(:,
0)
= getVar(sample,
2_IK, iweight)
525 call disp%show(
"var(:,0) ! reference")
528 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
529 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
532 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
533 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
535 call disp%show(
"call setVarMerged(varMerged, var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
536 call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
539 call disp%show(
"call setVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
540 call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
543 call disp%show(
"var(:,0) ! reference")
550 complex(TKG),
allocatable :: mean(:,:)
551 complex(TKG),
allocatable :: sample(:,:)
552 real(TKG),
allocatable :: var(:,:), varMerged(:)
555 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
558 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
559 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
562 lb(isam)
= ub(isam
- 1)
+ 1
569 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
571 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
572 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
573 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
574 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
575 call disp%show(
"call setResized(varMerged, ndim)")
577 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ndim, ub(nsam)), getUnifRand(-1., +1., ndim, ub(nsam)), TKG)")
581 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
585 call disp%show(
"var(:,0) = getVar(sample, 2_IK, iweight)")
586 var(:,
0)
= getVar(sample,
2_IK, iweight)
587 call disp%show(
"var(:,0) ! reference")
590 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
591 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
594 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
595 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
597 call disp%show(
"call setVarMerged(varMerged, var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
598 call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
601 call disp%show(
"call setVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
602 call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
605 call disp%show(
"var(:,0) ! reference")
612 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
613 call disp%show(
"!Compute the biased merged variance of a reliability weighted multivariate sample.")
614 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
618 real(TKG),
allocatable :: mean(:,:)
619 real(TKG),
allocatable :: sample(:,:)
620 real(TKG),
allocatable :: var(:,:), varMerged(:)
621 real(TKG),
allocatable :: rweight(:)
624 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
627 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
628 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
631 lb(isam)
= ub(isam
- 1)
+ 1
638 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
640 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
641 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
642 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
643 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
644 call disp%show(
"call setResized(varMerged, ndim)")
646 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
650 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
654 call disp%show(
"var(:,0) = getVar(sample, 2_IK, rweight)")
655 var(:,
0)
= getVar(sample,
2_IK, rweight)
656 call disp%show(
"var(:,0) ! reference")
659 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
660 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
663 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
664 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
666 call disp%show(
"call setVarMerged(varMerged, var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
667 call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
670 call disp%show(
"call setVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
671 call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
674 call disp%show(
"var(:,0) ! reference")
681 complex(TKG),
allocatable :: mean(:,:)
682 complex(TKG),
allocatable :: sample(:,:)
683 real(TKG),
allocatable :: var(:,:), varMerged(:)
684 real(TKG),
allocatable :: rweight(:)
687 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
690 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
691 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
694 lb(isam)
= ub(isam
- 1)
+ 1
701 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
703 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
704 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
705 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
706 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
707 call disp%show(
"call setResized(varMerged, ndim)")
709 call disp%show(
"sample = cmplx(getUnifRand(-1., +1., ndim, ub(nsam)), getUnifRand(-1., +1., ndim, ub(nsam)), TKG)")
713 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
717 call disp%show(
"var(:,0) = getVar(sample, 2_IK, rweight)")
718 var(:,
0)
= getVar(sample,
2_IK, rweight)
719 call disp%show(
"var(:,0) ! reference")
722 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
723 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
726 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
727 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
729 call disp%show(
"call setVarMerged(varMerged, var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
730 call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
733 call disp%show(
"call setVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
734 call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
737 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.282681763,
+0.458146751,
+0.750241160,
+0.940724730,
+0.566639543,
+0.155023754
23 var(isam)
= getVar(sample(lb(isam):ub(isam)))
24 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
26call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
29call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
38 lb(isam)
= ub(isam
- 1)
+ 1
47+0.616792977,
+0.222903907,
+0.774933517,
+0.932100832,
+0.511356890,
+0.324885726,
+0.974224746
52 var(isam)
= getVar(sample(lb(isam):ub(isam)))
53 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
55call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
58call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
67 lb(isam)
= ub(isam
- 1)
+ 1
76+0.628690124E-1,
+0.192238867,
+0.161434114,
+0.699594796,
+0.331201732,
+0.262139618,
+0.944855928,
+0.578415036,
+0.875757456,
+0.306986511,
+0.882413089
81 var(isam)
= getVar(sample(lb(isam):ub(isam)))
82 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
84call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
87call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
96 lb(isam)
= ub(isam
- 1)
+ 1
105+0.684865713E-1,
+0.222889483,
+0.927111328,
+0.343428493,
+0.806115270,
+0.395017862E-1,
+0.929943383,
+0.151185930
110 var(isam)
= getVar(sample(lb(isam):ub(isam)))
111 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
113call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
116call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
125 lb(isam)
= ub(isam
- 1)
+ 1
134+0.620729387,
+0.505263984,
+0.100430250,
+0.608765006,
+0.619660914,
+0.949781060,
+0.385886490,
+0.605916917,
+0.131013751
139 var(isam)
= getVar(sample(lb(isam):ub(isam)))
140 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
142call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
145call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
154 lb(isam)
= ub(isam
- 1)
+ 1
163(
+0.169601917,
-0.252396464E-1), (
+0.394556463,
-0.991977513), (
+0.132779360,
-0.284171999), (
+0.657557845,
-0.564487040), (
+0.288988829,
-0.403937459), (
+0.707922578,
-0.987175703), (
+0.464062154,
-0.663727582), (
+0.224713802,
-0.506084085), (
+0.886198521,
-0.507968962), (
+0.944580555,
-0.947833300), (
+0.495980680,
-0.714342713), (
+0.611722946,
-0.490783453)
168 var(isam)
= getVar(sample(lb(isam):ub(isam)))
169 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
171call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
174call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
183 lb(isam)
= ub(isam
- 1)
+ 1
192(
+0.472440183,
-0.234886527), (
+0.152583003,
-0.570800841), (
+0.431730568,
-0.846615672), (
+0.727401614,
-0.148321211), (
+0.810052812,
-0.618840754), (
+0.914186239E-2,
-0.509600461), (
+0.935771823,
-0.593167305), (
+0.597235978,
-0.578141093), (
+0.940371096,
-0.884296060), (
+0.989422202E-1,
-0.311418772E-1), (
+0.588447094,
-0.438360870)
197 var(isam)
= getVar(sample(lb(isam):ub(isam)))
198 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
200call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
203call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
212 lb(isam)
= ub(isam
- 1)
+ 1
221(
+0.345816255,
-0.303678036), (
+0.886424541,
-0.461009383), (
+0.251661718,
-0.986440718), (
+0.547977686,
-0.857630372E-1), (
+0.522475243E-1,
-0.501887977), (
+0.229573250E-2,
-0.614783347), (
+0.519573033,
-0.507852733)
226 var(isam)
= getVar(sample(lb(isam):ub(isam)))
227 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
229call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
232call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
241 lb(isam)
= ub(isam
- 1)
+ 1
250(
+0.274543166,
-0.556082845), (
+0.660316288,
-0.984699130), (
+0.959764123E-1,
-0.426200569), (
+0.768561840,
-0.367408693), (
+0.311403275E-1,
-0.204791427E-1), (
+0.931620955,
-0.879353285E-2), (
+0.483336031,
-0.251451313), (
+0.465781748,
-0.884848595)
255 var(isam)
= getVar(sample(lb(isam):ub(isam)))
256 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
258call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
261call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
270 lb(isam)
= ub(isam
- 1)
+ 1
279(
+0.686789572,
-0.118137598), (
+0.784673393,
-0.786661148), (
+0.506868184,
-0.270289004), (
+0.854646027,
-0.958091438), (
+0.125769973E-1,
-0.307230294), (
+0.789779425,
-0.302168727), (
+0.949921250,
-0.413799822), (
+0.724042237,
-0.279695094), (
+0.930103660E-1,
-0.577169001), (
+0.297805965,
-0.953446388), (
+0.311858714,
-0.156486988), (
+0.896656513E-1,
-0.975592434), (
+0.910786152,
-0.258129358)
284 var(isam)
= getVar(sample(lb(isam):ub(isam)))
285 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
287call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
290call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
304 lb(isam)
= ub(isam
- 1)
+ 1
313+0.347384691,
+0.588463366,
+0.781132758,
+0.501089036,
+0.321985483E-1,
+0.877985358,
+0.335842133
316+2,
+4,
+3,
+6,
+4,
+1,
+8
317var(
0)
= getVar(sample, iweight)
321 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
322 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
324call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
327call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
336 lb(isam)
= ub(isam
- 1)
+ 1
345+0.201636553E-2,
+0.670176923,
+0.998277843,
+0.622785985,
+0.888784349,
+0.625781059,
+0.116732538,
+0.767450988
348+1,
+2,
+2,
+8,
+9,
+1,
+9,
+7
349var(
0)
= getVar(sample, iweight)
353 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
354 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
356call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
359call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
368 lb(isam)
= ub(isam
- 1)
+ 1
377+0.475506961,
+0.199834108E-1,
+0.525703907,
+0.646607757,
+0.922822058,
+0.441325545
380+9,
+10,
+1,
+2,
+4,
+2
381var(
0)
= getVar(sample, iweight)
385 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
386 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
388call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
391call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
400 lb(isam)
= ub(isam
- 1)
+ 1
409+0.820667863,
+0.912171483,
+0.198023915E-1,
+0.439678431,
+0.375530899,
+0.788326263E-1,
+0.466294289,
+0.926025510,
+0.851589441E-1,
+0.730477273,
+0.596710980
412+10,
+7,
+9,
+5,
+2,
+2,
+4,
+9,
+2,
+7,
+5
413var(
0)
= getVar(sample, iweight)
417 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
418 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
420call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
423call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
432 lb(isam)
= ub(isam
- 1)
+ 1
441+0.121591091E-1,
+0.829657674,
+0.365109503,
+0.946821511,
+0.986628592,
+0.519187450,
+0.399219930,
+0.244624496,
+0.805534422,
+0.237654686,
+0.680320680,
+0.304630697
444+10,
+5,
+10,
+7,
+2,
+1,
+1,
+6,
+5,
+4,
+9,
+5
445var(
0)
= getVar(sample, iweight)
449 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
450 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
452call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
455call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
464 lb(isam)
= ub(isam
- 1)
+ 1
473(
+0.728536308,
-0.290056169), (
+0.280530989,
-0.483283699), (
+0.787321568,
-0.967009306), (
+0.829695821,
-0.602520764)
477var(
0)
= getVar(sample, iweight)
481 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
482 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
484call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
487call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
496 lb(isam)
= ub(isam
- 1)
+ 1
505(
+0.460407555,
-0.470099330), (
+0.215232611,
-0.678497672), (
+0.568859339,
-0.352218807), (
+0.140970647,
-0.105994463), (
+0.302153349,
-0.468982816), (
+0.960202098,
-0.525167942), (
+0.988228142,
-0.233498216), (
+0.518978477,
-0.810231030), (
+0.880361557,
-0.780780017), (
+0.823689997,
-0.632086337), (
+0.131901503,
-0.304824173)
508+8,
+8,
+4,
+7,
+8,
+4,
+2,
+9,
+6,
+1,
+4
509var(
0)
= getVar(sample, iweight)
513 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
514 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
516call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
519call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
528 lb(isam)
= ub(isam
- 1)
+ 1
537(
+0.845366180,
-0.746111989), (
+0.993096232E-1,
-0.645383596E-1), (
+0.989494085,
-0.441907108), (
+0.793754458E-1,
-0.808745325), (
+0.498732507,
-0.361979783), (
+0.263713598,
-0.359771371), (
+0.403559089,
-0.199794054), (
+0.664330661,
-0.364197969), (
+0.113191128,
-0.459438026)
540+6,
+2,
+3,
+5,
+5,
+2,
+4,
+9,
+9
541var(
0)
= getVar(sample, iweight)
545 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
546 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
548call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
551call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
560 lb(isam)
= ub(isam
- 1)
+ 1
569(
+0.681120157,
-0.561142266), (
+0.600067675,
-0.741244435), (
+0.685153008E-1,
-0.924964428), (
+0.905477643,
-0.878580391), (
+0.346492469,
-0.460470140), (
+0.837336302,
-0.440788269E-1), (
+0.678418398,
-0.112699866E-1), (
+0.254854977,
-0.881408811), (
+0.738232911,
-0.559942603), (
+0.570324421,
-0.901490748), (
+0.235462368,
-0.394087851), (
+0.855239868,
-0.552689850)
572+4,
+9,
+2,
+7,
+4,
+7,
+3,
+8,
+4,
+3,
+4,
+9
573var(
0)
= getVar(sample, iweight)
577 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
578 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
580call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
583call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
592 lb(isam)
= ub(isam
- 1)
+ 1
601(
+0.623735428,
-0.348751128), (
+0.141771138,
-0.579480410), (
+0.478169918,
-0.747342110), (
+0.696623206,
-0.304574430), (
+0.850297511,
-0.158303976E-1), (
+0.268961787,
-0.826235890), (
+0.421791315,
-0.218959808), (
+0.539953530,
-0.142040253), (
+0.501198530,
-0.444777429), (
+0.887419522,
-0.755490184), (
+0.679285288,
-0.228869796)
604+2,
+4,
+5,
+7,
+3,
+10,
+8,
+9,
+9,
+10,
+10
605var(
0)
= getVar(sample, iweight)
609 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
610 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
612call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
615call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
629 lb(isam)
= ub(isam
- 1)
+ 1
638+0.345471382,
+0.485727489,
+0.984534383,
+0.731319189,
+0.335763514,
+0.298995018,
+0.679663479,
+0.350297570,
+0.975026190,
+0.780901492,
+0.701982319,
+0.564911366,
+0.115896642
641+1.04257512,
+1.49429750,
+1.75459826,
+1.29751611,
+1.43655801,
+1.89791799,
+1.80559790,
+1.98860657,
+1.08174181,
+1.03304338,
+1.69656157,
+1.98976898,
+1.27701283
642var(
0)
= getVar(sample, rweight)
646 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
647 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
649call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
652call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
661 lb(isam)
= ub(isam
- 1)
+ 1
670+0.245393157,
+0.815335929,
+0.399940431,
+0.595761418,
+0.920845389,
+0.545531273
673+1.82141924,
+1.97131014,
+1.58715928,
+1.94621909,
+1.99396801,
+1.31281590
674var(
0)
= getVar(sample, rweight)
678 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
679 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
681call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
684call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
693 lb(isam)
= ub(isam
- 1)
+ 1
702+0.245913804,
+0.961104751,
+0.371187866,
+0.606853068,
+0.326335430E-2,
+0.593347907,
+0.484552979E-1,
+0.830431640,
+0.460434139,
+0.633784354
705+1.67454028,
+1.67554069,
+1.35289037,
+1.83354068,
+1.38854563,
+1.54515886,
+1.17616153,
+1.28903317,
+1.23424947,
+1.97806406
706var(
0)
= getVar(sample, rweight)
710 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
711 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
713call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
716call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
725 lb(isam)
= ub(isam
- 1)
+ 1
734+0.372254908,
+0.584946930,
+0.894830823,
+0.184308171,
+0.407962501,
+0.690055490
737+1.15169764,
+1.71129096,
+1.24037766,
+1.09339750,
+1.85883105,
+1.68955743
738var(
0)
= getVar(sample, rweight)
742 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
743 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
745call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
748call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
757 lb(isam)
= ub(isam
- 1)
+ 1
766+0.439015031,
+0.787458003,
+0.866388619,
+0.822199941,
+0.863702118,
+0.941337585,
+0.452802062,
+0.948480010,
+0.552669406
769+1.75023162,
+1.05153394,
+1.76741982,
+1.14104867,
+1.53707969,
+1.66175675,
+1.42000031,
+1.35965967,
+1.31342244
770var(
0)
= getVar(sample, rweight)
774 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
775 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
777call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
780call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
789 lb(isam)
= ub(isam
- 1)
+ 1
798(
+0.732771158E-1,
-0.363276601E-1), (
+0.423681021,
-0.384647727), (
+0.698369265,
-0.722870886), (
+0.888666511,
-0.409296513), (
+0.460444391,
-0.958284259)
801+1.83342385,
+1.65144396,
+1.53762484,
+1.82303226,
+1.22206521
802var(
0)
= getVar(sample, rweight)
806 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
807 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
809call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
812call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
821 lb(isam)
= ub(isam
- 1)
+ 1
830(
+0.115033925,
-0.477619410), (
+0.948059320,
-0.638965130), (
+0.256881893,
-0.571313798), (
+0.338890910,
-0.457587600), (
+0.429126322,
-0.853367925), (
+0.986338854E-1,
-0.423292696), (
+0.287814200,
-0.851023614), (
+0.359498382,
-0.873471677)
833+1.86694026,
+1.79731154,
+1.28275156,
+1.33540332,
+1.16399193,
+1.31400251,
+1.55825734,
+1.82298696
834var(
0)
= getVar(sample, rweight)
838 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
839 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
841call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
844call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
853 lb(isam)
= ub(isam
- 1)
+ 1
862(
+0.648760021,
-0.667374194), (
+0.544691563,
-0.913228512), (
+0.893469453E-1,
-0.996351361), (
+0.155209601,
-0.436864376), (
+0.148003876,
-0.527336597), (
+0.359869003,
-0.790384293), (
+0.330712914,
-0.365891874), (
+0.748180091,
-0.688499928), (
+0.168662131,
-0.428182244), (
+0.104938209,
-0.205892563), (
+0.931405127,
-0.669851661), (
+0.835194290,
-0.728791833)
865+1.98377371,
+1.22494948,
+1.45565784,
+1.07273340,
+1.23618960,
+1.75740647,
+1.95626235,
+1.29106581,
+1.00725722,
+1.06805372,
+1.56059408,
+1.44038749
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 setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
876call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
885 lb(isam)
= ub(isam
- 1)
+ 1
894(
+0.267716348,
-0.111458361), (
+0.212560236,
-0.540555537), (
+0.806218445,
-0.320773602), (
+0.138623834,
-0.656682968), (
+0.137818456E-1,
-0.407212257), (
+0.517320633E-1,
-0.899588168), (
+0.612607539,
-0.837301731), (
+0.468982339,
-0.549971819), (
+0.455741227,
-0.174955249), (
+0.420116067,
-0.343016148), (
+0.841679692,
-0.197088242)
897+1.86270618,
+1.89528823,
+1.81775546,
+1.39717209,
+1.38556528,
+1.64302206,
+1.88139200,
+1.48550129,
+1.32791328,
+1.86805677,
+1.95111132
898var(
0)
= getVar(sample, rweight)
902 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
903 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
905call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
908call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
917 lb(isam)
= ub(isam
- 1)
+ 1
926(
+0.581725895,
-0.196485758), (
+0.173325658,
-0.618625522), (
+0.159565389,
-0.522693634), (
+0.276407361,
-0.632627487), (
+0.764621377,
-0.837316155), (
+0.333611608,
-0.823956728E-1), (
+0.548387468,
-0.778472126), (
+0.823851168,
-0.336416721), (
+0.792722046,
-0.535879791)
929+1.97886717,
+1.56527865,
+1.32357621,
+1.37239552,
+1.76769733,
+1.59945011,
+1.95721519,
+1.23022389,
+1.34191442
930var(
0)
= getVar(sample, rweight)
934 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
935 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
937call setVarMerged(varMerged, var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
940call setVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
954 lb(isam)
= ub(isam
- 1)
+ 1
962call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
963call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
967-0.794386268,
-0.325577378,
+0.662713289,
+0.180385590,
-0.439923048,
-0.561428070,
-0.996729016,
+0.142635107E-1,
-0.157487512,
-0.967031837
968+0.190051913,
+0.947661400E-1,
-0.491064668,
+0.399760723,
+0.345922351,
-0.980438590,
-0.352344751,
-0.540076494E-1,
+0.332165956,
-0.534865141
969+0.999535441,
-0.878053188,
-0.390816092,
-0.744003057E-1,
+0.991106033,
+0.105976939,
-0.224736214,
+0.936365008,
+0.832479000E-1,
+0.137254000E-1
970+0.633215666,
+0.373368144,
+0.952402830,
+0.183950782,
-0.450929642,
+0.868240833,
-0.205165744,
+0.989269853,
+0.228426814,
-0.983087301
971var(:,
0)
= getVar(sample, dim)
973+0.252517849,
+0.194019303,
+0.361227870,
+0.380719155
975 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
976 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
978call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
980+0.252517849,
+0.194019288,
+0.361227870,
+0.380719155
981call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
983+0.252517849,
+0.194019288,
+0.361227870,
+0.380719155
985+0.252517849,
+0.194019303,
+0.361227870,
+0.380719155
990 lb(isam)
= ub(isam
- 1)
+ 1
998call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
999call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1003-0.938377142,
+0.155410767,
+0.140862226,
-0.298054576,
-0.197319269,
+0.177136660,
-0.966215134E-2,
+0.106565833,
+0.551277161,
+0.822041273,
+0.799075365,
-0.921103954E-1,
+0.973501205
1004-0.479531527,
-0.556624413,
+0.564124703,
+0.464079738,
+0.352094173E-1,
-0.609748721,
-0.814321518,
+0.104710698,
+0.845510721,
+0.176680088E-1,
+0.941410065E-1,
-0.547937512,
+0.668950677
1005var(:,
0)
= getVar(sample, dim)
1007+0.253881991,
+0.276060283
1009 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1010 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1012call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1014+0.253881991,
+0.276060283
1015call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1017+0.253881991,
+0.276060283
1019+0.253881991,
+0.276060283
1024 lb(isam)
= ub(isam
- 1)
+ 1
1032call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1033call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1037+0.262256742,
+0.326944947,
+0.917922258,
+0.538078070,
-0.210764289,
+0.823613048,
-0.299220443,
-0.170571446
1038var(:,
0)
= getVar(sample, dim)
1042 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1043 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1045call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1048call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1057 lb(isam)
= ub(isam
- 1)
+ 1
1065call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1066call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1070+0.643550515,
-0.891264558,
+0.447398901,
-0.870571613,
+0.927083969,
-0.677865624,
-0.132524729,
-0.825998783E-1,
-0.359449863,
+0.400723577,
+0.430679202
1071var(:,
0)
= getVar(sample, dim)
1075 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1076 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1078call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1081call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1090 lb(isam)
= ub(isam
- 1)
+ 1
1098call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1099call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1103+0.591102123,
+0.158055663,
-0.653327107,
+0.767853379,
+0.100939274E-1,
-0.740361810,
+0.951668024,
+0.518260241,
+0.815117359E-1,
-0.577678084,
-0.668644071
1104var(:,
0)
= getVar(sample, dim)
1108 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1109 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1111call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1114call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1123 lb(isam)
= ub(isam
- 1)
+ 1
1131call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1132call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1136(
-0.126476288E-1,
+0.912351608E-1), (
-0.256464481,
+0.613908529), (
+0.413810015,
+0.288720489), (
+0.204523206,
+0.715348601), (
+0.447022915E-1,
-0.995965838), (
+0.739053011,
-0.789644718), (
+0.133874893,
+0.537430048E-1), (
+0.105212927E-1,
+0.138975978), (
-0.571441174,
-0.597128034), (
+0.592389822,
-0.125735402), (
+0.315291286,
+0.409012198), (
-0.252358913,
+0.636641979), (
+0.883441925,
-0.331953049)
1137(
-0.502955794,
-0.464842200), (
+0.277886868,
-0.920887589), (
+0.887005091,
+0.352511406), (
-0.268300533,
+0.305666924E-1), (
-0.959437847,
+0.143795490), (
-0.867940426,
+0.295158148), (
-0.869491935,
-0.200284839), (
-0.164155722,
-0.831591129), (
+0.914741993,
-0.191862583E-1), (
-0.112508535,
+0.962959886), (
+0.612155199,
-0.238772631), (
-0.796514988,
+0.737625003), (
+0.532059550,
+0.678499818)
1138var(:,
0)
= getVar(sample, dim)
1140+0.442380220,
+0.742350042
1142 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1143 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1145call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1147+0.442380279,
+0.742350101
1148call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1150+0.442380279,
+0.742350101
1152+0.442380220,
+0.742350042
1157 lb(isam)
= ub(isam
- 1)
+ 1
1165call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1166call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1170(
+0.808130503,
-0.390991807), (
+0.338146329,
-0.618746161), (
+0.513154984,
-0.433777928), (
+0.115792871,
-0.444639325), (
-0.587761402,
-0.652993917), (
-0.539530277,
+0.969924569), (
-0.120231152,
+0.250126362), (
+0.977964401E-1,
+0.893638015), (
+0.239616513,
+0.757834792), (
-0.111284256,
-0.711120367), (
+0.290513515,
-0.188244820), (
-0.339172363,
-0.690890670)
1171(
+0.644578695,
+0.321047783), (
-0.735005021,
+0.856800079E-1), (
-0.232035279,
+0.433197021), (
+0.555115104,
+0.379326344E-1), (
+0.161392331,
-0.176887155), (
-0.320952535,
+0.506532431), (
+0.582489371,
-0.482062697), (
-0.345465660,
+0.221789360), (
-0.251729369,
-0.899287701), (
+0.630455971,
+0.305621266), (
-0.773149729,
-0.413843036), (
-0.493311524,
+0.381389499)
1172var(:,
0)
= getVar(sample, dim)
1174+0.545243382,
+0.438230515
1176 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1177 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1179call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1181+0.545243382,
+0.438230544
1182call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1184+0.545243382,
+0.438230544
1186+0.545243382,
+0.438230515
1191 lb(isam)
= ub(isam
- 1)
+ 1
1199call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1200call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1204(
-0.260365009E-1,
-0.481344819), (
+0.448194981,
-0.184910178), (
-0.572889686,
+0.477628112), (
-0.571240306,
+0.958866119), (
+0.556555271,
-0.372226238E-1), (
+0.969378710,
+0.769129515), (
-0.598190784,
-0.230672479), (
+0.936912179,
-0.196218729), (
-0.190650821,
+0.814921618)
1205var(:,
0)
= getVar(sample, dim)
1209 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1210 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1212call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1215call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1224 lb(isam)
= ub(isam
- 1)
+ 1
1232call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1233call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1237(
+0.548242569,
+0.817435503), (
-0.714928627,
+0.910810232E-1), (
+0.625148296,
-0.754720807), (
+0.516124845,
+0.329550385), (
+0.743469834,
+0.668153286), (
+0.605001569,
+0.373288989), (
+0.159842968,
-0.960388184E-1), (
+0.161159396,
+0.335026741), (
-0.358446360,
-0.178236485), (
+0.148166537,
+0.804740667), (
+0.306304574,
+0.520152569)
1238(
+0.581233978,
-0.576761484), (
+0.190872312,
+0.789755464), (
+0.489081144E-1,
-0.497144461), (
+0.645821333,
-0.113402486), (
+0.206929564,
-0.841662526), (
-0.871820450,
-0.791832328), (
-0.728684902,
+0.139729023), (
-0.869853258,
+0.293052077), (
-0.919104338,
-0.686519861), (
-0.400845766,
+0.998614788), (
-0.944448709E-1,
-0.147337794)
1239(
-0.765865088,
+0.999357939), (
-0.106465697,
-0.657376647), (
+0.200438976,
+0.125033259), (
+0.362173796,
+0.450037241), (
-0.760166049,
+0.633422136E-1), (
-0.885819197,
-0.740878224), (
-0.931785226,
-0.705699682), (
+0.620398521E-1,
-0.394546270), (
+0.688846469,
+0.603306651), (
-0.802195549,
+0.538636208), (
-0.325897217,
-0.636561036)
1240var(:,
0)
= getVar(sample, dim)
1242+0.385701269,
+0.674074352,
+0.651144922
1244 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1245 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1247call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1249+0.385701239,
+0.674074292,
+0.651144862
1250call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1252+0.385701239,
+0.674074292,
+0.651144862
1254+0.385701269,
+0.674074352,
+0.651144922
1259 lb(isam)
= ub(isam
- 1)
+ 1
1267call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1268call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1272(
+0.725114226,
-0.197496891), (
-0.710519433,
-0.430784702), (
+0.560350180,
-0.958535552), (
-0.543453932,
-0.204392076), (
+0.821977735,
+0.724628925), (
-0.467937350,
+0.312224865), (
+0.244568229,
+0.471392632)
1273(
-0.152471066E-1,
+0.249912262), (
+0.916773558,
+0.265466809), (
+0.818207145,
-0.819630027), (
-0.563442469,
-0.602471709), (
-0.327858925E-1,
-0.392822981), (
-0.450326204,
-0.927865505E-2), (
-0.467052698,
+0.502619624)
1274(
-0.808493733,
+0.584389329), (
-0.767815351,
+0.932148099), (
+0.144183397,
-0.133702278), (
+0.739225626,
-0.882296443), (
-0.756782293,
+0.406535029), (
+0.641397715,
+0.880883813), (
+0.807796717,
+0.821290135)
1275var(:,
0)
= getVar(sample, dim)
1277+0.650880337,
+0.532205462,
+0.870996356
1279 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1280 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1282call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1284+0.650880337,
+0.532205462,
+0.870996356
1285call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1287+0.650880337,
+0.532205462,
+0.870996356
1289+0.650880337,
+0.532205462,
+0.870996356
1299 lb(isam)
= ub(isam
- 1)
+ 1
1307call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1308call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1312-0.779642463,
+0.371027350,
+0.983049035,
+0.724993110,
+0.754500628E-1,
-0.395085812E-1,
+0.228537321,
+0.546303391
1313-0.514171004,
+0.443769574,
+0.689239621,
+0.910249829,
-0.611693978,
-0.966423988,
-0.661731243,
+0.298497200
1314-0.261814475,
-0.914427042E-1,
+0.331742883,
-0.883565307,
-0.164275289,
+0.419911146E-1,
+0.334653974,
-0.766766906
1317+8,
+6,
+9,
+2,
+3,
+4,
+3,
+5
1318var(:,
0)
= getVar(sample,
2_IK, iweight)
1320+0.368060052,
+0.396077842,
+0.149744704
1322 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1323 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1325call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1327+0.368060082,
+0.396077842,
+0.149744704
1328call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1330+0.368060082,
+0.396077842,
+0.149744704
1332+0.368060052,
+0.396077842,
+0.149744704
1337 lb(isam)
= ub(isam
- 1)
+ 1
1345call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1346call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1350+0.790132046,
+0.995968819,
-0.117847443,
-0.196304917
1354var(:,
0)
= getVar(sample,
2_IK, iweight)
1358 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1359 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1361call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1364call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1373 lb(isam)
= ub(isam
- 1)
+ 1
1381call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1382call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1386+0.880890489,
-0.225948930,
+0.695374489,
-0.715105891,
-0.650427461,
-0.210168004,
-0.793861270,
+0.989555955,
+0.521224380,
-0.521412969,
+0.255403876,
-0.745594501E-2,
-0.667250633
1387-0.952096462,
+0.294235349,
-0.127303600,
+0.655461192,
+0.755432963,
+0.754771471,
+0.330180645,
+0.184328914,
-0.461623788,
-0.647416949,
+0.923807383,
+0.682372689,
-0.703441501
1388-0.774365664E-1,
+0.482311606,
-0.224624634,
+0.532164693,
+0.253500819,
-0.967351794,
-0.606431961E-1,
-0.901088357,
-0.880862355,
-0.675177693,
+0.315366030,
+0.376101136,
+0.136048079
1391+6,
+10,
+3,
+6,
+1,
+2,
+2,
+4,
+8,
+6,
+6,
+7,
+7
1392var(:,
0)
= getVar(sample,
2_IK, iweight)
1394+0.333630919,
+0.396812975,
+0.291316450
1396 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1397 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1399call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1401+0.333630919,
+0.396813065,
+0.291316479
1402call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1404+0.333630919,
+0.396813065,
+0.291316479
1406+0.333630919,
+0.396812975,
+0.291316450
1411 lb(isam)
= ub(isam
- 1)
+ 1
1419call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1420call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1424-0.870981216,
-0.293064356,
-0.476996541,
+0.612924337,
-0.142254233,
+0.146249533E-1,
+0.483549595,
-0.191353083
1425+0.610053539E-1,
-0.745706081,
-0.313334703,
+0.648646593,
-0.636996031E-1,
+0.152570724,
-0.305884004,
-0.939720511
1426+0.205127716,
-0.928659797,
+0.393368483,
+0.923581481,
+0.225257874E-2,
-0.596172810,
-0.672949672,
+0.180585384
1429+6,
+3,
+10,
+8,
+6,
+4,
+4,
+10
1430var(:,
0)
= getVar(sample,
2_IK, iweight)
1432+0.208262056,
+0.261129141,
+0.273321450
1434 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1435 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1437call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1439+0.208262071,
+0.261129111,
+0.273321420
1440call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1442+0.208262071,
+0.261129111,
+0.273321420
1444+0.208262056,
+0.261129141,
+0.273321450
1449 lb(isam)
= ub(isam
- 1)
+ 1
1457call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1458call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1462+0.543098927,
+0.470154405,
-0.573689699,
+0.151789427,
+0.116292834,
+0.306769252,
+0.586919188,
+0.474532962,
+0.681506157
1465+3,
+3,
+8,
+6,
+2,
+2,
+5,
+10,
+9
1466var(:,
0)
= getVar(sample,
2_IK, iweight)
1470 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1471 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1473call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1476call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1485 lb(isam)
= ub(isam
- 1)
+ 1
1493call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1494call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1498(
-0.110805631,
+0.439788103E-1), (
-0.408427835,
+0.661906719), (
+0.763092518,
-0.279502869E-1), (
-0.907828569,
-0.599026680E-3), (
+0.925187469,
-0.242151618), (
-0.469572783,
-0.868839622), (
+0.841831684,
+0.544157147), (
+0.764955044,
+0.999729037), (
-0.830718398,
+0.848671794), (
-0.405682564,
+0.131775856), (
-0.347328067,
+0.254027843)
1499(
+0.805351973,
-0.249896288), (
+0.350865483,
-0.575051188), (
-0.677763820,
+0.253325582), (
+0.190366507E-1,
+0.110494018), (
+0.895068169,
+0.574787736), (
+0.285439491,
+0.768496037), (
-0.506423950,
+0.479465127), (
+0.452468157,
+0.881928563), (
-0.155337095,
-0.144932508), (
-0.838672400,
+0.782617092), (
+0.104426503,
+0.787577391)
1500(
+0.798611879,
+0.723740339), (
-0.636818051,
-0.709580898), (
-0.337017894,
-0.815467119), (
+0.905082941,
-0.958695173), (
-0.132172585,
-0.349843621), (
-0.648707986,
-0.303443193), (
+0.444703221,
-0.297963619), (
+0.683942914,
-0.484762311), (
+0.413931131,
+0.464705110), (
+0.926955462,
+0.775889754), (
-0.363178015,
-0.409186363)
1503+3,
+8,
+10,
+9,
+3,
+10,
+10,
+4,
+1,
+7,
+5
1504var(:,
0)
= getVar(sample,
2_IK, iweight)
1506+0.709448814,
+0.472461492,
+0.678329170
1508 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1509 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1511call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1513+0.709448814,
+0.472461462,
+0.678329110
1514call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1516+0.709448814,
+0.472461462,
+0.678329110
1518+0.709448814,
+0.472461492,
+0.678329170
1523 lb(isam)
= ub(isam
- 1)
+ 1
1531call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1532call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1536(
+0.575099587,
-0.471512079), (
+0.165558100,
-0.468282700E-1), (
-0.658837199,
+0.264920592), (
+0.451726437,
+0.967319965), (
+0.342851877E-1,
-0.605301261), (
-0.115012884,
-0.474712253), (
+0.317890763,
-0.746603370), (
+0.775372744,
-0.376402974), (
+0.577189803,
+0.929316282), (
-0.736213923,
-0.590994477)
1537(
+0.470925212,
-0.374642253), (
-0.166514754,
+0.990130663), (
+0.510759592,
+0.173046350), (
-0.484069705,
+0.586073756), (
-0.619891047,
+0.556636572), (
+0.611860633,
+0.454061270), (
-0.555001736,
+0.942466259E-1), (
+0.903034806,
-0.988619804), (
+0.890161157,
+0.940995455), (
-0.975872040,
+0.471503973)
1538(
+0.451868176,
-0.790648222), (
-0.728680253,
-0.599483848), (
-0.760388374E-2,
-0.214623213E-1), (
-0.425395846,
+0.865186930), (
-0.115696669,
-0.113111734), (
-0.212013364,
-0.849994540), (
+0.980087042,
+0.210381150), (
-0.389831781,
-0.978950143), (
+0.751599908,
+0.954835415), (
+0.231431603,
+0.429676414)
1539(
-0.851721883,
-0.968904972), (
+0.277927518,
+0.457118273), (
-0.167479038,
+0.226785779), (
+0.215460896,
+0.448988795), (
-0.153288841E-1,
-0.977933407E-1), (
-0.705790520E-2,
-0.917219281), (
+0.718338251,
-0.527290583), (
-0.595315814,
-0.355173588), (
+0.676217318,
+0.833532214), (
-0.888010144,
-0.279803038)
1542+10,
+4,
+1,
+1,
+6,
+6,
+8,
+7,
+2,
+10
1543var(:,
0)
= getVar(sample,
2_IK, iweight)
1545+0.404116631,
+0.823002040,
+0.614940166,
+0.591328502
1547 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1548 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1550call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1552+0.404116660,
+0.823002040,
+0.614940166,
+0.591328502
1553call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1555+0.404116660,
+0.823002040,
+0.614940166,
+0.591328502
1557+0.404116631,
+0.823002040,
+0.614940166,
+0.591328502
1562 lb(isam)
= ub(isam
- 1)
+ 1
1570call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1571call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1575(
-0.577630997E-1,
-0.467085600), (
+0.229264498,
-0.222207665), (
+0.896920800,
+0.581840277), (
-0.696803451,
-0.667212129), (
+0.356100321,
-0.512730718), (
-0.464267969,
-0.433839083), (
-0.269405484,
+0.604135036), (
-0.981544495,
+0.772737384), (
-0.383558869,
+0.820497155), (
-0.307947040,
-0.910818934), (
-0.198361635,
+0.624143124)
1576(
-0.360828519,
-0.514466643), (
+0.378947377,
-0.639093757), (
+0.587266088,
+0.482337713), (
+0.183734059,
-0.765872002E-2), (
+0.784380078,
+0.158514261), (
-0.862934589E-1,
-0.717438698), (
-0.610728383,
-0.822508335E-2), (
+0.824229479,
-0.541968346), (
-0.239370227,
+0.143245697), (
-0.950540900,
-0.675868988E-2), (
-0.759616375,
+0.329940557)
1579+5,
+3,
+9,
+5,
+7,
+5,
+5,
+1,
+5,
+9,
+3
1580var(:,
0)
= getVar(sample,
2_IK, iweight)
1582+0.665977597,
+0.512763381
1584 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1585 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1587call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1589+0.665977657,
+0.512763381
1590call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1592+0.665977657,
+0.512763381
1594+0.665977597,
+0.512763381
1599 lb(isam)
= ub(isam
- 1)
+ 1
1607call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1608call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1612(
-0.661615491,
+0.133234262), (
+0.615461707,
-0.850202799), (
+0.737609029,
-0.233716607), (
+0.835180759,
+0.689581156), (
+0.416536570,
+0.848124862), (
-0.387149096,
+0.902991295E-1), (
-0.431600571,
-0.740782976), (
-0.367417812,
+0.453714728), (
-0.771054387,
+0.427305102)
1613(
+0.645748377,
+0.166920424E-1), (
-0.424215794E-1,
-0.673459411), (
+0.171316266,
+0.828112125), (
-0.778745174,
+0.665221214E-1), (
+0.499007702,
+0.330409408), (
-0.949708223,
-0.432239652), (
+0.992193222E-1,
-0.925495744), (
-0.186007142,
-0.913597822), (
+0.357642651,
+0.451224208)
1616+9,
+1,
+7,
+5,
+6,
+5,
+10,
+5,
+4
1617var(:,
0)
= getVar(sample,
2_IK, iweight)
1619+0.630700648,
+0.638704479
1621 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1622 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1624call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1626+0.630700588,
+0.638704479
1627call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1629+0.630700588,
+0.638704479
1631+0.630700648,
+0.638704479
1636 lb(isam)
= ub(isam
- 1)
+ 1
1644call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1645call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1649(
-0.533411860,
-0.463769436), (
-0.419325590,
+0.961975336), (
+0.571228266,
-0.208961368), (
+0.468949318,
-0.235208035), (
+0.972064614,
+0.319884419), (
+0.421098590,
+0.627359748), (
-0.563890576,
-0.351084948), (
+0.652729869,
+0.930560827E-1), (
-0.895299792,
-0.917794704E-1), (
+0.228907585,
-0.995567203)
1650(
+0.147814989,
+0.318635464), (
+0.644891143,
+0.670288801), (
+0.967107773,
+0.634981751), (
+0.627987385,
+0.374973893), (
-0.834370852E-1,
+0.407241583), (
+0.437807202,
+0.906183481), (
+0.960109353,
+0.387710214), (
+0.213797331,
-0.320132852), (
-0.686890960,
-0.962248683), (
-0.908545136,
+0.853506327)
1651(
-0.964905500,
+0.780199647), (
-0.713003039,
-0.394178391), (
-0.419936061,
+0.179493308), (
-0.630067468,
+0.271220446), (
-0.901868463,
+0.133824468), (
-0.216668963,
+0.429481626), (
+0.978754640,
+0.908523679), (
-0.839982986,
-0.903547645), (
-0.133272886,
+0.731321454), (
+0.961488485,
+0.928944707)
1652(
-0.843516707,
+0.706558824), (
-0.407084465,
+0.292482615), (
-0.658698678,
+0.465146542), (
-0.910426259,
-0.161078095), (
+0.418271542,
+0.997216582), (
+0.527197719,
+0.391407013), (
+0.371357918,
-0.138603210), (
-0.264609814,
+0.761967421), (
-0.214399099,
-0.124849081E-1), (
-0.219334722,
+0.839722157E-2)
1655+2,
+3,
+3,
+3,
+9,
+8,
+10,
+3,
+9,
+5
1656var(:,
0)
= getVar(sample,
2_IK, iweight)
1658+0.697316885,
+0.798430622,
+0.756463230,
+0.383753330
1660 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1661 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1663call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1665+0.697316885,
+0.798430681,
+0.756463230,
+0.383753359
1666call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1668+0.697316885,
+0.798430681,
+0.756463230,
+0.383753359
1670+0.697316885,
+0.798430622,
+0.756463230,
+0.383753330
1680 lb(isam)
= ub(isam
- 1)
+ 1
1688call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1689call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1693-0.294872522E-1,
+0.536496997,
-0.259611011,
+0.925314546,
+0.271661162
1696+1.09982061,
+1.00587916,
+1.34108925,
+1.43254757,
+1.70560205
1697var(:,
0)
= getVar(sample,
2_IK, rweight)
1701 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1702 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1704call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1707call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1716 lb(isam)
= ub(isam
- 1)
+ 1
1724call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1725call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1729+0.420074224,
-0.892275095,
+0.426792622,
+0.376040459,
+0.320909858,
+0.618546009E-1,
+0.970496058,
-0.543833852
1732+1.57970023,
+1.44116843,
+1.31975746,
+1.38179147,
+1.39179289,
+1.61909580,
+1.54918575,
+1.42732441
1733var(:,
0)
= getVar(sample,
2_IK, rweight)
1737 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1738 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1740call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1743call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1752 lb(isam)
= ub(isam
- 1)
+ 1
1760call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1761call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1765+0.144684315E-1,
-0.536210418,
-0.258387208,
-0.239539146E-2,
-0.330352783,
+0.758076787,
+0.949331999,
+0.208339691E-1,
+0.683854461,
+0.550213337,
+0.201261401,
-0.459007025E-1,
+0.515185595
1766+0.778736472,
-0.917596817E-1,
+0.267517686,
-0.250741005,
-0.339260101E-1,
+0.525115848,
-0.230787158,
-0.903751612,
+0.979999900,
-0.249299645,
+0.684985399,
-0.758591890,
+0.836594224
1767+0.394351959,
+0.909280539,
-0.593260646,
-0.124391675,
+0.259449840,
-0.432466984,
+0.812274814,
+0.914189816,
-0.690861344,
-0.929305077,
-0.467175603,
-0.614592314,
+0.689827800
1768+0.348502398,
-0.968733430,
+0.823031425,
-0.679515123,
+0.109736443,
-0.933242202,
+0.596804857,
-0.963299990,
-0.980430484,
-0.952750802,
-0.137513876,
-0.601627827E-1,
+0.684206486
1769+0.270177007,
+0.329240441,
+0.990133762,
-0.927620530,
+0.809823275E-1,
+0.424442410,
+0.130565763,
+0.580305099,
-0.693318844,
-0.639566183E-1,
+0.418526053,
+0.404731631,
-0.409157276
1770-0.729846954E-1,
-0.956584215E-1,
-0.199806333,
-0.903074384,
+0.630893588,
+0.139374852,
-0.705280185,
-0.338845253E-1,
-0.311741352,
+0.424057245E-1,
-0.710073233,
-0.315279245,
-0.382762790
1773+1.82331169,
+1.96306002,
+1.72512794,
+1.79210806,
+1.89955807,
+1.12479019,
+1.32384932,
+1.26713276,
+1.61302137,
+1.38932955,
+1.58290017,
+1.87146449,
+1.78625298
1774var(:,
0)
= getVar(sample,
2_IK, rweight)
1776+0.189214170,
+0.337006986,
+0.409951091,
+0.448933899,
+0.268291026,
+0.157651469
1778 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1779 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1781call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1783+0.189214185,
+0.337007016,
+0.409951121,
+0.448933870,
+0.268290997,
+0.157651514
1784call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1786+0.189214185,
+0.337007016,
+0.409951121,
+0.448933870,
+0.268290997,
+0.157651514
1788+0.189214170,
+0.337006986,
+0.409951091,
+0.448933899,
+0.268291026,
+0.157651469
1793 lb(isam)
= ub(isam
- 1)
+ 1
1801call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1802call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1806+0.929083824,
+0.220651746,
-0.921595097,
-0.616780877,
-0.756878138,
+0.697091460,
-0.441668153
1807-0.753165483E-1,
+0.732868910,
-0.658041239E-1,
+0.801245809,
-0.230922699,
-0.145011902,
+0.565424323
1810+1.79460859,
+1.08830285,
+1.17355311,
+1.95332336,
+1.09231496,
+1.48970294,
+1.07883644
1811var(:,
0)
= getVar(sample,
2_IK, rweight)
1813+0.502466500,
+0.180449650
1815 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1816 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1818call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1820+0.502466559,
+0.180449665
1821call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1823+0.502466559,
+0.180449665
1825+0.502466500,
+0.180449650
1830 lb(isam)
= ub(isam
- 1)
+ 1
1838call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1839call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1843-0.987186193,
-0.218339920,
+0.429356575,
+0.136347055,
+0.817507148,
+0.541962266,
-0.132597804,
-0.758635402,
-0.193426609E-1
1846+1.31794667,
+1.63985586,
+1.52921987,
+1.84324276,
+1.55274868,
+1.05074430,
+1.69497228,
+1.90544558,
+1.02734423
1847var(:,
0)
= getVar(sample,
2_IK, rweight)
1851 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1852 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1854call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1857call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1866 lb(isam)
= ub(isam
- 1)
+ 1
1874call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1875call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1879(
+0.897367001E-1,
+0.373042822E-1), (
+0.934708595,
-0.476573348), (
-0.989857554,
-0.377285838), (
+0.194489837,
+0.339296699), (
-0.413433909,
+0.244913220), (
-0.181436419,
-0.851235747), (
-0.744381309,
-0.117990971E-1), (
-0.561195850,
+0.125545263E-1), (
-0.393722057E-1,
-0.810455799), (
-0.538261652,
-0.831589699), (
+0.447227955,
-0.219803929)
1880(
-0.829147696,
-0.148072720), (
+0.219634175,
-0.111380696), (
-0.766499758,
+0.982205629), (
+0.831144810,
+0.210317016), (
+0.872107506,
+0.118392706), (
-0.716125727,
-0.613348246), (
-0.308665633,
+0.285043240), (
-0.694833398,
-0.439177155), (
+0.719134092,
-0.163181663), (
+0.246385455,
+0.510793686), (
-0.255047083,
+0.972297072)
1881(
-0.957393169,
-0.128304243), (
-0.573877215,
+0.396553874), (
+0.314751148,
-0.282728672E-2), (
+0.768921137,
-0.179357529), (
+0.162391543,
-0.263881683), (
-0.784497261E-1,
+0.740525603), (
-0.733306408,
-0.756566525), (
+0.953020453,
+0.483871818), (
-0.440106869,
-0.979465723), (
-0.839006424,
-0.956128240), (
+0.126659870E-1,
-0.856366038)
1884+1.44226301,
+1.19237030,
+1.08518887,
+1.08148789,
+1.71804023,
+1.03792644,
+1.25479221,
+1.52199984,
+1.34539700,
+1.23475838,
+1.84177411
1885var(:,
0)
= getVar(sample,
2_IK, rweight)
1887+0.436836034,
+0.653117299,
+0.694399536
1889 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1890 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1892call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1894+0.436836064,
+0.653117359,
+0.694399655
1895call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1897+0.436836064,
+0.653117359,
+0.694399655
1899+0.436836034,
+0.653117299,
+0.694399536
1904 lb(isam)
= ub(isam
- 1)
+ 1
1912call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1913call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1917(
+0.891407251,
-0.303508401), (
-0.470538139E-1,
+0.883267045), (
-0.800948739,
+0.801616907), (
+0.898117542,
+0.298117518), (
+0.728914499,
-0.370856881), (
+0.180066705,
+0.231957912), (
+0.681900740,
-0.899252176), (
+0.237375021,
+0.170433640), (
-0.429615855,
+0.561780930), (
-0.630203962,
+0.558818817), (
-0.861097574,
+0.458370209)
1918(
+0.918472648,
-0.517752171), (
+0.508970737,
-0.787662983), (
-0.421175838,
+0.928377390), (
-0.509467721,
+0.975701690), (
+0.369640470,
+0.513610959), (
-0.167221904,
-0.500504136), (
+0.764301896,
+0.455288649), (
+0.459194064,
+0.152505755), (
-0.779030085,
+0.125599980), (
+0.733575583,
-0.347634196), (
-0.393433928,
+0.413806200)
1921+1.95097530,
+1.76496887,
+1.08978975,
+1.86465764,
+1.35181022,
+1.11587453,
+1.75375426,
+1.96317148,
+1.31884217,
+1.13044310,
+1.35277224
1922var(:,
0)
= getVar(sample,
2_IK, rweight)
1924+0.680343688,
+0.660446346
1926 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1927 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1929call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1931+0.680343747,
+0.660446405
1932call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1934+0.680343747,
+0.660446405
1936+0.680343688,
+0.660446346
1941 lb(isam)
= ub(isam
- 1)
+ 1
1949call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1950call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1954(
+0.542622209,
-0.695941329), (
-0.272163749,
+0.176460981), (
+0.834719539,
+0.395533204), (
-0.206228733,
-0.282016754), (
+0.344046235,
+0.908018947), (
+0.793401122,
-0.153975606), (
-0.114643574E-1,
+0.660843015), (
-0.585276365,
+0.931176782)
1955(
-0.900015235,
-0.835314393), (
+0.501311183,
+0.977442384), (
-0.924220562,
-0.333323717), (
+0.288808942,
+0.866656780), (
-0.867475510,
+0.883567691), (
+0.570622444,
+0.308579206), (
-0.925772548,
-0.805828214), (
+0.443373919E-1,
+0.942093372)
1958+1.70636415,
+1.79143858,
+1.09774244,
+1.42500544,
+1.98922646,
+1.47903538,
+1.96072054,
+1.52217817
1959var(:,
0)
= getVar(sample,
2_IK, rweight)
1961+0.534890652,
+1.00125861
1963 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1964 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1966call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1968+0.534890711,
+1.00125849
1969call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1971+0.534890711,
+1.00125849
1973+0.534890652,
+1.00125861
1978 lb(isam)
= ub(isam
- 1)
+ 1
1986call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1987call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1991(
+0.179611444E-1,
-0.557123303), (
+0.236685276,
+0.385071158), (
+0.526175976,
+0.412249565E-2), (
-0.308982849,
-0.638037920), (
-0.297726631,
-0.428469300), (
-0.453162670,
+0.853210568), (
-0.384537578,
-0.854053140), (
-0.389002085,
-0.772761106), (
-0.910327315,
-0.531966209), (
+0.897499323,
+0.131095290)
1994+1.16543996,
+1.51121092,
+1.10337925,
+1.04596114,
+1.75680780,
+1.55163360,
+1.32653379,
+1.60578036,
+1.68373585,
+1.63415122
1995var(:,
0)
= getVar(sample,
2_IK, rweight)
1999 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2000 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2002call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2005call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2014 lb(isam)
= ub(isam
- 1)
+ 1
2022call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2023call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
2027(
-0.342988253,
-0.356286168), (
-0.388448715,
-0.787526011), (
-0.970162868,
+0.244405985), (
-0.336713076,
+0.547979236), (
-0.892581463,
-0.895258307), (
-0.686578393,
-0.241750360), (
-0.779499650,
-0.794441700)
2028(
+0.687131763,
+0.215635896), (
+0.204660892,
-0.475227594), (
+0.601498604,
+0.304438353), (
-0.406438947,
-0.424983978), (
-0.196833730,
-0.402655244), (
-0.747307181,
+0.447738886), (
-0.974823475,
+0.923445106)
2031+1.38008273,
+1.43141520,
+1.28345358,
+1.25248182,
+1.21455884,
+1.58050847,
+1.83932710
2032var(:,
0)
= getVar(sample,
2_IK, rweight)
2034+0.306904405,
+0.639992476
2036 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2037 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2039call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2041+0.306904435,
+0.639992416
2042call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2044+0.306904435,
+0.639992416
2046+0.306904405,
+0.639992476