17 integer(IK) ,
parameter :: nsam
= 2
18 integer(IK) ,
allocatable :: iweight(:)
19 integer(IK) :: isam, ndim, lb(nsam), ub(nsam)
20 integer(IK) :: dim, itry, ntry
= 10
21 type(display_type) :: disp
25 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
26 call disp%show(
"!Compute the biased merged covariance of a multivariate sample.")
27 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
31 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), meanMerged(:), covMerged(:,:)
32 real(TKG),
allocatable :: sample(:,:)
35 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
38 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
39 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
42 lb(isam)
= ub(isam
- 1)
+ 1
49 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
51 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
52 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
53 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
54 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
55 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
57 call disp%show(
"call setResized(meanMerged, ndim)")
59 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
63 call disp%show(
"cov(:,:,0) = getCov(sample, dim)")
64 cov(:,:,
0)
= getCov(sample, dim)
65 call disp%show(
"cov(:,:,0) ! reference")
67 call disp%show(
"mean(:,0) = getMean(sample, dim)")
68 mean(:,
0)
= getMean(sample, dim)
69 call disp%show(
"mean(:,0) ! reference")
72 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim)")
73 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim)")
76 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
77 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
79 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG), uppDia)")
80 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
83 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG), lowDia)")
84 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
87 call disp%show(
"call setCovMeanMerged(cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(ub(1), TKG) / real(ub(2), TKG), uppDia)")
88 call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
91 call disp%show(
"cov(:,:,0) ! reference")
98 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
99 call disp%show(
"!Compute the biased merged covariance of a frequency weighted multivariate sample.")
100 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
104 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), meanMerged(:), covMerged(:,:)
105 real(TKG),
allocatable :: sample(:,:)
108 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
111 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
112 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
115 lb(isam)
= ub(isam
- 1)
+ 1
122 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
124 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
125 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
126 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
127 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
128 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
130 call disp%show(
"call setResized(meanMerged, ndim)")
132 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
136 call disp%show(
"iweight = getUnifRand(1, 10, size(sample, dim, IK))")
140 call disp%show(
"cov(:,:,0) = getCov(sample, dim, iweight)")
141 cov(:,:,
0)
= getCov(sample, dim, iweight)
142 call disp%show(
"cov(:,:,0) ! reference")
144 call disp%show(
"mean(:,0) = getMean(sample, dim, iweight)")
145 mean(:,
0)
= getMean(sample, dim, iweight)
146 call disp%show(
"mean(:,0) ! reference")
149 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))")
150 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))")
153 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
154 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
156 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), uppDia)")
157 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
160 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), lowDia)")
161 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
164 call disp%show(
"call setCovMeanMerged(cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(iweight(:ub(1))), TKG) / real(sum(iweight), TKG), uppDia)")
165 call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
168 call disp%show(
"cov(:,:,0) ! reference")
175 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
176 call disp%show(
"!Compute the biased merged covariance of a reliability weighted multivariate sample.")
177 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
181 real(TKG),
allocatable :: mean(:,:), cov(:,:,:), meanMerged(:), covMerged(:,:)
182 real(TKG),
allocatable :: sample(:,:)
183 real(TKG),
allocatable :: rweight(:)
186 call disp%show(
"dim = 2; lb(1) = 1; ub(1) = getUnifRand(2, 7)")
189 call disp%show(
" lb(isam) = ub(isam - 1) + 1")
190 call disp%show(
" ub(isam) = ub(isam - 1) + getUnifRand(2, 7)")
193 lb(isam)
= ub(isam
- 1)
+ 1
200 call disp%show(
"ndim = getUnifRand(1, minval(ub - lb + 1, 1))")
202 call disp%show(
"call setRebound(cov, [1_IK, 1_IK, 0_IK], [ndim, ndim, nsam])")
203 call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
204 call disp%show(
"call setRebound(mean, [1_IK, 0_IK], [ndim, nsam])")
205 call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
206 call disp%show(
"call setResized(covMerged, [ndim, ndim])")
208 call disp%show(
"call setResized(meanMerged, ndim)")
210 call disp%show(
"sample = getUnifRand(-1., +1., ndim, ub(nsam))")
214 call disp%show(
"rweight = getUnifRand(1., 2., size(sample, dim, IK))")
218 call disp%show(
"cov(:,:,0) = getCov(sample, 2_IK, rweight)")
219 cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
220 call disp%show(
"cov(:,:,0) ! reference")
222 call disp%show(
"mean(:,0) = getMean(sample, dim, rweight)")
223 mean(:,
0)
= getMean(sample, dim, rweight)
224 call disp%show(
"mean(:,0) ! reference")
227 call disp%show(
" cov(:,:,isam) = getCov(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
228 call disp%show(
" mean(:,isam) = getMean(sample(:,lb(isam):ub(isam)), 2_IK, rweight(lb(isam):ub(isam)))")
231 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
232 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
234 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), uppDia)")
235 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
238 call disp%show(
"call setCovMeanMerged(covMerged, meanMerged, cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), lowDia)")
239 call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
242 call disp%show(
"call setCovMeanMerged(cov(:,:,2), mean(:,2), cov(:,:,1), mean(:,1), real(sum(rweight(:ub(1))), TKG) / real(sum(rweight), TKG), uppDia)")
243 call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
246 call disp%show(
"cov(:,:,0) ! reference")
Generate minimally-spaced character, integer, real sequences or sequences at fixed intervals of size ...
Resize (shrink or expand) an input allocatable array of rank 1..3 to arbitrary new lower and upper bo...
Allocate or resize (shrink or expand) an input allocatable scalar string or array of rank 1....
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
This is a generic method of the derived type display_type with pass attribute.
This is a generic method of the derived type display_type with pass attribute.
Generate and return the (optionally unbiased) covariance matrix of a pair of (potentially weighted) t...
Generate and return the (weighted) mean of an input sample of nsam observations with ndim = 1 or 2 at...
This module contains procedures and generic interfaces for generating ranges of discrete character,...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for input/output (IO) or generic display operations on st...
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
This module contains classes and procedures for computing the first moment (i.e., the statistical mea...
Generate and return an object of type display_type.
9 lb(isam)
= ub(isam
- 1)
+ 1
17call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
18call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
23-0.983392239,
-0.906538248,
+0.566596746,
+0.285900354,
+0.919802070,
+0.801415443,
-0.394834876,
+0.815028310,
+0.655870080
24cov(:,:,
0)
= getCov(sample, dim)
31 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
32 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
34call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
37call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
40call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
49 lb(isam)
= ub(isam
- 1)
+ 1
57call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
58call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
63-0.476266146,
-0.145911813,
+0.230103374,
+0.302848935,
+0.429608107,
+0.946276426,
+0.772754312,
+0.711741447
64cov(:,:,
0)
= getCov(sample, dim)
71 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
72 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
74call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
77call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
80call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
89 lb(isam)
= ub(isam
- 1)
+ 1
97call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
98call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
103+0.762442112,
+0.728666782E-2,
+0.231051087,
+0.820343733,
-0.950583220E-1,
+0.120338082,
-0.888931751E-2
104+0.368424654E-1,
-0.805558801,
-0.838758945E-1,
-0.458779693,
+0.910453200,
+0.856227756,
+0.996220112
105+0.412815809,
-0.711351991,
-0.280172348,
-0.586521268,
+0.456707716,
-0.254994631,
-0.517393827
106cov(:,:,
0)
= getCov(sample, dim)
108+0.121280611,
-0.106701516,
+0.118461680E-1
109-0.106701516,
+0.445901543,
+0.124278843
110+0.118461680E-1,
+0.124278843,
+0.189563289
111mean(:,
0)
= getMean(sample, dim)
113+0.262502015,
+0.207361311,
-0.211558655
115 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
116 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
118call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
120+0.121280707,
-0.106701545,
+0.118461791E-1
121+0.00000000,
+0.445901543,
+0.124278821
122+0.00000000,
+0.00000000,
+0.189563334
123call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
125+0.121280707,
-0.106701545,
+0.118461791E-1
126-0.106701545,
+0.445901543,
+0.124278821
127+0.118461791E-1,
+0.124278821,
+0.189563334
128call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
130+0.121280707,
-0.106701545,
+0.118461791E-1
131-0.212676644,
+0.445901543,
+0.124278821
132-0.904808640E-1,
+0.117704809,
+0.189563334
134+0.121280611,
-0.106701516,
+0.118461680E-1
135-0.106701516,
+0.445901543,
+0.124278843
136+0.118461680E-1,
+0.124278843,
+0.189563289
141 lb(isam)
= ub(isam
- 1)
+ 1
149call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
150call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
155+0.998388290,
+0.230247140,
+0.619487762E-1,
-0.617218018E-2,
+0.607528329,
-0.175619245,
+0.369630814,
-0.142510176,
-0.411397934,
-0.274683237E-1,
+0.649638891
156cov(:,:,
0)
= getCov(sample, dim)
159mean(:,
0)
= getMean(sample, dim)
163 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
164 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
166call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
169call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
172call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
181 lb(isam)
= ub(isam
- 1)
+ 1
189call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
190call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
195+0.294937968,
+0.349892616,
+0.674491882,
-0.816441059,
-0.647302985,
-0.428435206
196cov(:,:,
0)
= getCov(sample, dim)
199mean(:,
0)
= getMean(sample, dim)
203 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
204 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
206call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
209call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
212call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
221 lb(isam)
= ub(isam
- 1)
+ 1
229call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
230call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
235+0.922245979,
-0.145381093,
-0.990380406,
+0.878242135,
-0.139081359
236cov(:,:,
0)
= getCov(sample, dim)
239mean(:,
0)
= getMean(sample, dim)
243 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
244 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
246call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
249call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
252call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
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,
0_IK], [ndim, nsam])
275-0.370319843,
-0.512055159,
-0.463492513,
+0.932714224,
-0.507686853,
-0.410154581,
+0.200248837
276-0.742630363,
-0.347755909,
+0.612773657,
+0.655895710,
-0.401390791E-1,
+0.340716124,
+0.818808079E-1
277cov(:,:,
0)
= getCov(sample, dim)
279+0.252505064,
+0.109776556
280+0.109776556,
+0.222515658
281mean(:,
0)
= getMean(sample, dim)
283-0.161535129,
+0.801058561E-1
285 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
286 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
288call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
290+0.252505124,
+0.109776571
291+0.307935338E-40,
+0.222515613
292call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
294+0.252505124,
+0.109776571
295+0.109776571,
+0.222515613
296call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
298+0.252505124,
+0.109776571
299+0.113237143,
+0.222515613
301+0.252505064,
+0.109776556
302+0.109776556,
+0.222515658
307 lb(isam)
= ub(isam
- 1)
+ 1
315call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
316call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
321+0.410302639,
-0.943424106,
+0.937252045E-1,
+0.418286800,
-0.195232153
322cov(:,:,
0)
= getCov(sample, dim)
325mean(:,
0)
= getMean(sample, dim)
329 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
330 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
332call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
335call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
338call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
347 lb(isam)
= ub(isam
- 1)
+ 1
355call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
356call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
361-0.959451199,
-0.289603949,
-0.944484591,
+0.637496710E-1,
+0.522669315,
-0.953766465
362+0.618827820,
+0.365375876,
-0.806772828,
+0.877617121,
+0.617314458,
-0.152364731
363cov(:,:,
0)
= getCov(sample, dim)
365+0.331727028,
+0.205851972
366+0.205851972,
+0.326128066
367mean(:,
0)
= getMean(sample, dim)
369-0.426814497,
+0.253332973
371 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
372 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
374call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
376+0.331726998,
+0.205851942
377+0.307935338E-40,
+0.326128095
378call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
380+0.331726998,
+0.205851942
381+0.205851942,
+0.326128095
382call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
384+0.331726998,
+0.205851942
385+0.365405500,
+0.326128095
387+0.331727028,
+0.205851972
388+0.205851972,
+0.326128066
393 lb(isam)
= ub(isam
- 1)
+ 1
401call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
402call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
407+0.673266053,
+0.837436795,
-0.464135408E-1,
+0.327112556,
+0.704258919,
-0.931052923,
+0.939259171,
-0.103528619,
-0.196568489,
+0.914087772,
+0.763027191
408+0.830568790,
-0.993834376,
-0.204558730,
+0.117491722,
-0.396186113E-1,
-0.548318982,
+0.644092560E-1,
-0.575642586E-1,
-0.384115100,
+0.937119126,
-0.563284516
409+0.692670345E-1,
-0.798385859,
+0.469762087,
+0.861128926,
+0.125003457,
+0.617935061,
-0.892008305,
+0.848717213,
-0.754354358,
-0.434735656,
-0.305904150E-1
410-0.547821879,
-0.689647555,
+0.261740804,
-0.999858379,
+0.903812289,
-0.384301186,
-0.429983735,
+0.866266370,
-0.947055221,
-0.982833385,
+0.407021046
411-0.867243886,
+0.338432431,
-0.280110002,
+0.750241518,
-0.786336899,
-0.574440956,
-0.652234674,
+0.628817081E-1,
+0.972767830,
+0.254431129,
+0.448332071
412cov(:,:,
0)
= getCov(sample, dim)
414+0.327882856,
+0.102108091,
-0.188759327,
-0.263094474E-1,
-0.162821263E-1
415+0.102108091,
+0.301954985,
+0.283485502E-1,
-0.749415532E-1,
-0.930879787E-1
416-0.188759327,
+0.283485502E-1,
+0.388752818,
+0.174273610,
-0.445051007E-1
417-0.263094474E-1,
-0.749415532E-1,
+0.174273610,
+0.471328914,
-0.163878277
418-0.162821263E-1,
-0.930879787E-1,
-0.445051007E-1,
-0.163878277,
+0.371596426
419mean(:,
0)
= getMean(sample, dim)
421+0.352807760,
-0.765186995E-1,
+0.743083563E-2,
-0.231150985,
-0.302981585E-1
423 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim)
424 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim)
426call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
428+0.327882707,
+0.102108099,
-0.188759327,
-0.263094380E-1,
-0.162821002E-1
429+0.205851942,
+0.301954955,
+0.283485651E-1,
-0.749415755E-1,
-0.930880532E-1
430+0.00000000,
+0.00000000,
+0.388752788,
+0.174273610,
-0.445050597E-1
431+0.00000000,
+0.00000000,
+0.00000000,
+0.471328855,
-0.163878277
432+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.371596485
433call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), lowDia)
435+0.327882707,
+0.102108099,
-0.188759327,
-0.263094380E-1,
-0.162821002E-1
436+0.102108099,
+0.301954955,
+0.283485651E-1,
-0.749415755E-1,
-0.930880532E-1
437-0.188759327,
+0.283485651E-1,
+0.388752788,
+0.174273610,
-0.445050597E-1
438-0.263094380E-1,
-0.749415755E-1,
+0.174273610,
+0.471328855,
-0.163878277
439-0.162821002E-1,
-0.930880532E-1,
-0.445050634E-1,
-0.163878277,
+0.371596485
440call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(ub(
1), TKG)
/ real(ub(
2), TKG), uppDia)
442+0.327882707,
+0.102108099,
-0.188759327,
-0.263094380E-1,
-0.162821002E-1
443+0.201111168,
+0.301954955,
+0.283485651E-1,
-0.749415755E-1,
-0.930880532E-1
444-0.270758241,
-0.100810155,
+0.388752788,
+0.174273610,
-0.445050597E-1
445-0.329711549E-1,
-0.131437525,
+0.308733165,
+0.471328855,
-0.163878277
446+0.301681366E-1,
-0.118925376E-1,
-0.871535167E-1,
-0.512859225E-1,
+0.371596485
448+0.327882856,
+0.102108091,
-0.188759327,
-0.263094474E-1,
-0.162821263E-1
449+0.102108091,
+0.301954985,
+0.283485502E-1,
-0.749415532E-1,
-0.930879787E-1
450-0.188759327,
+0.283485502E-1,
+0.388752818,
+0.174273610,
-0.445051007E-1
451-0.263094474E-1,
-0.749415532E-1,
+0.174273610,
+0.471328914,
-0.163878277
452-0.162821263E-1,
-0.930879787E-1,
-0.445051007E-1,
-0.163878277,
+0.371596426
462 lb(isam)
= ub(isam
- 1)
+ 1
470call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
471call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
476-0.872341275,
+0.439786077,
+0.664734125,
-0.315796137E-1,
+0.412347674,
+0.121695399,
+0.274711847,
+0.570886016
477+0.323488116,
-0.585140228,
-0.647064328,
+0.549629092,
+0.733035088,
+0.445165873,
+0.611363888,
-0.431539536
478+0.550215840,
-0.536445975,
-0.914926291,
+0.879891038,
+0.932551742,
-0.722442031,
-0.433287382,
+0.852670908
481+8,
+6,
+3,
+6,
+7,
+6,
+8,
+9
482cov(:,:,
0)
= getCov(sample, dim, iweight)
484+0.227531940,
-0.961641073E-1,
-0.615194812E-1
485-0.961641073E-1,
+0.263486385,
+0.687750205E-1
486-0.615194812E-1,
+0.687750205E-1,
+0.502598107
487mean(:,
0)
= getMean(sample, dim, iweight)
489+0.158811063,
+0.174395114,
+0.190916374
491 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
492 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
494call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
496+0.227531955,
-0.961641148E-1,
-0.615195073E-1
497+0.307935338E-40,
+0.263486415,
+0.687749982E-1
498+0.00000000,
-0.192902848,
+0.502598166
499call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
501+0.227531955,
-0.961641148E-1,
-0.615195073E-1
502-0.961641148E-1,
+0.263486415,
+0.687749982E-1
503-0.615195073E-1,
+0.687749982E-1,
+0.502598166
504call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
506+0.227531955,
-0.961641148E-1,
-0.615195073E-1
507-0.799326450E-1,
+0.263486415,
+0.687749982E-1
508+0.128076866,
-0.319137603,
+0.502598166
510+0.227531940,
-0.961641073E-1,
-0.615194812E-1
511-0.961641073E-1,
+0.263486385,
+0.687750205E-1
512-0.615194812E-1,
+0.687750205E-1,
+0.502598107
517 lb(isam)
= ub(isam
- 1)
+ 1
525call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
526call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
531+0.871731997,
+0.906964898,
-0.902264357,
-0.408882141,
+0.938417435,
+0.426919460,
-0.851143360,
-0.659917593,
-0.386589766
532+0.527657747,
-0.542777300,
+0.201543689,
+0.843955278,
-0.584905267,
-0.989788890,
+0.452311397,
-0.598481297,
-0.406501412
533+0.533214092,
+0.397715449,
+0.177038670,
-0.640622854,
+0.970257401,
+0.342923641,
-0.348046303,
+0.161872268,
+0.480622530
536+8,
+2,
+7,
+8,
+1,
+1,
+3,
+6,
+2
537cov(:,:,
0)
= getCov(sample, dim, iweight)
539+0.518954813,
+0.105255526E-1,
+0.182831869
540+0.105255526E-1,
+0.316870838,
-0.136866704
541+0.182831869,
-0.136866704,
+0.206183985
542mean(:,
0)
= getMean(sample, dim, iweight)
544-0.176839486,
+0.175697476,
+0.888671726E-1
546 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
547 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
549call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
551+0.518954694,
+0.105255451E-1,
+0.182831764
552-0.961641148E-1,
+0.316870809,
-0.136866733
553-0.615195073E-1,
+0.687749982E-1,
+0.206183940
554call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
556+0.518954694,
+0.105255451E-1,
+0.182831764
557+0.105255451E-1,
+0.316870809,
-0.136866733
558+0.182831764,
-0.136866733,
+0.206183940
559call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
561+0.518954694,
+0.105255451E-1,
+0.182831764
562-0.605103485E-1,
+0.316870809,
-0.136866733
563+0.832844526E-1,
-0.303807318,
+0.206183940
565+0.518954813,
+0.105255526E-1,
+0.182831869
566+0.105255526E-1,
+0.316870838,
-0.136866704
567+0.182831869,
-0.136866704,
+0.206183985
572 lb(isam)
= ub(isam
- 1)
+ 1
580call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
581call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
586-0.719815969,
+0.434856772,
-0.514058352,
+0.184169888,
+0.243343830,
-0.580325365,
+0.915624022,
-0.459622145E-1,
-0.802342415,
+0.198796272,
-0.453128815E-1
587-0.457105041,
-0.195128560,
+0.943547845,
-0.193790197E-1,
+0.144334912,
-0.268786192,
+0.531285524,
+0.960040808,
-0.706946969,
+0.150730610E-1,
+0.212450504
588+0.350733995,
+0.494531274,
-0.960105777,
-0.529353499,
-0.607080579,
+0.116647720,
+0.608706713,
+0.672276020,
+0.596568823,
-0.230395079,
+0.773772597
589-0.366394758,
-0.333453178,
-0.935181856,
+0.971811891,
-0.473067760E-1,
+0.823431373,
-0.553089738,
+0.940012217,
+0.684583187E-2,
+0.703319788,
+0.784388065
592+6,
+10,
+7,
+9,
+4,
+9,
+9,
+10,
+10,
+10,
+7
593cov(:,:,
0)
= getCov(sample, dim, iweight)
595+0.270221680,
+0.101914041,
+0.223563369E-1,
-0.224917065E-1
596+0.101914041,
+0.272711068,
-0.444898419E-1,
-0.125027336E-1
597+0.223563369E-1,
-0.444898419E-1,
+0.299578100,
+0.287952218E-1
598-0.224917065E-1,
-0.125027336E-1,
+0.287952218E-1,
+0.425186962
599mean(:,
0)
= getMean(sample, dim, iweight)
601-0.520045944E-1,
+0.971996114E-1,
+0.169951305,
+0.229708537
603 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
604 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
606call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
608+0.270221651,
+0.101914056,
+0.223563053E-1,
-0.224917587E-1
609+0.105255451E-1,
+0.272711009,
-0.444898307E-1,
-0.125027318E-1
610+0.182831764,
-0.136866733,
+0.299578130,
+0.287952051E-1
611+0.00000000,
+0.145355135,
-0.765186995E-1,
+0.425186902
612call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
614+0.270221651,
+0.101914056,
+0.223563053E-1,
-0.224917587E-1
615+0.101914056,
+0.272711009,
-0.444898307E-1,
-0.125027318E-1
616+0.223563053E-1,
-0.444898307E-1,
+0.299578130,
+0.287952051E-1
617-0.224917587E-1,
-0.125027318E-1,
+0.287952051E-1,
+0.425186902
618call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
620+0.270221651,
+0.101914056,
+0.223563053E-1,
-0.224917587E-1
621+0.207746193,
+0.272711009,
-0.444898307E-1,
-0.125027318E-1
622+0.269905431E-2,
+0.605359562E-1,
+0.299578130,
+0.287952051E-1
623-0.118380889,
+0.492995605E-1,
-0.562449619E-1,
+0.425186902
625+0.270221680,
+0.101914041,
+0.223563369E-1,
-0.224917065E-1
626+0.101914041,
+0.272711068,
-0.444898419E-1,
-0.125027336E-1
627+0.223563369E-1,
-0.444898419E-1,
+0.299578100,
+0.287952218E-1
628-0.224917065E-1,
-0.125027336E-1,
+0.287952218E-1,
+0.425186962
633 lb(isam)
= ub(isam
- 1)
+ 1
641call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
642call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
647+0.653579712,
-0.467783213,
-0.604078650,
-0.804414868,
-0.870658159,
-0.807228088,
+0.510380983,
-0.383740783,
+0.675415397,
+0.493945479,
-0.204651356E-1
650+8,
+9,
+4,
+8,
+3,
+2,
+7,
+8,
+1,
+7,
+8
651cov(:,:,
0)
= getCov(sample, dim, iweight)
654mean(:,
0)
= getMean(sample, dim, iweight)
658 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
659 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
661call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
664call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
667call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
676 lb(isam)
= ub(isam
- 1)
+ 1
684call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
685call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
690+0.567971230,
-0.926035285,
+0.879066229,
-0.123556852E-1,
+0.388380527,
+0.511034966,
-0.978226900,
+0.212329507,
+0.125435352
691+0.522053123,
-0.451764941,
+0.848307967,
+0.150354147,
-0.912380099,
-0.872391701,
-0.863469839E-1,
-0.189986944,
-0.715277195
692+0.137670159,
-0.813375473,
+0.347590446,
+0.969531536,
+0.100210547,
-0.365713000,
-0.763778090,
-0.679090381,
-0.292024493
695+9,
+4,
+5,
+1,
+8,
+6,
+8,
+7,
+2
696cov(:,:,
0)
= getCov(sample, dim, iweight)
698+0.405663967,
+0.759069994E-1,
+0.214887172
699+0.759069994E-1,
+0.358743310,
+0.930816457E-1
700+0.214887172,
+0.930816457E-1,
+0.207162485
701mean(:,
0)
= getMean(sample, dim, iweight)
703+0.117503822,
-0.174026325,
-0.242949724
705 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
706 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
708call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
710+0.405664027,
+0.759069696E-1,
+0.214887172
711+0.307935338E-40,
+0.358743191,
+0.930816159E-1
712+0.00000000,
-0.136866733,
+0.207162529
713call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
715+0.405664027,
+0.759069696E-1,
+0.214887172
716+0.759069771E-1,
+0.358743191,
+0.930816159E-1
717+0.214887172,
+0.930816159E-1,
+0.207162529
718call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
720+0.405664027,
+0.759069696E-1,
+0.214887172
721-0.169425681,
+0.358743191,
+0.930816159E-1
722+0.132891744,
-0.115795963,
+0.207162529
724+0.405663967,
+0.759069994E-1,
+0.214887172
725+0.759069994E-1,
+0.358743310,
+0.930816457E-1
726+0.214887172,
+0.930816457E-1,
+0.207162485
731 lb(isam)
= ub(isam
- 1)
+ 1
739call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
740call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
745-0.734853745E-2,
-0.715180397,
+0.350714684,
-0.124251366,
-0.147824407,
-0.818401694,
+0.807071567,
-0.686981916,
-0.659407377E-1
746+0.626269817,
+0.148324251,
+0.164194703,
-0.286100745,
+0.353755951E-1,
-0.219994903,
+0.274710178,
-0.650122285,
-0.862091541
747+0.450922728,
-0.259552598,
+0.153432369,
+0.936720610,
-0.798994303,
-0.578143239,
-0.677618980E-1,
+0.186652184,
-0.628881335
748-0.681552887E-1,
-0.739362001,
+0.229146242,
-0.758110285E-1,
+0.323850393,
+0.585678458,
-0.902093887,
+0.684251070,
-0.742348075
751+5,
+2,
+3,
+4,
+10,
+3,
+2,
+8,
+2
752cov(:,:,
0)
= getCov(sample, dim, iweight)
754+0.173242986,
+0.977050737E-1,
+0.195526704E-1,
-0.112603299
755+0.977050737E-1,
+0.184577450,
+0.372023461E-2,
-0.693374649E-1
756+0.195526704E-1,
+0.372023461E-2,
+0.329481274,
-0.238904469E-1
757-0.112603299,
-0.693374649E-1,
-0.238904469E-1,
+0.230867073
758mean(:,
0)
= getMean(sample, dim, iweight)
760-0.227154046,
-0.100148484,
-0.944039449E-1,
+0.147316888
762 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
763 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
765call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
767+0.173242956,
+0.977050886E-1,
+0.195526686E-1,
-0.112603292
768+0.759069771E-1,
+0.184577420,
+0.372027140E-2,
-0.693374872E-1
769+0.214887172,
+0.930816159E-1,
+0.329481363,
-0.238904692E-1
770+0.00000000,
-0.125027318E-1,
+0.287952051E-1,
+0.230867088
771call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
773+0.173242956,
+0.977050886E-1,
+0.195526686E-1,
-0.112603292
774+0.977050811E-1,
+0.184577420,
+0.372027140E-2,
-0.693374872E-1
775+0.195526686E-1,
+0.372027140E-2,
+0.329481363,
-0.238904692E-1
776-0.112603292,
-0.693374872E-1,
-0.238904692E-1,
+0.230867088
777call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
779+0.173242956,
+0.977050886E-1,
+0.195526686E-1,
-0.112603292
780+0.109069780,
+0.184577420,
+0.372027140E-2,
-0.693374872E-1
781-0.225117225E-1,
-0.205860306E-1,
+0.329481363,
-0.238904692E-1
782-0.323776186,
-0.874623284E-1,
+0.994289443E-1,
+0.230867088
784+0.173242986,
+0.977050737E-1,
+0.195526704E-1,
-0.112603299
785+0.977050737E-1,
+0.184577450,
+0.372023461E-2,
-0.693374649E-1
786+0.195526704E-1,
+0.372023461E-2,
+0.329481274,
-0.238904469E-1
787-0.112603299,
-0.693374649E-1,
-0.238904469E-1,
+0.230867073
792 lb(isam)
= ub(isam
- 1)
+ 1
800call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
801call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
806-0.684488416,
-0.759531975,
-0.127464533E-1,
-0.742536902,
-0.352996707,
-0.317959905,
-0.808501840,
+0.903011560E-1,
+0.377865553,
+0.791180611
807+0.992210746,
-0.855833292E-1,
+0.904832125,
-0.457458973,
+0.201971173,
-0.945067763,
-0.203007102,
+0.907687664,
-0.740822196,
-0.991106629
810+4,
+2,
+10,
+7,
+5,
+8,
+7,
+3,
+8,
+6
811cov(:,:,
0)
= getCov(sample, dim, iweight)
813+0.260188103,
-0.876784027E-1
814-0.876784027E-1,
+0.551104903
815mean(:,
0)
= getMean(sample, dim, iweight)
817-0.191825107,
-0.124635115
819 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
820 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
822call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
824+0.260188133,
-0.876784623E-1
825+0.977050811E-1,
+0.551104903
826call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
828+0.260188133,
-0.876784623E-1
829-0.876784623E-1,
+0.551104903
830call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
832+0.260188133,
-0.876784623E-1
833-0.187523782,
+0.551104903
835+0.260188103,
-0.876784027E-1
836-0.876784027E-1,
+0.551104903
841 lb(isam)
= ub(isam
- 1)
+ 1
849call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
850call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
855+0.565159321E-1,
-0.808794141,
-0.753669381,
-0.701167822,
-0.391706228,
-0.636185765,
-0.694110990
858+10,
+3,
+3,
+3,
+5,
+6,
+4
859cov(:,:,
0)
= getCov(sample, dim, iweight)
862mean(:,
0)
= getMean(sample, dim, iweight)
866 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
867 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
869call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
872call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
875call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
884 lb(isam)
= ub(isam
- 1)
+ 1
892call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
893call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
898+0.956273079E-1,
+0.669129610,
+0.911620498,
+0.222380161,
-0.843547702,
-0.329590321,
-0.329339266,
-0.580144048
899-0.234357595,
+0.613723636,
+0.272315264,
-0.175609231,
-0.727025032,
+0.619257689,
-0.450201750,
+0.233683825
902+1,
+1,
+1,
+3,
+6,
+3,
+2,
+2
903cov(:,:,
0)
= getCov(sample, dim, iweight)
905+0.279724836,
+0.147062108
906+0.147062108,
+0.264122218
907mean(:,
0)
= getMean(sample, dim, iweight)
909-0.290816098,
-0.148029447
911 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
912 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
914call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
916+0.279724866,
+0.147062108
917+0.307935338E-40,
+0.264122248
918call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
920+0.279724866,
+0.147062108
921+0.147062108,
+0.264122248
922call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
924+0.279724866,
+0.147062108
925-0.220523356E-2,
+0.264122248
927+0.279724836,
+0.147062108
928+0.147062108,
+0.264122218
933 lb(isam)
= ub(isam
- 1)
+ 1
941call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
942call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
947+0.123158336,
-0.346078873,
-0.694132805,
-0.367417574,
+0.769196153,
+0.770443678E-1,
+0.255964875
950+1,
+4,
+6,
+7,
+8,
+8,
+3
951cov(:,:,
0)
= getCov(sample, dim, iweight)
954mean(:,
0)
= getMean(sample, dim, iweight)
958 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
959 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)), dim, iweight(lb(isam):ub(isam)))
961call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
964call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), lowDia)
967call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(iweight(:ub(
1))), TKG)
/ real(
sum(iweight), TKG), uppDia)
981 lb(isam)
= ub(isam
- 1)
+ 1
989call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
990call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
995+0.901126862E-2,
+0.763417244,
+0.437683105,
+0.337855816E-1,
+0.980162740,
-0.672264338,
-0.126809120,
+0.309169769,
-0.149852991,
-0.210831523,
-0.498251319,
-0.843937874,
+0.993048191
998+1.64471579,
+1.16781831,
+1.39789736,
+1.39093733,
+1.43399119,
+1.63431549,
+1.41562998,
+1.88521135,
+1.53496337,
+1.41511810,
+1.24166846,
+1.51786137,
+1.07092404
999cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1002mean(:,
0)
= getMean(sample, dim, rweight)
1006 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1007 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1009call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1012call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1015call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1024 lb(isam)
= ub(isam
- 1)
+ 1
1032call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1033call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1038+0.937492847,
+0.549838424,
+0.444339871,
+0.861517429,
-0.308555484,
+0.156323910
1041+1.04534757,
+1.10362363,
+1.93200493,
+1.11134470,
+1.70671940,
+1.98849988
1042cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1045mean(:,
0)
= getMean(sample, dim, rweight)
1049 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1050 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1052call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1055call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1058call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1067 lb(isam)
= ub(isam
- 1)
+ 1
1075call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1076call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1081-0.736506224,
+0.759213686,
-0.779147387,
-0.152210593,
+0.975085735,
-0.186195374E-1,
+0.111787200,
-0.733381987,
-0.632111311,
+0.696297050
1082+0.127996206E-1,
+0.447777271,
-0.354464293,
-0.496356487E-1,
-0.936763883,
+0.384133935,
+0.759760618,
+0.267163515,
-0.425451279,
+0.563687921
1083+0.538259506,
+0.142455816,
+0.719458580,
+0.691131234,
-0.453511119,
-0.940337539,
+0.225146294,
+0.683842421,
+0.653871775,
-0.276131392
1086+1.97445488,
+1.50471163,
+1.42919898,
+1.68178129,
+1.27092385,
+1.85202765,
+1.72171569,
+1.11505628,
+1.64370000,
+1.48051727
1087cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1089+0.380958050,
+0.493932217E-1,
-0.216762081
1090+0.493932217E-1,
+0.229406312,
-0.568168201E-1
1091-0.216762081,
-0.568168201E-1,
+0.313119262
1092mean(:,
0)
= getMean(sample, dim, rweight)
1094-0.708165318E-1,
+0.874738097E-1,
+0.189223647
1096 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1097 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1099call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1101+0.380958080,
+0.493932180E-1,
-0.216762111
1102+0.307935338E-40,
+0.229406297,
-0.568168275E-1
1103+0.00000000,
+0.112103877E-43,
+0.313119233
1104call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1106+0.380958080,
+0.493932180E-1,
-0.216762111
1107+0.493932180E-1,
+0.229406297,
-0.568168275E-1
1108-0.216762111,
-0.568168275E-1,
+0.313119233
1109call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1111+0.380958080,
+0.493932180E-1,
-0.216762111
1112+0.209716454,
+0.229406297,
-0.568168275E-1
1113-0.293673217,
-0.148596227,
+0.313119233
1115+0.380958050,
+0.493932217E-1,
-0.216762081
1116+0.493932217E-1,
+0.229406312,
-0.568168201E-1
1117-0.216762081,
-0.568168201E-1,
+0.313119262
1122 lb(isam)
= ub(isam
- 1)
+ 1
1130call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1131call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1136+0.583632231,
-0.532859445,
+0.465836287,
+0.685912132,
-0.863176107,
-0.957157612E-1,
+0.302388906,
+0.888311863,
-0.182961583,
-0.445081353,
+0.198978186,
+0.790772200
1139+1.20540094,
+1.63233685,
+1.86911774,
+1.55541825,
+1.55842292,
+1.90102434,
+1.51036286,
+1.86065269,
+1.99727106,
+1.14192414,
+1.26569366,
+1.27255952
1140cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1143mean(:,
0)
= getMean(sample, dim, rweight)
1147 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1148 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1150call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1153call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1156call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1165 lb(isam)
= ub(isam
- 1)
+ 1
1173call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1174call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1179+0.839756727E-1,
+0.182847857,
+0.230667114,
-0.878359079,
+0.567645073,
+0.144941807,
+0.186116934,
-0.837831497E-1,
+0.924179554E-1
1180-0.159936428,
+0.270127058,
+0.442338228,
+0.937236547E-1,
+0.263673544,
+0.125091672,
+0.156951427,
+0.848851562,
-0.281935573
1181+0.121101379,
+0.667280316,
+0.388617039,
-0.954250813,
+0.832293034E-1,
-0.789279461,
+0.385317922,
+0.140124559E-1,
-0.419551373
1182+0.451750159,
-0.935984373,
+0.812960029,
-0.791021585,
-0.528596640E-1,
-0.944455743,
-0.517652392,
+0.110246778,
+0.245691776
1185+1.94029844,
+1.50409281,
+1.77924180,
+1.51067281,
+1.39584327,
+1.32783890,
+1.07654309,
+1.56558418,
+1.91140485
1186cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1188+0.130324572,
+0.524980761E-2,
+0.113309860,
+0.701396987E-1
1189+0.524980761E-2,
+0.106535293,
+0.541637056E-1,
+0.223275810E-2
1190+0.113309860,
+0.541637056E-1,
+0.259021789,
+0.967521518E-1
1191+0.701396987E-1,
+0.223275810E-2,
+0.967521518E-1,
+0.376871705
1192mean(:,
0)
= getMean(sample, dim, rweight)
1194+0.536773987E-1,
+0.179691672,
-0.577051714E-1,
-0.108676754
1196 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1197 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1199call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1201+0.130324572,
+0.524980016E-2,
+0.113309853,
+0.701397061E-1
1202+0.700649232E-44,
+0.106535286,
+0.541637018E-1,
+0.223276229E-2
1203+0.00000000,
+0.372027140E-2,
+0.259021789,
+0.967521369E-1
1204+0.00000000,
-0.693374872E-1,
-0.238904692E-1,
+0.376871616
1205call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1207+0.130324572,
+0.524980016E-2,
+0.113309853,
+0.701397061E-1
1208+0.524980016E-2,
+0.106535286,
+0.541637018E-1,
+0.223276229E-2
1209+0.113309853,
+0.541637018E-1,
+0.259021789,
+0.967521369E-1
1210+0.701397061E-1,
+0.223276252E-2,
+0.967521369E-1,
+0.376871616
1211call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1213+0.130324572,
+0.524980016E-2,
+0.113309853,
+0.701397061E-1
1214-0.205041841E-1,
+0.106535286,
+0.541637018E-1,
+0.223276229E-2
1215+0.185011253E-1,
+0.653388873E-1,
+0.259021789,
+0.967521369E-1
1216-0.121255601E-1,
+0.734357687E-3,
+0.484819412E-1,
+0.376871616
1218+0.130324572,
+0.524980761E-2,
+0.113309860,
+0.701396987E-1
1219+0.524980761E-2,
+0.106535293,
+0.541637056E-1,
+0.223275810E-2
1220+0.113309860,
+0.541637056E-1,
+0.259021789,
+0.967521518E-1
1221+0.701396987E-1,
+0.223275810E-2,
+0.967521518E-1,
+0.376871705
1226 lb(isam)
= ub(isam
- 1)
+ 1
1234call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1235call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1240-0.907207847,
-0.409345388,
-0.174140930E-2,
-0.540048003,
+0.678702831,
+0.533585548,
+0.923105955,
-0.608271241,
-0.657650948
1241-0.492187977,
+0.847451806,
+0.559273958,
-0.649985790,
-0.840529442,
+0.577440143,
+0.131551266,
+0.880162597,
+0.478424907
1244+1.07763088,
+1.03222203,
+1.15953946,
+1.87008238,
+1.23746681,
+1.01219761,
+1.04810524,
+1.08845937,
+1.30671573
1245cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1247+0.377780110,
-0.304243080E-1
1248-0.304243080E-1,
+0.408677459
1249mean(:,
0)
= getMean(sample, dim, rweight)
1251-0.146420389,
+0.962626189E-1
1253 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1254 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1256call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1258+0.377780169,
-0.304243006E-1
1259+0.524980016E-2,
+0.408677518
1260call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1262+0.377780169,
-0.304243006E-1
1263-0.304242931E-1,
+0.408677518
1264call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1266+0.377780169,
-0.304243006E-1
1267+0.491824234E-2,
+0.408677518
1269+0.377780110,
-0.304243080E-1
1270-0.304243080E-1,
+0.408677459
1275 lb(isam)
= ub(isam
- 1)
+ 1
1283call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1284call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1289-0.162543654,
+0.391954541,
+0.662426233,
-0.365625501,
+0.726747990,
-0.638747096
1290+0.687010765,
+0.799016237,
-0.860929489,
-0.789881587,
+0.746846437,
+0.371324539
1293+1.40831494,
+1.57480097,
+1.53521490,
+1.46696043,
+1.25989509,
+1.52375507
1294cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1296+0.269242853,
+0.103353746E-1
1297+0.103353746E-1,
+0.508305311
1298mean(:,
0)
= getMean(sample, dim, rweight)
1300+0.925172418E-1,
+0.142792404
1302 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1303 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1305call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1307+0.269242853,
+0.103353951E-1
1308-0.304242931E-1,
+0.508305311
1309call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1311+0.269242853,
+0.103353951E-1
1312+0.103353951E-1,
+0.508305311
1313call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1315+0.269242853,
+0.103353951E-1
1316+0.127041519,
+0.508305311
1318+0.269242853,
+0.103353746E-1
1319+0.103353746E-1,
+0.508305311
1324 lb(isam)
= ub(isam
- 1)
+ 1
1332call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1333call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1338+0.458455205,
-0.691320419,
+0.682012916,
+0.416854501,
-0.749879241,
-0.584628344,
-0.368337631,
+0.646575212
1341+1.72194004,
+1.22958946,
+1.79436946,
+1.83464265,
+1.80115485,
+1.98471534,
+1.95079637,
+1.45613647
1342cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1345mean(:,
0)
= getMean(sample, dim, rweight)
1349 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1350 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1352call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1355call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1358call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1367 lb(isam)
= ub(isam
- 1)
+ 1
1375call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1376call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1381+0.141153336,
+0.413580298,
+0.347508907,
-0.342417717,
+0.262140512,
+0.808933616,
-0.625546455,
-0.460398197E-1
1382-0.647318125,
-0.633061409,
+0.671309233,
-0.327589512E-1,
-0.433012247,
+0.198728561,
+0.848019719,
-0.295746446
1385+1.91358566,
+1.74581051,
+1.28563809,
+1.98341489,
+1.90950274,
+1.71698928,
+1.37907112,
+1.23091888
1386cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1388+0.174270958,
-0.610865094E-1
1389-0.610865094E-1,
+0.265407890
1390mean(:,
0)
= getMean(sample, dim, rweight)
1392+0.131402120,
-0.931263044E-1
1394 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1395 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1397call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1399+0.174270943,
-0.610865019E-1
1400+0.307935338E-40,
+0.265407890
1401call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1403+0.174270943,
-0.610865019E-1
1404-0.610865019E-1,
+0.265407890
1405call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1407+0.174270943,
-0.610865019E-1
1408-0.165171102,
+0.265407890
1410+0.174270958,
-0.610865094E-1
1411-0.610865094E-1,
+0.265407890
1416 lb(isam)
= ub(isam
- 1)
+ 1
1424call setRebound(cov, [
1_IK,
1_IK,
0_IK], [ndim, ndim, nsam])
1425call setRebound(mean, [
1_IK,
0_IK], [ndim, nsam])
1430+0.772733092,
-0.817067862,
-0.407605529,
+0.373222709,
-0.326901197,
+0.919963002,
+0.623701334,
-0.489149570,
+0.443457246,
+0.147727132
1431+0.656027436,
+0.260014057,
-0.371141911,
-0.570060730,
+0.347370505,
+0.853718519,
-0.352362990,
-0.440835953E-2,
-0.660302162,
+0.708751917
1432+0.796075106,
+0.390463591,
-0.400475025,
-0.634405136,
+0.764120221,
+0.519655824,
-0.290368795,
-0.171039104E-1,
+0.251866937,
+0.583487272
1435+1.04155898,
+1.73018765,
+1.70824099,
+1.05889463,
+1.19051361,
+1.37068951,
+1.58402085,
+1.47646117,
+1.14814258,
+1.93045866
1436cov(:,:,
0)
= getCov(sample,
2_IK, rweight)
1438+0.327239871,
+0.383769460E-1,
+0.211389828E-1
1439+0.383769460E-1,
+0.266241580,
+0.193893760
1440+0.211389828E-1,
+0.193893760,
+0.213558450
1441mean(:,
0)
= getMean(sample, dim, rweight)
1443+0.717712939E-1,
+0.107077792,
+0.189703077
1445 cov(:,:,isam)
= getCov(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1446 mean(:,isam)
= getMean(sample(:,lb(isam):ub(isam)),
2_IK, rweight(lb(isam):ub(isam)))
1448call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1450+0.327239841,
+0.383769386E-1,
+0.211389568E-1
1451-0.610865019E-1,
+0.266241550,
+0.193893746
1452+0.00000000,
+0.134451985,
+0.213558435
1453call setCovMeanMerged(covMerged, meanMerged, cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), lowDia)
1455+0.327239841,
+0.383769386E-1,
+0.211389568E-1
1456+0.383769386E-1,
+0.266241550,
+0.193893746
1457+0.211389568E-1,
+0.193893746,
+0.213558435
1458call setCovMeanMerged(cov(:,:,
2), mean(:,
2), cov(:,:,
1), mean(:,
1),
real(
sum(rweight(:ub(
1))), TKG)
/ real(
sum(rweight), TKG), uppDia)
1460+0.327239841,
+0.383769386E-1,
+0.211389568E-1
1461-0.308326874E-1,
+0.266241550,
+0.193893746
1462+0.625626892E-1,
+0.929234698E-1,
+0.213558435
1464+0.327239871,
+0.383769460E-1,
+0.211389828E-1
1465+0.383769460E-1,
+0.266241580,
+0.193893760
1466+0.211389828E-1,
+0.193893760,
+0.213558450