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(:,:,:), 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, 1_IK], [ndim, nsam])")
54 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
55 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
57 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
61 call disp%show(
"cov(:,:,0) = getCov(sample, dim)")
62 cov(:,:,
0)
= getCov(sample, dim)
63 call disp%show(
"cov(:,:,0) ! reference")
66 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim)")
67 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
70 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
71 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
73 call disp%show(
"covMerged = getCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
74 covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
77 call disp%show(
"covMerged = getCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG))")
78 covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
81 call disp%show(
"cov(:,:,0) ! reference")
88 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
89 call disp%show(
"!Compute the biased merged covariance of a frequency weighted multivariate sample.")
90 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
94 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), covMerged(:,:)
95 real(TKG),
allocatable :: sample(:,:)
98 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
101 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
102 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
105 lb(isam)
= ub(isam
- 1)
+ 1
112 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
114 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
115 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
116 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
117 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
118 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
120 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
124 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
128 call disp%show(
"cov(:,:,0) = getCov(sample, 2_IK, iweight)")
129 cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
130 call disp%show(
"cov(:,:,0) ! reference")
133 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
134 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, iweight(lb(isam):ub(isam)))")
137 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
138 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
140 call disp%show(
"covMerged = getCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
141 covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
144 call disp%show(
"covMerged = getCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG))")
145 covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
148 call disp%show(
"cov(:,:,0) ! reference")
155 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
156 call disp%show(
"!Compute the biased merged covariance of a reliability weighted multivariate sample.")
157 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
161 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), covMerged(:,:)
162 real(TKG),
allocatable :: sample(:,:)
163 real(TKG),
allocatable :: rweight(:)
166 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
169 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
170 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
173 lb(isam)
= ub(isam
- 1)
+ 1
180 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
182 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
183 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
184 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
185 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
186 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
188 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
192 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
196 call disp%show(
"cov(:,:,0) = getCov(sample, 2_IK, rweight)")
197 cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
198 call disp%show(
"cov(:,:,0) ! reference")
201 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
202 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
205 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
206 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
208 call disp%show(
"covMerged = getCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
209 covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
212 call disp%show(
"covMerged = getCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG))")
213 covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
216 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,
1_IK], [ndim, nsam])
22+0.516826034,
+0.792227030,
+0.328134894,
-0.440360904,
-0.657570362E-1
23cov(:,:,
0)
= getCov(sample, dim)
27 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
28 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
30covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
33covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
42 lb(isam)
= ub(isam
- 1)
+ 1
50call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
51call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
55-0.929594994,
+0.436816573,
+0.632729769,
-0.539220691,
-0.394586325E-1,
+0.920131207,
+0.127034664,
-0.955117226,
-0.688022375E-1,
+0.993142366
56cov(:,:,
0)
= getCov(sample, dim)
60 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
61 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
63covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
66covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
75 lb(isam)
= ub(isam
- 1)
+ 1
83call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
84call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
88+0.361089468,
-0.535390854,
-0.968912482,
+0.264635324,
+0.120630264E-1,
-0.544026017,
+0.192158103,
-0.520289779,
+0.350813031,
+0.291080236,
-0.880505919
89+0.230461359E-1,
-0.993716717,
+0.169888854,
-0.231464386,
+0.228006840E-1,
+0.904440999,
+0.623671889,
+0.202726364,
-0.415415168,
+0.859754682,
-0.943331957
90cov(:,:,
0)
= getCov(sample, dim)
92+0.241568461,
+0.738365054E-1
93+0.738365054E-1,
+0.374196976
95 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
96 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
98covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
100+0.241568431,
+0.738364905E-1
101+0.738364905E-1,
+0.374196887
102covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
104+0.241568431,
+0.738364905E-1
105+0.738364905E-1,
+0.374196887
107+0.241568461,
+0.738365054E-1
108+0.738365054E-1,
+0.374196976
113 lb(isam)
= ub(isam
- 1)
+ 1
121call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
122call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
126+0.568608642,
+0.869689822,
+0.628063679,
-0.522292852E-1,
-0.551448584,
-0.512825489
127+0.163981318,
-0.136163354,
+0.784253836,
+0.885677934,
-0.510733008,
+0.114840508
128cov(:,:,
0)
= getCov(sample, dim)
130+0.315596938,
+0.729629099E-1
131+0.729629099E-1,
+0.239412397
133 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
134 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
136covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
138+0.315596938,
+0.729629025E-1
139+0.729629025E-1,
+0.239412367
140covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
142+0.315596938,
+0.729629025E-1
143+0.729629025E-1,
+0.239412367
145+0.315596938,
+0.729629099E-1
146+0.729629099E-1,
+0.239412397
151 lb(isam)
= ub(isam
- 1)
+ 1
159call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
160call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
164-0.534685254,
-0.672600985,
-0.407761097,
-0.432819366,
+0.605178475,
+0.538234591,
+0.866307020E-1,
+0.102880359,
+0.457062840,
-0.804008722,
-0.598505855,
-0.173407555,
+0.516078472
165-0.796876550,
-0.429785609,
+0.338528395,
-0.896916628,
-0.511332273,
+0.690571547,
-0.848548889,
-0.212320924,
+0.161644340,
+0.696067929,
+0.913596034,
-0.505671740,
-0.359255791
166+0.197930813,
-0.555932641,
-0.174853206,
-0.202628613,
-0.555898070,
-0.216881752,
-0.492387295,
+0.386825919,
+0.472961426,
-0.590429783,
+0.975373507,
+0.333472133,
+0.550759554
167+0.718626618,
+0.314231277,
+0.543329477,
+0.456759453,
+0.812556386,
+0.734400988,
+0.160040140,
+0.821339369,
+0.929993272,
+0.505266190,
-0.156544566,
-0.173611283,
-0.507399321
168cov(:,:,
0)
= getCov(sample, dim)
170+0.241714865,
-0.289614219E-1,
+0.244316831E-1,
+0.287885591E-1
171-0.289614219E-1,
+0.364150763,
+0.677974597E-1,
+0.135139870E-1
172+0.244316831E-1,
+0.677974597E-1,
+0.238841236,
-0.762695372E-1
173+0.287885591E-1,
+0.135139870E-1,
-0.762695372E-1,
+0.184905261
175 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
176 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
178covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
180+0.241714865,
-0.289613716E-1,
+0.244316962E-1,
+0.287885480E-1
181-0.289613716E-1,
+0.364150882,
+0.677974373E-1,
+0.135139860E-1
182+0.244316962E-1,
+0.677974373E-1,
+0.238841221,
-0.762695372E-1
183+0.287885480E-1,
+0.135139860E-1,
-0.762695372E-1,
+0.184905261
184covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
186+0.241714865,
-0.289613716E-1,
+0.244316962E-1,
+0.287885480E-1
187-0.289613716E-1,
+0.364150882,
+0.677974373E-1,
+0.135139860E-1
188+0.244316962E-1,
+0.677974373E-1,
+0.238841221,
-0.762695372E-1
189+0.287885480E-1,
+0.135139860E-1,
-0.762695372E-1,
+0.184905261
191+0.241714865,
-0.289614219E-1,
+0.244316831E-1,
+0.287885591E-1
192-0.289614219E-1,
+0.364150763,
+0.677974597E-1,
+0.135139870E-1
193+0.244316831E-1,
+0.677974597E-1,
+0.238841236,
-0.762695372E-1
194+0.287885591E-1,
+0.135139870E-1,
-0.762695372E-1,
+0.184905261
199 lb(isam)
= ub(isam
- 1)
+ 1
207call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
208call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
212-0.609524846,
+0.903454900,
-0.992637634,
-0.367419004,
-0.395920634,
+0.878209352,
+0.877282619E-1,
+0.142601967,
-0.619130135E-1
213-0.225129008,
+0.984055519,
-0.917810440,
+0.593317866,
+0.360621452,
+0.929785371,
+0.132153511,
-0.274257064,
+0.233273268
214cov(:,:,
0)
= getCov(sample, dim)
216+0.360974699,
+0.270551324
217+0.270551324,
+0.332073838
219 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
220 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
222covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
224+0.360974669,
+0.270551264
225+0.270551264,
+0.332073718
226covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
228+0.360974669,
+0.270551264
229+0.270551264,
+0.332073718
231+0.360974699,
+0.270551324
232+0.270551324,
+0.332073838
237 lb(isam)
= ub(isam
- 1)
+ 1
245call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
246call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
250+0.365581155,
-0.731999874E-1,
-0.106020689,
+0.480873466,
-0.866870403,
-0.144724965,
+0.486977220,
-0.958377838,
+0.492727280,
-0.502603650
251+0.846296549E-1,
-0.704701543,
+0.119326115E-1,
+0.612583756,
+0.363431931,
-0.785944223,
+0.561885834E-1,
+0.121662021,
-0.946765304,
+0.375353932
252cov(:,:,
0)
= getCov(sample, dim)
254+0.273675501,
-0.636867508E-1
255-0.636867508E-1,
+0.261829853
257 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
258 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
260covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
262+0.273675442,
-0.636867881E-1
263-0.636867881E-1,
+0.261829883
264covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
266+0.273675442,
-0.636867881E-1
267-0.636867881E-1,
+0.261829883
269+0.273675501,
-0.636867508E-1
270-0.636867508E-1,
+0.261829853
275 lb(isam)
= ub(isam
- 1)
+ 1
283call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
284call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
288-0.349876881E-1,
-0.626652241E-1,
-0.591000199,
+0.294861317,
+0.597044230,
-0.914378405,
-0.837277174E-1,
+0.373550534,
-0.429550409,
-0.663632393,
+0.433650494,
-0.258920670
289-0.908768177E-2,
+0.549182296,
-0.112284660,
+0.657070637,
-0.665625811,
-0.785636902,
+0.635945082,
-0.436999202,
-0.299843550,
-0.986699581,
+0.663479686,
-0.250535011
290+0.660431743,
-0.825237751,
+0.444016695,
+0.427774310,
+0.234203815,
-0.325801373E-1,
-0.164177299,
-0.801920891E-1,
+0.848542571,
-0.327888131,
-0.937605739,
-0.472491264
291cov(:,:,
0)
= getCov(sample, dim)
293+0.209242567,
+0.112537101,
-0.342129432E-1
294+0.112537101,
+0.323154479,
-0.734827146E-1
295-0.342129432E-1,
-0.734827146E-1,
+0.292690605
297 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
298 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
300covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
302+0.209242582,
+0.112537108,
-0.342129506E-1
303+0.112537108,
+0.323154479,
-0.734827220E-1
304-0.342129506E-1,
-0.734827220E-1,
+0.292690665
305covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
307+0.209242582,
+0.112537108,
-0.342129506E-1
308+0.112537108,
+0.323154479,
-0.734827220E-1
309-0.342129506E-1,
-0.734827220E-1,
+0.292690665
311+0.209242567,
+0.112537101,
-0.342129432E-1
312+0.112537101,
+0.323154479,
-0.734827146E-1
313-0.342129432E-1,
-0.734827146E-1,
+0.292690605
318 lb(isam)
= ub(isam
- 1)
+ 1
326call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
327call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
331+0.424404979,
+0.785670638,
+0.388367772,
+0.858890414,
-0.148493886,
+0.726894140E-1,
+0.813536286,
+0.248022437,
+0.544670820E-1,
-0.638126850,
-0.220623732,
-0.635957003
332+0.773164868,
+0.796957731,
-0.723076344,
+0.990809798,
+0.232962012,
-0.180916667,
-0.681890249E-1,
+0.666097045,
-0.925797224,
-0.814512849,
+0.454158306,
-0.361469984
333-0.754204154,
-0.475651622,
+0.171577454,
+0.533966899,
-0.457488298,
+0.907611251,
+0.941787362,
-0.582931757,
+0.145285845,
+0.744918346,
-0.427311778,
-0.724592209
334-0.320654511,
+0.454595089,
+0.458431005,
+0.430775881,
-0.268006086,
-0.505381823E-1,
+0.928588629,
+0.728483081,
-0.914470553,
+0.148860455,
-0.463771820,
+0.725812316
335+0.641360521,
+0.543218255,
+0.381977201,
+0.208697319E-1,
-0.302205801,
+0.199525356E-1,
-0.134611130E-2,
+0.403347015,
+0.905563593,
+0.311652541,
+0.488184690,
+0.285637736
336cov(:,:,
0)
= getCov(sample, dim)
338+0.247134715,
+0.170436263,
+0.558982007E-1,
+0.772630274E-1,
-0.380213070E-2
339+0.170436263,
+0.422619432,
-0.162958264,
+0.583852939E-1,
-0.303879976E-1
340+0.558982007E-1,
-0.162958264,
+0.387754440,
+0.423029289E-1,
-0.670315623E-1
341+0.772630274E-1,
+0.583852939E-1,
+0.423029289E-1,
+0.290428221,
-0.646240339E-1
342-0.380213070E-2,
-0.303879976E-1,
-0.670315623E-1,
-0.646240339E-1,
+0.100378618
344 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
345 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
347covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
349+0.247134715,
+0.170436263,
+0.558982380E-1,
+0.772630274E-1,
-0.380212814E-2
350+0.170436263,
+0.422619313,
-0.162958220,
+0.583852828E-1,
-0.303879790E-1
351+0.558982380E-1,
-0.162958220,
+0.387754291,
+0.423028767E-1,
-0.670315847E-1
352+0.772630274E-1,
+0.583852828E-1,
+0.423028767E-1,
+0.290428281,
-0.646240488E-1
353-0.380212814E-2,
-0.303879790E-1,
-0.670315847E-1,
-0.646240488E-1,
+0.100378662
354covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
356+0.247134715,
+0.170436263,
+0.558982380E-1,
+0.772630274E-1,
-0.380212814E-2
357+0.170436263,
+0.422619313,
-0.162958220,
+0.583852828E-1,
-0.303879790E-1
358+0.558982380E-1,
-0.162958220,
+0.387754291,
+0.423028767E-1,
-0.670315847E-1
359+0.772630274E-1,
+0.583852828E-1,
+0.423028767E-1,
+0.290428281,
-0.646240488E-1
360-0.380212814E-2,
-0.303879790E-1,
-0.670315847E-1,
-0.646240488E-1,
+0.100378662
362+0.247134715,
+0.170436263,
+0.558982007E-1,
+0.772630274E-1,
-0.380213070E-2
363+0.170436263,
+0.422619432,
-0.162958264,
+0.583852939E-1,
-0.303879976E-1
364+0.558982007E-1,
-0.162958264,
+0.387754440,
+0.423029289E-1,
-0.670315623E-1
365+0.772630274E-1,
+0.583852939E-1,
+0.423029289E-1,
+0.290428221,
-0.646240339E-1
366-0.380213070E-2,
-0.303879976E-1,
-0.670315623E-1,
-0.646240339E-1,
+0.100378618
371 lb(isam)
= ub(isam
- 1)
+ 1
379call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
380call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
384+0.515513420,
-0.795791507,
-0.320049405,
-0.189167023,
-0.302109838,
+0.883353472,
-0.627524376,
+0.487047076,
-0.615162849E-1,
-0.470565796,
+0.543937206,
-0.271256089
385cov(:,:,
0)
= getCov(sample, dim)
389 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
390 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
392covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
395covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG))
409 lb(isam)
= ub(isam
- 1)
+ 1
417call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
418call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
422-0.197373867,
-0.153590918,
+0.845330715,
-0.190570354,
-0.973837495,
-0.570324659,
-0.456250310,
-0.760878444,
-0.542709589
423-0.406556606,
-0.383972764,
-0.432080865,
+0.657816529,
+0.954630733,
+0.608223319,
+0.843742132,
+0.764674664,
-0.775207281
426+3,
+2,
+9,
+9,
+4,
+1,
+3,
+1,
+8
427cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
429+0.342936337,
-0.156519756
430-0.156519756,
+0.446186632
432 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
433 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
435covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
437+0.342936277,
-0.156519666
438-0.156519666,
+0.446186662
439covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
441+0.342936277,
-0.156519666
442-0.156519666,
+0.446186662
444+0.342936337,
-0.156519756
445-0.156519756,
+0.446186632
450 lb(isam)
= ub(isam
- 1)
+ 1
458call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
459call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
463-0.384662628,
+0.264180303,
+0.100095510,
+0.913428068,
-0.277535796,
-0.534496546,
-0.880666971,
-0.863998413,
-0.291821718,
-0.999720454,
+0.670798779
464-0.105460286,
-0.751714826,
-0.794672489,
+0.556562185,
-0.973960280,
-0.960669875,
+0.726837754,
+0.622811913,
-0.120126009,
+0.690144539,
+0.343933940
465+0.705355644,
+0.634185076,
-0.231496096,
-0.968549728,
+0.924779296,
+0.926171899,
+0.108165264,
+0.832095504,
+0.675756454,
-0.378453255,
+0.681734085E-2
466-0.474578500,
+0.600169659,
-0.998362780,
-0.747842431,
-0.254931927,
+0.345475554,
-0.141309261,
+0.793461800,
+0.535066962,
+0.179133773,
+0.934324384
467+0.664455295,
-0.411200047,
+0.939453840,
-0.361445189,
+0.378231645,
-0.636002302,
-0.525096059,
+0.697192550,
+0.723150253,
-0.696919084,
+0.528118610E-1
470+3,
+4,
+6,
+10,
+4,
+5,
+5,
+2,
+8,
+9,
+2
471cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
473+0.430911988,
-0.385164954E-1,
-0.200423896,
-0.168991342,
+0.507045947E-1
474-0.385164954E-1,
+0.442107826,
-0.275711030,
-0.496802677E-3,
-0.170983419
475-0.200423896,
-0.275711030,
+0.452619851,
+0.235795990,
+0.142566517
476-0.168991342,
-0.496802677E-3,
+0.235795990,
+0.345864385,
-0.490205027E-1
477+0.507045947E-1,
-0.170983419,
+0.142566517,
-0.490205027E-1,
+0.380937964
479 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
480 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
482covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
484+0.430911988,
-0.385165215E-1,
-0.200423926,
-0.168991357,
+0.507046133E-1
485-0.385165215E-1,
+0.442107826,
-0.275711060,
-0.496804714E-3,
-0.170983404
486-0.200423926,
-0.275711060,
+0.452619910,
+0.235795975,
+0.142566472
487-0.168991357,
-0.496804714E-3,
+0.235795975,
+0.345864415,
-0.490205102E-1
488+0.507046133E-1,
-0.170983404,
+0.142566472,
-0.490205102E-1,
+0.380937964
489covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
491+0.430911988,
-0.385165215E-1,
-0.200423926,
-0.168991357,
+0.507046133E-1
492-0.385165215E-1,
+0.442107826,
-0.275711060,
-0.496804714E-3,
-0.170983404
493-0.200423926,
-0.275711060,
+0.452619910,
+0.235795975,
+0.142566472
494-0.168991357,
-0.496804714E-3,
+0.235795975,
+0.345864415,
-0.490205102E-1
495+0.507046133E-1,
-0.170983404,
+0.142566472,
-0.490205102E-1,
+0.380937964
497+0.430911988,
-0.385164954E-1,
-0.200423896,
-0.168991342,
+0.507045947E-1
498-0.385164954E-1,
+0.442107826,
-0.275711030,
-0.496802677E-3,
-0.170983419
499-0.200423896,
-0.275711030,
+0.452619851,
+0.235795990,
+0.142566517
500-0.168991342,
-0.496802677E-3,
+0.235795990,
+0.345864385,
-0.490205027E-1
501+0.507045947E-1,
-0.170983419,
+0.142566517,
-0.490205027E-1,
+0.380937964
506 lb(isam)
= ub(isam
- 1)
+ 1
514call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
515call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
519-0.296853900,
-0.379853964,
+0.870966554,
+0.126925588,
+0.780221224E-1,
-0.622567534
520+0.928081036,
+0.133993626E-1,
-0.363834620,
+0.561187863,
+0.687178612,
+0.110262394
523+2,
+5,
+5,
+6,
+3,
+9
524cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
526+0.268085241,
-0.531103946E-1
527-0.531103946E-1,
+0.146063551
529 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
530 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
532covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
534+0.268085212,
-0.531104133E-1
535-0.531104133E-1,
+0.146063656
536covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
538+0.268085212,
-0.531104133E-1
539-0.531104133E-1,
+0.146063656
541+0.268085241,
-0.531103946E-1
542-0.531103946E-1,
+0.146063551
547 lb(isam)
= ub(isam
- 1)
+ 1
555call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
556call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
560-0.529403329,
+0.820819020,
+0.507654667,
-0.804889202E-2,
+0.704899311,
-0.602318883,
+0.720745087,
+0.362958908,
-0.800908685,
-0.688358665
563+4,
+5,
+10,
+9,
+10,
+10,
+2,
+10,
+2,
+4
564cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
568 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
569 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
571covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
574covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
583 lb(isam)
= ub(isam
- 1)
+ 1
591call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
592call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
596-0.452725649,
-0.384719014,
-0.788704991,
-0.223403573,
+0.546823025,
+0.925979495,
+0.454782724,
+0.745723605,
-0.829160333,
-0.710417032E-1,
+0.296776175
597+0.507106304,
-0.965466380,
+0.688716888,
-0.790937424,
-0.313876271,
-0.975387692,
+0.271956563,
-0.691972375,
+0.256676555,
+0.690219283,
+0.616319060
598-0.337959409,
+0.128120303,
-0.649802446,
+0.649259925,
+0.413934350,
+0.543892026,
-0.140561342,
+0.119246483,
+0.930256128,
-0.431856155,
+0.112541199
601+1,
+8,
+5,
+9,
+3,
+3,
+2,
+4,
+10,
+9,
+10
602cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
604+0.282207668,
-0.686423779E-1,
-0.348747969E-1
605-0.686423779E-1,
+0.462637752,
-0.149972469
606-0.348747969E-1,
-0.149972469,
+0.242351234
608 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
609 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
611covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
613+0.282207549,
-0.686423257E-1,
-0.348748006E-1
614-0.686423257E-1,
+0.462637722,
-0.149972484
615-0.348748006E-1,
-0.149972484,
+0.242351264
616covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
618+0.282207549,
-0.686423257E-1,
-0.348748006E-1
619-0.686423257E-1,
+0.462637722,
-0.149972484
620-0.348748006E-1,
-0.149972484,
+0.242351264
622+0.282207668,
-0.686423779E-1,
-0.348747969E-1
623-0.686423779E-1,
+0.462637752,
-0.149972469
624-0.348747969E-1,
-0.149972469,
+0.242351234
629 lb(isam)
= ub(isam
- 1)
+ 1
637call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
638call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
642+0.380292654,
-0.982146263,
-0.219104171,
-0.135291338,
-0.215328932,
+0.342729926,
+0.449106336,
+0.941773057
643+0.977139354,
-0.985900640,
-0.570864797,
-0.101829410,
+0.989561081E-1,
-0.368391275E-1,
+0.221996188,
+0.205919862
644+0.177469492,
+0.279110074,
+0.181578159,
-0.536322594E-3,
+0.705205321,
+0.276971459,
+0.277076006,
+0.868615270
647+10,
+3,
+5,
+5,
+10,
+3,
+9,
+3
648cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
650+0.197466493,
+0.163436100,
+0.698185479E-2
651+0.163436100,
+0.287774801,
-0.498529291E-2
652+0.698185479E-2,
-0.498529291E-2,
+0.666568279E-1
654 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
655 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
657covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
659+0.197466537,
+0.163436085,
+0.698185898E-2
660+0.163436085,
+0.287774682,
-0.498528825E-2
661+0.698185898E-2,
-0.498528825E-2,
+0.666568205E-1
662covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
664+0.197466537,
+0.163436085,
+0.698185898E-2
665+0.163436085,
+0.287774682,
-0.498528825E-2
666+0.698185898E-2,
-0.498528825E-2,
+0.666568205E-1
668+0.197466493,
+0.163436100,
+0.698185479E-2
669+0.163436100,
+0.287774801,
-0.498529291E-2
670+0.698185479E-2,
-0.498529291E-2,
+0.666568279E-1
675 lb(isam)
= ub(isam
- 1)
+ 1
683call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
684call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
688+0.682350397E-1,
-0.602595210,
+0.156207919,
-0.538625360,
+0.998635054,
-0.239477992,
-0.199894071,
+0.371229649E-1,
-0.600934029E-3,
-0.694355965,
-0.770249844
691+5,
+5,
+4,
+5,
+4,
+2,
+9,
+8,
+10,
+2,
+8
692cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
696 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
697 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
699covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
702covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
711 lb(isam)
= ub(isam
- 1)
+ 1
719call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
720call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
724+0.368094802,
+0.337949157,
+0.402952194,
-0.881861567,
-0.992952585E-1,
+0.905221820,
+0.160231233,
-0.220838785,
-0.728072286
725-0.115045190,
-0.815880299E-1,
+0.759442091,
-0.148627639,
+0.154563546,
+0.899271965,
+0.382173538,
+0.759997129,
+0.426724195
728+8,
+8,
+4,
+9,
+2,
+9,
+6,
+6,
+4
729cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
731+0.345047832,
+0.107825622
732+0.107825622,
+0.175821573
734 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
735 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
737covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
739+0.345047832,
+0.107825622
740+0.107825622,
+0.175821558
741covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
743+0.345047832,
+0.107825622
744+0.107825622,
+0.175821558
746+0.345047832,
+0.107825622
747+0.107825622,
+0.175821573
752 lb(isam)
= ub(isam
- 1)
+ 1
760call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
761call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
765-0.818233490E-1,
-0.633889079,
-0.253235579,
+0.598838449,
+0.257411838
769cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
773 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
774 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
776covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
779covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
788 lb(isam)
= ub(isam
- 1)
+ 1
796call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
797call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
801-0.170882106,
-0.252037168,
-0.446684361E-1,
-0.131878853E-1,
-0.694365621,
-0.913623571,
-0.168187737
802+0.512237549,
+0.184028625,
+0.647021651,
-0.366358876,
+0.116049528,
+0.941858411,
-0.915260077
805+4,
+1,
+8,
+4,
+3,
+10,
+3
806cov(:,:,
0)
= getCov(sample,
2_IK, iweight)
808+0.148460373,
-0.116931848
809-0.116931848,
+0.342437923
811 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
812 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, iweight(lb(isam):ub(isam)))
814covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
816+0.148460373,
-0.116931841
817-0.116931841,
+0.342437953
818covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG))
820+0.148460373,
-0.116931841
821-0.116931841,
+0.342437953
823+0.148460373,
-0.116931848
824-0.116931848,
+0.342437923
834 lb(isam)
= ub(isam
- 1)
+ 1
842call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
843call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
847+0.360957503,
-0.739465952,
-0.761714697,
+0.958295822,
-0.159062266,
+0.927077293,
-0.633233786E-1
848-0.396288633,
-0.423802495,
+0.672270775,
+0.399061918,
-0.297021270,
-0.704549670,
+0.610934496E-1
849-0.169243693,
+0.266524911,
+0.288313627,
+0.390744209E-1,
+0.202130079,
-0.134053469,
+0.866860628
852+1.91824138,
+1.46134067,
+1.41583574,
+1.18216491,
+1.66033697,
+1.68307734,
+1.14657569
853cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
855+0.417825907,
-0.937820226E-1,
-0.108462065
856-0.937820226E-1,
+0.204837471,
+0.591445044E-1
857-0.108462065,
+0.591445044E-1,
+0.942178518E-1
859 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
860 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
862covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
864+0.417825997,
-0.937819853E-1,
-0.108462051
865-0.937819853E-1,
+0.204837486,
+0.591444857E-1
866-0.108462051,
+0.591444857E-1,
+0.942178518E-1
867covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
869+0.417825997,
-0.937819853E-1,
-0.108462051
870-0.937819853E-1,
+0.204837486,
+0.591444857E-1
871-0.108462051,
+0.591444857E-1,
+0.942178518E-1
873+0.417825907,
-0.937820226E-1,
-0.108462065
874-0.937820226E-1,
+0.204837471,
+0.591445044E-1
875-0.108462065,
+0.591445044E-1,
+0.942178518E-1
880 lb(isam)
= ub(isam
- 1)
+ 1
888call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
889call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
893+0.241536617,
+0.966992378E-1,
+0.677903891,
-0.175760627,
+0.161564469,
+0.111931086
896+1.39979434,
+1.51533961,
+1.57947898,
+1.67973089,
+1.95610476,
+1.14286280
897cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
901 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
902 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
904covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
907covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
916 lb(isam)
= ub(isam
- 1)
+ 1
924call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
925call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
929-0.275091529,
+0.294095993,
+0.250872612,
+0.567560077,
+0.501770258,
-0.346876025,
+0.428892970,
+0.861753225,
-0.974786043,
+0.809074640,
-0.530629635,
-0.335462451
930+0.423863530,
+0.972639322E-1,
+0.408923745,
+0.893060088,
+0.618352056,
-0.463502526,
+0.986799121,
-0.983544111,
-0.857935071,
-0.385099649,
+0.371596694,
-0.472127795
931+0.308703542,
-0.321539640E-1,
+0.142372847E-1,
-0.566674709,
+0.925955296,
-0.825324059E-1,
+0.103193045,
-0.779335976,
-0.735169053,
+0.551013947,
+0.461679220,
-0.670419097
932-0.492211699,
-0.188394785E-1,
+0.138161540,
+0.250766397,
-0.181661248,
+0.338420391,
+0.475397587,
-0.192470551E-1,
+0.847903252,
+0.280212522,
+0.498305559,
+0.599375725
933-0.931526661,
+0.703163385,
+0.376219392,
-0.298035264,
-0.919851780,
+0.448298454E-2,
-0.282403111,
-0.406238437,
+0.398800850,
-0.980658293,
+0.758873940,
+0.936976552
936+1.31860042,
+1.82277846,
+1.62192857,
+1.01879179,
+1.42741191,
+1.81362391,
+1.51937342,
+1.38301361,
+1.89903402,
+1.97993922,
+1.17469597,
+1.48918772
937cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
939+0.325481683,
+0.962510332E-1,
+0.952842236E-1,
-0.993375927E-1,
-0.206139103
940+0.962510332E-1,
+0.393903434,
+0.173604637,
-0.770341381E-1,
-0.726224557E-1
941+0.952842236E-1,
+0.173604637,
+0.274354756,
-0.849287510E-1,
-0.175725743
942-0.993375927E-1,
-0.770341381E-1,
-0.849287510E-1,
+0.122914165,
+0.120091721
943-0.206139103,
-0.726224557E-1,
-0.175725743,
+0.120091721,
+0.437241524
945 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
946 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
948covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
950+0.325481683,
+0.962509960E-1,
+0.952842012E-1,
-0.993376151E-1,
-0.206139058
951+0.962509960E-1,
+0.393903196,
+0.173604608,
-0.770341828E-1,
-0.726224408E-1
952+0.952842012E-1,
+0.173604608,
+0.274354726,
-0.849287659E-1,
-0.175725788
953-0.993376151E-1,
-0.770341828E-1,
-0.849287659E-1,
+0.122914158,
+0.120091811
954-0.206139058,
-0.726224408E-1,
-0.175725788,
+0.120091811,
+0.437241614
955covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
957+0.325481683,
+0.962509960E-1,
+0.952842012E-1,
-0.993376151E-1,
-0.206139058
958+0.962509960E-1,
+0.393903196,
+0.173604608,
-0.770341828E-1,
-0.726224408E-1
959+0.952842012E-1,
+0.173604608,
+0.274354726,
-0.849287659E-1,
-0.175725788
960-0.993376151E-1,
-0.770341828E-1,
-0.849287659E-1,
+0.122914158,
+0.120091811
961-0.206139058,
-0.726224408E-1,
-0.175725788,
+0.120091811,
+0.437241614
963+0.325481683,
+0.962510332E-1,
+0.952842236E-1,
-0.993375927E-1,
-0.206139103
964+0.962510332E-1,
+0.393903434,
+0.173604637,
-0.770341381E-1,
-0.726224557E-1
965+0.952842236E-1,
+0.173604637,
+0.274354756,
-0.849287510E-1,
-0.175725743
966-0.993375927E-1,
-0.770341381E-1,
-0.849287510E-1,
+0.122914165,
+0.120091721
967-0.206139103,
-0.726224557E-1,
-0.175725743,
+0.120091721,
+0.437241524
972 lb(isam)
= ub(isam
- 1)
+ 1
980call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
981call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
985-0.501985431,
-0.985388756,
-0.310525894E-1,
+0.490321159,
+0.100234747,
-0.250560045E-1,
+0.955755830,
+0.430842042,
+0.938860416,
+0.954124570
986+0.543662786,
-0.551708102,
+0.823263526,
-0.388222694,
+0.811121941,
-0.486875772E-1,
-0.180379510,
-0.840799212,
-0.377599359,
+0.914110184
987+0.168731928,
+0.126706004,
-0.110031009,
+0.958979726,
-0.949530721,
+0.806047559,
+0.114392042E-1,
+0.262134910,
-0.936820865,
-0.743309379
990+1.04637778,
+1.20557380,
+1.99812746,
+1.88331008,
+1.80594468,
+1.35627460,
+1.56584227,
+1.77880526,
+1.84208536,
+1.50244689
991cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
993+0.329413861,
-0.298964009E-1,
-0.106314972
994-0.298964009E-1,
+0.391267359,
-0.187977403
995-0.106314972,
-0.187977403,
+0.426224440
997 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
998 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1000covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1002+0.329413712,
-0.298963618E-1,
-0.106314927
1003-0.298963618E-1,
+0.391267389,
-0.187977389
1004-0.106314927,
-0.187977389,
+0.426224411
1005covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1007+0.329413712,
-0.298963618E-1,
-0.106314927
1008-0.298963618E-1,
+0.391267389,
-0.187977389
1009-0.106314927,
-0.187977389,
+0.426224411
1011+0.329413861,
-0.298964009E-1,
-0.106314972
1012-0.298964009E-1,
+0.391267359,
-0.187977403
1013-0.106314972,
-0.187977403,
+0.426224440
1018 lb(isam)
= ub(isam
- 1)
+ 1
1026call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1027call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1031+0.979285717,
-0.128044128,
+0.301945210,
+0.528178811,
-0.708107352,
+0.364736080,
-0.414368153,
-0.852784991,
+0.862962365,
-0.163995028E-1
1032-0.637146592,
+0.662392378E-1,
+0.865000248,
+0.953804970,
+0.839768767,
+0.282483578,
+0.898770571,
+0.206402183,
-0.852016926,
-0.859671116
1033+0.935975909,
-0.462703466,
-0.945620418,
+0.438893676,
+0.751902223,
+0.797420502,
+0.574779272,
+0.449635267,
+0.169126511,
+0.612174273
1034+0.505215049,
-0.478659153,
-0.633904934E-1,
+0.417311192E-1,
-0.539137363,
+0.555067658,
+0.886781216,
+0.982346058,
-0.460291386,
+0.332939386
1037+1.63534141,
+1.31228423,
+1.68069339,
+1.72839546,
+1.05993712,
+1.94147635,
+1.68227911,
+1.08322215,
+1.99006557,
+1.35536504
1038cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1040+0.325529844,
-0.185684830,
+0.972874754E-3,
-0.740441978E-1
1041-0.185684830,
+0.508236170,
-0.901889503E-1,
+0.482842401E-1
1042+0.972874754E-3,
-0.901889503E-1,
+0.327678591,
+0.140955418
1043-0.740441978E-1,
+0.482842401E-1,
+0.140955418,
+0.263842344
1045 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1046 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1048covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1050+0.325529903,
-0.185684890,
+0.972872600E-3,
-0.740442127E-1
1051-0.185684890,
+0.508236170,
-0.901889205E-1,
+0.482842550E-1
1052+0.972872600E-3,
-0.901889205E-1,
+0.327678651,
+0.140955433
1053-0.740442127E-1,
+0.482842550E-1,
+0.140955433,
+0.263842285
1054covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1056+0.325529903,
-0.185684890,
+0.972872600E-3,
-0.740442127E-1
1057-0.185684890,
+0.508236170,
-0.901889205E-1,
+0.482842550E-1
1058+0.972872600E-3,
-0.901889205E-1,
+0.327678651,
+0.140955433
1059-0.740442127E-1,
+0.482842550E-1,
+0.140955433,
+0.263842285
1061+0.325529844,
-0.185684830,
+0.972874754E-3,
-0.740441978E-1
1062-0.185684830,
+0.508236170,
-0.901889503E-1,
+0.482842401E-1
1063+0.972874754E-3,
-0.901889503E-1,
+0.327678591,
+0.140955418
1064-0.740441978E-1,
+0.482842401E-1,
+0.140955418,
+0.263842344
1069 lb(isam)
= ub(isam
- 1)
+ 1
1077call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1078call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1082-0.135654211,
+0.261668444,
+0.874604225,
+0.378290892,
-0.826758862,
+0.647881866,
-0.351764679,
+0.443034291,
+0.207801342,
-0.931807876,
+0.808244109
1085+1.69800603,
+1.46708083,
+1.31121588,
+1.08555400,
+1.87509310,
+1.28202796,
+1.35673285,
+1.10928679,
+1.08771443,
+1.47755742,
+1.41808116
1086cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1090 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1091 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1093covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1096covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1105 lb(isam)
= ub(isam
- 1)
+ 1
1113call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1114call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1118+0.147144079,
+0.823786020,
+0.364029646,
-0.510857105,
+0.781095386,
-0.190358996,
+0.902266502,
+0.613397360E-1
1119-0.638231754,
+0.772516131,
+0.289416313E-2,
-0.357980251,
-0.350979209,
+0.807772756,
+0.941998839,
-0.485590339
1120+0.545637965,
+0.239905119E-1,
-0.612406373,
+0.906431675,
+0.411424518,
-0.803205967,
-0.479002476,
+0.565589666E-1
1123+1.79489017,
+1.93938768,
+1.36774540,
+1.25689793,
+1.80379796,
+1.04966843,
+1.05916047,
+1.68959737
1124cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1126+0.212614149,
+0.116323739,
-0.574993715E-1
1127+0.116323739,
+0.358453304,
-0.207802683
1128-0.574993715E-1,
-0.207802683,
+0.272262603
1130 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1131 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1133covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1135+0.212614119,
+0.116323702,
-0.574993566E-1
1136+0.116323702,
+0.358453274,
-0.207802549
1137-0.574993566E-1,
-0.207802549,
+0.272262514
1138covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1140+0.212614119,
+0.116323702,
-0.574993566E-1
1141+0.116323702,
+0.358453274,
-0.207802549
1142-0.574993566E-1,
-0.207802549,
+0.272262514
1144+0.212614149,
+0.116323739,
-0.574993715E-1
1145+0.116323739,
+0.358453304,
-0.207802683
1146-0.574993715E-1,
-0.207802683,
+0.272262603
1151 lb(isam)
= ub(isam
- 1)
+ 1
1159call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1160call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1164+0.212894559,
+0.683618546,
+0.593861341E-1,
+0.463629127,
-0.756304860,
+0.343446732E-1,
-0.154616714,
-0.733219385,
-0.365341902,
-0.107770324,
+0.412575841,
-0.650551796,
-0.925751925E-1,
-0.619743824
1167+1.18728578,
+1.27279472,
+1.67459893,
+1.82329762,
+1.05707896,
+1.88751626,
+1.92389882,
+1.33283567,
+1.66768956,
+1.60123825,
+1.74192226,
+1.05816483,
+1.16486549,
+1.82178879
1168cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1172 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1173 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1175covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1178covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1187 lb(isam)
= ub(isam
- 1)
+ 1
1195call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1196call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1200-0.259560227,
-0.778334379,
-0.363828182,
+0.985556841,
-0.169667959,
+0.471833587,
+0.505243897,
-0.814437866E-1,
+0.252072573,
-0.876434326,
-0.777854919E-1,
+0.775524139,
-0.801629424
1201+0.572081327,
+0.174535513E-1,
-0.198926687,
-0.139424443,
-0.792771578,
-0.146719813,
-0.429209948,
+0.488437533,
-0.218845963,
-0.342070699,
-0.873639226,
-0.917390585,
+0.359002471
1202-0.964350581,
-0.381152272,
-0.543881416,
-0.717974782,
+0.318305492E-1,
-0.972483277,
-0.456343651,
-0.543370247,
+0.537569046,
-0.984147668,
+0.533637881,
+0.113233566,
-0.331536055
1203-0.712396026,
+0.138616800,
+0.996387005,
-0.660705447,
-0.448454618E-1,
-0.787701488,
-0.582031369,
+0.786630154,
-0.568483591,
-0.445982814,
-0.465743899,
-0.519416332E-1,
+0.893786788
1206+1.57856393,
+1.90027249,
+1.68442070,
+1.19526815,
+1.54564667,
+1.33161330,
+1.20656312,
+1.07392216,
+1.45212913,
+1.54400754,
+1.10952568,
+1.75179410,
+1.88956857
1207cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1209+0.348148704,
-0.111350976,
+0.587605573E-1,
-0.172193736
1210-0.111350976,
+0.217035726,
-0.124739900,
+0.701254010E-1
1211+0.587605573E-1,
-0.124739900,
+0.239982545,
+0.210398231E-1
1212-0.172193736,
+0.701254010E-1,
+0.210398231E-1,
+0.382413417
1214 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1215 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1217covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1219+0.348148704,
-0.111350961,
+0.587605797E-1,
-0.172193736
1220-0.111350961,
+0.217035696,
-0.124739885,
+0.701253414E-1
1221+0.587605797E-1,
-0.124739885,
+0.239982516,
+0.210398026E-1
1222-0.172193736,
+0.701253414E-1,
+0.210398026E-1,
+0.382413477
1223covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1225+0.348148704,
-0.111350961,
+0.587605797E-1,
-0.172193736
1226-0.111350961,
+0.217035696,
-0.124739885,
+0.701253414E-1
1227+0.587605797E-1,
-0.124739885,
+0.239982516,
+0.210398026E-1
1228-0.172193736,
+0.701253414E-1,
+0.210398026E-1,
+0.382413477
1230+0.348148704,
-0.111350976,
+0.587605573E-1,
-0.172193736
1231-0.111350976,
+0.217035726,
-0.124739900,
+0.701254010E-1
1232+0.587605573E-1,
-0.124739900,
+0.239982545,
+0.210398231E-1
1233-0.172193736,
+0.701254010E-1,
+0.210398231E-1,
+0.382413417
1238 lb(isam)
= ub(isam
- 1)
+ 1
1246call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1247call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1251-0.544495344,
-0.987998962,
+0.473976135E-1,
-0.510042906E-1,
+0.658811331,
-0.758180976,
-0.235369205,
+0.116745710,
+0.692964196,
+0.813944101,
-0.598946929
1252-0.353805780,
-0.989767671,
+0.469312668,
-0.466500640,
-0.830650330E-3,
-0.202172995E-1,
-0.273268580,
+0.621523976,
+0.156655550,
-0.832636833,
+0.976485014
1255+1.41977108,
+1.02416992,
+1.77285099,
+1.38984084,
+1.27823973,
+1.80695510,
+1.75225282,
+1.11161232,
+1.48210871,
+1.67843461,
+1.45311165
1256cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1258+0.334554911,
-0.196566731E-1
1259-0.196566731E-1,
+0.307370782
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)))
1264covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1266+0.334554851,
-0.196566507E-1
1267-0.196566507E-1,
+0.307370812
1268covMerged
= getCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG))
1270+0.334554851,
-0.196566507E-1
1271-0.196566507E-1,
+0.307370812
1273+0.334554911,
-0.196566731E-1
1274-0.196566731E-1,
+0.307370782