17 integer(IK) ,
parameter :: nsam
= 2
18 integer(IK) ,
allocatable :: iweight(:)
19 integer(IK) :: isam, ndim, lb(nsam), ub(nsam)
20 integer(IK) :: dim, itry, ntry
= 10
21 type(display_type) :: disp
25 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
26 call disp%show(
"!Compute the biased merged covariance of a multivariate sample.")
27 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
31 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), meanMerged(:), covMerged(:,:)
32 real(TKG),
allocatable :: sample(:,:)
35 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
38 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
39 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
42 lb(isam)
= ub(isam
- 1)
+ 1
49 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
51 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
52 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
53 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
54 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
55 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
57 call disp%show(
"call setResized(meanMerged, ndim)")
59 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
63 call disp%show(
"cov(:,:,0) = getCov(sample, dim)")
64 cov(:,:,
0)
= getCov(sample, dim)
65 call disp%show(
"cov(:,:,0) ! reference")
67 call disp%show(
"mean(:,0) = getMean(sample, dim)")
68 mean(:,
0)
= getMean(sample, dim)
69 call disp%show(
"mean(:,0) ! reference")
72 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim)")
73 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
76 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
77 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
79 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG), uppDia)")
80 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
83 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG), lowDia)")
84 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
87 call disp%show(
"call setCovMeanMerged(cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG), uppDia)")
88 call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
91 call disp%show(
"cov(:,:,0) ! reference")
98 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
99 call disp%show(
"!Compute the biased merged covariance of a frequency weighted multivariate sample.")
100 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
104 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), meanMerged(:), covMerged(:,:)
105 real(TKG),
allocatable :: sample(:,:)
108 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
111 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
112 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
115 lb(isam)
= ub(isam
- 1)
+ 1
122 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
124 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
125 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
126 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
127 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
128 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
130 call disp%show(
"call setResized(meanMerged, ndim)")
132 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
136 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
140 call disp%show(
"cov(:,:,0) = getCov(sample, dim, iweight)")
141 cov(:,:,
0)
= getCov(sample, dim, iweight)
142 call disp%show(
"cov(:,:,0) ! reference")
144 call disp%show(
"mean(:,0) = getMean(sample, dim, iweight)")
145 mean(:,
0)
= getMean(sample, dim, iweight)
146 call disp%show(
"mean(:,0) ! reference")
149 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))")
150 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))")
153 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
154 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
156 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), uppDia)")
157 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
160 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), lowDia)")
161 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
164 call disp%show(
"call setCovMeanMerged(cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), uppDia)")
165 call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
168 call disp%show(
"cov(:,:,0) ! reference")
175 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
176 call disp%show(
"!Compute the biased merged covariance of a reliability weighted multivariate sample.")
177 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
181 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), meanMerged(:), covMerged(:,:)
182 real(TKG),
allocatable :: sample(:,:)
183 real(TKG),
allocatable :: rweight(:)
186 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
189 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
190 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
193 lb(isam)
= ub(isam
- 1)
+ 1
200 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
202 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
203 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
204 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
205 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
206 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
208 call disp%show(
"call setResized(meanMerged, ndim)")
210 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
214 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
218 call disp%show(
"cov(:,:,0) = getCov(sample, 2_IK, rweight)")
219 cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
220 call disp%show(
"cov(:,:,0) ! reference")
222 call disp%show(
"mean(:,0) = getMean(sample, dim, rweight)")
223 mean(:,
0)
= getMean(sample, dim, rweight)
224 call disp%show(
"mean(:,0) ! reference")
227 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
228 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
231 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
232 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
234 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), uppDia)")
235 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
238 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), lowDia)")
239 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
242 call disp%show(
"call setCovMeanMerged(cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), uppDia)")
243 call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
246 call disp%show(
"cov(:,:,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 (optionally unbiased) covariance matrix of a pair of (potentially weighted) t...
Generate and return the (weighted) mean of an input sample of nsam observations with ndim = 1 or 2 at...
This module contains procedures and generic interfaces for generating ranges of discrete character,...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for input/output (IO) or generic display operations on st...
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
This module contains classes and procedures for computing the first moment (i.e., the statistical mea...
Generate and return an object of type display_type.
9 lb(isam)
= ub(isam
- 1)
+ 1
17call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
18call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
23+0.401394963,
+0.979276299,
+0.534193277,
-0.450837612,
-0.838486075,
+0.965392590E-2,
+0.893495798
24+0.513590097,
-0.936435103,
-0.955090642,
-0.972116470,
+0.523586631,
-0.497346282,
+0.180513144
25cov(:,:,
0)
= getCov(sample, dim)
27+0.396623105,
-0.853262097E-1
28-0.853262097E-1,
+0.413675517
31+0.218384385,
-0.306185544
33 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
34 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
36call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
38+0.396623135,
-0.853261948E-1
39+0.307542974E-40,
+0.413675517
40call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
42+0.396623135,
-0.853261948E-1
43-0.853261948E-1,
+0.413675517
44call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
46+0.396623135,
-0.853261948E-1
47+0.149780139,
+0.413675517
49+0.396623105,
-0.853262097E-1
50-0.853262097E-1,
+0.413675517
55 lb(isam)
= ub(isam
- 1)
+ 1
63call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
64call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
69-0.203351617,
+0.236961603,
+0.381265283,
+0.188914657,
-0.636925578,
+0.796635389,
-0.868143082,
-0.752728581
70+0.822335243,
-0.457805395,
+0.132377625,
+0.770981908,
+0.220284700,
-0.338171721E-1,
+0.530756593,
-0.149759889
71-0.495133638,
+0.865869164,
+0.375614882,
+0.918408394,
-0.147484064,
-0.844699740,
-0.100102186,
-0.708728790
72-0.293231487,
-0.166109800E-1,
+0.600808501,
-0.688944459,
+0.862112999,
-0.589067459,
-0.767288208,
-0.496743679
73cov(:,:,
0)
= getCov(sample, dim)
75+0.318405479,
-0.497720167E-1,
+0.811716318E-1,
+0.342581235E-2
76-0.497720167E-1,
+0.178811550,
+0.432243943E-2,
-0.611906350E-1
77+0.811716318E-1,
+0.432243943E-2,
+0.403089046,
+0.624237731E-1
78+0.342581235E-2,
-0.611906350E-1,
+0.624237731E-1,
+0.325805992
81-0.107171491,
+0.229419202,
-0.170319974E-1,
-0.173620597
83 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
84 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
86call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
88+0.318405479,
-0.497720242E-1,
+0.811716467E-1,
+0.342581421E-2
89-0.853261948E-1,
+0.178811550,
+0.432246178E-2,
-0.611906350E-1
90+0.00000000,
+0.00000000,
+0.403089046,
+0.624237657E-1
91+0.00000000,
+0.00000000,
+0.00000000,
+0.325805932
92call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
94+0.318405479,
-0.497720242E-1,
+0.811716467E-1,
+0.342581421E-2
95-0.497720242E-1,
+0.178811550,
+0.432246178E-2,
-0.611906350E-1
96+0.811716467E-1,
+0.432246178E-2,
+0.403089046,
+0.624237657E-1
97+0.342581421E-2,
-0.611906350E-1,
+0.624237657E-1,
+0.325805932
98call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
100+0.318405479,
-0.497720242E-1,
+0.811716467E-1,
+0.342581421E-2
101-0.769999474E-1,
+0.178811550,
+0.432246178E-2,
-0.611906350E-1
102-0.154122844,
+0.761472434E-1,
+0.403089046,
+0.624237657E-1
103-0.850856304E-1,
+0.439180434E-2,
+0.882764161E-1,
+0.325805932
105+0.318405479,
-0.497720167E-1,
+0.811716318E-1,
+0.342581235E-2
106-0.497720167E-1,
+0.178811550,
+0.432243943E-2,
-0.611906350E-1
107+0.811716318E-1,
+0.432243943E-2,
+0.403089046,
+0.624237731E-1
108+0.342581235E-2,
-0.611906350E-1,
+0.624237731E-1,
+0.325805992
113 lb(isam)
= ub(isam
- 1)
+ 1
121call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
122call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
127+0.624962926,
-0.127657890,
-0.970286727,
-0.968180299,
+0.668535471,
+0.501106858,
-0.853803635,
+0.305868983,
-0.810999513,
-0.363536358
128-0.563328266,
+0.157229662,
-0.184266567E-1,
+0.781577706,
-0.996863008,
-0.805929184,
+0.193685293E-1,
-0.669980168,
+0.466191769,
-0.538042665
129+0.680117488,
-0.558356285,
+0.272832036,
+0.855279326,
+0.558489203,
+0.770297170,
-0.625860453,
+0.632045031,
-0.495243192,
+0.996196985
130cov(:,:,
0)
= getCov(sample, dim)
132+0.419856757,
-0.301753849,
+0.154538617
133-0.301753849,
+0.308248430,
-0.151841283
134+0.154538617,
-0.151841283,
+0.356217146
135mean(:,
0)
= getMean(sample, dim)
137-0.199399024,
-0.216820240,
+0.308579773
139 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
140 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
142call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
144+0.419856817,
-0.301753879,
+0.154538602
145-0.497720242E-1,
+0.308248371,
-0.151841313
146+0.811716467E-1,
+0.432246178E-2,
+0.356217146
147call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
149+0.419856817,
-0.301753879,
+0.154538602
150-0.301753879,
+0.308248371,
-0.151841313
151+0.154538602,
-0.151841313,
+0.356217146
152call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
154+0.419856817,
-0.301753879,
+0.154538602
155-0.178762853,
+0.308248371,
-0.151841313
156+0.246388942,
-0.277714342,
+0.356217146
158+0.419856757,
-0.301753849,
+0.154538617
159-0.301753849,
+0.308248430,
-0.151841283
160+0.154538617,
-0.151841283,
+0.356217146
165 lb(isam)
= ub(isam
- 1)
+ 1
173call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
174call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
179+0.181960821,
-0.686401606,
+0.763722539,
+0.147343755,
-0.626788020,
+0.541644573,
+0.115449786,
+0.175318241,
-0.708830237
180+0.932233572,
+0.178673267,
+0.900178432,
-0.703949094,
-0.655135632,
-0.518699408,
+0.518129945,
+0.346210957,
+0.430200100E-1
181+0.765156150,
-0.833746195,
-0.107750058,
+0.464287281,
+0.309900045E-1,
-0.961233974,
-0.632779241,
-0.907546401,
-0.949697256
182cov(:,:,
0)
= getCov(sample, dim)
184+0.260105550,
+0.957356989E-1,
+0.627730861E-1
185+0.957356989E-1,
+0.352771342,
+0.331181437E-1
186+0.627730861E-1,
+0.331181437E-1,
+0.385390908
187mean(:,
0)
= getMean(sample, dim)
189-0.107311280E-1,
+0.115629114,
-0.348035514
191 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
192 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
194call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
196+0.260105610,
+0.957356989E-1,
+0.627730638E-1
197-0.301753879,
+0.352771223,
+0.331181213E-1
198+0.154538602,
-0.151841313,
+0.385390788
199call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
201+0.260105610,
+0.957356989E-1,
+0.627730638E-1
202+0.957356989E-1,
+0.352771223,
+0.331181213E-1
203+0.627730638E-1,
+0.331181213E-1,
+0.385390788
204call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
206+0.260105610,
+0.957356989E-1,
+0.627730638E-1
207+0.385766290E-1,
+0.352771223,
+0.331181213E-1
208-0.886560902E-1,
-0.773828030E-1,
+0.385390788
210+0.260105550,
+0.957356989E-1,
+0.627730861E-1
211+0.957356989E-1,
+0.352771342,
+0.331181437E-1
212+0.627730861E-1,
+0.331181437E-1,
+0.385390908
217 lb(isam)
= ub(isam
- 1)
+ 1
225call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
226call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
231+0.496601105,
+0.771960616,
-0.545367956,
+0.342359900,
+0.116581798,
-0.579060316,
-0.746165395,
-0.928222656,
-0.196714878
232cov(:,:,
0)
= getCov(sample, dim)
235mean(:,
0)
= getMean(sample, dim)
239 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
240 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
242call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
245call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
248call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
257 lb(isam)
= ub(isam
- 1)
+ 1
265call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
266call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
271-0.264775157,
+0.259539127,
+0.402713537,
-0.878988385,
-0.187849283,
-0.149602890,
+0.157950640
272-0.169064045,
-0.792035818,
+0.915624857,
+0.210077763,
-0.282754064,
-0.894224644,
-0.325002909
273cov(:,:,
0)
= getCov(sample, dim)
275+0.156066045,
+0.464977697E-2
276+0.464977697E-2,
+0.324015439
277mean(:,
0)
= getMean(sample, dim)
279-0.944303498E-1,
-0.191054136
281 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
282 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
284call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
286+0.156066030,
+0.464977697E-2
287+0.307542974E-40,
+0.324015379
288call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
290+0.156066030,
+0.464977697E-2
291+0.464977697E-2,
+0.324015379
292call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
294+0.156066030,
+0.464977697E-2
295+0.437665395E-1,
+0.324015379
297+0.156066045,
+0.464977697E-2
298+0.464977697E-2,
+0.324015439
303 lb(isam)
= ub(isam
- 1)
+ 1
311call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
312call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
317+0.680088997,
+0.856828094,
+0.618274808,
+0.935130596,
+0.841377854,
+0.744893074,
+0.379459620,
-0.134975195,
+0.127558351,
+0.221639037,
-0.991853237
318-0.470487833,
+0.567016482,
-0.510230660,
-0.450091243,
+0.720406771,
+0.488170505,
-0.115825653,
-0.432596207,
-0.758372426,
-0.263102651,
-0.746771097E-1
319cov(:,:,
0)
= getCov(sample, dim)
321+0.296680480,
+0.762073770E-1
322+0.762073770E-1,
+0.223634332
323mean(:,
0)
= getMean(sample, dim)
325+0.388947457,
-0.118162736
327 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
328 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
330call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
332+0.296680480,
+0.762073845E-1
333+0.464977697E-2,
+0.223634273
334call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
336+0.296680480,
+0.762073845E-1
337+0.762073845E-1,
+0.223634273
338call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
340+0.296680480,
+0.762073845E-1
341-0.799481496E-1,
+0.223634273
343+0.296680480,
+0.762073770E-1
344+0.762073770E-1,
+0.223634332
349 lb(isam)
= ub(isam
- 1)
+ 1
357call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
358call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
363+0.352897048,
-0.625674367,
+0.769135237,
-0.107919097,
+0.280082226,
-0.278752685,
-0.608349562,
+0.600519776,
+0.405497074,
+0.741304159,
-0.499057293,
-0.325697422,
-0.781980634
364-0.694901466,
+0.871395946,
+0.504343867,
+0.553126574,
-0.502991438,
-0.261566043,
-0.342728138,
+0.821214318,
+0.220376253E-1,
+0.654630184,
+0.924949408,
-0.322317839,
+0.324513316
365-0.978813648,
+0.257770300,
-0.579929829,
-0.538145423,
+0.966723800,
-0.719702721,
+0.267160177,
+0.563987613,
-0.191644073,
+0.923226237,
-0.569140911E-1,
+0.117556810,
-0.152100086
366cov(:,:,
0)
= getCov(sample, dim)
368+0.283554226,
+0.543091865E-2,
+0.359909013E-1
369+0.543091865E-2,
+0.300668955,
+0.670931935E-1
370+0.359909013E-1,
+0.670931935E-1,
+0.340043515
371mean(:,
0)
= getMean(sample, dim)
373-0.599965686E-2,
+0.196285114,
-0.929422583E-2
375 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
376 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
378call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
380+0.283554256,
+0.543093588E-2,
+0.359908566E-1
381+0.762073845E-1,
+0.300668985,
+0.670932829E-1
382+0.00000000,
+0.331181213E-1,
+0.340043724
383call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
385+0.283554256,
+0.543093588E-2,
+0.359908566E-1
386+0.543093588E-2,
+0.300668985,
+0.670932829E-1
387+0.359908566E-1,
+0.670932829E-1,
+0.340043724
388call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
390+0.283554256,
+0.543093588E-2,
+0.359908566E-1
391+0.533605441E-1,
+0.300668985,
+0.670932829E-1
392+0.171034560,
+0.727906153E-1,
+0.340043724
394+0.283554226,
+0.543091865E-2,
+0.359909013E-1
395+0.543091865E-2,
+0.300668955,
+0.670931935E-1
396+0.359909013E-1,
+0.670931935E-1,
+0.340043515
401 lb(isam)
= ub(isam
- 1)
+ 1
409call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
410call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
415-0.220445395E-1,
+0.101730824E-1,
+0.319839120,
-0.790473342,
-0.979961276,
+0.523224354,
+0.361259222,
-0.985275745,
+0.806435108,
+0.418053746
416+0.806536674E-1,
-0.528206944,
+0.816679001E-2,
-0.585064888,
-0.757138729,
+0.746034622,
-0.862583160,
+0.560268998,
-0.294501543,
+0.642369270
417-0.357433915,
+0.894447803,
+0.903111100,
-0.908998728,
+0.777324438E-1,
+0.873299837,
+0.632227063,
+0.111449003,
-0.929121614,
+0.984597325
418cov(:,:,
0)
= getCov(sample, dim)
420+0.387672901,
+0.724123493E-1,
+0.126333743
421+0.724123493E-1,
+0.321701169,
+0.127962053
422+0.126333743,
+0.127962053,
+0.506277025
423mean(:,
0)
= getMean(sample, dim)
425-0.338770263E-1,
-0.990001932E-1,
+0.228131041
427 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
428 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
430call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
432+0.387672842,
+0.724123493E-1,
+0.126333728
433+0.543093588E-2,
+0.321701199,
+0.127962068
434+0.359908566E-1,
+0.670932829E-1,
+0.506277204
435call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
437+0.387672842,
+0.724123493E-1,
+0.126333728
438+0.724123493E-1,
+0.321701199,
+0.127962068
439+0.126333728,
+0.127962068,
+0.506277204
440call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
442+0.387672842,
+0.724123493E-1,
+0.126333728
443+0.498162918E-1,
+0.321701199,
+0.127962068
444+0.199596100E-1,
+0.167634130,
+0.506277204
446+0.387672901,
+0.724123493E-1,
+0.126333743
447+0.724123493E-1,
+0.321701169,
+0.127962053
448+0.126333743,
+0.127962053,
+0.506277025
453 lb(isam)
= ub(isam
- 1)
+ 1
461call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
462call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
467+0.403508544,
+0.520493984,
-0.125038624,
+0.593883276,
+0.493101239,
-0.395837665,
-0.303887606,
+0.722774386,
+0.105382919
468cov(:,:,
0)
= getCov(sample, dim)
471mean(:,
0)
= getMean(sample, dim)
475 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
476 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
478call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
481call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
484call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
498 lb(isam)
= ub(isam
- 1)
+ 1
506call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
507call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
512+0.988867044,
+0.862279415,
+0.291793227,
+0.295010805,
+0.802412033,
+0.140203595,
+0.631700754,
+0.820457101,
+0.116583347
515+4,
+1,
+7,
+10,
+9,
+1,
+3,
+8,
+3
516cov(:,:,
0)
= getCov(sample, dim, iweight)
519mean(:,
0)
= getMean(sample, dim, iweight)
523 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
524 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
526call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
529call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
532call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
541 lb(isam)
= ub(isam
- 1)
+ 1
549call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
550call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
555-0.285781384,
-0.375538111,
+0.799287915,
-0.833766460,
-0.325526595,
+0.790790319E-1,
+0.122757912,
-0.928408980,
+0.960289836,
+0.735855579,
+0.305666804
556+0.518711209,
-0.551659346,
-0.326926589,
+0.836314917,
-0.141353846,
+0.624673367,
-0.825314760,
-0.127631426E-1,
+0.699443460,
-0.769469857,
+0.466434956E-1
559+10,
+3,
+6,
+4,
+3,
+6,
+9,
+5,
+7,
+6,
+9
560cov(:,:,
0)
= getCov(sample, dim, iweight)
562+0.329561979,
-0.641065612E-1
563-0.641065612E-1,
+0.331481755
564mean(:,
0)
= getMean(sample, dim, iweight)
566+0.107721582,
+0.212829579E-1
568 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
569 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
571call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
573+0.329561949,
-0.641065687E-1
574+0.307542974E-40,
+0.331481785
575call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
577+0.329561949,
-0.641065687E-1
578-0.641065687E-1,
+0.331481785
579call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
581+0.329561949,
-0.641065687E-1
582+0.703652129E-1,
+0.331481785
584+0.329561979,
-0.641065612E-1
585-0.641065612E-1,
+0.331481755
590 lb(isam)
= ub(isam
- 1)
+ 1
598call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
599call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
604-0.781103492,
-0.193438172,
-0.248659492,
-0.353055000,
-0.694214106,
-0.232177138,
+0.348670483,
-0.623793602,
-0.994449615,
-0.746173382
607+9,
+5,
+9,
+6,
+3,
+1,
+2,
+6,
+10,
+2
608cov(:,:,
0)
= getCov(sample, dim, iweight)
611mean(:,
0)
= getMean(sample, dim, iweight)
615 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
616 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
618call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
621call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
624call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
633 lb(isam)
= ub(isam
- 1)
+ 1
641call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
642call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
647+0.785812497,
-0.824784756,
+0.666745663,
+0.365077257,
-0.185629249,
+0.420706153
650+10,
+4,
+6,
+6,
+1,
+1
651cov(:,:,
0)
= getCov(sample, dim, iweight)
654mean(:,
0)
= getMean(sample, dim, iweight)
658 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
659 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
661call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
664call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
667call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
676 lb(isam)
= ub(isam
- 1)
+ 1
684call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
685call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
690+0.570910573,
-0.127544522,
+0.428713202,
-0.257845879,
-0.195810795,
-0.793071866
693+4,
+4,
+7,
+10,
+3,
+5
694cov(:,:,
0)
= getCov(sample, dim, iweight)
697mean(:,
0)
= getMean(sample, dim, iweight)
701 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
702 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
704call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
707call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
710call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
719 lb(isam)
= ub(isam
- 1)
+ 1
727call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
728call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
733+0.340297818,
-0.314716816,
+0.447024822,
-0.253781915,
+0.531564593,
-0.360658050,
-0.883807182,
-0.783993602,
+0.303907156
736+1,
+6,
+5,
+3,
+6,
+2,
+4,
+7,
+10
737cov(:,:,
0)
= getCov(sample, dim, iweight)
740mean(:,
0)
= getMean(sample, dim, iweight)
744 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
745 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
747call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
750call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
753call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
762 lb(isam)
= ub(isam
- 1)
+ 1
770call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
771call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
776+0.526983738,
+0.734173536,
-0.750278234E-1,
+0.138120294,
-0.859922290,
-0.789736509,
-0.469662428,
+0.126173973,
-0.336406827,
-0.912330151E-1
777+0.595440626,
+0.790780544,
-0.110277057,
+0.260060906,
-0.881801963,
-0.463621616E-1,
+0.955803394,
-0.987929702,
-0.481496334,
-0.146150231
778-0.400498033,
+0.673115015,
-0.780857801E-1,
+0.560301185,
+0.240870476,
+0.764454842,
+0.842976093,
-0.204507709,
+0.836968541,
-0.985610485
781+6,
+10,
+8,
+7,
+8,
+6,
+2,
+7,
+10,
+2
782cov(:,:,
0)
= getCov(sample, dim, iweight)
784+0.277378827,
+0.209649324,
-0.492883772E-1
785+0.209649324,
+0.391249955,
+0.596357696E-1
786-0.492883772E-1,
+0.596357696E-1,
+0.239725173
787mean(:,
0)
= getMean(sample, dim, iweight)
789-0.659118518E-1,
-0.761377886E-1,
+0.315032393
791 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
792 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
794call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
796+0.277378827,
+0.209649414,
-0.492884852E-1
797+0.307542974E-40,
+0.391250044,
+0.596356615E-1
798+0.00000000,
+0.140129846E-44,
+0.239725292
799call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
801+0.277378827,
+0.209649414,
-0.492884852E-1
802+0.209649414,
+0.391250044,
+0.596356615E-1
803-0.492884852E-1,
+0.596356615E-1,
+0.239725292
804call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
806+0.277378827,
+0.209649414,
-0.492884852E-1
807-0.885012969E-1,
+0.391250044,
+0.596356615E-1
808-0.117471650,
+0.118432455,
+0.239725292
810+0.277378827,
+0.209649324,
-0.492883772E-1
811+0.209649324,
+0.391249955,
+0.596357696E-1
812-0.492883772E-1,
+0.596357696E-1,
+0.239725173
817 lb(isam)
= ub(isam
- 1)
+ 1
825call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
826call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
831-0.842973709,
-0.172415376,
+0.549086452,
+0.757463813,
-0.449250698,
-0.330590129,
-0.171786427
832-0.987430811E-1,
-0.832129240,
-0.561766267,
+0.797159672,
-0.800980210,
-0.812227726,
+0.451050162
833+0.260770321E-1,
+0.187863469,
+0.721985102,
+0.766978621,
-0.580928326,
-0.888453722,
+0.660766840
836+2,
+10,
+2,
+5,
+6,
+1,
+1
837cov(:,:,
0)
= getCov(sample, dim, iweight)
839+0.238151833,
+0.209738970,
+0.196257070
840+0.209738970,
+0.412076622,
+0.218169644
841+0.196257070,
+0.218169644,
+0.264301926
842mean(:,
0)
= getMean(sample, dim, iweight)
844-0.637959167E-1,
-0.400873035,
+0.129496112
846 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
847 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
849call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
851+0.238151759,
+0.209738970,
+0.196257055
852+0.209649414,
+0.412076652,
+0.218169659
853-0.492884852E-1,
+0.596356615E-1,
+0.264301926
854call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
856+0.238151759,
+0.209738970,
+0.196257055
857+0.209738970,
+0.412076652,
+0.218169659
858+0.196257055,
+0.218169659,
+0.264301926
859call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
861+0.238151759,
+0.209738970,
+0.196257055
862+0.417133272,
+0.412076652,
+0.218169659
863+0.357336581,
+0.529201031,
+0.264301926
865+0.238151833,
+0.209738970,
+0.196257070
866+0.209738970,
+0.412076622,
+0.218169644
867+0.196257070,
+0.218169644,
+0.264301926
872 lb(isam)
= ub(isam
- 1)
+ 1
880call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
881call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
886-0.899042368,
+0.415559888,
+0.658540368,
-0.353461504E-1,
-0.461444497,
+0.355355263,
+0.992943525,
-0.881807804,
-0.572038889E-1,
+0.793961287
887+0.375169516E-1,
+0.421079159,
+0.912170768,
+0.649054885,
+0.140879393,
-0.322211623,
+0.166822314,
-0.508674026,
-0.735155463,
+0.518001437
888+0.701062679E-1,
-0.610464692,
+0.537359357,
+0.243846416,
-0.413061261,
+0.626007318,
-0.272456288,
+0.724389315,
-0.332440138E-1,
-0.903913736
891+1,
+3,
+6,
+5,
+8,
+7,
+6,
+2,
+7,
+9
892cov(:,:,
0)
= getCov(sample, dim, iweight)
894+0.298760593,
+0.118387401,
-0.699651837E-1
895+0.118387401,
+0.261046827,
-0.596563146E-1
896-0.699651837E-1,
-0.596563146E-1,
+0.290581763
897mean(:,
0)
= getMean(sample, dim, iweight)
899+0.256617576,
+0.155372545,
-0.887820199E-1
901 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
902 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
904call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
906+0.298760474,
+0.118387401,
-0.699651986E-1
907+0.209738970,
+0.261046827,
-0.596563220E-1
908+0.196257055,
+0.218169659,
+0.290581733
909call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
911+0.298760474,
+0.118387401,
-0.699651986E-1
912+0.118387401,
+0.261046827,
-0.596563220E-1
913-0.699651986E-1,
-0.596563220E-1,
+0.290581733
914call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
916+0.298760474,
+0.118387401,
-0.699651986E-1
917+0.931684598E-1,
+0.261046827,
-0.596563220E-1
918-0.112460896,
-0.134293810,
+0.290581733
920+0.298760593,
+0.118387401,
-0.699651837E-1
921+0.118387401,
+0.261046827,
-0.596563146E-1
922-0.699651837E-1,
-0.596563146E-1,
+0.290581763
927 lb(isam)
= ub(isam
- 1)
+ 1
935call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
936call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
941+0.888599515,
-0.585449457,
+0.648303151,
-0.978484273,
+0.484613538,
-0.912868023,
-0.652709126,
+0.463062048
942+0.143523693,
+0.640499353,
+0.379489541,
-0.746734619,
+0.141961813,
+0.387100220,
+0.328666449,
-0.565282702
945+1,
+3,
+3,
+5,
+6,
+1,
+1,
+5
946cov(:,:,
0)
= getCov(sample, dim, iweight)
948+0.464220256,
+0.767516866E-1
949+0.767516866E-1,
+0.252795905
950mean(:,
0)
= getMean(sample, dim, iweight)
952-0.631384831E-2,
-0.715623572E-1
954 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
955 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
957call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
959+0.464220226,
+0.767517015E-1
960+0.118387401,
+0.252795935
961call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
963+0.464220226,
+0.767517015E-1
964+0.767517015E-1,
+0.252795935
965call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
967+0.464220226,
+0.767517015E-1
968-0.138533592,
+0.252795935
970+0.464220256,
+0.767516866E-1
971+0.767516866E-1,
+0.252795905
981 lb(isam)
= ub(isam
- 1)
+ 1
989call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
990call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
995-0.571417451,
+0.931302309,
+0.977839231,
-0.977189541,
+0.124298334E-1,
-0.929996610,
-0.237989902,
-0.665398240
996-0.721635222,
-0.580464959,
+0.696216345,
+0.908992887,
+0.697550774,
+0.947196960,
-0.743318677,
+0.271038175
999+1.74241102,
+1.07246280,
+1.34209347,
+1.74612665,
+1.54444861,
+1.87357593,
+1.24614286,
+1.18447518
1000cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1002+0.493098676,
-0.114245765
1003-0.114245765,
+0.498654991
1004mean(:,
0)
= getMean(sample, dim, rweight)
1006-0.272193044,
+0.245787114
1008 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1009 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1011call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1013+0.493098617,
-0.114245795
1014+0.307542974E-40,
+0.498655081
1015call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1017+0.493098617,
-0.114245795
1018-0.114245795,
+0.498655081
1019call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1021+0.493098617,
-0.114245795
1022-0.625649914E-1,
+0.498655081
1024+0.493098676,
-0.114245765
1025-0.114245765,
+0.498654991
1030 lb(isam)
= ub(isam
- 1)
+ 1
1038call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1039call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1044-0.100600958,
+0.465979815,
+0.937115550,
+0.947117805E-2,
-0.884447813,
+0.688210726E-1,
-0.206218958E-1
1045-0.982659101,
-0.230036139,
-0.986904025,
-0.227110505,
+0.249925733,
+0.933954835,
-0.398386002
1048+1.35797143,
+1.48018575,
+1.08264220,
+1.86524141,
+1.07807255,
+1.06541657,
+1.98629522
1049cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1051+0.210854188,
-0.101127788
1052-0.101127788,
+0.314686835
1053mean(:,
0)
= getMean(sample, dim, rweight)
1055+0.669851825E-1,
-0.271669000
1057 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1058 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1060call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1062+0.210854143,
-0.101127759
1063-0.114245795,
+0.314686865
1064call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1066+0.210854143,
-0.101127759
1067-0.101127759,
+0.314686865
1068call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1070+0.210854143,
-0.101127759
1071-0.250287931E-1,
+0.314686865
1073+0.210854188,
-0.101127788
1074-0.101127788,
+0.314686835
1079 lb(isam)
= ub(isam
- 1)
+ 1
1087call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1088call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1093+0.828327656,
+0.342382431,
-0.813050389,
+0.176590443,
+0.405714512,
-0.922640324,
+0.829164505
1096+1.14408064,
+1.77327693,
+1.16854942,
+1.93906152,
+1.79191124,
+1.28518343,
+1.46284866
1097cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1100mean(:,
0)
= getMean(sample, dim, rweight)
1104 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1105 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1107call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1110call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1113call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1122 lb(isam)
= ub(isam
- 1)
+ 1
1130call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1131call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1136+0.888656020,
+0.824515104,
+0.657945395,
-0.835144401,
-0.136049867,
+0.254159451
1137+0.439628601,
+0.224233150,
-0.341976762,
-0.504526258,
+0.470494628,
+0.399492979E-1
1138+0.936748981E-1,
+0.357176065E-1,
-0.518700480,
+0.749618888,
+0.589431524,
+0.201505542
1141+1.11516237,
+1.84511876,
+1.37388158,
+1.87129498,
+1.12207687,
+1.71832609
1142cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1144+0.404780746,
+0.130590633,
-0.214089751
1145+0.130590633,
+0.132172212,
-0.106344931E-1
1146-0.214089751,
-0.106344931E-1,
+0.166609764
1147mean(:,
0)
= getMean(sample, dim, rweight)
1149+0.236299470,
+0.957544707E-2,
+0.206517741
1151 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1152 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1154call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1156+0.404780746,
+0.130590633,
-0.214089751
1157+0.307542974E-40,
+0.132172182,
-0.106344968E-1
1158+0.00000000,
+0.140203595,
+0.166609764
1159call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1161+0.404780746,
+0.130590633,
-0.214089751
1162+0.130590633,
+0.132172182,
-0.106344968E-1
1163-0.214089751,
-0.106344968E-1,
+0.166609764
1164call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1166+0.404780746,
+0.130590633,
-0.214089751
1167+0.135784522,
+0.132172182,
-0.106344968E-1
1168-0.110218026,
-0.434923284E-1,
+0.166609764
1170+0.404780746,
+0.130590633,
-0.214089751
1171+0.130590633,
+0.132172212,
-0.106344931E-1
1172-0.214089751,
-0.106344931E-1,
+0.166609764
1177 lb(isam)
= ub(isam
- 1)
+ 1
1185call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1186call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1191+0.832262516,
-0.127016664,
+0.675074220,
-0.692545176E-1,
-0.312504053,
-0.147016168,
-0.656745791,
-0.496482849,
-0.112376213,
-0.515277267,
+0.206839681,
+0.527799249,
+0.666803122E-1
1192-0.189897180,
+0.834278584,
+0.197322369E-1,
+0.380127549,
+0.828397870,
+0.599107385,
-0.388283610,
-0.889161468,
+0.149561167E-1,
+0.280150890,
+0.230778813,
+0.975948691,
-0.463267922
1193-0.672663331,
-0.887493372,
+0.931710243,
-0.195828080,
-0.598317385,
-0.292548776,
+0.232463956,
-0.939492345,
-0.344697952,
-0.279079795,
-0.455149651,
+0.930317283,
+0.960207701
1196+1.17317104,
+1.08978057,
+1.03113854,
+1.65376806,
+1.15174007,
+1.87458849,
+1.07473564,
+1.15723777,
+1.65951502,
+1.98242521,
+1.76555133,
+1.12642741,
+1.52484775
1197cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1199+0.171400830,
+0.239375029E-1,
+0.849001184E-1
1200+0.239375029E-1,
+0.252791524,
-0.173558444E-1
1201+0.849001184E-1,
-0.173558444E-1,
+0.368183792
1202mean(:,
0)
= getMean(sample, dim, rweight)
1204-0.352020413E-1,
+0.183240026,
-0.142958015
1206 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1207 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1209call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1211+0.171400696,
+0.239375569E-1,
+0.849001408E-1
1212+0.130590633,
+0.252791554,
-0.173558611E-1
1213-0.214089751,
-0.106344968E-1,
+0.368183732
1214call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1216+0.171400696,
+0.239375569E-1,
+0.849001408E-1
1217+0.239375588E-1,
+0.252791554,
-0.173558611E-1
1218+0.849001408E-1,
-0.173558593E-1,
+0.368183732
1219call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1221+0.171400696,
+0.239375569E-1,
+0.849001408E-1
1222+0.107305847,
+0.252791554,
-0.173558611E-1
1223+0.114116281,
+0.839109495E-1,
+0.368183732
1225+0.171400830,
+0.239375029E-1,
+0.849001184E-1
1226+0.239375029E-1,
+0.252791524,
-0.173558444E-1
1227+0.849001184E-1,
-0.173558444E-1,
+0.368183792
1232 lb(isam)
= ub(isam
- 1)
+ 1
1240call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1241call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1246+0.404199958,
+0.385584474,
+0.962077379,
+0.989631295,
+0.906902075,
+0.369826794,
+0.214111447,
+0.719902635,
+0.955797791,
+0.859641910,
+0.374192119
1247+0.916485667,
+0.685104728,
-0.771127343,
-0.129209995,
+0.691563368,
-0.682598472,
-0.353188515E-1,
-0.270760059,
+0.834522605,
+0.212278366,
-0.663943291E-1
1248-0.448723078,
+0.341016531,
-0.316381097,
+0.161957860,
-0.454055309,
-0.789866090,
-0.953052044E-1,
-0.617180824,
-0.371762156,
+0.395909309,
-0.237938404
1251+1.25719142,
+1.88753462,
+1.78973460,
+1.63784313,
+1.07069731,
+1.45528889,
+1.52669919,
+1.99930871,
+1.52823269,
+1.76109016,
+1.27821636
1252cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1254+0.798297375E-1,
-0.112736505E-1,
+0.807016809E-2
1255-0.112736505E-1,
+0.312828720,
+0.651230067E-1
1256+0.807016809E-2,
+0.651230067E-1,
+0.144341275
1257mean(:,
0)
= getMean(sample, dim, rweight)
1259+0.657695055,
+0.913084820E-1,
-0.198439956
1261 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1262 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1264call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1266+0.798297450E-1,
-0.112736579E-1,
+0.807017367E-2
1267+0.239375588E-1,
+0.312828660,
+0.651229993E-1
1268+0.849001408E-1,
-0.173558593E-1,
+0.144341305
1269call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1271+0.798297450E-1,
-0.112736579E-1,
+0.807017367E-2
1272-0.112736579E-1,
+0.312828660,
+0.651229993E-1
1273+0.807017367E-2,
+0.651229993E-1,
+0.144341305
1274call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1276+0.798297450E-1,
-0.112736579E-1,
+0.807017367E-2
1277+0.901869610E-1,
+0.312828660,
+0.651229993E-1
1278+0.167856105E-1,
+0.656217933E-1,
+0.144341305
1280+0.798297375E-1,
-0.112736505E-1,
+0.807016809E-2
1281-0.112736505E-1,
+0.312828720,
+0.651230067E-1
1282+0.807016809E-2,
+0.651230067E-1,
+0.144341275
1287 lb(isam)
= ub(isam
- 1)
+ 1
1295call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1296call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1301-0.863177776E-1,
+0.153170824,
-0.867692471,
+0.265655756,
+0.258558869,
+0.170848370E-1,
+0.649551392,
-0.763564944,
-0.666157007,
-0.504313231,
+0.824491262,
+0.474032640,
+0.132736564
1302+0.615460396,
+0.999525785E-1,
+0.774540186,
+0.651717782,
+0.994927049,
-0.955973148,
+0.117524505,
+0.375386596,
-0.455659151,
-0.404311299,
+0.775979519,
+0.570208192,
+0.464315414
1303+0.677419186,
-0.516756177,
-0.764681458,
+0.816069722,
-0.223775864,
-0.983417034E-2,
-0.323994160,
+0.495145798,
+0.383114815,
+0.322707653,
+0.425106883,
+0.324522257E-1,
-0.161771059
1304+0.717651963,
+0.416641474,
-0.685375571,
+0.877829790,
-0.771937609,
-0.419826031,
+0.795935512,
-0.852182627,
+0.297023416,
-0.382172465,
-0.691636205,
-0.291558981,
+0.744650722
1305-0.110471249E-1,
-0.543419719,
-0.533681154,
+0.893100619,
-0.272954702,
-0.570667624,
-0.508905411,
-0.869824409,
-0.593479633,
-0.367645383,
-0.568282604E-1,
+0.868855476,
-0.526344299
1306-0.669735193,
-0.908916593,
+0.354645371,
+0.552516103,
+0.761848688,
-0.593866110,
+0.405617833,
-0.139012694,
+0.181772470,
+0.539301634,
-0.897852182E-1,
-0.131992459,
-0.417092800
1309+1.58328199,
+1.16560173,
+1.15800667,
+1.47889853,
+1.47997761,
+1.28861511,
+1.32251430,
+1.47728014,
+1.29625607,
+1.83890510,
+1.22764826,
+1.77725387,
+1.29353333
1310cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1312+0.262752771,
+0.880049840E-1,
-0.105997641E-1,
+0.883777365E-1,
+0.142275095,
-0.193362106E-1
1313+0.880049840E-1,
+0.301021129,
-0.737055438E-2,
-0.137093151E-1,
+0.120418891,
+0.496521890E-1
1314-0.105997641E-1,
-0.737055438E-2,
+0.198423773,
+0.504211709E-1,
+0.874293372E-1,
+0.475427136E-2
1315+0.883777365E-1,
-0.137093151E-1,
+0.504211709E-1,
+0.410423517,
+0.821109265E-1,
-0.693959594E-1
1316+0.142275095,
+0.120418891,
+0.874293372E-1,
+0.821109265E-1,
+0.299637049,
+0.404402576E-1
1317-0.193362106E-1,
+0.496521890E-1,
+0.475427136E-2,
-0.693959594E-1,
+0.404402576E-1,
+0.254095823
1318mean(:,
0)
= getMean(sample, dim, rweight)
1320-0.108147785E-1,
+0.279249042,
+0.120246388,
-0.263803434E-1,
-0.201071024,
+0.784909632E-2
1322 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1323 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1325call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1327+0.262752801,
+0.880050063E-1,
-0.105997575E-1,
+0.883777142E-1,
+0.142275095,
-0.193362311E-1
1328-0.112736579E-1,
+0.301021159,
-0.737057161E-2,
-0.137093123E-1,
+0.120418914,
+0.496522300E-1
1329+0.807017367E-2,
+0.651229993E-1,
+0.198423848,
+0.504211709E-1,
+0.874293670E-1,
+0.475430535E-2
1330+0.00000000,
+0.00000000,
+0.00000000,
+0.410423577,
+0.821109340E-1,
-0.693959594E-1
1331+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.299637049,
+0.404402614E-1
1332+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.254095703
1333call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1335+0.262752801,
+0.880050063E-1,
-0.105997575E-1,
+0.883777142E-1,
+0.142275095,
-0.193362311E-1
1336+0.880050063E-1,
+0.301021159,
-0.737057161E-2,
-0.137093123E-1,
+0.120418914,
+0.496522300E-1
1337-0.105997585E-1,
-0.737057161E-2,
+0.198423848,
+0.504211709E-1,
+0.874293670E-1,
+0.475430535E-2
1338+0.883777142E-1,
-0.137093123E-1,
+0.504211709E-1,
+0.410423577,
+0.821109340E-1,
-0.693959594E-1
1339+0.142275095,
+0.120418914,
+0.874293670E-1,
+0.821109340E-1,
+0.299637049,
+0.404402614E-1
1340-0.193362311E-1,
+0.496522300E-1,
+0.475430535E-2,
-0.693959594E-1,
+0.404402614E-1,
+0.254095703
1341call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1343+0.262752801,
+0.880050063E-1,
-0.105997575E-1,
+0.883777142E-1,
+0.142275095,
-0.193362311E-1
1344+0.209802270,
+0.301021159,
-0.737057161E-2,
-0.137093123E-1,
+0.120418914,
+0.496522300E-1
1345-0.605140515E-1,
-0.326006003E-1,
+0.198423848,
+0.504211709E-1,
+0.874293670E-1,
+0.475430535E-2
1346+0.482197665E-2,
-0.531651005E-1,
-0.840576440E-1,
+0.410423577,
+0.821109340E-1,
-0.693959594E-1
1347+0.249660119,
+0.116730630,
-0.608433075E-1,
-0.240780152E-1,
+0.299637049,
+0.404402614E-1
1348-0.853117481E-1,
-0.121517599,
+0.345167629E-1,
-0.404614992E-1,
-0.235077329E-1,
+0.254095703
1350+0.262752771,
+0.880049840E-1,
-0.105997641E-1,
+0.883777365E-1,
+0.142275095,
-0.193362106E-1
1351+0.880049840E-1,
+0.301021129,
-0.737055438E-2,
-0.137093151E-1,
+0.120418891,
+0.496521890E-1
1352-0.105997641E-1,
-0.737055438E-2,
+0.198423773,
+0.504211709E-1,
+0.874293372E-1,
+0.475427136E-2
1353+0.883777365E-1,
-0.137093151E-1,
+0.504211709E-1,
+0.410423517,
+0.821109265E-1,
-0.693959594E-1
1354+0.142275095,
+0.120418891,
+0.874293372E-1,
+0.821109265E-1,
+0.299637049,
+0.404402576E-1
1355-0.193362106E-1,
+0.496521890E-1,
+0.475427136E-2,
-0.693959594E-1,
+0.404402576E-1,
+0.254095823
1360 lb(isam)
= ub(isam
- 1)
+ 1
1368call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1369call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1374+0.457989335,
+0.810496807E-1,
-0.654122233,
+0.952951312,
-0.175797582,
-0.288653135,
+0.610257149,
-0.605512977,
+0.164504528
1375-0.379562378E-3,
+0.612921000,
-0.233135819,
+0.323578358,
+0.360003829,
+0.951012969,
-0.248130918,
+0.459004998,
+0.377686262
1378+1.99604058,
+1.03367066,
+1.79838407,
+1.52714586,
+1.87341881,
+1.16860247,
+1.90648484,
+1.74027216,
+1.75829124
1379cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1381+0.274544597,
-0.367375016E-1
1382-0.367375016E-1,
+0.125712201
1383mean(:,
0)
= getMean(sample, dim, rweight)
1385+0.681750923E-1,
+0.235321149
1387 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1388 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1390call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1392+0.274544597,
-0.367375202E-1
1393+0.880050063E-1,
+0.125712141
1394call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1396+0.274544597,
-0.367375202E-1
1397-0.367375202E-1,
+0.125712141
1398call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1400+0.274544597,
-0.367375202E-1
1401-0.305830371E-1,
+0.125712141
1403+0.274544597,
-0.367375016E-1
1404-0.367375016E-1,
+0.125712201
1409 lb(isam)
= ub(isam
- 1)
+ 1
1417call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1418call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1423-0.165395021,
+0.555638313,
-0.414738774,
+0.955900669,
+0.472225070,
+0.573829889,
-0.314169765,
+0.488227963,
+0.341623902
1426+1.10429144,
+1.66944027,
+1.59734917,
+1.46626925,
+1.16047907,
+1.60012817,
+1.47653294,
+1.06352961,
+1.09223926
1427cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1430mean(:,
0)
= getMean(sample, dim, rweight)
1434 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1435 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1437call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1440call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1443call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1452 lb(isam)
= ub(isam
- 1)
+ 1
1460call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1461call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1466+0.661407113,
+0.921583056,
+0.328138232,
-0.931953192,
-0.989997864,
-0.397196651,
+0.530970097E-1,
+0.857881069,
+0.103220701,
-0.173329234,
+0.921910048,
-0.611423254E-1
1467+0.162328124,
-0.912275314,
+0.334692001E-1,
-0.149823785,
+0.845754743,
+0.403283119,
+0.381153345,
+0.251645207,
+0.405446649,
-0.707323909,
+0.597941875E-2,
+0.683746815
1468-0.376863718,
-0.151410222,
-0.118592024,
-0.125443101,
-0.994588852,
-0.908079505,
-0.937530160,
-0.672951102,
-0.983699679,
-0.565156698,
+0.949061751,
+0.775063634
1469+0.313614607,
+0.163513303,
+0.915796041,
+0.971271992E-1,
-0.107908368,
-0.381834149,
+0.244665146,
-0.621290803,
-0.405495882,
-0.538779974,
-0.101119280E-1,
+0.921919703
1472+1.83676660,
+1.25321245,
+1.51819849,
+1.94060791,
+1.59117222,
+1.13043523,
+1.51070869,
+1.50350833,
+1.81314301,
+1.87763572,
+1.20856524,
+1.23270035
1473cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1475+0.411606401,
-0.870113969E-1,
+0.112263426,
+0.210920442E-1
1476-0.870113969E-1,
+0.238117605,
-0.557373241E-1,
+0.350888185E-1
1477+0.112263426,
-0.557373241E-1,
+0.342847049,
+0.151569068
1478+0.210920442E-1,
+0.350888185E-1,
+0.151569068,
+0.229862705
1479mean(:,
0)
= getMean(sample, dim, rweight)
1481+0.708988383E-1,
+0.104679950,
-0.384704679,
+0.309142582E-1
1483 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1484 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1486call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1488+0.411606431,
-0.870114118E-1,
+0.112263419,
+0.210920833E-1
1489+0.307542974E-40,
+0.238117635,
-0.557373203E-1,
+0.350888222E-1
1490+0.00000000,
+0.731962323E-1,
+0.342847049,
+0.151569054
1491+0.00000000,
+0.353250235,
+0.704036851E-3,
+0.229862720
1492call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1494+0.411606431,
-0.870114118E-1,
+0.112263419,
+0.210920833E-1
1495-0.870114118E-1,
+0.238117635,
-0.557373203E-1,
+0.350888222E-1
1496+0.112263419,
-0.557373203E-1,
+0.342847049,
+0.151569054
1497+0.210920833E-1,
+0.350888222E-1,
+0.151569054,
+0.229862720
1498call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1500+0.411606431,
-0.870114118E-1,
+0.112263419,
+0.210920833E-1
1501+0.528442077E-1,
+0.238117635,
-0.557373203E-1,
+0.350888222E-1
1502+0.830376446E-1,
+0.758205354E-1,
+0.342847049,
+0.151569054
1503-0.440159403E-1,
+0.142612323,
+0.318926066,
+0.229862720
1505+0.411606401,
-0.870113969E-1,
+0.112263426,
+0.210920442E-1
1506-0.870113969E-1,
+0.238117605,
-0.557373241E-1,
+0.350888185E-1
1507+0.112263426,
-0.557373241E-1,
+0.342847049,
+0.151569068
1508+0.210920442E-1,
+0.350888185E-1,
+0.151569068,
+0.229862705