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.164151967,
+0.699272633,
+0.862817764,
+0.376528263,
+0.534596741,
+0.294695497E-1,
+0.889389098,
+0.715736091,
+0.475658476,
+0.611594260
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.434075296,
+0.330543101,
+0.773452222,
+0.254093409,
+0.986161292,
+0.615804732
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.183619320,
+0.253177941,
+0.813970685,
+0.469569981,
+0.426874161E-1,
+0.665751100E-1,
+0.445553541,
+0.119808137
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.681556582,
+0.781108320,
+0.188434124E-1,
+0.719868779,
+0.230368555,
+0.232569277,
+0.380671501
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.717068911,
+0.940654099,
+0.183019102,
+0.757230997,
+0.475270569,
+0.127272129,
+0.813057840,
+0.927254856,
+0.639357269,
+0.753042042,
+0.326402128,
+0.532595575,
+0.513611674
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.879619598,
-0.386423469E-1), (
+0.315917373,
-0.232920110), (
+0.723547399,
-0.756328940), (
+0.802050054,
-0.297957420), (
+0.857516170,
-0.642530799), (
+0.641659975,
-0.441883683), (
+0.117792904,
-0.161728621)
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.624335170,
-0.132989883E-1), (
+0.987618148,
-0.566642046), (
+0.285093844,
-0.234333992), (
+0.233627439,
-0.244343936), (
+0.474358678,
-0.670908153), (
+0.639385581,
-0.400373220), (
+0.645031869,
-0.478297472E-1), (
+0.774315715,
-0.427086115), (
+0.609079778,
-0.974845648), (
+0.878086686E-1,
-0.735144198), (
+0.911273599,
-0.124685645), (
+0.444563866,
-0.395017862E-2), (
+0.581902266,
-0.701951981E-2)
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.359709203,
-0.447478652), (
+0.350123107,
-0.605803728E-2), (
+0.715538919,
-0.527207315), (
+0.337770581E-1,
-0.777663052), (
+0.302694142,
-0.891901553), (
+0.265041113,
-0.134768605)
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.294067085,
-0.681287169), (
+0.259248435,
-0.364429355E-1), (
+0.681817532E-1,
-0.639215946), (
+0.333932161,
-0.428721488)
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.529882908,
-0.356199145), (
+0.259251833,
-0.933011949), (
+0.320232272,
-0.267866015), (
+0.985427797,
-0.399567246)
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.972808599,
+0.112576008,
+0.745032430E-1,
+0.421082973E-2,
+0.930948377,
+0.307939649,
+0.235144258,
+0.886963248,
+0.804508686,
+0.664887667,
+0.701334357,
+0.855493009
316+10,
+8,
+7,
+3,
+9,
+9,
+10,
+4,
+6,
+7,
+6,
+6
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.887729228,
+0.636436760,
+0.470908344,
+0.320635200,
+0.819119155,
+0.535648823,
+0.703419447,
+0.265076697,
+0.748179793,
+0.847605288
348+6,
+8,
+7,
+1,
+4,
+7,
+10,
+10,
+9,
+2
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.239746094,
+0.202776194E-1,
+0.877468109,
+0.448650301,
+0.306854367,
+0.809815407,
+0.934142530,
+0.937414169,
+0.670443952,
+0.415828824,
+0.182000399E-1,
+0.150106728
380+7,
+10,
+5,
+4,
+5,
+3,
+9,
+10,
+2,
+8,
+8,
+7
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.750622749E-1,
+0.621466041,
+0.694541991,
+0.788949668,
+0.811946869,
+0.792292118,
+0.114213109,
+0.983950377,
+0.857075453E-1,
+0.788395405E-1
412+7,
+3,
+3,
+2,
+7,
+7,
+1,
+8,
+9,
+10
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.850398779,
+0.534661055,
+0.107200682,
+0.988739252,
+0.327584565,
+0.226667047
444+5,
+4,
+1,
+10,
+5,
+2
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.753573120,
-0.205585778), (
+0.588616014,
-0.349650562), (
+0.758927584,
-0.719285011E-1), (
+0.519736707,
-0.957044363), (
+0.171725452,
-0.674497306), (
+0.591776192,
-0.467490256), (
+0.988852501,
-0.184851646), (
+0.734142959,
-0.841705441), (
+0.288699746,
-0.757220507)
476+8,
+1,
+4,
+7,
+1,
+10,
+9,
+8,
+6
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.828162432,
-0.540539920), (
+0.862169683,
-0.111547709E-1), (
+0.494877040,
-0.217702150), (
+0.944724858,
-0.791101575), (
+0.260692835E-2,
-0.587025881E-1), (
+0.366137445,
-0.828595340), (
+0.719009876,
-0.146787107), (
+0.413095415,
-0.503180265), (
+0.726852417E-1,
-0.372728705)
508+6,
+7,
+1,
+8,
+1,
+4,
+3,
+2,
+1
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.389677763,
-0.630508661E-1), (
+0.207754016,
-0.319278359), (
+0.232211649,
-0.636393428), (
+0.163720429,
-0.785705566), (
+0.862738192,
-0.505026460), (
+0.315562487,
-0.887426078), (
+0.561940551,
-0.252704620), (
+0.444091499,
-0.446763635E-1), (
+0.367281854,
-0.108294845), (
+0.944713354E-1,
-0.822873354), (
+0.326592028,
-0.391823649), (
+0.279402971,
-0.146980762)
540+5,
+8,
+9,
+9,
+8,
+6,
+4,
+2,
+8,
+2,
+7,
+7
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.496831298,
-0.943834424), (
+0.239272177,
-0.851670563), (
+0.259641171,
-0.214396656), (
+0.844989777,
-0.201076686), (
+0.379852235,
-0.759978116), (
+0.546731710,
-0.559752047), (
+0.961523235,
-0.259603262E-1)
572+8,
+2,
+6,
+1,
+1,
+2,
+6
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.205968797,
-0.937522531), (
+0.779543877,
-0.250976980), (
+0.277280211,
-0.796820819), (
+0.722867966,
-0.931649923), (
+0.459817529,
-0.106618345), (
+0.731259227,
-0.623011410), (
+0.544086158,
-0.988216400E-1), (
+0.365303397,
-0.888077796), (
+0.605010152,
-0.193560719E-1), (
+0.190733075,
-0.425789535), (
+0.588288009,
-0.860783219)
604+6,
+9,
+8,
+1,
+3,
+4,
+2,
+8,
+10,
+7,
+6
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.522112548,
+0.745396018E-1,
+0.848525703,
+0.216225803,
+0.734398484,
+0.972545028,
+0.393302679,
+0.138874829,
+0.719871163
641+1.59913039,
+1.24323881,
+1.34636045,
+1.86488664,
+1.37167716,
+1.39705443,
+1.18799472,
+1.43345761,
+1.27433681
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.289219141,
+0.287060022,
+0.970878839,
+0.551008344,
+0.729447126,
+0.899676502,
+0.837872207
673+1.62280560,
+1.41476130,
+1.79379821,
+1.90319538,
+1.42796981,
+1.43536758,
+1.81797791
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.413454175,
+0.831113577,
+0.534334242,
+0.816870451,
+0.680686951,
+0.289376438,
+0.229618847,
+0.237056673,
+0.285225511E-1,
+0.480395913,
+0.249648213
705+1.82679033,
+1.66666186,
+1.05207324,
+1.13198483,
+1.35660696,
+1.98397875,
+1.27012563,
+1.31911206,
+1.28169096,
+1.06015491,
+1.44365478
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.930402756,
+0.907879710,
+0.993824005E-1,
+0.561941504,
+0.491538048E-1,
+0.411397696,
+0.367166162,
+0.961572289,
+0.924229622E-3
737+1.66278017,
+1.89248300,
+1.75872231,
+1.06546760,
+1.42512679,
+1.60996890,
+1.78678441,
+1.58781862,
+1.02365780
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.521515906,
+0.283242226,
+0.489792824,
+0.577539802E-1,
+0.495089650,
+0.768620431,
+0.854355156,
+0.256272912,
+0.542074740,
+0.959119678,
+0.162405491
769+1.37470174,
+1.54892147,
+1.73006225,
+1.83361077,
+1.52070498,
+1.53563273,
+1.97179508,
+1.68564868,
+1.56750321,
+1.75964260,
+1.28962326
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.187669277,
-0.294457674), (
+0.131798446,
-0.977987409), (
+0.417990625,
-0.681630909), (
+0.195452213,
-0.912569404), (
+0.684241712,
-0.982518971), (
+0.657462597,
-0.902616739), (
+0.430532873,
-0.209922612)
801+1.17799830,
+1.59848356,
+1.26334071,
+1.64681363,
+1.75805664,
+1.11584365,
+1.50780797
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.257757664,
-0.653854191), (
+0.467664778,
-0.808455586), (
+0.115035176E-1,
-0.624481857), (
+0.229232669,
-0.688477397), (
+0.429523230,
-0.307994723), (
+0.977592111,
-0.574360251), (
+0.269044757,
-0.100793242E-1), (
+0.753784180,
-0.449637890), (
+0.269126534,
-0.674464285), (
+0.433862627,
-0.610498667), (
+0.652330577,
-0.761833847), (
+0.622174978,
-0.484607458)
833+1.07039714,
+1.67559683,
+1.44247675,
+1.48888659,
+1.09602714,
+1.40373707,
+1.45310903,
+1.62354350,
+1.68173528,
+1.29635501,
+1.94903636,
+1.27637291
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.144823611,
-0.353288710), (
+0.840760231,
-0.141086400), (
+0.390878677,
-0.581986129), (
+0.783770502,
-0.263178945), (
+0.731932521E-1,
-0.816559553), (
+0.408051968,
-0.871779025), (
+0.439551771,
-0.803559959), (
+0.303816915,
-0.804044962)
865+1.83297598,
+1.30447459,
+1.00364125,
+1.86158073,
+1.60139298,
+1.48901224,
+1.72583795,
+1.22281647
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.809294760,
-0.540455997), (
+0.138245285,
-0.147022605), (
+0.938309848,
-0.708187222E-1), (
+0.725421965,
-0.907560110), (
+0.216003120,
-0.364010930E-1)
897+1.27158403,
+1.16790915,
+1.96443295,
+1.18513393,
+1.37741137
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.102448106,
-0.441954315), (
+0.417125225E-1,
-0.115974545E-1), (
+0.359860063E-1,
-0.324419618), (
+0.929557323,
-0.703358173), (
+0.776709795,
-0.928337038)
929+1.25677538,
+1.03377628,
+1.97580612,
+1.08447456,
+1.59794366
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.902274489,
-0.648801923,
+0.145853758,
-0.249810338,
+0.917064190,
+0.510713816,
-0.962661624,
+0.375113130,
-0.338170648,
+0.776476741,
-0.987750292
968+0.249561548,
-0.489351034,
-0.220516443,
+0.153056979,
+0.222944021E-1,
+0.232488394,
+0.760874510,
+0.293239236,
-0.925330043,
+0.114853740,
-0.964010239
969+0.816024661,
+0.461724997E-1,
-0.819743872E-1,
+0.259068012E-1,
+0.164117813E-1,
+0.216339231,
-0.377901912,
+0.735760689,
-0.210621357,
+0.264455557,
+0.403519869
970var(:,
0)
= getVar(sample, dim)
972+0.455610245,
+0.257973731,
+0.124659672
974 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
975 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
977call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
979+0.455610305,
+0.257973731,
+0.124659665
980call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
982+0.455610305,
+0.257973731,
+0.124659665
984+0.455610245,
+0.257973731,
+0.124659672
989 lb(isam)
= ub(isam
- 1)
+ 1
997call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
998call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1002-0.499877453,
-0.409510136E-1,
+0.640338659,
+0.903482795,
+0.445844412,
-0.943887234,
+0.566644311,
-0.288667440
1003-0.338787079,
-0.669931889,
-0.287321806E-1,
+0.545727015E-1,
+0.615508199,
-0.469040871,
+0.284605980,
-0.746615648
1004var(:,
0)
= getVar(sample, dim)
1006+0.361920327,
+0.199242145
1008 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1009 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1011call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1013+0.361920357,
+0.199242160
1014call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1016+0.361920357,
+0.199242160
1018+0.361920327,
+0.199242145
1023 lb(isam)
= ub(isam
- 1)
+ 1
1031call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1032call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1036+0.788258314E-1,
+0.606618047,
-0.561826229E-1,
+0.178701878,
+0.845884919,
+0.984233141
1037+0.525810957,
-0.680465698E-1,
-0.624785066,
+0.494570374,
+0.886281610,
-0.326054931
1038var(:,
0)
= getVar(sample, dim)
1040+0.155602321,
+0.279418826
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))
1047+0.155602306,
+0.279418826
1048call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1050+0.155602306,
+0.279418826
1052+0.155602321,
+0.279418826
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.323566437,
-0.561246753,
-0.921070218,
-0.576605082,
-0.196160913,
+0.966912389,
+0.590512037,
-0.265647769,
+0.426393032,
+0.527081728,
-0.592560768
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.615378737,
+0.337416172,
-0.408806086,
-0.981052637,
+0.788019776,
-0.842159867
1104+0.426098228,
+0.692504525,
+0.589112997,
+0.869750619,
-0.195138454E-1,
+0.726882935
1105var(:,
0)
= getVar(sample, dim)
1107+0.485353678,
+0.825041234E-1
1109 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1110 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1112call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1114+0.485353708,
+0.825041234E-1
1115call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1117+0.485353708,
+0.825041234E-1
1119+0.485353678,
+0.825041234E-1
1124 lb(isam)
= ub(isam
- 1)
+ 1
1132call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1133call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1137(
-0.898863435,
-0.490910411), (
+0.531323195,
-0.810796022), (
-0.507146120,
-0.634126306), (
+0.477120876,
-0.859382629), (
-0.312059522,
+0.434140325), (
+0.768672824,
-0.405973792), (
+0.709519625,
+0.343047380E-1), (
-0.919806957E-2,
-0.501530647), (
+0.206574202,
+0.383299470), (
-0.313534975,
+0.618358254)
1138(
+0.922214270,
-0.722887874), (
-0.720769763,
-0.591884851E-1), (
-0.565511823,
-0.443517804), (
+0.113196850,
-0.296497107), (
+0.287744403,
+0.379383564E-1), (
-0.206608176,
+0.644248009), (
-0.811513186,
-0.735806465), (
-0.991337299E-1,
-0.285777688), (
+0.583809853,
-0.547064424), (
-0.463418126,
+0.835097194)
1139(
+0.439243317,
+0.507998109), (
+0.903160334,
-0.502853155), (
-0.226024508,
+0.944223046), (
+0.585753918E-1,
-0.388604164), (
-0.710730672,
-0.215749502), (
-0.686340690,
+0.721068263), (
+0.410027504,
+0.900212049), (
+0.488995552,
+0.959570408), (
-0.468291402,
-0.312693119), (
-0.956358671,
+0.838008165)
1140var(:,
0)
= getVar(sample, dim)
1142+0.554112434,
+0.555925667,
+0.698470771
1144 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1145 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1147call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1149+0.554112434,
+0.555925667,
+0.698470771
1150call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1152+0.554112434,
+0.555925667,
+0.698470771
1154+0.554112434,
+0.555925667,
+0.698470771
1159 lb(isam)
= ub(isam
- 1)
+ 1
1167call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1168call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1172(
-0.476198077,
+0.982520342), (
+0.131975055,
-0.308599353), (
+0.397111177E-1,
+0.812063694), (
+0.925534487,
+0.819882274), (
+0.493365526,
-0.481618643), (
+0.769308805E-1,
-0.115126967), (
+0.415607214,
+0.775631070), (
-0.527027845,
-0.893667459), (
-0.462110877,
+0.350874305), (
+0.616801262,
-0.757657766)
1173(
-0.252541542,
+0.291156769E-1), (
-0.427683711,
+0.288555980), (
+0.165092945E-1,
+0.160862923), (
-0.924117327,
-0.349455237), (
-0.534555793,
+0.938708663), (
+0.109989643E-1,
-0.989046574), (
+0.909367561,
-0.851778150), (
-0.193476677E-1,
+0.834021568E-1), (
+0.582243204E-1,
-0.649415135), (
+0.469630957E-1,
+0.515188813)
1174(
-0.909083366,
-0.413929105), (
-0.527043343E-1,
+0.539823771E-1), (
-0.314494371E-1,
+0.893720627), (
-0.397528410,
-0.654806137), (
-0.106633425,
+0.418690085), (
-0.894980311,
+0.225264668), (
-0.665152550,
+0.247690678E-1), (
-0.963680029,
-0.664603710), (
+0.588335276,
-0.933631063), (
+0.704765677,
-0.181089640)
1175var(:,
0)
= getVar(sample, dim)
1177+0.683834910,
+0.553842247,
+0.609315455
1179 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1180 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1182call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1184+0.683834910,
+0.553842187,
+0.609315455
1185call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1187+0.683834910,
+0.553842187,
+0.609315455
1189+0.683834910,
+0.553842247,
+0.609315455
1194 lb(isam)
= ub(isam
- 1)
+ 1
1202call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1203call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1207(
+0.278764367,
+0.698744893), (
-0.236132145E-1,
+0.257015705), (
+0.269762754,
-0.875359535), (
+0.831024885,
-0.182772875E-1), (
-0.113333106,
+0.788612127), (
-0.532136440,
+0.887642264), (
-0.344892502,
+0.937377930), (
+0.403262019,
-0.349161148), (
+0.303972006,
+0.373735785), (
-0.651593685,
+0.477085114)
1208(
-0.108979940,
-0.651787877), (
-0.281770349,
-0.733716726), (
-0.986692309,
+0.316802025), (
+0.502418518,
+0.227830887), (
-0.954000354,
-0.689566135E-2), (
+0.166068554,
-0.827320814E-1), (
-0.418071628,
+0.650168657), (
+0.320088029,
-0.777449012), (
-0.381816745,
+0.956978798E-1), (
+0.169273496,
+0.634987593)
1209(
-0.736998320E-1,
+0.410800457), (
-0.517692566E-1,
+0.134818673), (
-0.808278918,
+0.511107564), (
+0.468002677,
-0.112135172), (
-0.897903204,
+0.375828743E-1), (
+0.509718418,
-0.208347321), (
+0.304595232E-1,
+0.508501172), (
-0.550749183,
+0.151139736), (
-0.456794858,
+0.991290450), (
+0.806081533,
-0.955107808)
1210(
-0.866036534,
-0.169871092), (
-0.774098396,
-0.303337693), (
-0.913855076,
+0.206735134E-1), (
+0.971655846E-1,
+0.564742804), (
+0.103779435,
-0.784039617), (
+0.998983622,
-0.224384785), (
+0.366091967,
+0.332816958), (
+0.424306512,
-0.831377506E-2), (
+0.105937004,
+0.627314568), (
+0.500897527,
+0.527920842)
1211var(:,
0)
= getVar(sample, dim)
1213+0.500744760,
+0.486831188,
+0.547000527,
+0.563286126
1215 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1216 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1218call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1220+0.500744760,
+0.486831188,
+0.547000527,
+0.563286126
1221call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1223+0.500744760,
+0.486831188,
+0.547000527,
+0.563286126
1225+0.500744760,
+0.486831188,
+0.547000527,
+0.563286126
1230 lb(isam)
= ub(isam
- 1)
+ 1
1238call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1239call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1243(
-0.642212510,
-0.355322480), (
-0.562590718,
-0.614030123), (
+0.898363829,
+0.523419380E-1), (
+0.429342985,
+0.896628499), (
+0.841554403E-1,
-0.827000380), (
-0.433498502,
-0.868036389), (
-0.424673915,
-0.323034286), (
+0.987298727,
+0.363167763), (
+0.367053151,
+0.626986146), (
+0.654906869,
-0.490441561), (
+0.624285221,
-0.337656498), (
-0.491318703E-1,
-0.451584458)
1244(
+0.653537154,
+0.478377700), (
-0.132556796,
-0.908329487E-1), (
-0.153341413,
+0.246960282), (
+0.106026292,
-0.826819420), (
+0.293275237,
+0.441260338E-1), (
+0.453273535,
-0.722075582), (
-0.278672218,
+0.658713460), (
-0.583637357,
-0.470158458), (
+0.570943475,
+0.362503052), (
-0.543933868,
-0.918868065), (
-0.662287712,
-0.425020099), (
+0.249499321,
-0.949183702)
1245(
-0.376461387,
-0.739163518), (
-0.981158853,
+0.231155753), (
-0.552382946,
+0.273451805), (
-0.864680171,
-0.447574019), (
+0.505847573,
-0.768286109), (
-0.554777384E-1,
-0.537810564), (
-0.843819618,
+0.813514709), (
+0.834091425,
-0.687093019), (
-0.765749216,
-0.322658539), (
-0.931522965,
-0.444607139), (
-0.305900812,
-0.637520790), (
+0.643431067,
-0.796520829)
1246var(:,
0)
= getVar(sample, dim)
1248+0.599861920,
+0.496708930,
+0.625316739
1250 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1251 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1253call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1255+0.599861920,
+0.496708930,
+0.625316679
1256call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1258+0.599861920,
+0.496708930,
+0.625316679
1260+0.599861920,
+0.496708930,
+0.625316739
1265 lb(isam)
= ub(isam
- 1)
+ 1
1273call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1274call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1278(
-0.643377662,
+0.164959192), (
-0.942224979,
-0.485823870), (
-0.851301670,
-0.937904477), (
+0.290958762,
+0.888182044), (
-0.188605905,
+0.656538486)
1279var(:,
0)
= getVar(sample, dim)
1283 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1284 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1286call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1289call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1303 lb(isam)
= ub(isam
- 1)
+ 1
1311call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1312call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1316+0.705856204,
+0.809519410,
+0.434411526,
+0.927069902
1317+0.337166667,
+0.843157053,
+0.500416756E-2,
-0.819877625
1321var(:,
0)
= getVar(sample,
2_IK, iweight)
1323+0.343582705E-1,
+0.285232842
1325 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1326 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1328call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1330+0.343582742E-1,
+0.285232842
1331call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1333+0.343582742E-1,
+0.285232842
1335+0.343582705E-1,
+0.285232842
1340 lb(isam)
= ub(isam
- 1)
+ 1
1348call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1349call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1353+0.415035725,
+0.702162504,
+0.234350920,
+0.251109600,
+0.978246212,
-0.125225663,
+0.788504958,
-0.314499736
1356+3,
+3,
+6,
+2,
+3,
+10,
+6,
+2
1357var(:,
0)
= getVar(sample,
2_IK, iweight)
1361 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1362 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1364call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1367call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1376 lb(isam)
= ub(isam
- 1)
+ 1
1384call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1385call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1389-0.354228020,
+0.884099007E-1,
+0.845666528,
+0.429026008,
-0.278605938,
-0.906816125,
+0.899421811,
-0.508985162,
-0.189560056,
-0.343139172,
+0.784993291
1392+2,
+1,
+6,
+6,
+3,
+1,
+2,
+10,
+8,
+8,
+9
1393var(:,
0)
= getVar(sample,
2_IK, iweight)
1397 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1398 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1400call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1403call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1412 lb(isam)
= ub(isam
- 1)
+ 1
1420call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1421call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1425-0.658561945,
+0.179705620,
+0.278295755,
+0.317425489,
-0.154611588
1429var(:,
0)
= getVar(sample,
2_IK, iweight)
1433 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1434 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1436call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1439call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1448 lb(isam)
= ub(isam
- 1)
+ 1
1456call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1457call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1461-0.741678834,
-0.274475217,
-0.197572112,
+0.302025795,
+0.368824124,
-0.834445834,
+0.544283390E-1,
-0.577114582
1464+4,
+3,
+3,
+9,
+6,
+4,
+6,
+5
1465var(:,
0)
= getVar(sample,
2_IK, iweight)
1469 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1470 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1472call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1475call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1484 lb(isam)
= ub(isam
- 1)
+ 1
1492call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1493call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1497(
+0.666962147,
-0.961538076), (
+0.706293583,
-0.441780448), (
-0.796339393,
-0.575214148), (
+0.967309833,
-0.177862644E-1), (
+0.747689009E-1,
+0.160133600), (
-0.139149666,
+0.177966237), (
+0.447077513,
+0.793694496), (
-0.821147561,
-0.908249497), (
+0.761312246,
+0.711736441), (
+0.571289897,
+0.556646943)
1498(
-0.132532358,
+0.421752214), (
-0.363661766,
+0.360602021), (
+0.328309298,
+0.415349007E-2), (
+0.387532711,
-0.820196986), (
-0.222299576,
+0.394693851), (
+0.976669669,
+0.872541547), (
-0.651199937,
-0.795539260), (
+0.325256109,
-0.494964123E-1), (
+0.973410606E-1,
-0.652403951), (
-0.949118137E-1,
+0.748809934)
1499(
-0.650966525,
+0.197395444), (
+0.973906994,
+0.477362156), (
-0.238532782,
-0.295203686), (
+0.162520528,
+0.657770038), (
-0.138929009,
+0.727949739), (
-0.511267900,
+0.713053107), (
+0.203129292,
+0.773863077), (
-0.665916562,
-0.989996076), (
-0.934052229,
-0.475013375), (
+0.830846667,
-0.926733494)
1500(
-0.210882306,
-0.909214854), (
+0.299995899,
-0.405943632), (
+0.230109453,
+0.812988400), (
-0.430759788,
+0.390757561), (
+0.825081468,
+0.565663576), (
+0.795649409,
+0.811537147), (
+0.118882895,
-0.599910617), (
+0.635424256,
+0.913573384), (
-0.177166104,
+0.941385627), (
+0.600104213,
-0.216313601)
1503+10,
+4,
+2,
+10,
+7,
+6,
+7,
+6,
+5,
+4
1504var(:,
0)
= getVar(sample,
2_IK, iweight)
1506+0.706464827,
+0.581372440,
+0.709974051,
+0.687819958
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.706464887,
+0.581372499,
+0.709974051,
+0.687820077
1514call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1516+0.706464887,
+0.581372499,
+0.709974051,
+0.687820077
1518+0.706464827,
+0.581372440,
+0.709974051,
+0.687819958
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.245508194,
+0.855350137), (
-0.228887796,
-0.817185760), (
+0.728183150,
+0.239063263), (
+0.683184385,
+0.375074267), (
+0.581895113,
-0.130779743), (
-0.817774296,
+0.730833769), (
-0.143849850,
-0.148228765), (
-0.418767333,
+0.826942086), (
+0.928836346,
-0.591916919), (
-0.578984261,
-0.625543594)
1537(
-0.135410190,
-0.624632359), (
-0.523886323,
+0.923788071), (
-0.231252909E-1,
-0.772210360E-1), (
-0.117218971,
-0.604552150), (
+0.291371107,
-0.253922939E-1), (
+0.406780720,
-0.623906374), (
-0.853809834,
-0.425007820), (
+0.692703605,
+0.481491446), (
-0.680694580,
-0.575056434), (
+0.819445848E-1,
-0.926185012)
1540+4,
+6,
+8,
+7,
+10,
+6,
+5,
+3,
+7,
+10
1541var(:,
0)
= getVar(sample,
2_IK, iweight)
1543+0.677568495,
+0.451545179
1545 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1546 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1548call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1550+0.677568495,
+0.451545179
1551call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1553+0.677568495,
+0.451545179
1555+0.677568495,
+0.451545179
1560 lb(isam)
= ub(isam
- 1)
+ 1
1568call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1569call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1573(
+0.106148362,
+0.910557151), (
+0.378223419,
+0.434238672), (
-0.618231773,
-0.588397026), (
+0.982396483,
-0.301270604), (
-0.487957478,
+0.552514195), (
+0.935700774,
-0.354831100), (
-0.773203254,
+0.264267921E-1), (
-0.395789146E-1,
+0.236356258E-2), (
-0.556547523,
+0.582742333)
1576+5,
+1,
+8,
+9,
+1,
+10,
+1,
+2,
+5
1577var(:,
0)
= getVar(sample,
2_IK, iweight)
1581 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1582 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1584call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1587call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1596 lb(isam)
= ub(isam
- 1)
+ 1
1604call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1605call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1609(
-0.542162061,
+0.190398693E-1), (
-0.285102844,
-0.459957123E-1), (
-0.877113342E-1,
+0.783468246), (
+0.930641532,
-0.167688847), (
-0.783588886E-1,
-0.917996287), (
+0.819467068,
-0.301821232), (
-0.129559755,
-0.802323461), (
-0.758509278,
-0.810205936)
1612+7,
+4,
+8,
+5,
+5,
+4,
+9,
+3
1613var(:,
0)
= getVar(sample,
2_IK, iweight)
1617 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1618 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1620call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1623call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1632 lb(isam)
= ub(isam
- 1)
+ 1
1640call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1641call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1645(
-0.926922679,
-0.298818350E-1), (
-0.243337274,
+0.745286942E-1), (
-0.280132413,
-0.182285786), (
+0.755516768,
-0.802497745), (
+0.118303895,
+0.560160875), (
-0.950950503,
-0.429252744), (
-0.554244518E-1,
+0.524238229), (
+0.770309091,
+0.249984622), (
+0.607469916,
+0.616116524), (
-0.398754001,
+0.540415406), (
-0.574567437,
-0.278787255), (
+0.107298613,
-0.965413809), (
+0.110342622,
+0.416818857E-1)
1646(
-0.159470677,
+0.931211710E-1), (
-0.698794603,
+0.891720057E-1), (
+0.460260749,
+0.388891220), (
+0.516233444E-1,
+0.875191450), (
-0.528159380,
+0.766205907), (
-0.646242499,
+0.493923426E-1), (
-0.737263799,
-0.848336577), (
-0.145791888,
+0.556642175), (
+0.805462122,
-0.323294759), (
-0.973893166,
-0.555580139), (
+0.375540376,
-0.310533404), (
-0.149048328,
+0.442517161), (
+0.266939521,
+0.240548134)
1649+1,
+3,
+9,
+4,
+5,
+10,
+10,
+6,
+7,
+4,
+1,
+9,
+2
1650var(:,
0)
= getVar(sample,
2_IK, iweight)
1652+0.591235518,
+0.570302665
1654 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1655 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1657call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1659+0.591235518,
+0.570302665
1660call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1662+0.591235518,
+0.570302665
1664+0.591235518,
+0.570302665
1674 lb(isam)
= ub(isam
- 1)
+ 1
1682call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1683call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1687+0.460029602,
-0.293368578,
-0.373516560,
-0.962636828,
+0.375983477
1688+0.930958867,
+0.256311297,
+0.736812711,
+0.138683319,
+0.487774849
1691+1.01017797,
+1.62548566,
+1.40532470,
+1.49459112,
+1.90453887
1692var(:,
0)
= getVar(sample,
2_IK, rweight)
1694+0.267559111,
+0.743860304E-1
1696 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1697 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1699call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1701+0.267559052,
+0.743860304E-1
1702call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1704+0.267559052,
+0.743860304E-1
1706+0.267559111,
+0.743860304E-1
1711 lb(isam)
= ub(isam
- 1)
+ 1
1719call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1720call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1724+0.614552379,
-0.423786402,
-0.229943037,
+0.696812987,
-0.525009632E-2,
-0.733222723,
-0.185187101,
+0.466558933,
-0.722688913,
-0.403021455,
+0.717701316
1725+0.433647037,
-0.598072052,
-0.236315727,
-0.953127384,
-0.476506710,
-0.639668465,
-0.527256846,
-0.715587497,
-0.324638128,
+0.554684401E-1,
-0.615184903
1726+0.110324264,
+0.934567928,
-0.818084836,
+0.529495358,
-0.194543481,
-0.141818523,
+0.246729612,
+0.331275821,
-0.857547402,
+0.961508036,
-0.407112360
1729+1.77659607,
+1.96177340,
+1.05248678,
+1.07430673,
+1.14621592,
+1.55237627,
+1.47137332,
+1.55097818,
+1.49967313,
+1.00736797,
+1.24252176
1730var(:,
0)
= getVar(sample,
2_IK, rweight)
1732+0.286406636,
+0.144501895,
+0.338863969
1734 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1735 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1737call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1739+0.286406636,
+0.144501910,
+0.338864028
1740call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1742+0.286406636,
+0.144501910,
+0.338864028
1744+0.286406636,
+0.144501895,
+0.338863969
1749 lb(isam)
= ub(isam
- 1)
+ 1
1757call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1758call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1762+0.961881995,
+0.891190052,
-0.572649479,
-0.596723557E-1,
+0.998573303,
-0.878392100
1763-0.326285481,
-0.361691713,
+0.106955647,
-0.224904656,
-0.234822989,
-0.245733142
1764-0.541052938,
-0.131974101,
-0.803843379,
-0.576705217,
+0.884593606,
-0.139188766E-1
1767+1.82439256,
+1.16331935,
+1.01992369,
+1.45991325,
+1.39180076,
+1.22219563
1768var(:,
0)
= getVar(sample,
2_IK, rweight)
1770+0.562728822,
+0.189643465E-1,
+0.306391031
1772 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1773 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1775call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1777+0.562728822,
+0.189643428E-1,
+0.306390971
1778call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1780+0.562728822,
+0.189643428E-1,
+0.306390971
1782+0.562728822,
+0.189643465E-1,
+0.306391031
1787 lb(isam)
= ub(isam
- 1)
+ 1
1795call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1796call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1800-0.655491114,
-0.816720009,
-0.714972854,
+0.192758560,
+0.663073182
1801+0.732982516,
-0.912226558,
-0.194191337,
+0.576962471,
-0.823698878
1804+1.15768600,
+1.77890849,
+1.56300688,
+1.28153014,
+1.49125910
1805var(:,
0)
= getVar(sample,
2_IK, rweight)
1807+0.355272084,
+0.448565066
1809 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1810 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1812call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1814+0.355272084,
+0.448565036
1815call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1817+0.355272084,
+0.448565036
1819+0.355272084,
+0.448565066
1824 lb(isam)
= ub(isam
- 1)
+ 1
1832call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1833call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1837+0.808606744,
+0.220362544,
+0.801791430,
-0.757616401,
-0.679708362,
+0.911671281,
-0.647293568,
-0.262669921,
+0.786936045,
+0.309132457,
-0.845189333
1838+0.172615051E-2,
+0.589989424E-1,
+0.167469859,
+0.561426044,
-0.298482895,
-0.735116720,
-0.654819131,
+0.248589396,
+0.869643688E-1,
+0.595354557,
+0.348871350
1839+0.250681400,
+0.728107452,
+0.918661118,
+0.674109340,
+0.154230356,
+0.596910715,
-0.519001603,
+0.778382659,
+0.352741480E-1,
-0.463998318E-2,
-0.673311830
1842+1.28772223,
+1.69019306,
+1.27230287,
+1.46205974,
+1.63165152,
+1.42136574,
+1.16030455,
+1.74902666,
+1.87387681,
+1.00522423,
+1.69825566
1843var(:,
0)
= getVar(sample,
2_IK, rweight)
1845+0.465581715,
+0.158381343,
+0.256403744
1847 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1848 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1850call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1852+0.465581685,
+0.158381343,
+0.256403714
1853call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1855+0.465581685,
+0.158381343,
+0.256403714
1857+0.465581715,
+0.158381343,
+0.256403744
1862 lb(isam)
= ub(isam
- 1)
+ 1
1870call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1871call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1875(
-0.261031628,
-0.914617062), (
+0.409027338E-1,
-0.744035244), (
-0.721457601,
+0.585199475), (
-0.716470003,
-0.490597486E-1), (
-0.798190951,
+0.701974392), (
-0.992403269,
+0.839413047), (
+0.899656534,
-0.973448753E-1), (
+0.317848682,
-0.997044325), (
+0.637656808,
-0.521872759)
1878+1.67616272,
+1.33454847,
+1.48370922,
+1.68066669,
+1.84223580,
+1.74137402,
+1.14436841,
+1.76129508,
+1.78658557
1879var(:,
0)
= getVar(sample,
2_IK, rweight)
1883 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1884 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1886call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1889call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1898 lb(isam)
= ub(isam
- 1)
+ 1
1906call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1907call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1911(
-0.593937635E-1,
-0.617805839), (
+0.252352953E-1,
-0.486315131), (
-0.919135332,
+0.405859947E-1), (
-0.629978776,
-0.240323663), (
+0.448878407,
-0.718716145), (
-0.495951176E-1,
-0.970683813), (
+0.640217185,
-0.747477770)
1912(
+0.416382909,
-0.603934407), (
-0.380209088,
-0.375291348), (
+0.216463089,
-0.332289577), (
-0.114862084,
+0.683399439), (
-0.761282444E-1,
-0.953667283), (
+0.362300634,
+0.996231675), (
+0.427666426,
+0.954665184)
1915+1.60993636,
+1.99905765,
+1.19328082,
+1.36388254,
+1.98567915,
+1.46525574,
+1.32392824
1916var(:,
0)
= getVar(sample,
2_IK, rweight)
1918+0.312767416,
+0.648638785
1920 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1921 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1923call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1925+0.312767446,
+0.648638785
1926call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1928+0.312767446,
+0.648638785
1930+0.312767416,
+0.648638785
1935 lb(isam)
= ub(isam
- 1)
+ 1
1943call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1944call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1948(
+0.695660114E-1,
-0.623174906E-1), (
+0.972343206,
-0.553104758), (
+0.711107731,
+0.517208099), (
-0.896051288,
+0.727480292), (
-0.714633346,
+0.843813419), (
+0.151866317,
-0.440363050), (
+0.112660170,
-0.499837041), (
+0.680850625,
+0.773581266)
1951+1.05890477,
+1.78127253,
+1.63326573,
+1.54850364,
+1.37721264,
+1.80129933,
+1.64056325,
+1.73799264
1952var(:,
0)
= getVar(sample,
2_IK, rweight)
1956 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1957 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1959call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1962call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1971 lb(isam)
= ub(isam
- 1)
+ 1
1979call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1980call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1984(
-0.169258118E-1,
+0.951236725), (
-0.974802494,
+0.643455744), (
+0.784508109,
-0.166407466), (
+0.137954831,
-0.281824827), (
+0.553798914,
-0.654775381), (
-0.208607912,
+0.722259283), (
-0.491883397,
-0.414582253), (
+0.307537675,
+0.605530024)
1987+1.54709601,
+1.67517757,
+1.62471104,
+1.34143257,
+1.62061715,
+1.45655537,
+1.69151711,
+1.62985349
1988var(:,
0)
= getVar(sample,
2_IK, rweight)
1992 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1993 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1995call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1998call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2007 lb(isam)
= ub(isam
- 1)
+ 1
2015call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
2016call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
2020(
-0.378898263,
-0.270808458), (
-0.627832174,
-0.321482301), (
+0.441128135,
-0.945881248), (
-0.309360027,
-0.674990296), (
+0.108477712,
-0.660212517), (
-0.485677481,
-0.420334458), (
+0.584257126,
+0.291362047), (
-0.302208543,
-0.687203050), (
-0.508160233,
+0.882365823)
2021(
+0.957346797,
-0.228343368), (
-0.483841538,
-0.329884887), (
+0.407219410,
+0.964086652), (
-0.603586316,
+0.722148418E-1), (
+0.527539730,
+0.616263151), (
-0.197870851,
-0.160362840), (
+0.224440098E-1,
-0.829549313), (
-0.560439825E-1,
-0.664375663), (
-0.445518970,
-0.263706446)
2022(
-0.273316979,
-0.936767817), (
-0.580408454,
-0.836801529E-2), (
+0.346042395,
+0.349774480), (
-0.201066017,
+0.473378897), (
+0.476075888,
+0.987301707), (
-0.473621607,
+0.879487634), (
+0.982446909,
-0.424958110), (
+0.877273202,
+0.946827650), (
+0.599318027,
+0.547075629)
2025+1.34842634,
+1.66167831,
+1.85626805,
+1.43410659,
+1.33140016,
+1.84539771,
+1.59292006,
+1.14236569,
+1.31127191
2026var(:,
0)
= getVar(sample,
2_IK, rweight)
2028+0.459440827,
+0.542389929,
+0.674313724
2030 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2031 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
2033call setVarMerged(varMerged, var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2035+0.459440857,
+0.542389929,
+0.674313784
2036call setVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
2038+0.459440857,
+0.542389929,
+0.674313784
2040+0.459440827,
+0.542389929,
+0.674313724