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(
"call setCovMerged(covMerged, cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG), uppDia)")
74 call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
77 call disp%show(
"call setCovMerged(covMerged, cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG), lowDia)")
78 call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
81 call disp%show(
"call setCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(ub(1), TKG) / real(ub(2), TKG), uppDia)")
82 call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
85 call disp%show(
"cov(:,:,0) ! reference")
92 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
93 call disp%show(
"!Compute the biased merged covariance of a frequency weighted multivariate sample.")
94 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
98 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), covMerged(:,:)
99 real(TKG),
allocatable :: sample(:,:)
102 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
105 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
106 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
109 lb(isam)
= ub(isam
- 1)
+ 1
116 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
118 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
119 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
120 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
121 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
122 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
124 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
128 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
132 call disp%show(
"cov(:,:,0) = getCov(sample, dim, iweight)")
133 cov(:,:,
0)
= getCov(sample, dim, iweight)
134 call disp%show(
"cov(:,:,0) ! reference")
137 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))")
138 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))")
141 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
142 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
144 call disp%show(
"call setCovMerged(covMerged, cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), uppDia)")
145 call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
148 call disp%show(
"call setCovMerged(covMerged, cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), lowDia)")
149 call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
152 call disp%show(
"call setCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), uppDia)")
153 call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
156 call disp%show(
"cov(:,:,0) ! reference")
163 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
164 call disp%show(
"!Compute the biased merged covariance of a reliability weighted multivariate sample.")
165 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
169 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), covMerged(:,:)
170 real(TKG),
allocatable :: sample(:,:)
171 real(TKG),
allocatable :: rweight(:)
174 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
177 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
178 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
181 lb(isam)
= ub(isam
- 1)
+ 1
188 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
190 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
191 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
192 call disp%show(
"call setRebound(mean, [1_IK, 1_IK], [ndim, nsam])")
193 call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
194 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
196 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
200 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
204 call disp%show(
"cov(:,:,0) = getCov(sample, dim, rweight)")
205 cov(:,:,
0)
= getCov(sample, dim, rweight)
206 call disp%show(
"cov(:,:,0) ! reference")
209 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))")
210 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))")
213 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
214 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
216 call disp%show(
"call setCovMerged(covMerged, cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), uppDia)")
217 call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
220 call disp%show(
"call setCovMerged(covMerged, cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), lowDia)")
221 call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
224 call disp%show(
"call setCovMerged(cov(:,:,2), cov(:,:,1), mean(:,1) - mean(:,2), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), lowDia)")
225 call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
228 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.730689526,
+0.174534321E-1,
+0.183695674,
+0.683780670,
-0.559766293,
+0.301911831,
+0.747577667,
+0.578095913E-1
23+0.477809548,
-0.452203751,
+0.482423067,
+0.439996362,
-0.726205826,
+0.746573091,
-0.845496297,
-0.195052624E-1
24-0.131714344E-1,
-0.830839396,
-0.650135636,
+0.553827286,
-0.140310526E-1,
+0.615285873,
+0.131123543,
-0.559346437
25cov(:,:,
0)
= getCov(sample, dim)
27+0.242581904,
+0.276127458E-2,
+0.726257265E-1
28+0.276127458E-2,
+0.332221717,
+0.848705173E-1
29+0.726257265E-1,
+0.848705173E-1,
+0.256889105
31 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
32 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
34call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
36+0.242581859,
+0.276129507E-2,
+0.726257265E-1
37+0.700649232E-44,
+0.332221687,
+0.848705173E-1
38+0.00000000,
+0.229218613E-17,
+0.256889105
39call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
41+0.242581859,
+0.276129507E-2,
+0.726257265E-1
42+0.276129507E-2,
+0.332221687,
+0.848705173E-1
43+0.726257265E-1,
+0.848705173E-1,
+0.256889105
44call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
46+0.242581859,
+0.276129507E-2,
+0.726257265E-1
47+0.285798013E-1,
+0.332221687,
+0.848705173E-1
48+0.589051247E-1,
+0.101532042,
+0.256889105
50+0.242581904,
+0.276127458E-2,
+0.726257265E-1
51+0.276127458E-2,
+0.332221717,
+0.848705173E-1
52+0.726257265E-1,
+0.848705173E-1,
+0.256889105
57 lb(isam)
= ub(isam
- 1)
+ 1
65call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
66call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
70-0.946599126,
-0.990906954,
+0.977254272,
-0.387032866
71+0.384290814,
-0.765278220,
-0.679953933,
+0.679383636
72cov(:,:,
0)
= getCov(sample, dim)
74+0.632243276,
-0.165349573
75-0.165349573,
+0.405208290
77 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
78 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
80call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
82+0.632243216,
-0.165349558
83+0.276129507E-2,
+0.405208290
84call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
86+0.632243216,
-0.165349558
87-0.165349558,
+0.405208290
88call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
90+0.632243216,
-0.165349558
91-0.463631690,
+0.405208290
93+0.632243276,
-0.165349573
94-0.165349573,
+0.405208290
99 lb(isam)
= ub(isam
- 1)
+ 1
107call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
108call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
112-0.462055445,
-0.787469745,
-0.821069121,
+0.628360391,
-0.953035951,
-0.536895037,
+0.833436251E-1,
+0.980564356E-1,
-0.853996992,
-0.755120277
113cov(:,:,
0)
= getCov(sample, dim)
117 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
118 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
120call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
123call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
126call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
135 lb(isam)
= ub(isam
- 1)
+ 1
143call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
144call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
148-0.379431725,
-0.836435914,
-0.912897825,
+0.529910326,
+0.443402410,
-0.593599200,
-0.311179161,
+0.589435935,
+0.177296877
149-0.901885629,
+0.936819196,
-0.620117545,
-0.370423198,
+0.496287465,
+0.361477137E-1,
+0.895470738,
+0.405763268,
+0.410896540
150cov(:,:,
0)
= getCov(sample, dim)
152+0.310727090,
+0.384052582E-1
153+0.384052582E-1,
+0.379015505
155 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
156 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
158call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
160+0.310727060,
+0.384052694E-1
161+0.306926403E-40,
+0.379015714
162call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
164+0.310727060,
+0.384052694E-1
165+0.384052694E-1,
+0.379015714
166call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
168+0.310727060,
+0.384052694E-1
169-0.375374369E-1,
+0.379015714
171+0.310727090,
+0.384052582E-1
172+0.384052582E-1,
+0.379015505
177 lb(isam)
= ub(isam
- 1)
+ 1
185call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
186call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
190+0.374810696E-1,
+0.411958218,
-0.492504120,
-0.254053473,
+0.359842300,
-0.803278089,
-0.631706715,
+0.978754759E-1
191+0.911535025,
-0.708054304,
-0.819551587,
-0.214745879,
+0.444084287,
-0.847451687,
+0.599342346,
+0.779874086
192+0.434346080,
+0.815168619,
-0.727621913,
-0.450674295E-1,
+0.637009263,
-0.814595938,
+0.837425828,
-0.515295625
193cov(:,:,
0)
= getCov(sample, dim)
195+0.182323232,
+0.952538848E-1,
+0.140626520
196+0.952538848E-1,
+0.491273522,
+0.185841322
197+0.140626520,
+0.185841322,
+0.421564698
199 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
200 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
202call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
204+0.182323217,
+0.952538773E-1,
+0.140626505
205+0.384052694E-1,
+0.491273493,
+0.185841322
206+0.00000000,
-0.536895037,
+0.421564728
207call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
209+0.182323217,
+0.952538773E-1,
+0.140626505
210+0.952538773E-1,
+0.491273493,
+0.185841322
211+0.140626505,
+0.185841322,
+0.421564728
212call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
214+0.182323217,
+0.952538773E-1,
+0.140626505
215+0.156050190,
+0.491273493,
+0.185841322
216+0.680148974E-1,
+0.213558272,
+0.421564728
218+0.182323232,
+0.952538848E-1,
+0.140626520
219+0.952538848E-1,
+0.491273522,
+0.185841322
220+0.140626520,
+0.185841322,
+0.421564698
225 lb(isam)
= ub(isam
- 1)
+ 1
233call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
234call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
238+0.364428043,
-0.515398860,
+0.652534008,
-0.669225335,
-0.665750146,
-0.984044671,
-0.577415824
239cov(:,:,
0)
= getCov(sample, dim)
243 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
244 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
246call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
249call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
252call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
261 lb(isam)
= ub(isam
- 1)
+ 1
269call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
270call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
274-0.346350789,
+0.476465344,
-0.276171088,
+0.197638154,
+0.397109389
275cov(:,:,
0)
= getCov(sample, dim)
279 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
280 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
282call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
285call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
288call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
297 lb(isam)
= ub(isam
- 1)
+ 1
305call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
306call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
310-0.859471560E-1,
-0.539392233E-1,
+0.373067975,
-0.930675507,
+0.681902170,
+0.910912275,
+0.373710155,
-0.789715886,
+0.948826909,
+0.799502254,
+0.781024218,
+0.891802073,
+0.666453838E-1,
-0.902711034
311-0.828424335,
+0.524303436,
-0.131525040,
-0.627688050,
-0.260316730,
+0.430219531,
-0.911553264,
-0.646403670,
-0.803606153,
-0.187167764,
-0.556549788,
-0.831809998,
+0.767171383,
+0.171775818E-1
312+0.175639510,
+0.520481825,
-0.715081453,
-0.355767012,
+0.539922595,
-0.254513741,
+0.178725600,
-0.904102325E-1,
+0.352688670,
+0.256346583,
+0.454981327,
+0.583103538,
-0.531791806,
-0.446020246
313+0.115087390,
-0.164973497,
+0.244309187,
-0.192922950,
+0.712866783E-1,
-0.646193981,
-0.452042103,
-0.457638502E-1,
-0.818601370,
-0.877954125,
-0.377118587,
-0.947598934,
+0.339508891,
-0.437502861
314cov(:,:,
0)
= getCov(sample, dim)
316+0.440503895,
-0.146509930E-1,
+0.137967303,
-0.120467961
317-0.146509930E-1,
+0.284001112,
-0.858512968E-1,
+0.653083324E-1
318+0.137967303,
-0.858512968E-1,
+0.180216253,
-0.763353407E-1
319-0.120467961,
+0.653083324E-1,
-0.763353407E-1,
+0.164279029
321 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
322 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
324call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
326+0.440503955,
-0.146510061E-1,
+0.137967288,
-0.120467931
327+0.00000000,
+0.284001112,
-0.858513117E-1,
+0.653083324E-1
328+0.00000000,
+0.00000000,
+0.180216268,
-0.763353258E-1
329+0.00000000,
+0.00000000,
+0.00000000,
+0.164279044
330call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
332+0.440503955,
-0.146510061E-1,
+0.137967288,
-0.120467931
333-0.146510061E-1,
+0.284001112,
-0.858513117E-1,
+0.653083324E-1
334+0.137967288,
-0.858513117E-1,
+0.180216268,
-0.763353258E-1
335-0.120467931,
+0.653083324E-1,
-0.763353258E-1,
+0.164279044
336call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
338+0.440503955,
-0.146510061E-1,
+0.137967288,
-0.120467931
339-0.138241649,
+0.284001112,
-0.858513117E-1,
+0.653083324E-1
340+0.243580267,
-0.177357689,
+0.180216268,
-0.763353258E-1
341-0.193252400,
+0.155635208,
-0.136775777,
+0.164279044
343+0.440503895,
-0.146509930E-1,
+0.137967303,
-0.120467961
344-0.146509930E-1,
+0.284001112,
-0.858512968E-1,
+0.653083324E-1
345+0.137967303,
-0.858512968E-1,
+0.180216253,
-0.763353407E-1
346-0.120467961,
+0.653083324E-1,
-0.763353407E-1,
+0.164279029
351 lb(isam)
= ub(isam
- 1)
+ 1
359call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
360call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
364-0.658514977,
+0.568775535,
+0.121245623,
+0.843169928,
-0.488507986,
+0.115345001,
-0.210928202
365+0.851544261,
+0.754020572,
+0.840786695E-1,
+0.820228815,
-0.905382276,
+0.559113622,
-0.138116002
366cov(:,:,
0)
= getCov(sample, dim)
368+0.252450913,
+0.145961121
369+0.145961121,
+0.362690657
371 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
372 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
374call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
376+0.252450913,
+0.145961136
377-0.146510061E-1,
+0.362690657
378call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
380+0.252450913,
+0.145961136
381+0.145961136,
+0.362690657
382call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
384+0.252450913,
+0.145961136
385+0.568718351E-1,
+0.362690657
387+0.252450913,
+0.145961121
388+0.145961121,
+0.362690657
393 lb(isam)
= ub(isam
- 1)
+ 1
401call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
402call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
406+0.233716965E-1,
-0.427911639,
-0.718568563E-1,
+0.105924726,
+0.804329872,
-0.567386031,
+0.342302322
407+0.142705083,
+0.978607059,
+0.813134789,
-0.265102863,
-0.338086009,
+0.600153804,
-0.785182476
408cov(:,:,
0)
= getCov(sample, dim)
410+0.182836547,
-0.202477112
411-0.202477112,
+0.373258680
413 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
414 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
416call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
418+0.182836577,
-0.202477127
419+0.145961136,
+0.373258680
420call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
422+0.182836577,
-0.202477127
423-0.202477127,
+0.373258680
424call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
426+0.182836577,
-0.202477127
427-0.194157600,
+0.373258680
429+0.182836547,
-0.202477112
430-0.202477112,
+0.373258680
440 lb(isam)
= ub(isam
- 1)
+ 1
448call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
449call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
453+0.663856387,
+0.370827556,
-0.847152233,
-0.812152386,
-0.829524755,
+0.984963179,
+0.806908131,
-0.675281048,
-0.852650881,
-0.291886568,
-0.352904320,
+0.327542543
454+0.928904533,
-0.202757239,
-0.771390915,
+0.441682816,
-0.559954405,
-0.913024783,
+0.421713829,
-0.983847380E-1,
+0.350212097,
-0.725432396,
-0.642382145,
-0.732728243E-1
457+7,
+6,
+7,
+9,
+2,
+1,
+7,
+10,
+8,
+4,
+7,
+2
458cov(:,:,
0)
= getCov(sample, dim, iweight)
460+0.413921237,
+0.114881843
461+0.114881843,
+0.299755603
463 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
464 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
466call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
468+0.413921297,
+0.114881814
469+0.306926403E-40,
+0.299755484
470call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
472+0.413921297,
+0.114881814
473+0.114881814,
+0.299755484
474call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
476+0.413921297,
+0.114881814
477+0.309136230E-1,
+0.299755484
479+0.413921237,
+0.114881843
480+0.114881843,
+0.299755603
485 lb(isam)
= ub(isam
- 1)
+ 1
493call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
494call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
498-0.407470942,
-0.713191986,
+0.470277667,
+0.375116706,
+0.545417070E-1,
+0.243642926,
+0.625978470,
-0.471876025,
-0.357685685
501+5,
+10,
+6,
+2,
+10,
+4,
+1,
+10,
+2
502cov(:,:,
0)
= getCov(sample, dim, iweight)
506 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
507 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
509call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
512call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
515call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
524 lb(isam)
= ub(isam
- 1)
+ 1
532call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
533call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
537-0.700062752,
+0.183719397,
-0.225917935,
-0.130042434,
+0.173640966,
-0.456657410,
-0.664977551,
+0.388418317
538-0.832275152,
-0.783651352,
+0.533238053,
+0.136787891E-1,
+0.270397305,
+0.264082909,
-0.162866950,
+0.855092525
541+4,
+1,
+8,
+2,
+2,
+10,
+8,
+4
542cov(:,:,
0)
= getCov(sample, dim, iweight)
544+0.118133374,
+0.104674026
545+0.104674026,
+0.227452010
547 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
548 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
550call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
552+0.118133366,
+0.104673967
553+0.306926403E-40,
+0.227452025
554call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
556+0.118133366,
+0.104673967
557+0.104673967,
+0.227452025
558call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
560+0.118133366,
+0.104673967
561+0.106482521,
+0.227452025
563+0.118133374,
+0.104674026
564+0.104674026,
+0.227452010
569 lb(isam)
= ub(isam
- 1)
+ 1
577call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
578call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
582-0.504294038,
+0.469903946E-1,
+0.690548897,
+0.807892323,
-0.184420228,
+0.149602890E-1,
-0.276815057,
+0.196006060
585+1,
+4,
+7,
+6,
+6,
+6,
+6,
+10
586cov(:,:,
0)
= getCov(sample, dim, iweight)
590 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
591 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
593call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
596call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
599call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
608 lb(isam)
= ub(isam
- 1)
+ 1
616call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
617call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
621-0.182912111,
+0.121653080E-2,
-0.447187424,
-0.544398546,
+0.342447400,
+0.698184967E-2,
+0.182927251
622+0.770518064,
-0.689107895,
-0.174644470,
+0.484755874,
+0.726112604,
-0.731326580,
+0.174746037
623-0.981148481,
-0.731261611,
-0.839406490,
-0.645702958,
-0.536799312,
+0.631785750,
-0.418595314
626+4,
+3,
+9,
+5,
+7,
+1,
+2
627cov(:,:,
0)
= getCov(sample, dim, iweight)
629+0.115784690,
+0.520413034E-1,
+0.402940586E-1
630+0.520413034E-1,
+0.262709767,
-0.217508189E-1
631+0.402940586E-1,
-0.217508189E-1,
+0.840230733E-1
633 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
634 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
636call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
638+0.115784682,
+0.520412959E-1,
+0.402940437E-1
639+0.306926403E-40,
+0.262709796,
-0.217507929E-1
640+0.00000000,
+0.840779079E-44,
+0.840230435E-1
641call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
643+0.115784682,
+0.520412959E-1,
+0.402940437E-1
644+0.520412959E-1,
+0.262709796,
-0.217507929E-1
645+0.402940437E-1,
-0.217507929E-1,
+0.840230435E-1
646call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
648+0.115784682,
+0.520412959E-1,
+0.402940437E-1
649+0.497264005E-1,
+0.262709796,
-0.217507929E-1
650-0.337774232E-1,
-0.147378668,
+0.840230435E-1
652+0.115784690,
+0.520413034E-1,
+0.402940586E-1
653+0.520413034E-1,
+0.262709767,
-0.217508189E-1
654+0.402940586E-1,
-0.217508189E-1,
+0.840230733E-1
659 lb(isam)
= ub(isam
- 1)
+ 1
667call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
668call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
672-0.802071691,
+0.828577161,
-0.238057494,
-0.982773066,
-0.387693882,
+0.437667966,
+0.105431795,
+0.973488092,
+0.410717010,
+0.643182755,
-0.126897931
675+7,
+1,
+4,
+1,
+8,
+5,
+2,
+10,
+7,
+9,
+6
676cov(:,:,
0)
= getCov(sample, dim, iweight)
680 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
681 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
683call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
686call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
689call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
698 lb(isam)
= ub(isam
- 1)
+ 1
706call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
707call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
711+0.869203568,
-0.175490379,
+0.161351442,
+0.192652464,
+0.454697609,
-0.140027046,
+0.888376713
714+3,
+3,
+3,
+7,
+7,
+3,
+4
715cov(:,:,
0)
= getCov(sample, dim, iweight)
719 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
720 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
722call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
725call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
728call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
737 lb(isam)
= ub(isam
- 1)
+ 1
745call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
746call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
750-0.312891006E-1,
-0.643205047,
+0.785211802,
+0.927758455,
-0.985159755,
+0.881144166,
+0.162091374,
+0.515345335E-1,
+0.675864816,
+0.680277228
751-0.465568423,
-0.680073023,
+0.890941620,
+0.264859319,
+0.829611778,
+0.123730898,
+0.706754684,
-0.701587677,
+0.966189981,
-0.619274139
754+7,
+9,
+7,
+10,
+9,
+4,
+5,
+5,
+9,
+3
755cov(:,:,
0)
= getCov(sample, dim, iweight)
757+0.467610121,
+0.100294508
758+0.100294508,
+0.436529011
760 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
761 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
763call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
765+0.467610151,
+0.100294486
766+0.306926403E-40,
+0.436528951
767call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
769+0.467610151,
+0.100294486
770+0.100294486,
+0.436528951
771call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
773+0.467610151,
+0.100294486
774+0.971634537E-1,
+0.436528951
776+0.467610121,
+0.100294508
777+0.100294508,
+0.436529011
782 lb(isam)
= ub(isam
- 1)
+ 1
790call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
791call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
795+0.668229699,
-0.351782441,
-0.265810490E-1,
-0.441505671,
-0.972831249,
+0.875248909E-1,
+0.550184011,
+0.733505964
796+0.907635450,
-0.304288387,
-0.650396466,
-0.265776157,
-0.367276669E-1,
-0.729820967,
-0.492727995,
+0.612545013E-1
799+6,
+6,
+1,
+7,
+8,
+4,
+1,
+7
800cov(:,:,
0)
= getCov(sample, dim, iweight)
802+0.405942261,
+0.125824735
803+0.125824735,
+0.217718601
805 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
806 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
808call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
810+0.405942172,
+0.125824720
811+0.100294486,
+0.217718601
812call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
814+0.405942172,
+0.125824720
815+0.125824720,
+0.217718601
816call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
818+0.405942172,
+0.125824720
819-0.332863000E-2,
+0.217718601
821+0.405942261,
+0.125824735
822+0.125824735,
+0.217718601
827 lb(isam)
= ub(isam
- 1)
+ 1
835call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
836call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
840+0.388753414E-1,
-0.322981477,
-0.633426547,
-0.999809146,
-0.717864394
841-0.377255678,
-0.483701944,
+0.146544456,
-0.766735077E-1,
-0.565060496
845cov(:,:,
0)
= getCov(sample, dim, iweight)
847+0.124290712,
-0.338921919E-1
848-0.338921919E-1,
+0.708793178E-1
850 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
851 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
853call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
855+0.124290705,
-0.338921882E-1
856+0.125824720,
+0.708793178E-1
857call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
859+0.124290705,
-0.338921882E-1
860-0.338921882E-1,
+0.708793178E-1
861call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
863+0.124290705,
-0.338921882E-1
864-0.343291759E-1,
+0.708793178E-1
866+0.124290712,
-0.338921919E-1
867-0.338921919E-1,
+0.708793178E-1
877 lb(isam)
= ub(isam
- 1)
+ 1
885call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
886call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
890+0.870274186,
-0.991506100,
+0.368181825,
-0.547364950,
+0.999474883,
-0.154445529
891+0.511259794,
+0.408414721,
+0.820556760,
+0.940234423,
+0.733433366,
-0.939947009
894+1.03623295,
+1.26095629,
+1.63207519,
+1.67217135,
+1.34095991,
+1.18148041
895cov(:,:,
0)
= getCov(sample, dim, rweight)
897+0.502747834,
+0.620557740E-1
898+0.620557740E-1,
+0.371575356
900 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
901 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
903call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
905+0.502747834,
+0.620558113E-1
906+0.306926403E-40,
+0.371575415
907call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
909+0.502747834,
+0.620558113E-1
910+0.620558113E-1,
+0.371575415
911call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
913+0.502747834,
+0.480807334
914+0.620558113E-1,
+0.371575415
916+0.502747834,
+0.620557740E-1
917+0.620557740E-1,
+0.371575356
922 lb(isam)
= ub(isam
- 1)
+ 1
930call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
931call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
935+0.691048265,
+0.740772843,
-0.840807080,
-0.682481170,
-0.177993774,
+0.467064857
938+1.15253329,
+1.55860722,
+1.55412745,
+1.76245379,
+1.66519260,
+1.57567096
939cov(:,:,
0)
= getCov(sample, dim, rweight)
943 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
944 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
946call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
949call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
952call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
961 lb(isam)
= ub(isam
- 1)
+ 1
969call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
970call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
974+0.398148179,
-0.527587295,
-0.583539009,
+0.976370931,
+0.653176665,
-0.795626044,
+0.250953317
975-0.837774038,
-0.343427896,
+0.202022076,
-0.647808313,
+0.799802780,
-0.298324823,
+0.417882442
978+1.27997518,
+1.20901752,
+1.17209601,
+1.52410865,
+1.29827309,
+1.04914808,
+1.79405916
979cov(:,:,
0)
= getCov(sample, dim, rweight)
981+0.384199560,
-0.108779976E-1
982-0.108779976E-1,
+0.311271399
984 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
985 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
987call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
989+0.384199560,
-0.108780172E-1
990+0.306926403E-40,
+0.311271489
991call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
993+0.384199560,
-0.108780172E-1
994-0.108780172E-1,
+0.311271489
995call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
997+0.384199560,
+0.174528942
998-0.108780172E-1,
+0.311271489
1000+0.384199560,
-0.108779976E-1
1001-0.108779976E-1,
+0.311271399
1006 lb(isam)
= ub(isam
- 1)
+ 1
1014call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1015call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1019-0.960645795,
-0.627111554,
+0.796313286,
+0.771232009,
-0.603791952,
+0.174230933,
+0.487685919,
-0.849461555E-2
1022+1.80196762,
+1.04828525,
+1.67848074,
+1.49191523,
+1.41728687,
+1.75760233,
+1.41154253,
+1.15368915
1023cov(:,:,
0)
= getCov(sample, dim, rweight)
1027 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1028 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1030call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1033call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1036call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1045 lb(isam)
= ub(isam
- 1)
+ 1
1053call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1054call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1058+0.974749804,
-0.615403652,
-0.426107287,
-0.308581352,
-0.417941689,
+0.298369527,
-0.159184217,
+0.648955822,
+0.950057149,
+0.478833795,
+0.916819692
1061+1.86715579,
+1.94725442,
+1.97527015,
+1.89408898,
+1.03421712,
+1.84572613,
+1.73083353,
+1.41008782,
+1.79386687,
+1.02597880,
+1.48380542
1062cov(:,:,
0)
= getCov(sample, dim, rweight)
1066 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1067 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1069call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1072call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1075call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1084 lb(isam)
= ub(isam
- 1)
+ 1
1092call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1093call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1097+0.418940544,
-0.238521576,
+0.577281237,
-0.907191515,
+0.229027987,
-0.419914961,
+0.664691925E-1,
+0.370019555,
+0.940028429
1098+0.336863637,
-0.952550173E-1,
+0.309643149,
-0.930635571,
+0.242483974,
+0.281979561,
+0.719059706E-1,
+0.711679459,
-0.552076578
1101+1.27995753,
+1.51105499,
+1.35377562,
+1.80416691,
+1.46107578,
+1.49078465,
+1.66627026,
+1.00531673,
+1.08269715
1102cov(:,:,
0)
= getCov(sample, dim, rweight)
1104+0.289658338,
+0.125951514
1105+0.125951514,
+0.229427144
1107 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1108 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1110call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1112+0.289658368,
+0.125951514
1113+0.306926403E-40,
+0.229427144
1114call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1116+0.289658368,
+0.125951514
1117+0.125951514,
+0.229427144
1118call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1120+0.289658368,
-0.109702945
1121+0.125951514,
+0.229427144
1123+0.289658338,
+0.125951514
1124+0.125951514,
+0.229427144
1129 lb(isam)
= ub(isam
- 1)
+ 1
1137call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1138call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1142-0.116342306E-1,
+0.122706294,
-0.108284235,
-0.845807314,
+0.435151815,
+0.247197628,
+0.976198673,
+0.488863587,
+0.270626068,
-0.777511120,
-0.244551063,
+0.132064700,
+0.472542763
1143-0.246724486,
-0.301316977,
+0.627304196,
-0.510698676,
-0.692843914,
-0.265240669E-1,
+0.625134349,
+0.279241681,
-0.310252547,
+0.229469538,
+0.637716413,
-0.969700933,
-0.936755300
1144+0.944617867,
+0.142844796,
+0.239712238,
+0.274327993E-1,
-0.611779690,
-0.647225857,
-0.930213690,
+0.360064149,
-0.620508194,
-0.534427166E-2,
-0.294285178,
+0.717393756,
-0.443050265
1145+0.592299819,
+0.553686380,
-0.252180696,
+0.519427180,
+0.879992962,
-0.513559222,
-0.963591814,
-0.139496446,
-0.697885871,
+0.527045250,
+0.909265280E-1,
+0.525159001,
+0.547252655
1148+1.56536078,
+1.55579054,
+1.60155630,
+1.32733440,
+1.66204166,
+1.08505321,
+1.27863312,
+1.58737767,
+1.92929327,
+1.82076478,
+1.32020164,
+1.03257108,
+1.37034059
1149cov(:,:,
0)
= getCov(sample, dim, rweight)
1151+0.235387415,
-0.194309521E-1,
-0.962856039E-1,
-0.113151498
1152-0.194309521E-1,
+0.282119334,
-0.147216897E-1,
-0.165110275
1153-0.962856039E-1,
-0.147216897E-1,
+0.283279181,
+0.132138297
1154-0.113151498,
-0.165110275,
+0.132138297,
+0.318178713
1156 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1157 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1159call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1161+0.235387519,
-0.194309540E-1,
-0.962856039E-1,
-0.113151543
1162+0.125951514,
+0.282119364,
-0.147216981E-1,
-0.165110230
1163+0.00000000,
-0.858513117E-1,
+0.283279151,
+0.132138208
1164+0.00000000,
+0.653083324E-1,
-0.763353258E-1,
+0.318178564
1165call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1167+0.235387519,
-0.194309540E-1,
-0.962856039E-1,
-0.113151543
1168-0.194309540E-1,
+0.282119364,
-0.147216981E-1,
-0.165110230
1169-0.962856039E-1,
-0.147216981E-1,
+0.283279151,
+0.132138208
1170-0.113151543,
-0.165110230,
+0.132138208,
+0.318178564
1171call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1173+0.235387519,
-0.363548845E-1,
-0.940104872E-1,
-0.189470693
1174-0.194309540E-1,
+0.282119364,
-0.666179433E-1,
-0.131700113
1175-0.962856039E-1,
-0.147216981E-1,
+0.283279151,
+0.177193880
1176-0.113151543,
-0.165110230,
+0.132138208,
+0.318178564
1178+0.235387415,
-0.194309521E-1,
-0.962856039E-1,
-0.113151498
1179-0.194309521E-1,
+0.282119334,
-0.147216897E-1,
-0.165110275
1180-0.962856039E-1,
-0.147216897E-1,
+0.283279181,
+0.132138297
1181-0.113151498,
-0.165110275,
+0.132138297,
+0.318178713
1186 lb(isam)
= ub(isam
- 1)
+ 1
1194call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1195call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1199-0.695487022,
-0.232330918,
-0.518701434,
-0.469264507,
-0.903761029,
-0.573341250,
-0.690165997
1202+1.69031453,
+1.28638148,
+1.11450469,
+1.47378707,
+1.89527369,
+1.28675365,
+1.86891818
1203cov(:,:,
0)
= getCov(sample, dim, rweight)
1207 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1208 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1210call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1213call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1216call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1225 lb(isam)
= ub(isam
- 1)
+ 1
1233call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1234call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1238+0.813830495,
-0.697535753,
+0.531675339,
-0.663321972,
-0.224798441,
-0.639760375,
-0.740989566,
-0.628053784,
+0.370065093,
-0.474365711,
+0.643649459
1239+0.641063213,
+0.108627081E-1,
+0.620541096,
+0.226497650E-1,
+0.727605820E-2,
+0.940457702,
+0.477403045,
-0.963360071,
+0.332710505,
-0.612947583,
+0.322023869
1240+0.463461876,
-0.269112587E-1,
+0.536152124,
+0.935295582,
-0.501708150,
+0.165737152,
-0.618223429,
-0.639381528,
-0.693918228,
+0.472055316,
+0.761637092
1241+0.976543784,
-0.595544338,
-0.214109421E-1,
-0.295390129,
+0.558245182E-2,
-0.423438191,
-0.586519599,
-0.966666698,
+0.111133933,
-0.193339229,
-0.206089258
1244+1.33977342,
+1.56714368,
+1.69980264,
+1.48273349,
+1.15075803,
+1.30249929,
+1.80330729,
+1.73553944,
+1.04150343,
+1.63476944,
+1.40844989
1245cov(:,:,
0)
= getCov(sample, dim, rweight)
1247+0.343572348,
+0.133891404,
+0.120037317,
+0.212991893
1248+0.133891404,
+0.306821913,
+0.724486113E-1,
+0.126839072
1249+0.120037317,
+0.724486113E-1,
+0.331439912,
+0.106509455
1250+0.212991893,
+0.126839072,
+0.106509455,
+0.225562170
1252 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1253 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1255call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1257+0.343572110,
+0.133891329,
+0.120037332,
+0.212991804
1258+0.700649232E-44,
+0.306821913,
+0.724486262E-1,
+0.126839086
1259+0.00000000,
-0.147216981E-1,
+0.331439883,
+0.106509343
1260+0.00000000,
-0.165110230,
+0.132138208,
+0.225562274
1261call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1263+0.343572110,
+0.133891329,
+0.120037332,
+0.212991804
1264+0.133891329,
+0.306821913,
+0.724486262E-1,
+0.126839086
1265+0.120037332,
+0.724486262E-1,
+0.331439883,
+0.106509343
1266+0.212991804,
+0.126839086,
+0.106509343,
+0.225562274
1267call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1269+0.343572110,
+0.583240241E-1,
+0.112679943,
+0.106323555
1270+0.133891329,
+0.306821913,
+0.396544114E-1,
+0.807652250E-1
1271+0.120037332,
+0.724486262E-1,
+0.331439883,
+0.628527701E-1
1272+0.212991804,
+0.126839086,
+0.106509343,
+0.225562274
1274+0.343572348,
+0.133891404,
+0.120037317,
+0.212991893
1275+0.133891404,
+0.306821913,
+0.724486113E-1,
+0.126839072
1276+0.120037317,
+0.724486113E-1,
+0.331439912,
+0.106509455
1277+0.212991893,
+0.126839072,
+0.106509455,
+0.225562170
1282 lb(isam)
= ub(isam
- 1)
+ 1
1290call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1291call setRebound(mean, [
1_IK,
1_IK], [ndim, nsam])
1295+0.582108974,
+0.301891685,
-0.126112819,
+0.113933206,
+0.951265812,
-0.782470226,
+0.366373420
1298+1.24149013,
+1.79682541,
+1.76411724,
+1.70669603,
+1.27832305,
+1.02916932,
+1.75521207
1299cov(:,:,
0)
= getCov(sample, dim, rweight)
1303 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1304 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, rweight(lb(isam):ub(isam)))
1306call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1309call setCovMerged(covMerged, cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1312call setCovMerged(cov(:,:,
2), cov(:,:,
1), mean(:,
1)
- mean(:,
2),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)