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.450917125,
+0.647481143,
+0.579088211,
+0.919775665,
+0.116213977,
+0.626923084,
+0.252830148,
+0.812254667,
+0.764650822,
+0.402873576,
+0.611293614,
+0.274069428,
+0.698515773
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.597962201,
+0.783627629,
+0.778012991,
+0.839513838,
+0.665288568E-1,
+0.574622750E-1,
+0.287638009
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.522541583,
+0.871412575,
+0.480284035,
+0.971569538,
+0.653817952,
+0.366228819E-1,
+0.211793125,
+0.485796511
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.701328456,
+0.717610896,
+0.981206298,
+0.235992670,
+0.513632417,
+0.684459209,
+0.128934205,
+0.514072537,
+0.752650261,
+0.330836773
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.641783416,
+0.943830848,
+0.969381273,
+0.338236094E-1,
+0.402671814,
+0.632590532,
+0.337394714,
+0.189478397E-1
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.926578522,
+0.492532432,
+0.584347427,
+0.594964683,
+0.367180228,
+0.699692965E-2,
+0.717345178,
+0.390479803,
+0.529461861
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.338360667,
+0.140468717,
+0.724098504
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.310075819,
+0.423070192E-1,
+0.579325616,
+0.569845974,
+0.194519162,
+0.624657273,
+0.625143647E-1,
+0.435734630,
+0.460889161,
+0.947877705,
+0.317531228E-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.701725006,
+0.641202450,
+0.524549246,
+0.442632616,
+0.168878436,
+0.173875213,
+0.876675248,
+0.943157434,
+0.165064931E-1
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.526680708,
+0.296144009,
+0.329744816E-2,
+0.148903370,
+0.488615394
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.762919188E-1,
+0.999307156,
+0.328115582,
+0.229218841,
+0.201227903,
+0.557934225,
+0.675776422,
+0.897380888,
+0.985947430
306+8,
+1,
+5,
+3,
+9,
+3,
+3,
+8,
+9
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.222939193,
+0.285459101,
+0.305185556,
+0.301230967,
+0.495533168,
+0.634538293,
+0.813965201,
+0.473072708
337+1,
+7,
+2,
+2,
+10,
+3,
+8,
+6
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.148988485,
+0.875638485,
+0.153747797E-1,
+0.803034306,
+0.470391035,
+0.685037196,
+0.842763603
368+9,
+1,
+2,
+2,
+6,
+1,
+9
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.826362610,
+0.792909563,
+0.241511822,
+0.954819739,
+0.388998926,
+0.910269499,
+0.960105479
399+5,
+1,
+1,
+6,
+2,
+3,
+7
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.280467212,
+0.837615788,
+0.817623377
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.234112203,
+0.263211846,
+0.483875990,
+0.937600076,
+0.293082535,
+0.730139375,
+0.170546353
461+2,
+3,
+4,
+9,
+10,
+6,
+7
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.372631550,
+0.318311393,
+0.158778369,
+0.263324797,
+0.575749218,
+0.199925721,
+0.569555163,
+0.341664970,
+0.218133628,
+0.784770250E-1,
+0.595189154,
+0.561443925,
+0.283682048,
+0.743308663E-1
492+7,
+9,
+10,
+7,
+10,
+3,
+1,
+6,
+8,
+5,
+1,
+5,
+10,
+6
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.964238048,
+0.451485991,
+0.515642166,
+0.842446625,
+0.720165610,
+0.527779520,
+0.763204515,
+0.412703097
523+5,
+1,
+3,
+1,
+1,
+10,
+8,
+3
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.818193436,
+0.719758928,
+0.379511297,
+0.931423962,
+0.775915027,
+0.996842086,
+0.627090812,
+0.884779632,
+0.324178994,
+0.762720704
554+8,
+10,
+10,
+1,
+5,
+8,
+6,
+10,
+10,
+4
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.152686119,
+0.208874285,
+0.506744385E-1,
+0.617766142,
+0.680997014,
+0.365606904,
+0.453252792E-1,
+0.304934323,
+0.695051908,
+0.777131855,
+0.759560645
585+3,
+1,
+2,
+10,
+8,
+8,
+9,
+2,
+10,
+1,
+4
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.219523549,
+0.530630350E-1,
+0.591719151,
+0.554355145,
+0.855650306E-1,
+0.390973747,
+0.686549723,
+0.876149595
621+1.98177838,
+1.93096590,
+1.66918647,
+1.06585813,
+1.89013302,
+1.56982970,
+1.27188182,
+1.20056570
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.198383570,
+0.804668903,
+0.708835721E-1,
+0.592108309,
+0.514955163,
+0.497140467,
+0.950913489,
+0.890209377,
+0.194909513
652+1.54535902,
+1.55188823,
+1.93018126,
+1.36929965,
+1.75334907,
+1.84446430,
+1.91420722,
+1.23851562,
+1.71659327
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.336120367,
+0.197943091,
+0.902010739,
+0.243265152,
+0.763902426,
+0.316251516
683+1.50913501,
+1.32110560,
+1.80416369,
+1.94302225,
+1.03593516,
+1.43990326
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.202634156,
+0.247774780,
+0.121375084,
+0.754495919,
+0.339238048E-1,
+0.470085084,
+0.981331944
714+1.59052038,
+1.59817326,
+1.90501332,
+1.15679741,
+1.14449167,
+1.57645488,
+1.53720415
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.790228426,
+0.782313526,
+0.158574581E-1,
+0.461433649,
+0.472879410E-2,
+0.663071275
745+1.52090764,
+1.49448633,
+1.40140331,
+1.67496121,
+1.09953547,
+1.62354589
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.134169936,
+0.822881997,
+0.938389301,
+0.721104383,
+0.373824298,
+0.668917716,
+0.153163910
776+1.52584398,
+1.71873546,
+1.60745549,
+1.77679586,
+1.63114715,
+1.35978794,
+1.74717808
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.212155700,
+0.303383946,
+0.101725161,
+0.863918900,
+0.623905778,
+0.472624063,
+0.570023000,
+0.147542298
807+1.75116062,
+1.94414186,
+1.92987323,
+1.17084336,
+1.24150038,
+1.29195631,
+1.89036167,
+1.71409440
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.925814509,
+0.818524837,
+0.663796067,
+0.374792159,
+0.661942065,
+0.123193860E-1,
+0.794063985,
+0.938345492,
+0.216213584,
+0.858668685
838+1.06552732,
+1.54301429,
+1.46615005,
+1.15717053,
+1.21852922,
+1.94094324,
+1.47867775,
+1.58476901,
+1.02172172,
+1.46851611
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.302460909,
+0.918179691,
+0.773340881,
+0.581627488E-1,
+0.236778378,
+0.314706564,
+0.416233838,
+0.892307341
869+1.12129414,
+1.42488420,
+1.04970360,
+1.67097211,
+1.34633660,
+1.71282721,
+1.61148119,
+1.92837989
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.282088637,
+0.648042798,
+0.295889616,
+0.151668191,
+0.713053942,
+0.620828569,
+0.523952007,
+0.509549379,
+0.329731762
900+1.43753850,
+1.78597784,
+1.91361976,
+1.66635203,
+1.79687297,
+1.18335116,
+1.43050218,
+1.29144895,
+1.57350898
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.792492628E-1,
-0.814815283,
-0.527313471,
-0.769526958,
-0.952939153,
+0.106447458,
-0.987016201
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.981962085,
-0.651926517,
-0.735267639
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.863795877,
-0.360506058,
+0.924755096,
+0.939900517,
+0.948748112,
+0.819427133
999-0.364199758,
+0.722237945,
-0.718837142,
+0.304505348,
-0.886546612,
+0.672012687
1000mean(:,
0)
= getMean(sample, dim)
1002+0.689353466,
-0.451379232E-1
1004 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1006call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1008+0.689353466,
-0.451379195E-1
1009call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1011+0.689353466,
-0.451379195E-1
1013+0.689353466,
-0.451379232E-1
1018 lb(isam)
= ub(isam
- 1)
+ 1
1026call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1030-0.830671430,
-0.171727538,
-0.417463899,
-0.731832266,
+0.671187639,
-0.894712806,
+0.508095384
1031mean(:,
0)
= getMean(sample, dim)
1035 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1037call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1040call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1049 lb(isam)
= ub(isam
- 1)
+ 1
1057call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1061-0.650840878,
-0.889423847,
+0.224046469,
-0.217802882,
-0.441380978,
-0.530807972
1062mean(:,
0)
= getMean(sample, dim)
1066 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1068call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1071call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1080 lb(isam)
= ub(isam
- 1)
+ 1
1088call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1092-0.348052502,
+0.466094136,
-0.484915733,
-0.759090066,
-0.997612596
1093+0.772763848,
-0.393468142E-1,
+0.243729353,
-0.548236370E-1,
-0.484283328
1094mean(:,
0)
= getMean(sample, dim)
1096-0.424715340,
+0.876078829E-1
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))
1102-0.424715370,
+0.876078829E-1
1103call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1105-0.424715370,
+0.876078829E-1
1107-0.424715340,
+0.876078829E-1
1112 lb(isam)
= ub(isam
- 1)
+ 1
1120call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1124+0.396198750,
-0.762101054,
+0.969689369,
+0.204674721,
+0.795172453,
-0.854734540,
+0.329230547,
-0.234186053,
-0.215780854,
+0.314898849
1125-0.185340166,
-0.274718761,
-0.306277514,
+0.301264405,
-0.267217159E-1,
+0.811058521,
+0.658574343,
-0.489765525,
-0.956857443,
-0.833419800
1126-0.399742723,
+0.608554840,
-0.359128714E-1,
+0.669092059,
-0.329301000,
+0.851661325,
-0.249085426,
+0.185821414,
+0.263315082,
+0.261289597
1127-0.332344055,
-0.396055222,
+0.368460655,
-0.634359360,
-0.381998777,
+0.901958227,
+0.560599923,
-0.559601068,
-0.351397157,
+0.666177034
1128-0.135547638,
+0.815477252,
+0.891591311E-1,
+0.271032333,
-0.780617595,
+0.978300214,
-0.767898560E-2,
-0.629989862,
+0.253776312,
-0.507682443
1129mean(:,
0)
= getMean(sample, dim)
1131+0.943062231E-1,
-0.130220369,
+0.182569236,
-0.158559810E-1,
+0.346228741E-1
1133 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1135call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1137+0.943062156E-1,
-0.130220369,
+0.182569236,
-0.158559754E-1,
+0.346228704E-1
1138call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1140+0.943062156E-1,
-0.130220369,
+0.182569236,
-0.158559754E-1,
+0.346228704E-1
1142+0.943062231E-1,
-0.130220369,
+0.182569236,
-0.158559810E-1,
+0.346228741E-1
1147 lb(isam)
= ub(isam
- 1)
+ 1
1155call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1159-0.438054562,
+0.977368116,
+0.262950063,
-0.437987804,
+0.531090498E-1,
+0.335477114,
-0.878040314,
-0.496594548
1160-0.887449384,
-0.255854487,
+0.749666452,
-0.704798341,
+0.632519722E-1,
+0.500726104,
+0.627201676,
-0.743068576
1161mean(:,
0)
= getMean(sample, dim)
1163-0.777216107E-1,
-0.812905729E-1
1165 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1167call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1169-0.777216107E-1,
-0.812905729E-1
1170call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1172-0.777216107E-1,
-0.812905729E-1
1174-0.777216107E-1,
-0.812905729E-1
1179 lb(isam)
= ub(isam
- 1)
+ 1
1187call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1191+0.279461145,
+0.356027722,
+0.428345919,
-0.444567561,
-0.917892456,
-0.291271448,
-0.481053829,
-0.439249873,
-0.753668666,
+0.574824929,
+0.199058056
1192-0.777057290,
-0.511084795,
+0.236571193,
-0.738955975,
+0.980914474,
+0.650825620,
-0.168111682,
+0.216277838,
-0.635682106,
+0.419185996,
+0.214123487
1193+0.944226980E-1,
-0.260097980,
-0.148499846,
-0.759849072,
-0.103393078,
-0.822131991,
+0.996027231,
+0.794895291,
+0.807016015,
-0.417787433,
-0.657907724E-1
1194-0.917671561,
-0.996422768E-2,
-0.553924561,
+0.124514103E-2,
+0.434498429,
+0.597463131,
+0.996397018,
+0.294252753,
-0.905705690,
+0.755450368,
+0.667023301
1195mean(:,
0)
= getMean(sample, dim)
1197-0.135453284,
-0.102721127E-1,
+0.104373693E-1,
+0.123551287
1199 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1201call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1203-0.135453284,
-0.102721229E-1,
+0.104373619E-1,
+0.123551294
1204call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1206-0.135453284,
-0.102721229E-1,
+0.104373619E-1,
+0.123551294
1208-0.135453284,
-0.102721127E-1,
+0.104373693E-1,
+0.123551287
1213 lb(isam)
= ub(isam
- 1)
+ 1
1221call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1225-0.684290051,
+0.782229781,
+0.281648278,
-0.645682335,
-0.726508856,
+0.632212996,
+0.282883763
1226mean(:,
0)
= getMean(sample, dim)
1230 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
1232call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1235call setMeanMerged(mean(:,
2), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
1249 lb(isam)
= ub(isam
- 1)
+ 1
1257call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1261+0.482810855,
+0.937024593
1265mean(:,
0)
= getMean(sample,
2_IK, iweight)
1269 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1271call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1274call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1283 lb(isam)
= ub(isam
- 1)
+ 1
1291call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1295-0.576741457,
-0.139735341,
-0.573289156,
+0.326125622E-1,
-0.897829294,
+0.283338547,
-0.511916757
1298+4,
+6,
+1,
+2,
+6,
+7,
+10
1299mean(:,
0)
= getMean(sample,
2_IK, iweight)
1303 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1305call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1308call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1317 lb(isam)
= ub(isam
- 1)
+ 1
1325call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1329-0.841204166,
+0.985905647,
-0.919012427,
+0.287443757
1333mean(:,
0)
= getMean(sample,
2_IK, iweight)
1337 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1339call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1342call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1351 lb(isam)
= ub(isam
- 1)
+ 1
1359call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1363-0.644922614,
+0.565256953,
-0.828196526,
+0.732427835,
+0.527273297,
-0.257063031
1366+2,
+10,
+8,
+10,
+3,
+1
1367mean(:,
0)
= getMean(sample,
2_IK, iweight)
1371 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1373call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1376call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1385 lb(isam)
= ub(isam
- 1)
+ 1
1393call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1397+0.484885097,
+0.110787153E-1,
-0.757265687,
+0.447370410,
-0.790321827,
-0.423355579,
+0.737119436,
-0.442779064E-2,
-0.373531342,
+0.897284985,
+0.957667589
1400+6,
+10,
+6,
+1,
+2,
+10,
+10,
+9,
+8,
+1,
+7
1401mean(:,
0)
= getMean(sample,
2_IK, iweight)
1405 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1407call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1410call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1419 lb(isam)
= ub(isam
- 1)
+ 1
1427call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1431-0.331791997,
+0.163888931E-2,
+0.600777745,
-0.913177252,
+0.160655975E-1,
-0.749979854
1434+2,
+1,
+8,
+7,
+3,
+1
1435mean(:,
0)
= getMean(sample,
2_IK, iweight)
1439 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1441call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1444call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1453 lb(isam)
= ub(isam
- 1)
+ 1
1461call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1465-0.487235069,
+0.339401484,
-0.905389786,
-0.504239202,
+0.351963758,
+0.431632638,
+0.205532193,
-0.898764491,
-0.145130992,
+0.679128885,
+0.553917170
1466+0.788933158,
-0.761665106E-1,
+0.873067021,
+0.236579061,
-0.158794761,
-0.958598852E-1,
+0.153317928,
+0.974743128,
-0.469558358,
+0.768145561,
-0.303305864
1467-0.758252621,
-0.166301608,
+0.885466933,
-0.781515956,
+0.444861531,
-0.961319327,
+0.109414339,
-0.751347542E-1,
+0.428465486,
-0.640737653,
-0.348250866E-1
1468+0.340287924,
+0.560660362,
+0.126956105,
+0.890405178,
-0.931557894,
+0.238565683,
+0.883358121,
-0.776199460,
-0.485056162,
-0.979454756,
-0.906894207
1471+7,
+7,
+1,
+5,
+10,
+3,
+6,
+5,
+6,
+1,
+3
1472mean(:,
0)
= getMean(sample,
2_IK, iweight)
1474-0.266188998E-1,
+0.148227811,
-0.107832253,
-0.538071059E-1
1476 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1478call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1480-0.266189016E-1,
+0.148227811,
-0.107832253,
-0.538071245E-1
1481call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1483-0.266189016E-1,
+0.148227811,
-0.107832253,
-0.538071245E-1
1485-0.266188998E-1,
+0.148227811,
-0.107832253,
-0.538071059E-1
1490 lb(isam)
= ub(isam
- 1)
+ 1
1498call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1502-0.557179093,
-0.528724432,
-0.443182349,
-0.707323909,
+0.916378140,
+0.123061895,
-0.873217344,
-0.624743938,
+0.630951285,
+0.595556021,
+0.743969679E-1,
-0.967899799,
+0.823269606,
+0.250709176
1503+0.674735785,
+0.910176396,
+0.586703658,
-0.718786836,
+0.308602452,
+0.693881512,
-0.855715632,
+0.352848887,
-0.156742334,
+0.991984129,
+0.867736578,
-0.141610742,
+0.650489926,
+0.775660276
1504+0.300323248,
+0.238235950,
+0.433225513,
-0.355881214,
+0.819830060,
+0.811681986,
+0.479061723,
-0.901137233,
+0.111945748,
+0.397492528,
-0.279687285,
+0.882044435,
-0.918330908,
+0.407699585
1505-0.393614769,
+0.168538332,
-0.803650618E-1,
+0.872630239,
-0.610213757,
+0.281566262,
+0.812224746,
+0.147985220,
-0.317014456,
-0.440087318,
+0.190901160,
-0.300594568E-1,
-0.749153733,
+0.874656796
1508+2,
+2,
+1,
+2,
+8,
+7,
+10,
+9,
+3,
+1,
+2,
+6,
+4,
+1
1509mean(:,
0)
= getMean(sample,
2_IK, iweight)
1511-0.169283718,
+0.156010613,
+0.205516815,
+0.766632408E-1
1513 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1515call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1517-0.169283733,
+0.156010598,
+0.205516815,
+0.766632408E-1
1518call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1520-0.169283733,
+0.156010598,
+0.205516815,
+0.766632408E-1
1522-0.169283718,
+0.156010613,
+0.205516815,
+0.766632408E-1
1527 lb(isam)
= ub(isam
- 1)
+ 1
1535call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1539+0.876397729,
+0.836336255,
+0.976139665,
-0.693070889E-2,
-0.476712465,
-0.982724071
1540-0.257339001,
+0.708578825E-1,
+0.947160125,
-0.254223228,
-0.348530769,
-0.551466703
1543+7,
+6,
+4,
+8,
+5,
+1
1544mean(:,
0)
= getMean(sample,
2_IK, iweight)
1546+0.375342816,
-0.617900491E-1
1548 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1550call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1552+0.375342846,
-0.617900528E-1
1553call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1555+0.375342846,
-0.617900528E-1
1557+0.375342816,
-0.617900491E-1
1562 lb(isam)
= ub(isam
- 1)
+ 1
1570call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1574-0.646990180,
-0.668978691E-1,
+0.583444357,
+0.650681853,
-0.343840957
1578mean(:,
0)
= getMean(sample,
2_IK, iweight)
1582 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
1584call setMeanMerged(meanMerged, mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1587call setMeanMerged(mean(:,
2), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
1601 lb(isam)
= ub(isam
- 1)
+ 1
1609call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1613+0.490081310E-2,
-0.618352771,
-0.959584713E-1,
+0.265206099E-1,
+0.385996222,
-0.763190150,
+0.227920651,
+0.899359584,
+0.703910232,
-0.857494235,
-0.138871908
1614-0.589534402,
+0.805036664,
+0.458170056,
-0.920660496E-1,
-0.412250161,
+0.467877626,
-0.306022406,
+0.731227398E-1,
+0.243213415,
-0.406769276,
+0.761890411
1615-0.390493512,
-0.837315321,
+0.380594492,
-0.259415269,
+0.214852452,
-0.805881739,
+0.429843307,
-0.336398005,
+0.998563766E-1,
+0.351866484E-1,
+0.867916346
1616-0.613202572,
+0.735453367E-1,
+0.917518139E-1,
-0.957813740,
-0.390845656,
+0.644869804,
+0.616778135,
-0.392547965,
+0.766606331,
+0.914812326,
+0.414704919
1617+0.212546587E-1,
-0.720114708E-1,
+0.494683504,
+0.113169074,
+0.904285908E-1,
-0.741189837,
-0.825785637,
+0.310242414,
-0.897150755,
+0.789320827,
+0.988754988