15 integer(IK) ,
parameter :: nsam
= 2
16 integer(IK) ,
allocatable :: iweight(:)
17 integer(IK) :: isam, ndim, lb(nsam), ub(nsam)
18 integer(IK) :: dim, itry, ntry
= 10
19 type(display_type) :: disp
23 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
24 call disp%show(
"!Compute the merged mean of a univariate sample.")
25 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
29 real(TKG) :: mean(
0:nsam), meanMerged
30 real(TKG),
allocatable :: sample(:)
33 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(1, 7)")
36 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
37 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(1, 7)")
40 lb(isam)
= ub(isam
- 1)
+ 1
47 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
51 call disp%show(
"mean(0) = getMean(sample)")
53 call disp%show(
"mean(0) ! reference")
56 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)))")
59 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
61 call disp%show(
"call setMeanMerged(meanMerged, mean(2), mean(1), ub(1) / real(ub(2), TKG))")
62 call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
65 call disp%show(
"call setMeanMerged(mean(2), mean(1), ub(1) / real(ub(2), TKG))")
66 call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
69 call disp%show(
"mean(0) ! reference")
76 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
77 call disp%show(
"!Compute the merged mean of a frequency weighted univariate sample.")
78 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
82 real(TKG) :: mean(
0:nsam), meanMerged
83 real(TKG),
allocatable :: sample(:)
86 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(1, 7)")
89 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
90 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(1, 7)")
93 lb(isam)
= ub(isam
- 1)
+ 1
100 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
104 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, 1, IK))")
108 call disp%show(
"mean(0) = getMean(sample, iweight)")
109 mean(
0)
= getMean(sample, iweight)
110 call disp%show(
"mean(0) ! reference")
113 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))")
116 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
118 call disp%show(
"call setMeanMerged(meanMerged, mean(2), mean(1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
119 call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
122 call disp%show(
"call setMeanMerged(mean(2), mean(1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
123 call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
126 call disp%show(
"mean(0) ! reference")
133 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
134 call disp%show(
"!Compute the merged mean of a reliability weighted univariate sample.")
135 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
139 real(TKG) :: mean(
0:nsam), meanMerged
140 real(TKG),
allocatable :: sample(:)
141 real(TKG),
allocatable :: rweight(:)
144 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(1, 7)")
147 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
148 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(1, 7)")
151 lb(isam)
= ub(isam
- 1)
+ 1
158 call disp%show(
"sample = getUnifRand(0., 1., ub(nsam))")
162 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, 1, IK))")
166 call disp%show(
"mean(0) = getMean(sample, rweight)")
167 mean(
0)
= getMean(sample, rweight)
168 call disp%show(
"mean(0) ! reference")
171 call disp%show(
" mean(isam) = getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))")
174 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
176 call disp%show(
"call setMeanMerged(meanMerged, mean(2), mean(1), real(sum(rweight(:ub(1))) / sum(rweight), TKG))")
177 call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
180 call disp%show(
"call setMeanMerged(mean(2), mean(1), real(sum(rweight(:ub(1))) / sum(rweight), TKG))")
181 call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
184 call disp%show(
"mean(0) ! reference")
191 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
192 call disp%show(
"!Compute the merged mean of a multivariate sample.")
193 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
197 real(TKG),
allocatable :: mean(:,:), meanMerged(:)
198 real(TKG),
allocatable :: sample(:,:)
201 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(1, 7)")
204 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
205 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(1, 7)")
208 lb(isam)
= ub(isam
- 1)
+ 1
215 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
217 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
218 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
219 call disp%show(
"call setResized(meanMerged, ndim)")
221 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
225 call disp%show(
"mean(:,0) = getMean(sample, dim)")
226 mean(:,
0)
= getMean(sample, dim)
227 call disp%show(
"mean(:,0) ! reference")
230 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
233 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
235 call disp%show(
"call setMeanMerged(meanMerged, mean(:,2), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG))")
236 call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
239 call disp%show(
"call setMeanMerged(mean(:,2), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG))")
240 call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
243 call disp%show(
"mean(:,0) ! reference")
250 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
251 call disp%show(
"!Compute the merged mean of a frequency weighted multivariate sample.")
252 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
256 real(TKG),
allocatable :: mean(:,:), meanMerged(:)
257 real(TKG),
allocatable :: sample(:,:)
260 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(1, 7)")
263 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
264 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(1, 7)")
267 lb(isam)
= ub(isam
- 1)
+ 1
274 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
276 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
277 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
278 call disp%show(
"call setResized(meanMerged, ndim)")
280 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
284 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
288 call disp%show(
"mean(:,0) = getMean(sample, 2_IK, iweight)")
289 mean(:,
0)
= getMean(sample,
2_IK, iweight)
290 call disp%show(
"mean(:,0) ! reference")
293 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
296 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
298 call disp%show(
"call setMeanMerged(meanMerged, mean(:,2), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
299 call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
302 call disp%show(
"call setMeanMerged(mean(:,2), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
303 call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
306 call disp%show(
"mean(:,0) ! reference")
313 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
314 call disp%show(
"!Compute the merged mean of a reliability weighted multivariate sample.")
315 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
319 real(TKG),
allocatable :: mean(:,:), meanMerged(:)
320 real(TKG),
allocatable :: sample(:,:)
321 real(TKG),
allocatable :: rweight(:)
324 call disp%show(
"lb(1) = 1; ub(1) = getUnifRand(1, 7)")
327 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
328 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(1, 7)")
331 lb(isam)
= ub(isam
- 1)
+ 1
338 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
340 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
341 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
342 call disp%show(
"call setResized(meanMerged, ndim)")
344 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
348 call disp%show(
"rweight = getUnifRand(1, 10, size(sample, dim, IK))")
352 call disp%show(
"mean(:,0) = getMean(sample, dim, rweight)")
353 mean(:,
0)
= getMean(sample, dim, rweight)
354 call disp%show(
"mean(:,0) ! reference")
357 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))")
360 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
362 call disp%show(
"call setMeanMerged(meanMerged, mean(:,2), mean(:,1), real(sum(rweight(:ub(1))) / sum(rweight), TKG))")
363 call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
366 call disp%show(
"call setMeanMerged(mean(:,2), mean(:,1), real(sum(rweight(:ub(1))) / sum(rweight), TKG))")
367 call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
370 call disp%show(
"mean(:,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...
Return the (weighted) merged mean of a sample resulting from the merger of two separate (weighted) sa...
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.
Generate and return an object of type display_type.
9 lb(isam)
= ub(isam
- 1)
+ 1
18+0.970227480,
+0.889174044,
+0.400529444,
+0.108696580,
+0.360504448,
+0.611867964,
+0.340796709E-1,
+0.398242354,
+0.731139064
23 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
25call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
37 lb(isam)
= ub(isam
- 1)
+ 1
46+0.988277555,
+0.352771282,
+0.193285286
51 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
53call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
65 lb(isam)
= ub(isam
- 1)
+ 1
74+0.434534550,
+0.703356326,
+0.483815253,
+0.922264576,
+0.775594234,
+0.500287414E-1,
+0.699615538,
+0.463437378,
+0.469291031
79 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
81call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
93 lb(isam)
= ub(isam
- 1)
+ 1
102+0.743840456,
+0.924445868,
+0.574701309
107 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
109call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
112call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
121 lb(isam)
= ub(isam
- 1)
+ 1
130+0.502027273E-1,
+0.385735393
135 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
137call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
140call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
149 lb(isam)
= ub(isam
- 1)
+ 1
158+0.523915768,
+0.731165648,
+0.849807918,
+0.303806782,
+0.393536031,
+0.887881219,
+0.285871685,
+0.907886326,
+0.840399981,
+0.298650146,
+0.518002212,
+0.303899825,
+0.847452641
163 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
165call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
168call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
177 lb(isam)
= ub(isam
- 1)
+ 1
186+0.559693992,
+0.321461499,
+0.423072040,
+0.928239048,
+0.113605142,
+0.384397566,
+0.105644166,
+0.674464464
191 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
193call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
196call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
205 lb(isam)
= ub(isam
- 1)
+ 1
214+0.638076544,
+0.192362010,
+0.272083640,
+0.726456702,
+0.531597555,
+0.298390865,
+0.611184180,
+0.677358985,
+0.449197292E-1
219 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
221call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
224call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
233 lb(isam)
= ub(isam
- 1)
+ 1
242+0.199613988,
+0.305162668,
+0.552734494,
+0.722108781,
+0.113323748,
+0.889423549,
+0.475358486,
+0.518952370,
+0.725527227
247 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
249call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
252call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
261 lb(isam)
= ub(isam
- 1)
+ 1
270+0.451303363,
+0.809018612E-1,
+0.106862247,
+0.504584491,
+0.960011005,
+0.395677865
275 mean(isam)
= getMean(sample(lb(isam):ub(isam)))
277call setMeanMerged(meanMerged, mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
280call setMeanMerged(mean(
2), mean(
1), ub(
1)
/ real(ub(
2), TKG))
294 lb(isam)
= ub(isam
- 1)
+ 1
303+0.931036353,
+0.360781252,
+0.658607483E-2,
+0.993314683,
+0.382083535,
+0.617320836,
+0.659363270E-1
306+6,
+10,
+4,
+7,
+9,
+8,
+2
307mean(
0)
= getMean(sample, iweight)
311 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
313call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
316call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
325 lb(isam)
= ub(isam
- 1)
+ 1
334+0.196274102,
+0.473127067,
+0.356733918,
+0.475377560,
+0.721541822,
+0.683517039,
+0.504666090,
+0.650321841,
+0.106697500,
+0.506007314
337+7,
+8,
+7,
+10,
+9,
+3,
+10,
+3,
+8,
+3
338mean(
0)
= getMean(sample, iweight)
342 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
344call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
347call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
356 lb(isam)
= ub(isam
- 1)
+ 1
365+0.666209996,
+0.372503519,
+0.644750595,
+0.935250759,
+0.491415143
369mean(
0)
= getMean(sample, iweight)
373 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
375call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
378call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
387 lb(isam)
= ub(isam
- 1)
+ 1
396+0.625750959,
+0.610696912,
+0.474483311,
+0.465015531,
+0.335899472,
+0.572187483,
+0.389681578,
+0.793304324,
+0.810409427
399+4,
+7,
+3,
+4,
+7,
+9,
+2,
+8,
+4
400mean(
0)
= getMean(sample, iweight)
404 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
406call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
409call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
418 lb(isam)
= ub(isam
- 1)
+ 1
427+0.775321066,
+0.794384539,
+0.403132200,
+0.215430498,
+0.133064210,
+0.293928564,
+0.561876714,
+0.731466770,
+0.475014567,
+0.935007513
430+2,
+4,
+4,
+4,
+10,
+6,
+6,
+10,
+4,
+8
431mean(
0)
= getMean(sample, iweight)
435 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
437call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
440call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
449 lb(isam)
= ub(isam
- 1)
+ 1
458+0.165893912,
+0.464924753,
+0.119281590,
+0.745239377,
+0.130114675,
+0.415790260
461+7,
+3,
+6,
+10,
+8,
+1
462mean(
0)
= getMean(sample, iweight)
466 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
468call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
471call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
480 lb(isam)
= ub(isam
- 1)
+ 1
489+0.558915734E-1,
+0.837118387,
+0.350727201,
+0.850502491,
+0.317502201,
+0.982917964,
+0.354779959E-1,
+0.796154797,
+0.596565247,
+0.725060821
492+6,
+4,
+5,
+1,
+3,
+4,
+6,
+2,
+8,
+2
493mean(
0)
= getMean(sample, iweight)
497 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
499call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
502call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
511 lb(isam)
= ub(isam
- 1)
+ 1
520+0.784939110,
+0.192044199,
+0.957277656,
+0.725010097,
+0.363101900,
+0.509126782E-1,
+0.925499678,
+0.814981043,
+0.230785012
523+2,
+8,
+6,
+1,
+1,
+6,
+4,
+1,
+5
524mean(
0)
= getMean(sample, iweight)
528 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
530call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
533call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
542 lb(isam)
= ub(isam
- 1)
+ 1
551+0.449719250,
+0.450491786,
+0.481030107,
+0.553094089,
+0.118318498,
+0.138391972,
+0.626205504
554+3,
+7,
+6,
+4,
+6,
+8,
+7
555mean(
0)
= getMean(sample, iweight)
559 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
561call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
564call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
573 lb(isam)
= ub(isam
- 1)
+ 1
582+0.464345932,
+0.421162665,
+0.924247861,
+0.947289884,
+0.879936099,
+0.106137991E-1,
+0.312648833,
+0.540764749
585+2,
+5,
+4,
+4,
+4,
+6,
+10,
+6
586mean(
0)
= getMean(sample, iweight)
590 mean(isam)
= getMean(sample(lb(isam):ub(isam)), iweight(lb(isam):ub(isam)))
592call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
595call setMeanMerged(mean(
2), mean(
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
609 lb(isam)
= ub(isam
- 1)
+ 1
618+0.336420298,
+0.900482416,
+0.269200981,
+0.801517010,
+0.172469258,
+0.463568211,
+0.422131181,
+0.408791304,
+0.407567441,
+0.581609607E-1
621+1.52702785,
+1.63808441,
+1.13411498,
+1.19323611,
+1.16349483,
+1.13752007,
+1.67987084,
+1.05819869,
+1.07211399,
+1.47440314
622mean(
0)
= getMean(sample, rweight)
626 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
628call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
631call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
640 lb(isam)
= ub(isam
- 1)
+ 1
649+0.536305487,
+0.286720932,
+0.871705592,
+0.769569814,
+0.206260383,
+0.863130748,
+0.832658350,
+0.107159734,
+0.559514821
652+1.82568407,
+1.86250114,
+1.23210704,
+1.05506516,
+1.40587664,
+1.62315893,
+1.33776939,
+1.55490136,
+1.99593031
653mean(
0)
= getMean(sample, rweight)
657 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
659call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
662call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
671 lb(isam)
= ub(isam
- 1)
+ 1
680+0.389442444,
+0.455424368,
+0.606934428,
+0.814193249
683+1.82544494,
+1.11103487,
+1.85694730,
+1.53554177
684mean(
0)
= getMean(sample, rweight)
688 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
690call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
693call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
702 lb(isam)
= ub(isam
- 1)
+ 1
711+0.847686410,
+0.987932384,
+0.198018849,
+0.307114780,
+0.762261212,
+0.389709711,
+0.909324944,
+0.977158546E-1,
+0.604008555,
+0.956929982,
+0.456336141E-1,
+0.382911325,
+0.955774307
714+1.42221522,
+1.76652336,
+1.58145356,
+1.16595054,
+1.66225612,
+1.79025745,
+1.23157299,
+1.34393477,
+1.75656569,
+1.40479469,
+1.89837289,
+1.52191210,
+1.41057038
715mean(
0)
= getMean(sample, rweight)
719 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
721call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
724call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
733 lb(isam)
= ub(isam
- 1)
+ 1
742+0.850012362,
+0.779312849,
+0.328576744,
+0.336791992,
+0.523952425,
+0.868522763,
+0.298385799,
+0.240834951E-1,
+0.172997057,
+0.719876945,
+0.482980490,
+0.681391358E-1
745+1.24348593,
+1.90975738,
+1.46687090,
+1.51497221,
+1.89248204,
+1.73911059,
+1.69809842,
+1.05052841,
+1.60351896,
+1.87277389,
+1.34884048,
+1.04988599
746mean(
0)
= getMean(sample, rweight)
750 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
752call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
755call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
764 lb(isam)
= ub(isam
- 1)
+ 1
773+0.876530707,
+0.157267094,
+0.883412004,
+0.154815853,
+0.122767150,
+0.352350056,
+0.487301290,
+0.369765401,
+0.905390680,
+0.301889539,
+0.915843129,
+0.521158874,
+0.933727980
776+1.31678510,
+1.94649076,
+1.80073905,
+1.52297151,
+1.16900313,
+1.00780749,
+1.23268461,
+1.50199580,
+1.00478733,
+1.70807457,
+1.21136546,
+1.25116885,
+1.21510315
777mean(
0)
= getMean(sample, rweight)
781 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
783call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
786call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
795 lb(isam)
= ub(isam
- 1)
+ 1
804+0.746999264,
+0.772125304,
+0.595147133,
+0.710202754,
+0.314260244,
+0.920885563,
+0.668401062,
+0.432263553
807+1.89723039,
+1.12448907,
+1.19223714,
+1.42174125,
+1.12921214,
+1.54256248,
+1.05259728,
+1.96224689
808mean(
0)
= getMean(sample, rweight)
812 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
814call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
817call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
826 lb(isam)
= ub(isam
- 1)
+ 1
835+0.939530909,
+0.367635071,
+0.766378045,
+0.414034128,
+0.259239018,
+0.915690482,
+0.369608164
838+1.69386923,
+1.45650268,
+1.58792543,
+1.20705438,
+1.69190228,
+1.56793857,
+1.01857913
839mean(
0)
= getMean(sample, rweight)
843 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
845call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
848call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
857 lb(isam)
= ub(isam
- 1)
+ 1
866+0.834842503,
+0.587104321,
+0.409085393,
+0.392835975,
+0.117182434,
+0.500619292,
+0.995307744,
+0.477947831
869+1.95303488,
+1.15331650,
+1.87951970,
+1.28134584,
+1.88879061,
+1.29159021,
+1.71719456,
+1.38875723
870mean(
0)
= getMean(sample, rweight)
874 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
876call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
879call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
888 lb(isam)
= ub(isam
- 1)
+ 1
897+0.209422886,
+0.320493460,
+0.593413830,
+0.574087381,
+0.248882771E-1,
+0.564098537
900+1.65988326,
+1.71720719,
+1.86698461,
+1.05854750,
+1.91033769,
+1.99657559
901mean(
0)
= getMean(sample, rweight)
905 mean(isam)
= getMean(sample(lb(isam):ub(isam)), rweight(lb(isam):ub(isam)))
907call setMeanMerged(meanMerged, mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
910call setMeanMerged(mean(
2), mean(
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
924 lb(isam)
= ub(isam
- 1)
+ 1
932call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
936+0.886028290,
+0.453812718,
-0.889998317,
-0.224975824
937mean(:,
0)
= getMean(sample, dim)
941 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
943call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
946call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
955 lb(isam)
= ub(isam
- 1)
+ 1
963call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
967+0.307964563,
+0.926474452,
+0.833592534,
-0.374875188,
+0.880347013,
+0.229835510E-3,
-0.374148250,
+0.428916812
968mean(:,
0)
= getMean(sample, dim)
972 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
974call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
977call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
986 lb(isam)
= ub(isam
- 1)
+ 1
994call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
998-0.214012623,
-0.400241256,
-0.444279909,
-0.655349255,
-0.230083466,
-0.128688693,
+0.103243113,
+0.799262524
999mean(:,
0)
= getMean(sample, dim)
1003 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1005call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1008call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1017 lb(isam)
= ub(isam
- 1)
+ 1
1025call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1029-0.961787224,
-0.559598684,
+0.768479586,
+0.691591859,
+0.405498743E-1,
+0.618446350
1030mean(:,
0)
= getMean(sample, dim)
1034 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1036call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1039call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1048 lb(isam)
= ub(isam
- 1)
+ 1
1056call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1060+0.636579394,
+0.233046174,
+0.716617703,
-0.461002350,
+0.853260994,
-0.512781858,
-0.700733304,
-0.631747961,
+0.554563046
1061+0.727544308,
+0.820618272,
-0.868051648,
-0.429713845,
-0.367620945,
+0.701023936,
-0.350947022,
-0.755479336,
-0.446876287E-1
1062-0.608101010,
-0.955167890,
+0.978379488,
-0.452326655,
+0.998688698,
-0.533783436E-1,
+0.919018388,
-0.635264516,
-0.449159265
1063mean(:,
0)
= getMean(sample, dim)
1065+0.764224306E-1,
-0.630348772E-1,
-0.285901222E-1
1067 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1069call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1071+0.764224529E-1,
-0.630348772E-1,
-0.285901204E-1
1072call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1074+0.764224529E-1,
-0.630348772E-1,
-0.285901204E-1
1076+0.764224306E-1,
-0.630348772E-1,
-0.285901222E-1
1081 lb(isam)
= ub(isam
- 1)
+ 1
1089call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1093+0.453999996,
-0.831409335,
-0.462193251,
-0.726779819
1094mean(:,
0)
= getMean(sample, dim)
1098 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1100call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1103call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1112 lb(isam)
= ub(isam
- 1)
+ 1
1120call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1124+0.191084504,
-0.496219158,
+0.360538960,
+0.476084352,
-0.906051159
1125-0.944288492,
+0.436568618,
+0.694154501,
+0.498358011E-1,
+0.589161873
1126mean(:,
0)
= getMean(sample, dim)
1128-0.749125034E-1,
+0.165086463
1130 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1132call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1134-0.749125034E-1,
+0.165086463
1135call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1137-0.749125034E-1,
+0.165086463
1139-0.749125034E-1,
+0.165086463
1144 lb(isam)
= ub(isam
- 1)
+ 1
1152call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1156+0.835810900E-1,
-0.682339191,
-0.553030729,
+0.464090824,
+0.968884826,
+0.670804024
1157mean(:,
0)
= getMean(sample, dim)
1161 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1163call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1166call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1175 lb(isam)
= ub(isam
- 1)
+ 1
1183call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1187+0.367890120,
+0.609076381,
+0.401162028,
-0.771420240,
+0.933075070,
-0.700249553,
-0.906868339,
-0.815922618,
+0.802740455,
-0.490402818,
-0.417932272
1188+0.738120317,
+0.189231277,
-0.340724111,
-0.359696507,
+0.565462232,
+0.604000092,
+0.469835281,
-0.866739750E-1,
-0.849789977,
-0.405657291E-2,
-0.722428322
1189-0.663097858,
-0.258183479,
+0.294872403,
+0.384969354,
+0.363638520,
-0.845884442,
-0.388226986,
-0.509362459,
-0.681086779E-1,
+0.399967074,
-0.389128923
1190mean(:,
0)
= getMean(sample, dim)
1192-0.898956209E-1,
+0.184799768E-1,
-0.152595043
1194 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1196call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1198-0.898955911E-1,
+0.184799843E-1,
-0.152595028
1199call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1201-0.898955911E-1,
+0.184799843E-1,
-0.152595028
1203-0.898956209E-1,
+0.184799768E-1,
-0.152595043
1208 lb(isam)
= ub(isam
- 1)
+ 1
1216call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1220-0.908590198,
-0.302977085,
-0.435143232
1221mean(:,
0)
= getMean(sample, dim)
1225 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1227call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1230call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1244 lb(isam)
= ub(isam
- 1)
+ 1
1252call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1256-0.384712934,
+0.549722552,
-0.137264729,
+0.115150094,
-0.100257397E-1,
+0.840308785,
-0.550587416,
-0.964916825,
+0.822875261,
-0.500732899,
-0.533808827,
+0.936962247
1259+5,
+9,
+2,
+10,
+5,
+9,
+5,
+3,
+4,
+2,
+6,
+7
1260mean(:,
0)
= getMean(sample,
2_IK, iweight)
1264 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1266call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1269call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1278 lb(isam)
= ub(isam
- 1)
+ 1
1286call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1290+0.279239058,
+0.994216681,
-0.575693250,
+0.922137499E-1,
-0.759993434
1291-0.427760720,
-0.801674724,
+0.427672267,
+0.767664909E-1,
-0.371416569
1295mean(:,
0)
= getMean(sample,
2_IK, iweight)
1297+0.101925634,
-0.226913929
1299 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1301call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1303+0.101925634,
-0.226913929
1304call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1306+0.101925634,
-0.226913929
1308+0.101925634,
-0.226913929
1313 lb(isam)
= ub(isam
- 1)
+ 1
1321call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1325+0.730616331,
-0.132999539,
+0.476830125,
-0.251829743,
+0.729217172,
+0.875183105,
+0.172075510
1326-0.682505846,
+0.933519244,
-0.287114024,
-0.702607870,
+0.840710402,
-0.154170752,
+0.926141858
1327+0.404745817,
+0.138210416,
+0.761262894,
-0.433515072,
+0.669943213,
+0.527891040,
+0.612375617
1330+9,
+8,
+2,
+2,
+4,
+2,
+3
1331mean(:,
0)
= getMean(sample,
2_IK, iweight)
1333+0.371500462,
+0.172636136,
+0.365885794
1335 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1337call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1339+0.371500432,
+0.172636122,
+0.365885794
1340call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1342+0.371500432,
+0.172636122,
+0.365885794
1344+0.371500462,
+0.172636136,
+0.365885794
1349 lb(isam)
= ub(isam
- 1)
+ 1
1357call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1361-0.645502090,
+0.669400930,
+0.343166232,
-0.371553659,
+0.558206797,
+0.305296898,
+0.789406538,
+0.235946059
1364+6,
+1,
+9,
+3,
+4,
+6,
+4,
+6
1365mean(:,
0)
= getMean(sample,
2_IK, iweight)
1369 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1371call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1374call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1383 lb(isam)
= ub(isam
- 1)
+ 1
1391call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1395+0.595678687,
-0.202252388,
+0.264178634,
+0.703871250,
-0.252048373
1399mean(:,
0)
= getMean(sample,
2_IK, iweight)
1403 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1405call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1408call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1417 lb(isam)
= ub(isam
- 1)
+ 1
1425call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1429-0.783409595,
-0.403064728,
+0.469662547,
-0.770673156,
+0.333247185,
+0.380345225,
+0.717551708
1432+7,
+3,
+4,
+1,
+10,
+10,
+6
1433mean(:,
0)
= getMean(sample,
2_IK, iweight)
1437 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1439call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1442call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1451 lb(isam)
= ub(isam
- 1)
+ 1
1459call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1463+0.168335438E-1,
+0.140553713,
-0.199712873,
+0.934781551,
-0.214332819,
+0.518136501,
+0.138473868
1464-0.638781548,
-0.714817882,
+0.895185471E-1,
-0.236177444E-1,
-0.561852098,
+0.222158670,
+0.131550312
1467+5,
+3,
+1,
+1,
+1,
+6,
+10
1468mean(:,
0)
= getMean(sample,
2_IK, iweight)
1470+0.204448983,
-0.117994718
1472 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1474call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1476+0.204448998,
-0.117994718
1477call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1479+0.204448998,
-0.117994718
1481+0.204448983,
-0.117994718
1486 lb(isam)
= ub(isam
- 1)
+ 1
1494call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1498+0.668281317E-1,
+0.803951263,
-0.436257243,
-0.256437182,
-0.922491908,
+0.282671809,
-0.533916473,
-0.124479532,
+0.271448731,
+0.451467514
1499+0.926629543,
-0.446412802,
-0.170309901,
-0.568793416,
-0.936034203,
+0.437713504,
-0.681233048,
-0.723099232,
+0.909014106,
-0.826038122E-1
1502+10,
+2,
+8,
+2,
+10,
+1,
+3,
+4,
+1,
+5
1503mean(:,
0)
= getMean(sample,
2_IK, iweight)
1505-0.222606018,
-0.162811413
1507 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1509call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1511-0.222606018,
-0.162811443
1512call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1514-0.222606018,
-0.162811443
1516-0.222606018,
-0.162811413
1521 lb(isam)
= ub(isam
- 1)
+ 1
1529call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1533+0.280038595,
+0.882654428,
+0.306001306,
+0.278992653E-1,
-0.862076402,
-0.309554338E-1,
-0.718109369,
-0.290211439,
-0.396250844,
-0.553926229,
+0.684447765,
-0.910490274
1536+5,
+8,
+9,
+8,
+8,
+5,
+7,
+7,
+6,
+7,
+3,
+6
1537mean(:,
0)
= getMean(sample,
2_IK, iweight)
1541 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1543call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1546call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1555 lb(isam)
= ub(isam
- 1)
+ 1
1563call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1567-0.527770281,
-0.955617428E-2,
+0.224102855,
+0.468953490,
+0.914654732,
-0.933868885E-1,
-0.558397412
1568-0.498374701,
+0.854009986,
+0.767431259E-1,
+0.936418176,
+0.796172738,
+0.176037550E-1,
+0.467400074
1571+6,
+9,
+3,
+10,
+7,
+5,
+8
1572mean(:,
0)
= getMean(sample,
2_IK, iweight)
1574+0.745351166E-1,
+0.493555874
1576 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1578call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1580+0.745351240E-1,
+0.493555844
1581call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1583+0.745351240E-1,
+0.493555844
1585+0.745351166E-1,
+0.493555874
1595 lb(isam)
= ub(isam
- 1)
+ 1
1603call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1607+0.800471544,
+0.140030146,
+0.504118919,
+0.608778596,
+0.364939332,
+0.792155981,
-0.663981557,
+0.506578684,
+0.644219875,
+0.904301763,
+0.311916947,
-0.754499316
1608-0.660890818,
+0.115471721,
+0.656057358,
-0.441813707,
-0.218908787,
-0.567436576,
+0.327087879,
+0.777155280,
+0.534326434,
+0.395199895,
+0.894017220,
-0.138114810
1609-0.742075562,
-0.189094543,
+0.891225100,
+0.959707975,
-0.786428452,
+0.627806187E-1,
+0.923795581,
-0.439095497E-1,
+0.679832339,
+0.613041639,
-0.507868528E-1,
-0.554790616
1610-0.423260093,
+0.428040028E-1,
-0.466392398,
-0.171620369,
+0.509369850,
+0.376764536,
-0.272278428,
-0.643194556,
-0.311586261,
+0.736231804E-1,
-0.832336307,
+0.790560603
1611+0.870278120,
+0.606098533,
+0.980907917,
+0.840212822,
-0.985554576,
-0.686565638,
-0.172982216,
+0.364545584E-1,
-0.954738140,
+0.489278674,
+0.535527706,
-0.369416118
1614+4.00000000,
+9.00000000,
+10.0000000,
+6.00000000,
+2.00000000,
+10.0000000,
+8.00000000,
+5.00000000,
+5.00000000,
+9.00000000,
+9.00000000,
+8.00000000
1615mean(:,
0)
= getMean(sample, dim, rweight)
1617+0.319524378,
+0.186651036,
+0.238202900,
-0.113778017,
+0.179320514
1619 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1621call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1623+0.319524378,
+0.186651006,
+0.238202810,
-0.113778003,
+0.179320499
1624call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1626+0.319524378,
+0.186651006,
+0.238202810,
-0.113778003,
+0.179320499
1628+0.319524378,
+0.186651036,
+0.238202900,
-0.113778017,
+0.179320514
1633 lb(isam)
= ub(isam
- 1)
+ 1
1641call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1645-0.503708124E-1,
+0.597526789,
+0.169150472,
+0.614407063,
-0.893617749
1648+4.00000000,
+8.00000000,
+9.00000000,
+1.00000000,
+5.00000000
1649mean(:,
0)
= getMean(sample, dim, rweight)
1653 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1655call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1658call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1667 lb(isam)
= ub(isam
- 1)
+ 1
1675call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1679+0.190246105,
-0.607671618,
+0.561362743,
-0.114269495,
+0.839279413,
+0.430876255,
+0.992139101,
-0.282889247,
+0.723344088E-1,
-0.988511562,
-0.454624891E-1,
+0.149570584,
+0.339024067E-1
1680+0.615563273,
-0.338718891,
-0.645978808,
-0.167361379,
+0.767328382,
+0.139996886,
+0.739565730,
+0.783224463,
-0.904332995,
+0.784183979,
+0.820348024,
-0.710049033,
+0.687514424
1681+0.204392195,
+0.935039520,
+0.618576288,
+0.794567704,
+0.462361574,
-0.517475605E-1,
+0.537747145E-1,
-0.694197774,
+0.891687989,
+0.696514249,
-0.946390033,
-0.865288258,
+0.324153781
1682+0.815904737,
+0.375840306,
-0.755804777E-1,
-0.968709588,
+0.730776310,
+0.837359190,
-0.927855849,
+0.693623185,
+0.728870034,
-0.109955907,
-0.629399538,
+0.817830086,
+0.666112304
1685+7.00000000,
+2.00000000,
+6.00000000,
+1.00000000,
+10.0000000,
+2.00000000,
+8.00000000,
+3.00000000,
+8.00000000,
+6.00000000,
+8.00000000,
+4.00000000,
+6.00000000
1686mean(:,
0)
= getMean(sample, dim, rweight)
1688+0.208434954,
+0.297547340,
+0.181687504,
+0.226554215
1690 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1692call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1694+0.208434939,
+0.297547340,
+0.181687504,
+0.226554215
1695call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1697+0.208434939,
+0.297547340,
+0.181687504,
+0.226554215
1699+0.208434954,
+0.297547340,
+0.181687504,
+0.226554215
1704 lb(isam)
= ub(isam
- 1)
+ 1
1712call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1716-0.921881676,
-0.398877025,
+0.370139599,
+0.702141523,
+0.120090127,
-0.752707124,
+0.725308895,
+0.978895426
1719+10.0000000,
+1.00000000,
+2.00000000,
+4.00000000,
+9.00000000,
+8.00000000,
+10.0000000,
+9.00000000
1720mean(:,
0)
= getMean(sample, dim, rweight)
1724 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1726call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1729call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1738 lb(isam)
= ub(isam
- 1)
+ 1
1746call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1750+0.234567642,
+0.312389135E-1,
-0.326205850,
-0.980954170,
-0.703580379,
-0.539213657,
-0.102962613,
-0.677458525
1753+5.00000000,
+4.00000000,
+1.00000000,
+10.0000000,
+5.00000000,
+9.00000000,
+5.00000000,
+6.00000000
1754mean(:,
0)
= getMean(sample, dim, rweight)
1758 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1760call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1763call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1772 lb(isam)
= ub(isam
- 1)
+ 1
1780call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1784-0.753197193,
+0.172114730,
+0.854082704,
+0.575641036,
+0.174626112,
-0.459991217,
-0.211122274,
+0.437294126,
-0.767268538,
+0.496149063E-2,
-0.348613262,
-0.847686291,
+0.812623978
1785+0.411628127,
+0.276562333,
+0.935594201,
+0.346701145E-1,
+0.276611090,
+0.495509148,
-0.803940415,
+0.611291170,
-0.729983568,
+0.119412065,
+0.936248541,
-0.183007598,
+0.776155591
1786-0.329596519,
+0.261848211,
+0.929041743,
-0.492847443,
+0.461137056,
+0.728097796,
+0.905899763,
-0.339400411,
-0.184212685,
-0.268014669,
+0.298787355E-1,
-0.435762405,
+0.909954667
1787+0.905327320,
-0.163976312,
-0.164536834,
+0.650094748,
+0.144749880E-1,
-0.764482379,
-0.767233610,
+0.967568874,
+0.788559556,
+0.142902136E-1,
+0.885273576,
-0.643710375,
-0.503027439
1788+0.846656203,
-0.382355452E-1,
-0.571891546,
+0.134841919,
-0.131641984,
-0.699191689,
+0.105100155,
+0.783560276E-1,
+0.401389241,
+0.941381097,
-0.987400174,
+0.383214235,
-0.843565702
1789+0.841886044,
+0.233924389E-1,
+0.848460913,
-0.457740545,
-0.128737092,
+0.648412108,
-0.692012429,
-0.231123924,
+0.350592852,
-0.860209465,
+0.380976915,
+0.340142846,
-0.316105008
1792+2.00000000,
+10.0000000,
+6.00000000,
+2.00000000,
+6.00000000,
+6.00000000,
+4.00000000,
+3.00000000,
+9.00000000,
+7.00000000,
+10.0000000,
+3.00000000,
+6.00000000
1793mean(:,
0)
= getMean(sample, dim, rweight)
1795-0.375459380E-1,
+0.273770928,
+0.232511654,
+0.934877768E-1,
-0.131943583,
+0.786087438E-1
1797 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1799call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1801-0.375459380E-1,
+0.273770928,
+0.232511654,
+0.934877843E-1,
-0.131943583,
+0.786087438E-1
1802call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1804-0.375459380E-1,
+0.273770928,
+0.232511654,
+0.934877843E-1,
-0.131943583,
+0.786087438E-1
1806-0.375459380E-1,
+0.273770928,
+0.232511654,
+0.934877768E-1,
-0.131943583,
+0.786087438E-1
1811 lb(isam)
= ub(isam
- 1)
+ 1
1819call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1823+0.672057867E-1,
-0.878316402,
-0.956241608,
+0.925262809,
+0.848703861,
-0.332255125,
+0.144010544,
-0.515627861E-2,
-0.229692221
1826+10.0000000,
+6.00000000,
+9.00000000,
+9.00000000,
+7.00000000,
+1.00000000,
+7.00000000,
+1.00000000,
+3.00000000
1827mean(:,
0)
= getMean(sample, dim, rweight)
1831 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1833call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1836call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1845 lb(isam)
= ub(isam
- 1)
+ 1
1853call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1857+0.133295655,
-0.931665063,
-0.268453121,
+0.452363372,
-0.309066534,
-0.750119686,
-0.393433928
1860+10.0000000,
+7.00000000,
+4.00000000,
+1.00000000,
+3.00000000,
+6.00000000,
+2.00000000
1861mean(:,
0)
= getMean(sample, dim, rweight)
1865 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1867call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1870call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1879 lb(isam)
= ub(isam
- 1)
+ 1
1887call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1891+0.687152743,
-0.565898180,
-0.154854059
1894+10.0000000,
+1.00000000,
+5.00000000
1895mean(:,
0)
= getMean(sample, dim, rweight)
1899 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1901call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1904call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1913 lb(isam)
= ub(isam
- 1)
+ 1
1921call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1925+0.726318836,
-0.640252709,
+0.172882915,
-0.439074516,
-0.703207850,
-0.568240166
1928+5.00000000,
+10.0000000,
+6.00000000,
+8.00000000,
+5.00000000,
+7.00000000
1929mean(:,
0)
= getMean(sample, dim, rweight)
1933 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1935call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))
1938call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(rweight(:ub(
1)))
/ sum(rweight), TKG))