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
= 10
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), var(
0:nsam), varMerged
31 real(TKG),
allocatable :: sample(:)
34 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
37 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
38 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
41 lb(isam)
= ub(isam
- 1)
+ 1
48 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
52 call disp%show(
"var(0) = getVar(sample)")
57 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)))")
58 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)))")
61 var(isam)
= getVar(sample(lb(isam):ub(isam)))
62 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
64 call disp%show(
"varMerged = getVarMerged(var(2), var(1), mean(1) - mean(2), ub(1) / real(ub(2), TKG))")
65 varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
75 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
76 call disp%show(
"!Compute the biased merged variance of a frequency weighted univariate sample.")
77 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
81 real(TKG) :: mean(
1:nsam), var(
0:nsam), varMerged
82 real(TKG),
allocatable :: sample(:)
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 = getUnifRand(0., 1., ub(nsam))")
103 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, 1, IK))")
107 call disp%show(
"var(0) = getVar(sample, iweight)")
108 var(
0)
= getVar(sample, iweight)
109 call disp%show(
"var(0) ! reference")
112 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
113 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
116 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
117 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
119 call disp%show(
"varMerged = getVarMerged(var(2), var(1), mean(1) - mean(2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
120 varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
123 call disp%show(
"var(0) ! reference")
130 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
131 call disp%show(
"!Compute the biased merged variance of a reliability weighted univariate sample.")
132 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
136 real(TKG) :: mean(
1:nsam), var(
0:nsam), varMerged
137 real(TKG),
allocatable :: sample(:), rweight(:)
140 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(2, 7)")
143 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
144 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
147 lb(isam)
= ub(isam
- 1)
+ 1
154 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
158 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, 1, IK))")
162 call disp%show(
"var(0) = getVar(sample, rweight)")
163 var(
0)
= getVar(sample, rweight)
164 call disp%show(
"var(0) ! reference")
167 call disp%show(
" var(isam) = getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
168 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
171 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
172 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
174 call disp%show(
"varMerged = getVarMerged(var(2), var(1), mean(1) - mean(2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
175 varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
178 call disp%show(
"var(0) ! reference")
185 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
186 call disp%show(
"!Compute the biased merged variance of a multivariate sample.")
187 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
191 real(TKG),
allocatable :: mean(:,:), var(:,:), varMerged(:)
192 real(TKG),
allocatable :: sample(:,:)
195 call disp%show(
"dim = 2; 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(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
211 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
212 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
213 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
214 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
215 call disp%show(
"call setResized(varMerged, ndim)")
217 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
221 call disp%show(
"var(:,0) = getVar(sample, dim)")
222 var(:,
0)
= getVar(sample, dim)
223 call disp%show(
"var(:,0) ! reference")
226 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), dim)")
227 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
230 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
231 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
233 call disp%show(
"varMerged = getVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
234 varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
237 call disp%show(
"var(:,0) ! reference")
244 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
245 call disp%show(
"!Compute the biased merged variance of a frequency weighted multivariate sample.")
246 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
250 real(TKG),
allocatable :: mean(:,:), var(:,:), varMerged(:)
251 real(TKG),
allocatable :: sample(:,:)
254 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
257 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
258 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
261 lb(isam)
= ub(isam
- 1)
+ 1
268 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
270 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
271 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
272 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
273 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
274 call disp%show(
"call setResized(varMerged, ndim)")
276 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
280 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
284 call disp%show(
"var(:,0) = getVar(sample, 2_IK, iweight)")
285 var(:,
0)
= getVar(sample,
2_IK, iweight)
286 call disp%show(
"var(:,0) ! reference")
289 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
290 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
293 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
294 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
296 call disp%show(
"varMerged = getVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
297 varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
300 call disp%show(
"var(:,0) ! reference")
307 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
308 call disp%show(
"!Compute the biased merged variance of a reliability weighted multivariate sample.")
309 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
313 real(TKG),
allocatable :: sample(:,:), mean(:,:)
314 real(TKG),
allocatable :: rweight(:), var(:,:), varMerged(:)
317 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
320 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
321 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
324 lb(isam)
= ub(isam
- 1)
+ 1
331 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
333 call disp%show(
"call setRebound(var, [1_IK, 0_IK], [ndim, nsam])")
334 call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
335 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
336 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
337 call disp%show(
"call setResized(varMerged, ndim)")
339 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
343 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
347 call disp%show(
"var(:,0) = getVar(sample, 2_IK, rweight)")
348 var(:,
0)
= getVar(sample,
2_IK, rweight)
349 call disp%show(
"var(:,0) ! reference")
352 call disp%show(
" var(:,isam) = getVar(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
353 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
356 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
357 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
359 call disp%show(
"varMerged = getVarMerged(var(:,2), var(:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
360 varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
363 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 RK
The default real kind in the ParaMonte library: real64 in Fortran, c_double in C-Fortran Interoperati...
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.365969598,
+0.539367974,
+0.860647857,
+0.669404268,
+0.639917970,
+0.934130669,
+0.761378646,
+0.769686580,
+0.910623789,
+0.107624352,
+0.204128623E-1,
+0.404817104
23 var(isam)
= getVar(sample(lb(isam):ub(isam)))
24 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
26varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
35 lb(isam)
= ub(isam
- 1)
+ 1
44+0.891571641E-1,
+0.317646980,
+0.591190279,
+0.989820838,
+0.288057327E-2,
+0.923178732,
+0.856376529,
+0.651625335,
+0.543845952,
+0.681804121,
+0.401710033
49 var(isam)
= getVar(sample(lb(isam):ub(isam)))
50 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
52varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
61 lb(isam)
= ub(isam
- 1)
+ 1
70+0.549304366,
+0.507123470E-1,
+0.278680503,
+0.427578390,
+0.692629993,
+0.272925079,
+0.277589381,
+0.598483503,
+0.275083840,
+0.243517518
75 var(isam)
= getVar(sample(lb(isam):ub(isam)))
76 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
78varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
87 lb(isam)
= ub(isam
- 1)
+ 1
96+0.635099828,
+0.994973958,
+0.706785679,
+0.517038107,
+0.955167711,
+0.644648612,
+0.403783441,
+0.643406510,
+0.512327194,
+0.206679702E-1,
+0.994967818E-1,
+0.920399427E-1,
+0.567415178
101 var(isam)
= getVar(sample(lb(isam):ub(isam)))
102 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
104varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
113 lb(isam)
= ub(isam
- 1)
+ 1
122+0.127256453,
+0.323118567E-1,
+0.867960989,
+0.121276140,
+0.780102551,
+0.505891681,
+0.692693591,
+0.235909045,
+0.668296993,
+0.820031404
127 var(isam)
= getVar(sample(lb(isam):ub(isam)))
128 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
130varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
139 lb(isam)
= ub(isam
- 1)
+ 1
148+0.138990223,
+0.164124370E-1,
+0.659975231,
+0.816543758,
+0.326628506,
+0.369476438,
+0.718675852E-1
153 var(isam)
= getVar(sample(lb(isam):ub(isam)))
154 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
156varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
165 lb(isam)
= ub(isam
- 1)
+ 1
174+0.747892261,
+0.292134821,
+0.224555314,
+0.347658992,
+0.881615758
179 var(isam)
= getVar(sample(lb(isam):ub(isam)))
180 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
182varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
191 lb(isam)
= ub(isam
- 1)
+ 1
200+0.388197541,
+0.239570796,
+0.423648715,
+0.863466620,
+0.826491952,
+0.997288525,
+0.358664215,
+0.289916992E-3
205 var(isam)
= getVar(sample(lb(isam):ub(isam)))
206 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
208varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
217 lb(isam)
= ub(isam
- 1)
+ 1
226+0.578138411,
+0.826375186,
+0.332468510,
+0.586797655,
+0.663269103,
+0.832929969,
+0.983767748,
+0.220828533,
+0.234472811
231 var(isam)
= getVar(sample(lb(isam):ub(isam)))
232 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
234varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
243 lb(isam)
= ub(isam
- 1)
+ 1
252+0.671956599,
+0.629860103,
+0.177184701,
+0.887980044,
+0.434041619,
+0.608365059,
+0.628340721,
+0.341608107,
+0.694782138E-1,
+0.421007276E-1
257 var(isam)
= getVar(sample(lb(isam):ub(isam)))
258 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
260varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2), ub(
1)
/ real(ub(
2), TKG))
274 lb(isam)
= ub(isam
- 1)
+ 1
283+0.192155838,
+0.215730846,
+0.244340599,
+0.760361016,
+0.366439998,
+0.792468250
286+1,
+7,
+4,
+1,
+2,
+9
287var(
0)
= getVar(sample, iweight)
291 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
292 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
294varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
303 lb(isam)
= ub(isam
- 1)
+ 1
312+0.185890496,
+0.705393016,
+0.708951175,
+0.382265210,
+0.243085206,
+0.628276885,
+0.576271594,
+0.778601766E-1
315+1,
+2,
+2,
+2,
+2,
+8,
+1,
+5
316var(
0)
= getVar(sample, iweight)
320 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
321 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
323varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
332 lb(isam)
= ub(isam
- 1)
+ 1
341+0.415813923E-2,
+0.174619198,
+0.167224348,
+0.997489572
345var(
0)
= getVar(sample, iweight)
349 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
350 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
352varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
361 lb(isam)
= ub(isam
- 1)
+ 1
370+0.768472910,
+0.586231112,
+0.418836236,
+0.253783464,
+0.916766822,
+0.117777646,
+0.934052646,
+0.425538063,
+0.554310918
373+2,
+2,
+3,
+2,
+8,
+1,
+3,
+1,
+6
374var(
0)
= getVar(sample, iweight)
378 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
379 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
381varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
390 lb(isam)
= ub(isam
- 1)
+ 1
399+0.403336704,
+0.437636077,
+0.112331212,
+0.430452704,
+0.875146091,
+0.120087147,
+0.333021164,
+0.831938624,
+0.318392038
402+10,
+8,
+2,
+4,
+7,
+6,
+2,
+9,
+1
403var(
0)
= getVar(sample, iweight)
407 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
408 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
410varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
419 lb(isam)
= ub(isam
- 1)
+ 1
428+0.499932528,
+0.195486009,
+0.106252372,
+0.968675613E-1,
+0.156602442
432var(
0)
= getVar(sample, iweight)
436 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
437 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
439varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
448 lb(isam)
= ub(isam
- 1)
+ 1
457+0.473103821,
+0.804166555,
+0.371435285,
+0.185356438,
+0.802984834,
+0.267572820,
+0.237235248,
+0.869830072,
+0.413538218,
+0.836687207
460+3,
+4,
+1,
+4,
+8,
+9,
+1,
+10,
+10,
+5
461var(
0)
= getVar(sample, iweight)
465 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
466 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
468varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
477 lb(isam)
= ub(isam
- 1)
+ 1
486+0.517180562E-1,
+0.344385505E-1,
+0.442240477,
+0.325714469,
+0.817954063,
+0.824210465,
+0.619122744
489+3,
+9,
+2,
+1,
+1,
+10,
+9
490var(
0)
= getVar(sample, iweight)
494 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
495 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
497varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
506 lb(isam)
= ub(isam
- 1)
+ 1
515+0.302187026,
+0.403615534,
+0.636334836,
+0.439223647,
+0.608977139,
+0.551796615
518+7,
+1,
+8,
+1,
+5,
+6
519var(
0)
= getVar(sample, iweight)
523 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
524 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
526varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
535 lb(isam)
= ub(isam
- 1)
+ 1
544+0.630253553,
+0.915598571,
+0.574551821E-1,
+0.878818750,
+0.508635521,
+0.545359790,
+0.559843779E-1,
+0.117216945,
+0.880172908,
+0.529791057,
+0.784381747
547+8,
+3,
+6,
+1,
+8,
+7,
+1,
+5,
+9,
+1,
+9
548var(
0)
= getVar(sample, iweight)
552 var(isam)
= getVar(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
553 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
555varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
569 lb(isam)
= ub(isam
- 1)
+ 1
578+0.911826611,
+0.678420842,
+0.241578221,
+0.291455209,
+0.881421745,
+0.462833643,
+0.529112577,
+0.102979422,
+0.205452263
581+1.12894595,
+1.98632443,
+1.10174203,
+1.31219566,
+1.29628932,
+1.80648267,
+1.92742968,
+1.41125298,
+1.57952511
582var(
0)
= getVar(sample, rweight)
586 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
587 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
589varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
598 lb(isam)
= ub(isam
- 1)
+ 1
607+0.892997205,
+0.973314822,
+0.266717613,
+0.949708164,
+0.277200639,
+0.874832392,
+0.977554679,
+0.665105820,
+0.217965484
610+1.70789027,
+1.55572355,
+1.11108840,
+1.58028388,
+1.53585505,
+1.16611719,
+1.23605156,
+1.70818591,
+1.33596683
611var(
0)
= getVar(sample, rweight)
615 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
616 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
618varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
627 lb(isam)
= ub(isam
- 1)
+ 1
636+0.762252808E-1,
+0.726244450E-1,
+0.747029185,
+0.969027877,
+0.935334861,
+0.989501536,
+0.640519679,
+0.570739210,
+0.515088022,
+0.572620869
639+1.41976547,
+1.06699586,
+1.46102047,
+1.67751622,
+1.08413863,
+1.15341854,
+1.15464854,
+1.84910917,
+1.43810034,
+1.45353150
640var(
0)
= getVar(sample, rweight)
644 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
645 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
647varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
656 lb(isam)
= ub(isam
- 1)
+ 1
665+0.178883016,
+0.933498085,
+0.399465799,
+0.527028441E-1,
+0.970194817
668+1.79133403,
+1.94495845,
+1.75885487,
+1.83700788,
+1.50504875
669var(
0)
= getVar(sample, rweight)
673 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
674 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
676varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
685 lb(isam)
= ub(isam
- 1)
+ 1
694+0.365253389,
+0.311100841,
+0.563895702,
+0.103803277,
+0.832665324
697+1.27340698,
+1.69329238,
+1.98719263,
+1.65814424,
+1.58595431
698var(
0)
= getVar(sample, rweight)
702 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
703 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
705varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
714 lb(isam)
= ub(isam
- 1)
+ 1
723+0.909076452,
+0.518222034,
+0.436950326,
+0.654068232,
+0.522867918,
+0.790539801,
+0.350945711,
+0.834322810,
+0.449866116,
+0.610065579
726+1.92141294,
+1.56948042,
+1.43730164,
+1.53060627,
+1.59464717,
+1.34180748,
+1.91020441,
+1.82805824,
+1.76202059,
+1.00744009
727var(
0)
= getVar(sample, rweight)
731 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
732 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
734varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
743 lb(isam)
= ub(isam
- 1)
+ 1
752+0.658622622,
+0.334244370,
+0.122356057,
+0.504886746,
+0.857854366
755+1.53603125,
+1.08443141,
+1.44010758,
+1.10636520,
+1.67062020
756var(
0)
= getVar(sample, rweight)
760 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
761 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
763varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
772 lb(isam)
= ub(isam
- 1)
+ 1
781+0.932872117,
+0.641629338,
+0.877119541,
+0.446920037,
+0.709850550,
+0.757925570,
+0.830348372,
+0.669314921
784+1.14050269,
+1.08093762,
+1.51856613,
+1.80405986,
+1.21332240,
+1.72392678,
+1.32793570,
+1.48838687
785var(
0)
= getVar(sample, rweight)
789 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
790 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
792varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
801 lb(isam)
= ub(isam
- 1)
+ 1
810+0.167209864,
+0.598109961E-1,
+0.134212732,
+0.333488286,
+0.864756107E-1,
+0.670296252,
+0.977488756E-1,
+0.655621290E-2,
+0.124229550,
+0.411810398
813+1.62893629,
+1.45038533,
+1.80866027,
+1.38803196,
+1.90581834,
+1.95609832,
+1.21783710,
+1.99536443,
+1.49530077,
+1.21245432
814var(
0)
= getVar(sample, rweight)
818 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
819 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
821varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
830 lb(isam)
= ub(isam
- 1)
+ 1
839+0.120754421,
+0.594324768,
+0.387402773,
+0.882316947,
+0.362297416,
+0.400876224,
+0.964011550,
+0.608533800,
+0.575663090,
+0.269784451,
+0.296675801,
+0.120380878
842+1.58186078,
+1.94900036,
+1.31332421,
+1.22539473,
+1.75035024,
+1.06165195,
+1.20171714,
+1.13796735,
+1.01178527,
+1.96979940,
+1.91439486,
+1.14375329
843var(
0)
= getVar(sample, rweight)
847 var(isam)
= getVar(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
848 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
850varMerged
= getVarMerged(var(
2), var(
1), mean(
1)
- mean(
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
864 lb(isam)
= ub(isam
- 1)
+ 1
872call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
873call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
877+0.822533369,
-0.201804757,
+0.185117960,
+0.926655412,
+0.898510456,
+0.854184985,
-0.198894858,
+0.246800184,
+0.257098556,
-0.657779813,
+0.649523735E-1
878-0.269520998,
-0.244022012,
-0.495082617,
-0.118978500,
+0.675124764,
-0.604276538,
+0.977416039E-1,
-0.935224771,
-0.769565582,
-0.327815652,
-0.323675990
879var(:,
0)
= getVar(sample, dim)
881+0.256479889,
+0.172895908
883 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
884 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
886varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
888+0.256479919,
+0.172895879
890+0.256479889,
+0.172895908
895 lb(isam)
= ub(isam
- 1)
+ 1
903call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
904call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
908+0.784901023,
-0.886325717,
+0.452937007,
-0.742119551E-1,
-0.546293736,
+0.917783022,
+0.463622808,
+0.296071410,
-0.133958697,
+0.855198860
909-0.178349018E-2,
-0.318868160,
+0.616620779,
+0.812220573E-2,
+0.662694097,
+0.828704715,
+0.459209561,
-0.170265436E-1,
-0.973475456,
+0.563909769
910var(:,
0)
= getVar(sample, dim)
912+0.335140675,
+0.275049657
914 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
915 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
917varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
919+0.335140675,
+0.275049627
921+0.335140675,
+0.275049657
926 lb(isam)
= ub(isam
- 1)
+ 1
934call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
935call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
939-0.771070242,
+0.114753842,
-0.316683888,
-0.300210595,
+0.403383136,
+0.710779905,
-0.787927747,
-0.677094460E-1,
-0.589005709,
-0.976577401,
-0.160291910
940+0.437226176,
+0.736133933,
-0.754325032,
-0.116185308,
-0.636162400,
-0.482977629E-1,
+0.594291925,
-0.546994328,
+0.868991137,
+0.338348985,
-0.172033668
941+0.221798420E-1,
+0.999202132,
+0.744336605,
+0.401337624,
+0.854566097E-1,
+0.143545985,
-0.813827395,
-0.670041323,
+0.722679019,
+0.192085505,
-0.812834501E-1
942-0.917914629,
-0.830164075,
-0.597445130,
+0.185401440,
+0.901916027,
-0.490061164,
-0.965566635E-1,
+0.884055257,
+0.676346183,
-0.291362405,
-0.378295302
943var(:,
0)
= getVar(sample, dim)
945+0.248637274,
+0.293594152,
+0.285629094,
+0.397294104
947 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
948 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
950varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
952+0.248637319,
+0.293594152,
+0.285629094,
+0.397294074
954+0.248637274,
+0.293594152,
+0.285629094,
+0.397294104
959 lb(isam)
= ub(isam
- 1)
+ 1
967call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
968call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
972+0.907283187,
+0.895410895,
+0.629059315,
-0.562860847,
+0.839654565,
-0.659306049,
-0.726124048,
-0.473100781,
+0.676214695E-1,
+0.378771424,
-0.253709197,
+0.327527642
973+0.317063570,
-0.187385082E-2,
-0.681246638,
+0.315501690,
+0.625920057,
-0.330458879E-1,
-0.609752417,
-0.734395981,
+0.196926594,
-0.786827326,
+0.898350477,
+0.845467329
974-0.897233605,
+0.404576778,
-0.795149207,
+0.408176303,
+0.347421646,
-0.169090033E-1,
+0.286790609,
+0.283926129,
+0.659218550,
-0.810095549,
-0.459394097,
+0.146639347
975-0.554223537,
+0.256314039,
-0.358138800,
-0.459490180,
+0.668642998,
+0.654584765,
+0.774138689,
-0.129163146,
-0.451373935,
+0.702569485,
-0.544707775,
+0.555702567
976-0.137287259,
-0.135267258,
-0.466312408,
-0.619010925E-1,
-0.324093103E-1,
-0.110737085,
-0.806897879E-1,
-0.446021676,
-0.330206275,
-0.725614309,
-0.888796568,
-0.493605852
977var(:,
0)
= getVar(sample, dim)
979+0.365954876,
+0.344797492,
+0.279876530,
+0.283589393,
+0.727569461E-1
981 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
982 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
984varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
986+0.365954936,
+0.344797522,
+0.279876560,
+0.283589393,
+0.727569610E-1
988+0.365954876,
+0.344797492,
+0.279876530,
+0.283589393,
+0.727569461E-1
993 lb(isam)
= ub(isam
- 1)
+ 1
1001call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1002call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1006+0.837407947,
-0.291133165,
+0.730831504,
+0.181254148E-1,
+0.238177776,
+0.837677240,
-0.703911781E-1,
+0.943988085,
-0.955191851E-1,
-0.292454958E-1
1007-0.874317527,
-0.965081692,
+0.829295158,
+0.738942623E-1,
-0.363408327,
+0.474955320,
+0.742375135,
-0.513280034,
+0.399335742,
+0.596296668
1008-0.911384463,
-0.165361166E-1,
-0.188089252,
-0.654661417,
-0.640970349,
-0.624086976,
+0.595846415,
+0.240228176,
-0.295981288,
+0.752123594E-1
1009var(:,
0)
= getVar(sample, dim)
1011+0.201154396,
+0.406026363,
+0.201534152
1013 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1014 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1016varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1018+0.201154366,
+0.406026363,
+0.201534137
1020+0.201154396,
+0.406026363,
+0.201534152
1025 lb(isam)
= ub(isam
- 1)
+ 1
1033call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1034call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1038+0.793565989,
+0.682521343,
-0.213110447E-1,
-0.220548153,
-0.605292916,
-0.656153202,
+0.276506782,
-0.641230226
1039-0.783159852,
+0.339449286,
+0.499422073,
+0.153867483,
+0.136586905,
-0.619032264,
+0.438970208,
+0.796170235E-1
1040var(:,
0)
= getVar(sample, dim)
1042+0.301253051,
+0.199375749
1044 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1045 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1047varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1049+0.301253080,
+0.199375764
1051+0.301253051,
+0.199375749
1056 lb(isam)
= ub(isam
- 1)
+ 1
1064call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1065call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1069-0.233726263,
+0.243409276,
-0.850091696,
+0.694321275,
+0.735538960,
+0.726227641,
-0.250964403,
+0.449663281
1070-0.129350901,
-0.768071175,
-0.770327449,
-0.187208772,
-0.795598865,
-0.979716659,
-0.741080523,
-0.546659231E-1
1071var(:,
0)
= getVar(sample, dim)
1073+0.295693815,
+0.116427571
1075 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1076 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1078varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1080+0.295693845,
+0.116427571
1082+0.295693815,
+0.116427571
1087 lb(isam)
= ub(isam
- 1)
+ 1
1095call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1096call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1100+0.157640338,
+0.685888290,
+0.667116046,
+0.782219768,
-0.533158422,
-0.829219818E-1,
+0.226912737,
+0.129403710,
-0.117244959,
-0.660314322,
-0.733315945E-1
1101var(:,
0)
= getVar(sample, dim)
1105 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1106 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1108varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1117 lb(isam)
= ub(isam
- 1)
+ 1
1125call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1126call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1130-0.837792277,
-0.716931462,
+0.793094635E-1,
+0.314378858,
+0.848489285,
+0.830800533E-1,
-0.549773335,
-0.464801788,
+0.171472073,
-0.786471963,
+0.265978813,
+0.950421095E-1,
+0.701242685E-1
1131-0.341052294,
-0.187232494,
+0.490494967,
-0.937614918,
-0.882568955,
+0.594423175,
+0.566866159,
-0.971008658,
-0.857906342E-1,
+0.902328134,
+0.215884447,
-0.122589350,
+0.603598356
1132+0.222824574,
+0.678692818,
-0.593200922,
-0.993439794,
-0.269551039,
-0.537260652,
+0.904311776,
+0.699999213,
-0.149301410,
+0.531911850,
-0.298003197,
-0.383806109,
+0.461785674
1133var(:,
0)
= getVar(sample, dim)
1135+0.241687492,
+0.377940953,
+0.328225732
1137 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1138 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1140varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1142+0.241687477,
+0.377941012,
+0.328225791
1144+0.241687492,
+0.377940953,
+0.328225732
1149 lb(isam)
= ub(isam
- 1)
+ 1
1157call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1158call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1162-0.887034059,
+0.971585512E-1,
+0.449972272,
+0.589379072E-1,
-0.770032287,
+0.403831482,
+0.688744545,
-0.487589121,
+0.463336825
1163-0.136800647,
+0.818308592,
+0.817033172,
+0.200883627,
-0.232460737,
+0.384944201,
+0.569100022,
-0.966175199,
+0.660362840
1164-0.400850534,
+0.482622504,
+0.847465158,
-0.753340483,
-0.932233572,
-0.960196733,
+0.634031653,
-0.497366428,
-0.455193639
1165var(:,
0)
= getVar(sample, dim)
1167+0.298334181,
+0.310532451,
+0.429640919
1169 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)), dim)
1170 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1172varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1174+0.298334152,
+0.310532451,
+0.429640949
1176+0.298334181,
+0.310532451,
+0.429640919
1186 lb(isam)
= ub(isam
- 1)
+ 1
1194call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1195call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1199-0.346815348,
-0.975849390,
+0.961728454,
-0.498436689E-1,
+0.330421448,
-0.333267450E-1,
+0.206635833,
-0.883876920
1200+0.141075373,
-0.623292327,
-0.831227303,
-0.602367163,
-0.938077092,
-0.193790674,
+0.887881994,
+0.650165439
1201+0.169270158,
+0.229446888E-1,
+0.480728030,
+0.586300850,
-0.562178493,
-0.660140872,
+0.419149041,
-0.325677276
1202+0.161718726,
+0.531843901E-1,
+0.882487416,
+0.228562117,
-0.490270257,
+0.275135636,
+0.947010517E-1,
-0.879087925
1205+7,
+7,
+3,
+1,
+8,
+4,
+10,
+7
1206var(:,
0)
= getVar(sample,
2_IK, iweight)
1208+0.333731890,
+0.495526791,
+0.168917492,
+0.213474572
1210 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1211 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1213varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1215+0.333731949,
+0.495526850,
+0.168917507,
+0.213474587
1217+0.333731890,
+0.495526791,
+0.168917492,
+0.213474572
1222 lb(isam)
= ub(isam
- 1)
+ 1
1230call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1231call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1235+0.996523380,
-0.337075353,
+0.660981417,
-0.592251658,
-0.744774818,
+0.103395224,
-0.325147510,
-0.838988304,
-0.149564028,
+0.992459655,
+0.167228341,
-0.941172719,
-0.679996967
1238+6,
+2,
+10,
+9,
+7,
+8,
+10,
+8,
+10,
+10,
+7,
+5,
+10
1239var(:,
0)
= getVar(sample,
2_IK, iweight)
1243 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1244 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1246varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1255 lb(isam)
= ub(isam
- 1)
+ 1
1263call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1264call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1268+0.184363127,
+0.348393202,
+0.526417613,
-0.357596278,
-0.811357856,
-0.642525315,
-0.211865783,
-0.503559589,
-0.298101902
1269+0.260335207E-1,
-0.941589832,
+0.606511831,
-0.795395136,
+0.473697662,
+0.523055553,
+0.454682112E-1,
+0.762981892,
+0.870145679
1272+7,
+7,
+6,
+4,
+8,
+2,
+3,
+6,
+8
1273var(:,
0)
= getVar(sample,
2_IK, iweight)
1275+0.202948958,
+0.405113071
1277 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1278 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1280varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1282+0.202948943,
+0.405113101
1284+0.202948958,
+0.405113071
1289 lb(isam)
= ub(isam
- 1)
+ 1
1297call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1298call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1302-0.884651184,
+0.244111896,
+0.590127826,
+0.292065501,
-0.900294185,
-0.539452672,
-0.871702194,
+0.317164063,
-0.292907476
1303+0.633534789,
+0.783125281,
-0.768205166,
-0.899906754,
-0.896088719,
+0.257887483,
-0.428955317,
+0.415202737,
-0.791176558
1304-0.688297391,
-0.197516441,
+0.167887926,
-0.289689898,
-0.192551613,
+0.760492921,
-0.347185135E-1,
+0.156588912,
-0.152474523
1307+2,
+4,
+5,
+8,
+5,
+9,
+9,
+6,
+4
1308var(:,
0)
= getVar(sample,
2_IK, iweight)
1310+0.306273907,
+0.366800100,
+0.143454388
1312 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1313 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1315varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1317+0.306273907,
+0.366800100,
+0.143454403
1319+0.306273907,
+0.366800100,
+0.143454388
1324 lb(isam)
= ub(isam
- 1)
+ 1
1332call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1333call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1337+0.792176247,
+0.247698784,
-0.120775580,
+0.415964007,
-0.284881592E-1,
-0.494365692E-1,
+0.899451017,
-0.954319715,
+0.636990070,
+0.730323792E-1,
+0.781589150,
+0.504717350
1338+0.194655299,
+0.965791941E-1,
-0.586743951,
+0.704725504,
-0.369609475,
+0.629681468,
+0.258773565,
+0.272131324,
+0.912307739,
+0.865791559,
-0.803885102,
-0.183542252
1339+0.442858934,
-0.847361207,
+0.661986351,
-0.212518334,
+0.836201668,
+0.806870103,
-0.761285543,
-0.938312531,
+0.723911881,
-0.670963526,
-0.238749623,
+0.248911381E-1
1340-0.661587596,
-0.655667543,
+0.649488091,
+0.144449711,
+0.920401812E-1,
+0.966109753,
+0.367500305,
-0.228299975,
-0.774057388,
+0.896820307,
+0.503402114,
+0.103186369E-1
1343+10,
+8,
+7,
+1,
+4,
+10,
+10,
+6,
+6,
+4,
+8,
+8
1344var(:,
0)
= getVar(sample,
2_IK, iweight)
1346+0.261162996,
+0.260256201,
+0.446453691,
+0.363874108
1348 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1349 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1351varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1353+0.261162996,
+0.260256201,
+0.446453750,
+0.363874108
1355+0.261162996,
+0.260256201,
+0.446453691,
+0.363874108
1360 lb(isam)
= ub(isam
- 1)
+ 1
1368call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1369call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1373+0.140013695E-1,
-0.776290894E-1,
-0.472523689,
-0.684229851,
-0.373419404,
+0.521303415E-1,
-0.971325397,
-0.307337046,
+0.638170123,
-0.417646766,
+0.301136494
1374-0.374989510,
-0.211510181,
-0.640881062E-2,
+0.395775914,
+0.134289265E-2,
+0.605209351,
-0.868660927,
+0.923526287,
-0.989169359,
-0.165797114,
+0.176223159
1377+9,
+5,
+10,
+5,
+3,
+5,
+1,
+4,
+8,
+8,
+3
1378var(:,
0)
= getVar(sample,
2_IK, iweight)
1380+0.164171442,
+0.261275947
1382 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1383 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1385varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1387+0.164171472,
+0.261275977
1389+0.164171442,
+0.261275947
1394 lb(isam)
= ub(isam
- 1)
+ 1
1402call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1403call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1407-0.461739540,
+0.360735059,
-0.173306227,
+0.136063457,
-0.729244232,
-0.446571589
1408-0.643176913,
+0.262316346,
+0.762135625,
-0.713714957,
-0.547412157,
-0.291866064
1411+1,
+9,
+9,
+10,
+9,
+2
1412var(:,
0)
= getVar(sample,
2_IK, iweight)
1414+0.160794541,
+0.345188975
1416 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1417 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1419varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1421+0.160794526,
+0.345188975
1423+0.160794541,
+0.345188975
1428 lb(isam)
= ub(isam
- 1)
+ 1
1436call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1437call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1441+0.735610485,
-0.993516803,
-0.554462671E-1,
-0.911352158,
+0.218448043,
-0.296621442,
+0.770020485,
-0.482940316
1444+8,
+10,
+6,
+9,
+7,
+10,
+9,
+5
1445var(:,
0)
= getVar(sample,
2_IK, iweight)
1449 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1450 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1452varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1461 lb(isam)
= ub(isam
- 1)
+ 1
1469call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1470call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1474+0.613574386,
-0.222642422,
-0.313583851,
+0.890709400,
-0.515742302,
-0.363742828,
+0.778588891,
-0.461034775,
+0.129873395,
-0.482453704
1477+6,
+5,
+4,
+3,
+10,
+7,
+1,
+3,
+2,
+3
1478var(:,
0)
= getVar(sample,
2_IK, iweight)
1482 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1483 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1485varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1494 lb(isam)
= ub(isam
- 1)
+ 1
1502call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1503call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1507-0.954358459,
-0.950636029,
+0.291834116,
+0.994037390,
-0.121437311,
+0.614416599E-2,
-0.804266214,
+0.487909436,
+0.456443667,
+0.105119705,
+0.631949186,
+0.988863826,
+0.329134822
1508-0.366969109E-1,
+0.583245873,
+0.327589750,
-0.303538680,
+0.260551572,
-0.442986488,
+0.698370099,
-0.194320440,
-0.909412265,
+0.226153851,
+0.299042821,
-0.789839268,
+0.526757479
1509+0.765581608,
-0.422755241,
+0.608958125,
+0.462303162,
-0.205655098E-1,
+0.290399790,
+0.831194997,
+0.488145351E-1,
-0.119636059E-1,
+0.669911742,
+0.657019138,
+0.759059668,
+0.344942570
1512+10,
+9,
+8,
+4,
+6,
+4,
+9,
+5,
+9,
+8,
+9,
+6,
+2
1513var(:,
0)
= getVar(sample,
2_IK, iweight)
1515+0.456371874,
+0.256696731,
+0.164911643
1517 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1518 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1520varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1522+0.456371874,
+0.256696761,
+0.164911643
1524+0.456371874,
+0.256696731,
+0.164911643
1534 lb(isam)
= ub(isam
- 1)
+ 1
1542call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1543call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1547+0.660363078,
-0.966759086,
-0.860092759,
-0.356497884,
-0.644602895
1548-0.533545017E-1,
+0.950456858,
-0.800315857,
+0.460979104,
+0.588412404
1551+1.57147098,
+1.79876924,
+1.14372087,
+1.13252592,
+1.82317734
1552var(:,
0)
= getVar(sample,
2_IK, rweight)
1554+0.359834194,
+0.337686747
1556 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1557 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1559varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1561+0.359834164,
+0.337686747
1563+0.359834194,
+0.337686747
1568 lb(isam)
= ub(isam
- 1)
+ 1
1576call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1577call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1581-0.791878700,
-0.491848469,
-0.960410953,
-0.100202799,
+0.893252730,
-0.587533832,
+0.475235939,
+0.610859752,
+0.706519961,
+0.591389418
1582-0.348970175,
-0.442520380,
+0.872310042,
+0.584486842,
+0.967864990E-1,
-0.393078208,
+0.565749526,
+0.876973033,
+0.873039484,
+0.733671665
1583+0.205409646,
-0.975620985,
+0.444479227,
-0.456634045,
+0.832700610,
-0.122840166,
+0.894105792,
+0.269454122,
+0.644260287,
-0.979103208
1584+0.172781944,
-0.797153831,
+0.922139406,
-0.428922176E-1,
+0.444456697,
+0.151862741,
-0.842796326,
+0.828892708,
+0.411669493,
+0.163495660
1587+1.24288142,
+1.26046181,
+1.58228040,
+1.67941308,
+1.19372642,
+1.60810089,
+1.29749918,
+1.79418826,
+1.83854413,
+1.86654139
1588var(:,
0)
= getVar(sample,
2_IK, rweight)
1590+0.425889641,
+0.272481233,
+0.419787347,
+0.289969563
1592 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1593 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1595varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1597+0.425889671,
+0.272481233,
+0.419787407,
+0.289969563
1599+0.425889641,
+0.272481233,
+0.419787347,
+0.289969563
1604 lb(isam)
= ub(isam
- 1)
+ 1
1612call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1613call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1617+0.693080068,
+0.380644560,
+0.389378667,
+0.397407651,
+0.901105642,
+0.882100344,
-0.511245847,
+0.842516065,
-0.100411177,
-0.207977295E-1,
-0.707914472,
-0.862921238,
+0.537495375
1618-0.346648693E-1,
-0.215473890,
+0.720063448E-1,
+0.316514850,
+0.403852701,
+0.975278497,
+0.628647685,
+0.580743670,
+0.246021152,
-0.453594923,
+0.974860907,
+0.999689102,
-0.191987038
1619+0.362644196,
+0.591808677,
+0.515258908,
+0.386660457,
-0.227061629,
-0.841854453,
-0.110002756,
+0.876263976,
+0.429481030,
-0.447105527,
+0.650120497,
+0.709533453,
-0.454079509
1620-0.473112345,
+0.651614308,
+0.766399860,
+0.262941003,
-0.192683578,
-0.271325111,
-0.696439862,
+0.258036613,
-0.573599100,
-0.905174017E-1,
+0.292499542,
+0.708933115,
+0.633116961
1621-0.168622732E-1,
+0.764720440E-1,
-0.535539389,
+0.501987219,
+0.962664723,
-0.365248442,
-0.859807730E-1,
+0.870487928,
+0.155248761,
-0.447557092,
+0.807462215,
-0.913901329E-1,
-0.595916629
1622+0.250450134,
-0.635687232,
-0.719767094,
+0.777232051,
-0.712450147,
+0.800201893,
-0.236097932,
+0.810490489,
-0.561327577,
+0.430908442,
+0.934778214,
-0.778853416,
-0.880653858
1625+1.09385085,
+1.36129546,
+1.17108357,
+1.61249471,
+1.46207738,
+1.89953399,
+1.02423453,
+1.90906286,
+1.98339987,
+1.90777624,
+1.63946795,
+1.79679871,
+1.91663694
1626var(:,
0)
= getVar(sample,
2_IK, rweight)
1628+0.345990241,
+0.232519254,
+0.299046457,
+0.224219635,
+0.279615343,
+0.503155470
1630 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1631 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1633varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1635+0.345990270,
+0.232519254,
+0.299046457,
+0.224219665,
+0.279615343,
+0.503155529
1637+0.345990241,
+0.232519254,
+0.299046457,
+0.224219635,
+0.279615343,
+0.503155470
1642 lb(isam)
= ub(isam
- 1)
+ 1
1650call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1651call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1655+0.756149054,
-0.774503589,
-0.671933651,
+0.986950755,
-0.499761939,
-0.419576168E-1,
-0.480232477,
+0.862828612
1658+1.37810230,
+1.20350397,
+1.84911561,
+1.26827312,
+1.88262701,
+1.24688554,
+1.93623257,
+1.51983833
1659var(:,
0)
= getVar(sample,
2_IK, rweight)
1663 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1664 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1666varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1675 lb(isam)
= ub(isam
- 1)
+ 1
1683call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1684call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1688-0.151960611,
-0.758148432,
+0.898804069,
+0.413609624,
+0.152255177,
+0.475582719,
+0.910844803,
-0.148987293,
+0.669309855,
-0.724211097,
-0.650645018
1689-0.479493976,
-0.699717879,
-0.154049277,
-0.858289242,
-0.916954398,
-0.304045558,
+0.641827703,
-0.954472065,
-0.527627110,
-0.466414332,
+0.712178707
1690-0.275203347,
+0.989192605,
-0.126069307,
-0.276707411E-1,
+0.214401841,
-0.615044832E-1,
+0.192118645,
-0.587381244,
+0.208611965,
-0.176807046,
-0.832059741
1693+1.04653037,
+1.30666733,
+1.85078406,
+1.30741811,
+1.10350037,
+1.88958347,
+1.50808024,
+1.12076640,
+1.56512475,
+1.25214911,
+1.48790658
1694var(:,
0)
= getVar(sample,
2_IK, rweight)
1696+0.376059651,
+0.295622885,
+0.194241762
1698 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1699 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1701varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1703+0.376059651,
+0.295622855,
+0.194241747
1705+0.376059651,
+0.295622885,
+0.194241762
1710 lb(isam)
= ub(isam
- 1)
+ 1
1718call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1719call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1723-0.299291253,
+0.862094164,
+0.214979172,
-0.552385569,
-0.726424932,
+0.252487063,
+0.326475739,
+0.711191654,
-0.774406672,
-0.234422445,
-0.180139065
1724-0.552764893,
-0.253407121,
-0.150469661,
+0.203090191,
+0.221945047,
+0.262630105,
-0.433197021,
+0.864746094,
+0.104330301,
+0.674487352,
+0.899980068E-1
1725+0.730610251,
-0.427738070,
+0.668025970,
+0.531019568,
-0.730203867,
-0.137626171,
-0.700907826,
+0.761593103,
+0.606769919,
+0.607911348,
+0.846730828
1728+1.71857429,
+1.43184614,
+1.12418556,
+1.39916444,
+1.68876648,
+1.52994049,
+1.75960684,
+1.10115457,
+1.34660971,
+1.89873075,
+1.71198916
1729var(:,
0)
= getVar(sample,
2_IK, rweight)
1731+0.258907527,
+0.169161752,
+0.365339160
1733 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1734 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1736varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1738+0.258907527,
+0.169161752,
+0.365339160
1740+0.258907527,
+0.169161752,
+0.365339160
1745 lb(isam)
= ub(isam
- 1)
+ 1
1753call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1754call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1758+0.727317214,
-0.444257975,
-0.845570683,
-0.235957503,
-0.861949801,
-0.228896856,
+0.879729629,
+0.530003190,
+0.947062969E-1,
+0.260903120
1761+1.16687036,
+1.07370543,
+1.71065629,
+1.92600822,
+1.54404330,
+1.66079140,
+1.08190048,
+1.65655899,
+1.21597373,
+1.73632741
1762var(:,
0)
= getVar(sample,
2_IK, rweight)
1766 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1767 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1769varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1778 lb(isam)
= ub(isam
- 1)
+ 1
1786call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1787call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1791+0.190700412,
-0.513439417,
-0.997154713,
-0.503218651,
+0.957273960,
+0.287846088,
+0.357133389
1794+1.94418788,
+1.43793130,
+1.50021744,
+1.48280120,
+1.04479861,
+1.82410073,
+1.26461220
1795var(:,
0)
= getVar(sample,
2_IK, rweight)
1799 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1800 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1802varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1811 lb(isam)
= ub(isam
- 1)
+ 1
1819call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1820call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1824+0.845301390,
-0.305573821,
+0.674729347E-1,
+0.892303944,
+0.279717922,
-0.799050927
1825-0.924136758,
+0.278034925,
-0.953253865,
+0.399394274,
-0.262610793,
-0.783197761
1828+1.44209433,
+1.96993470,
+1.90179718,
+1.97132242,
+1.40143275,
+1.01933360
1829var(:,
0)
= getVar(sample,
2_IK, rweight)
1831+0.320065200,
+0.333336234
1833 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1834 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1836varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1838+0.320065200,
+0.333336264
1840+0.320065200,
+0.333336234
1845 lb(isam)
= ub(isam
- 1)
+ 1
1853call setRebound(var, [
1_IK,
0_IK], [ndim, nsam])
1854call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1858+0.865042210E-1,
+0.537455678,
-0.533723950,
+0.667766571,
+0.107939363,
+0.363443971,
+0.951228142E-1,
-0.645814896,
-0.645524263
1859+0.337412357E-1,
-0.675078392,
-0.374974608,
+0.435822606,
-0.226102591,
-0.514562845,
+0.380314350,
+0.183895230,
-0.222754359
1862+1.75644445,
+1.22698855,
+1.81325841,
+1.67027819,
+1.14698446,
+1.45705974,
+1.58294988,
+1.00524092,
+1.64840508
1863var(:,
0)
= getVar(sample,
2_IK, rweight)
1865+0.222054631,
+0.134773955
1867 var(:,isam)
= getVar(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1868 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1870varMerged
= getVarMerged(var(:,
2), var(:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1872+0.222054631,
+0.134773955
1874+0.222054631,
+0.134773955