19 integer(IK) :: itry, ntry
= 10
20 type(display_type) :: disp
21 character(:),
allocatable :: format
25 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
26 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
27 call disp%show(
"!Convert correlation matrix and standard deviation to covariance matrix.")
28 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
29 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
37 real(TKG),
allocatable :: cov(:,:), cor(:,:), std(:)
38 format = getFormat(mold
= [
0._TKG], ed
= SK_
"es", signed
= .true._LK)
41 call disp%show(
"ndim = getUnifRand(1, 7)")
45 call disp%show(
"std = getUnifRand(1, 10, ndim)")
48 call disp%show( std ,
format = format )
49 call disp%show(
"cor = getCovRand(1., ndim)")
52 call disp%show( cor ,
format = format )
53 call disp%show(
"cov = getFilled(0._TKG, ndim, ndim + 1)")
55 call disp%show(
"call setMatCopy(cov(1:ndim, 1:ndim), rdpack, cor, rdpack, lowDia)")
56 call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
58 call disp%show( cov ,
format = format )
59 call disp%show(
"call setCov(cov(1:ndim, 2:ndim+1), uppDia, cov(1:ndim, 1:ndim), lowDia, std) ! convert lower-triangle correlation matrix to upper-triangle covariance matrix.")
60 call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
62 call disp%show( cov ,
format = format )
64 call disp%show(
"cov = getFilled(0._TKG, ndim, ndim + 1)")
66 call disp%show(
"call setMatCopy(cov(1:ndim, 2:ndim+1), rdpack, cor, rdpack, uppDia)")
67 call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
69 call disp%show( cov ,
format = format )
70 call disp%show(
"call setCov(cov(1:ndim, 1:ndim), lowDia, cov(1:ndim, 2:ndim+1), uppDia, std) ! convert upper-triangle correlation matrix to lower-triangle covariance matrix.")
71 call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
73 call disp%show( cov ,
format = format )
79 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
80 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
81 call disp%show(
"!Compute the covariance matrix of a 2-D sample.")
82 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
83 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
88 real(TKG),
allocatable :: sample(:,:), samShifted(:,:), cov(:,:), mean(:)
89 integer(IK) :: ndim, nsam, dim
90 call disp%show(
"ndim = 2; nsam = 10; dim = 2")
91 ndim
= 2; nsam
= 10; dim
= 2
92 call disp%show(
"sample = reshape(getUnifRand(1, 20, ndim * nsam), shape = [ndim, nsam], order = [2, 1])")
93 sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
96 call disp%show(
"mean = getMean(sample, dim)")
100 call disp%show(
"samShifted = getShifted(sample, dim, -mean)")
104 call disp%show(
"cov = getFilled(0., ndim, ndim)")
106 call disp%show(
"call setCov(cov, uppDia, mean, sample, dim)")
107 call setCov(cov, uppDia, mean, sample, dim)
110 call disp%show(
"cov = getFilled(0., ndim, ndim)")
112 call disp%show(
"call setCov(cov, uppDia, samShifted, dim) ! same result as above.")
113 call setCov(cov, uppDia, samShifted, dim)
117 call disp%show(
"cov = getFilled(0., ndim, ndim)")
119 call disp%show(
"call setCov(cov, lowDia, mean, sample, dim)")
120 call setCov(cov, lowDia, mean, sample, dim)
123 call disp%show(
"cov = getFilled(0., ndim, ndim)")
125 call disp%show(
"call setCov(cov, lowDia, samShifted, dim) ! same result as above.")
126 call setCov(cov, lowDia, samShifted, dim)
130 call disp%show(
"Compute the sample covariance along the first dimension.", deliml
= SK_
'''')
134 call disp%show(
"cov = getFilled(0., ndim, ndim)")
136 call disp%show(
"call setCov(cov, uppDia, mean, transpose(sample), dim)")
137 call setCov(cov, uppDia, mean,
transpose(sample), dim)
140 call disp%show(
"cov = getFilled(0., ndim, ndim)")
142 call disp%show(
"call setCov(cov, uppDia, transpose(samShifted), dim) ! same result as above.")
143 call setCov(cov, uppDia,
transpose(samShifted), dim)
147 call disp%show(
"cov = getFilled(0., ndim, ndim)")
149 call disp%show(
"call setCov(cov, lowDia, mean, transpose(sample), dim)")
150 call setCov(cov, lowDia, mean,
transpose(sample), dim)
153 call disp%show(
"cov = getFilled(0., ndim, ndim)")
155 call disp%show(
"call setCov(cov, lowDia, transpose(samShifted), dim) ! same result as above.")
156 call setCov(cov, lowDia,
transpose(samShifted), dim)
160 call disp%show(
"Compute the full sample covariance for a pair of time series.", deliml
= SK_
'''')
162 call disp%show(
"cov = getFilled(0., ndim, ndim)")
164 call disp%show(
"call setCov(cov, mean, sample(1,:), sample(2,:))")
165 call setCov(cov, mean, sample(
1,:), sample(
2,:))
168 call disp%show(
"cov = getFilled(0., ndim, ndim)")
170 call disp%show(
"call setCov(cov, samShifted(1,:), samShifted(2,:)) ! same result as above.")
171 call setCov(cov, samShifted(
1,:), samShifted(
2,:))
178 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
179 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
180 call disp%show(
"!Compute the biased covariance matrix of a weighted 2-D sample.")
181 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
182 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
188 integer(IK) :: iweisum
189 real(TKG),
allocatable :: rweight(:)
190 integer(IK),
allocatable :: iweight(:)
191 real(TKG),
allocatable :: sample(:,:), samShifted(:,:), cov(:,:), mean(:)
192 integer(IK) :: ndim, nsam, dim
193 call disp%show(
"ndim = 2; nsam = 10; dim = 2")
194 ndim
= 2; nsam
= 10; dim
= 2
195 call disp%show(
"sample = reshape(getUnifRand(1, 20, ndim * nsam), shape = [ndim, nsam], order = [2, 1])")
196 sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
199 call disp%show(
"call setResized(mean, ndim)")
201 call disp%show(
"iweight = getUnifRand(1, 10, nsam) ! integer-valued weights.")
205 call disp%show(
"call setMean(mean, sample, dim, iweight, iweisum)")
206 call setMean(mean, sample, dim, iweight, iweisum)
211 call disp%show(
"rweight = iweight ! or real-valued weights.")
215 call disp%show(
"call setMean(mean, sample, dim, rweight, rweisum)")
216 call setMean(mean, sample, dim, rweight, rweisum)
223 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
224 call disp%show(
"!Compute the covariance matrix integer weights.")
225 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
228 call disp%show(
"samShifted = getShifted(sample, dim, -mean)")
232 call disp%show(
"cov = getFilled(0., ndim, ndim)")
234 call disp%show(
"call setCov(cov, uppDia, mean, sample, dim, iweight, iweisum)")
235 call setCov(cov, uppDia, mean, sample, dim, iweight, iweisum)
238 call disp%show(
"cov = getFilled(0., ndim, ndim)")
240 call disp%show(
"call setCov(cov, uppDia, samShifted, dim, iweight, sum(iweight)) ! same result as above.")
241 call setCov(cov, uppDia, samShifted, dim, iweight,
sum(iweight))
245 call disp%show(
"cov = getFilled(0., ndim, ndim)")
247 call disp%show(
"call setCov(cov, lowDia, mean, sample, dim, iweight, iweisum)")
248 call setCov(cov, lowDia, mean, sample, dim, iweight, iweisum)
251 call disp%show(
"cov = getFilled(0., ndim, ndim)")
253 call disp%show(
"call setCov(cov, lowDia, samShifted, dim, iweight, sum(iweight)) ! same result as above.")
254 call setCov(cov, lowDia, samShifted, dim, iweight,
sum(iweight))
258 call disp%show(
"Compute the sample covariance along the first dimension.", deliml
= SK_
'''')
262 call disp%show(
"cov = getFilled(0., ndim, ndim)")
264 call disp%show(
"call setCov(cov, uppDia, mean, transpose(sample), dim, iweight, iweisum)")
265 call setCov(cov, uppDia, mean,
transpose(sample), dim, iweight, iweisum)
268 call disp%show(
"cov = getFilled(0., ndim, ndim)")
270 call disp%show(
"call setCov(cov, uppDia, transpose(samShifted), dim, iweight, sum(iweight)) ! same result as above.")
271 call setCov(cov, uppDia,
transpose(samShifted), dim, iweight,
sum(iweight))
275 call disp%show(
"cov = getFilled(0., ndim, ndim)")
277 call disp%show(
"call setCov(cov, lowDia, mean, transpose(sample), dim, iweight, iweisum)")
278 call setCov(cov, lowDia, mean,
transpose(sample), dim, iweight, iweisum)
281 call disp%show(
"cov = getFilled(0., ndim, ndim)")
283 call disp%show(
"call setCov(cov, lowDia, transpose(samShifted), dim, iweight, sum(iweight)) ! same result as above.")
284 call setCov(cov, lowDia,
transpose(samShifted), dim, iweight,
sum(iweight))
288 call disp%show(
"Compute the full sample covariance for a pair of time series.", deliml
= SK_
'''')
290 call disp%show(
"cov = getFilled(0., ndim, ndim)")
292 call disp%show(
"call setCov(cov, mean, sample(1,:), sample(2,:), iweight, iweisum)")
293 call setCov(cov, mean, sample(
1,:), sample(
2,:), iweight, iweisum)
296 call disp%show(
"cov = getFilled(0., ndim, ndim)")
298 call disp%show(
"call setCov(cov, samShifted(1,:), samShifted(2,:), iweight, sum(iweight)) ! same result as above.")
299 call setCov(cov, samShifted(
1,:), samShifted(
2,:), iweight,
sum(iweight))
305 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
306 call disp%show(
"!Compute the covariance matrix real weights.")
307 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
312 call disp%show(
"samShifted = getShifted(sample, dim, -mean)")
316 call disp%show(
"cov = getFilled(0., ndim, ndim)")
318 call disp%show(
"call setCov(cov, uppDia, mean, sample, dim, rweight, rweisum)")
319 call setCov(cov, uppDia, mean, sample, dim, rweight, rweisum)
322 call disp%show(
"cov = getFilled(0., ndim, ndim)")
324 call disp%show(
"call setCov(cov, uppDia, samShifted, dim, rweight, sum(rweight)) ! same result as above.")
325 call setCov(cov, uppDia, samShifted, dim, rweight,
sum(rweight))
329 call disp%show(
"cov = getFilled(0., ndim, ndim)")
331 call disp%show(
"call setCov(cov, lowDia, mean, sample, dim, rweight, rweisum)")
332 call setCov(cov, lowDia, mean, sample, dim, rweight, rweisum)
335 call disp%show(
"cov = getFilled(0., ndim, ndim)")
337 call disp%show(
"call setCov(cov, lowDia, samShifted, dim, rweight, sum(rweight)) ! same result as above.")
338 call setCov(cov, lowDia, samShifted, dim, rweight,
sum(rweight))
342 call disp%show(
"Compute the sample covariance along the first dimension.", deliml
= SK_
'''')
346 call disp%show(
"cov = getFilled(0., ndim, ndim)")
348 call disp%show(
"call setCov(cov, uppDia, mean, transpose(sample), dim, rweight, rweisum)")
349 call setCov(cov, uppDia, mean,
transpose(sample), dim, rweight, rweisum)
352 call disp%show(
"cov = getFilled(0., ndim, ndim)")
354 call disp%show(
"call setCov(cov, uppDia, transpose(samShifted), dim, rweight, sum(rweight)) ! same result as above.")
355 call setCov(cov, uppDia,
transpose(samShifted), dim, rweight,
sum(rweight))
359 call disp%show(
"cov = getFilled(0., ndim, ndim)")
361 call disp%show(
"call setCov(cov, lowDia, mean, transpose(sample), dim, rweight, rweisum)")
362 call setCov(cov, lowDia, mean,
transpose(sample), dim, rweight, rweisum)
365 call disp%show(
"cov = getFilled(0., ndim, ndim)")
367 call disp%show(
"call setCov(cov, lowDia, transpose(samShifted), dim, rweight, sum(rweight)) ! same result as above.")
368 call setCov(cov, lowDia,
transpose(samShifted), dim, rweight,
sum(rweight))
372 call disp%show(
"Compute the full sample covariance for a pair of time series.", deliml
= SK_
'''')
374 call disp%show(
"cov = getFilled(0., ndim, ndim)")
376 call disp%show(
"call setCov(cov, mean, sample(1,:), sample(2,:), rweight, rweisum)")
377 call setCov(cov, mean, sample(
1,:), sample(
2,:), rweight, rweisum)
380 call disp%show(
"cov = getFilled(0., ndim, ndim)")
382 call disp%show(
"call setCov(cov, samShifted(1,:), samShifted(2,:), rweight, sum(rweight)) ! same result as above.")
383 call setCov(cov, samShifted(
1,:), samShifted(
2,:), rweight,
sum(rweight))
Generate and return an array of the specified rank and shape of arbitrary intrinsic type and kind wit...
Allocate or resize (shrink or expand) an input allocatable scalar string or array of rank 1....
Generate count evenly spaced points over the interval [x1, x2] if x1 < x2, or [x2,...
Generate and return a random positive-definite (correlation or covariance) matrix using the Gram meth...
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.
Copy a desired subset of the input source matrix of arbitrary shape (:) or (:,:) to the target subset...
Generate and return the (weighted) mean of an input sample of nsam observations with ndim = 1 or 2 at...
Return the (weighted) mean of a pair of time series or of an input sample of nsam observations with n...
Generate a sample of shape (nsam), or (ndim, nsam) or (nsam, ndim) that is shifted by the specified i...
This module contains procedures and generic interfaces for convenient allocation and filling of array...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains procedures and generic interfaces for generating arrays with linear or logarithm...
This module contains classes and procedures for generating random matrices distributed on the space o...
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 LK
The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
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 procedures and generic interfaces relevant to copying (diagonal or upper/lower t...
This module contains classes and procedures for computing the first moment (i.e., the statistical mea...
This module contains classes and procedures for shifting univariate or multivariate samples by arbitr...
Generate and return an object of type display_type.
14+8.000000E+00,
+1.000000E+01,
+6.000000E+00,
+7.000000E+00,
+7.000000E+00,
+2.000000E+00
17+1.000000E+00,
-6.146030E-01,
+6.040904E-01,
+8.290617E-03,
-3.783300E-01,
+1.096828E-01
18-6.146030E-01,
+1.000000E+00,
-9.506434E-01,
-4.666466E-01,
-1.907399E-01,
-5.286536E-01
19+6.040904E-01,
-9.506434E-01,
+1.000000E+00,
+6.326929E-01,
+1.303558E-01,
+4.052165E-01
20+8.290617E-03,
-4.666466E-01,
+6.326929E-01,
+1.000000E+00,
+4.495044E-02,
-3.255494E-02
21-3.783300E-01,
-1.907399E-01,
+1.303558E-01,
+4.495044E-02,
+1.000000E+00,
+2.005435E-01
22+1.096828E-01,
-5.286536E-01,
+4.052165E-01,
-3.255494E-02,
+2.005435E-01,
+1.000000E+00
24call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
26+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
27-6.146030E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
28+6.040904E-01,
-9.506434E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
29+8.290617E-03,
-4.666466E-01,
+6.326929E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
30-3.783300E-01,
-1.907399E-01,
+1.303558E-01,
+4.495044E-02,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
31+1.096828E-01,
-5.286536E-01,
+4.052165E-01,
-3.255494E-02,
+2.005435E-01,
+1.000000E+00,
+0.000000E+00
32call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
34+1.000000E+00,
+6.400000E+01,
-4.916824E+01,
+2.899634E+01,
+4.642745E-01,
-2.118648E+01,
+1.754925E+00
35-6.146030E-01,
+1.000000E+00,
+1.000000E+02,
-5.703860E+01,
-3.266526E+01,
-1.335180E+01,
-1.057307E+01
36+6.040904E-01,
-9.506434E-01,
+1.000000E+00,
+3.600000E+01,
+2.657310E+01,
+5.474944E+00,
+4.862597E+00
37+8.290617E-03,
-4.666466E-01,
+6.326929E-01,
+1.000000E+00,
+4.900000E+01,
+2.202572E+00,
-4.557692E-01
38-3.783300E-01,
-1.907399E-01,
+1.303558E-01,
+4.495044E-02,
+1.000000E+00,
+4.900000E+01,
+2.807609E+00
39+1.096828E-01,
-5.286536E-01,
+4.052165E-01,
-3.255494E-02,
+2.005435E-01,
+1.000000E+00,
+4.000000E+00
42call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
44+0.000000E+00,
+1.000000E+00,
-6.146030E-01,
+6.040904E-01,
+8.290617E-03,
-3.783300E-01,
+1.096828E-01
45+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-9.506434E-01,
-4.666466E-01,
-1.907399E-01,
-5.286536E-01
46+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+6.326929E-01,
+1.303558E-01,
+4.052165E-01
47+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+4.495044E-02,
-3.255494E-02
48+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+2.005435E-01
49+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
50call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
52+6.400000E+01,
+1.000000E+00,
-6.146030E-01,
+6.040904E-01,
+8.290617E-03,
-3.783300E-01,
+1.096828E-01
53-4.916824E+01,
+1.000000E+02,
+1.000000E+00,
-9.506434E-01,
-4.666466E-01,
-1.907399E-01,
-5.286536E-01
54+2.899634E+01,
-5.703860E+01,
+3.600000E+01,
+1.000000E+00,
+6.326929E-01,
+1.303558E-01,
+4.052165E-01
55+4.642745E-01,
-3.266526E+01,
+2.657310E+01,
+4.900000E+01,
+1.000000E+00,
+4.495044E-02,
-3.255494E-02
56-2.118648E+01,
-1.335180E+01,
+5.474944E+00,
+2.202572E+00,
+4.900000E+01,
+1.000000E+00,
+2.005435E-01
57+1.754925E+00,
-1.057307E+01,
+4.862597E+00,
-4.557692E-01,
+2.807609E+00,
+4.000000E+00,
+1.000000E+00
65+2.000000E+00,
+1.000000E+00,
+1.000000E+01,
+5.000000E+00,
+8.000000E+00
68+1.000000E+00,
+8.009845E-01,
-2.793278E-01,
+3.857585E-01,
-5.876297E-01
69+8.009845E-01,
+1.000000E+00,
+2.067818E-02,
+4.578368E-01,
-6.384791E-01
70-2.793278E-01,
+2.067818E-02,
+1.000000E+00,
+4.020308E-01,
+4.063008E-01
71+3.857585E-01,
+4.578368E-01,
+4.020308E-01,
+1.000000E+00,
-4.943134E-01
72-5.876297E-01,
-6.384791E-01,
+4.063008E-01,
-4.943134E-01,
+1.000000E+00
74call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
76+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
77+8.009845E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
78-2.793278E-01,
+2.067818E-02,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
79+3.857585E-01,
+4.578368E-01,
+4.020308E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
80-5.876297E-01,
-6.384791E-01,
+4.063008E-01,
-4.943134E-01,
+1.000000E+00,
+0.000000E+00
81call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
83+1.000000E+00,
+4.000000E+00,
+1.601969E+00,
-5.586555E+00,
+3.857584E+00,
-9.402076E+00
84+8.009845E-01,
+1.000000E+00,
+1.000000E+00,
+2.067818E-01,
+2.289184E+00,
-5.107833E+00
85-2.793278E-01,
+2.067818E-02,
+1.000000E+00,
+1.000000E+02,
+2.010154E+01,
+3.250406E+01
86+3.857585E-01,
+4.578368E-01,
+4.020308E-01,
+1.000000E+00,
+2.500000E+01,
-1.977254E+01
87-5.876297E-01,
-6.384791E-01,
+4.063008E-01,
-4.943134E-01,
+1.000000E+00,
+6.400000E+01
90call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
92+0.000000E+00,
+1.000000E+00,
+8.009845E-01,
-2.793278E-01,
+3.857585E-01,
-5.876297E-01
93+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+2.067818E-02,
+4.578368E-01,
-6.384791E-01
94+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+4.020308E-01,
+4.063008E-01
95+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-4.943134E-01
96+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
97call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
99+4.000000E+00,
+1.000000E+00,
+8.009845E-01,
-2.793278E-01,
+3.857585E-01,
-5.876297E-01
100+1.601969E+00,
+1.000000E+00,
+1.000000E+00,
+2.067818E-02,
+4.578368E-01,
-6.384791E-01
101-5.586555E+00,
+2.067818E-01,
+1.000000E+02,
+1.000000E+00,
+4.020308E-01,
+4.063008E-01
102+3.857584E+00,
+2.289184E+00,
+2.010154E+01,
+2.500000E+01,
+1.000000E+00,
-4.943134E-01
103-9.402076E+00,
-5.107833E+00,
+3.250406E+01,
-1.977254E+01,
+6.400000E+01,
+1.000000E+00
116call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
118+1.000000E+00,
+0.000000E+00
119call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
121+1.000000E+00,
+4.900000E+01
124call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
126+0.000000E+00,
+1.000000E+00
127call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
129+4.900000E+01,
+1.000000E+00
137+3.000000E+00,
+9.000000E+00,
+2.000000E+00,
+9.000000E+00,
+5.000000E+00,
+8.000000E+00
140+1.000000E+00,
+2.053825E-01,
+3.161243E-01,
-1.263403E-03,
+2.475277E-01,
-3.331965E-01
141+2.053825E-01,
+1.000000E+00,
+7.257905E-01,
-1.875798E-02,
+4.492319E-01,
+4.462033E-01
142+3.161243E-01,
+7.257905E-01,
+1.000000E+00,
-4.131553E-01,
-3.410673E-02,
+5.687057E-01
143-1.263403E-03,
-1.875798E-02,
-4.131553E-01,
+1.000000E+00,
-1.294024E-01,
-1.136455E-01
144+2.475277E-01,
+4.492319E-01,
-3.410673E-02,
-1.294024E-01,
+1.000000E+00,
-1.098759E-01
145-3.331965E-01,
+4.462033E-01,
+5.687057E-01,
-1.136455E-01,
-1.098759E-01,
+1.000000E+00
147call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
149+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
150+2.053825E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
151+3.161243E-01,
+7.257905E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
152-1.263403E-03,
-1.875798E-02,
-4.131553E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
153+2.475277E-01,
+4.492319E-01,
-3.410673E-02,
-1.294024E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
154-3.331965E-01,
+4.462033E-01,
+5.687057E-01,
-1.136455E-01,
-1.098759E-01,
+1.000000E+00,
+0.000000E+00
155call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
157+1.000000E+00,
+9.000000E+00,
+5.545328E+00,
+1.896746E+00,
-3.411189E-02,
+3.712915E+00,
-7.996715E+00
158+2.053825E-01,
+1.000000E+00,
+8.100000E+01,
+1.306423E+01,
-1.519396E+00,
+2.021544E+01,
+3.212664E+01
159+3.161243E-01,
+7.257905E-01,
+1.000000E+00,
+4.000000E+00,
-7.436795E+00,
-3.410673E-01,
+9.099292E+00
160-1.263403E-03,
-1.875798E-02,
-4.131553E-01,
+1.000000E+00,
+8.100000E+01,
-5.823108E+00,
-8.182474E+00
161+2.475277E-01,
+4.492319E-01,
-3.410673E-02,
-1.294024E-01,
+1.000000E+00,
+2.500000E+01,
-4.395036E+00
162-3.331965E-01,
+4.462033E-01,
+5.687057E-01,
-1.136455E-01,
-1.098759E-01,
+1.000000E+00,
+6.400000E+01
165call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
167+0.000000E+00,
+1.000000E+00,
+2.053825E-01,
+3.161243E-01,
-1.263403E-03,
+2.475277E-01,
-3.331965E-01
168+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+7.257905E-01,
-1.875798E-02,
+4.492319E-01,
+4.462033E-01
169+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-4.131553E-01,
-3.410673E-02,
+5.687057E-01
170+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-1.294024E-01,
-1.136455E-01
171+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-1.098759E-01
172+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
173call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
175+9.000000E+00,
+1.000000E+00,
+2.053825E-01,
+3.161243E-01,
-1.263403E-03,
+2.475277E-01,
-3.331965E-01
176+5.545328E+00,
+8.100000E+01,
+1.000000E+00,
+7.257905E-01,
-1.875798E-02,
+4.492319E-01,
+4.462033E-01
177+1.896746E+00,
+1.306423E+01,
+4.000000E+00,
+1.000000E+00,
-4.131553E-01,
-3.410673E-02,
+5.687057E-01
178-3.411189E-02,
-1.519396E+00,
-7.436795E+00,
+8.100000E+01,
+1.000000E+00,
-1.294024E-01,
-1.136455E-01
179+3.712915E+00,
+2.021544E+01,
-3.410673E-01,
-5.823108E+00,
+2.500000E+01,
+1.000000E+00,
-1.098759E-01
180-7.996715E+00,
+3.212664E+01,
+9.099292E+00,
-8.182474E+00,
-4.395036E+00,
+6.400000E+01,
+1.000000E+00
188+9.000000E+00,
+4.000000E+00,
+9.000000E+00,
+6.000000E+00,
+1.000000E+01,
+7.000000E+00
191+1.000000E+00,
+8.902969E-01,
+7.304340E-01,
+5.993081E-01,
-7.672507E-01,
-6.349715E-01
192+8.902969E-01,
+1.000000E+00,
+3.396693E-01,
+3.598257E-01,
-6.172136E-01,
-5.159491E-01
193+7.304340E-01,
+3.396693E-01,
+1.000000E+00,
+6.754510E-01,
-6.397194E-01,
-5.528387E-01
194+5.993081E-01,
+3.598257E-01,
+6.754510E-01,
+1.000000E+00,
-9.078873E-01,
-1.553995E-01
195-7.672507E-01,
-6.172136E-01,
-6.397194E-01,
-9.078873E-01,
+1.000000E+00,
+1.096865E-01
196-6.349715E-01,
-5.159491E-01,
-5.528387E-01,
-1.553995E-01,
+1.096865E-01,
+1.000000E+00
198call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
200+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
201+8.902969E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
202+7.304340E-01,
+3.396693E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
203+5.993081E-01,
+3.598257E-01,
+6.754510E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
204-7.672507E-01,
-6.172136E-01,
-6.397194E-01,
-9.078873E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
205-6.349715E-01,
-5.159491E-01,
-5.528387E-01,
-1.553995E-01,
+1.096865E-01,
+1.000000E+00,
+0.000000E+00
206call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
208+1.000000E+00,
+8.100000E+01,
+3.205069E+01,
+5.916515E+01,
+3.236264E+01,
-6.905257E+01,
-4.000320E+01
209+8.902969E-01,
+1.000000E+00,
+1.600000E+01,
+1.222809E+01,
+8.635817E+00,
-2.468855E+01,
-1.444658E+01
210+7.304340E-01,
+3.396693E-01,
+1.000000E+00,
+8.100000E+01,
+3.647436E+01,
-5.757475E+01,
-3.482884E+01
211+5.993081E-01,
+3.598257E-01,
+6.754510E-01,
+1.000000E+00,
+3.600000E+01,
-5.447324E+01,
-6.526778E+00
212-7.672507E-01,
-6.172136E-01,
-6.397194E-01,
-9.078873E-01,
+1.000000E+00,
+1.000000E+02,
+7.678058E+00
213-6.349715E-01,
-5.159491E-01,
-5.528387E-01,
-1.553995E-01,
+1.096865E-01,
+1.000000E+00,
+4.900000E+01
216call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
218+0.000000E+00,
+1.000000E+00,
+8.902969E-01,
+7.304340E-01,
+5.993081E-01,
-7.672507E-01,
-6.349715E-01
219+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+3.396693E-01,
+3.598257E-01,
-6.172136E-01,
-5.159491E-01
220+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+6.754510E-01,
-6.397194E-01,
-5.528387E-01
221+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-9.078873E-01,
-1.553995E-01
222+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+1.096865E-01
223+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
224call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
226+8.100000E+01,
+1.000000E+00,
+8.902969E-01,
+7.304340E-01,
+5.993081E-01,
-7.672507E-01,
-6.349715E-01
227+3.205069E+01,
+1.600000E+01,
+1.000000E+00,
+3.396693E-01,
+3.598257E-01,
-6.172136E-01,
-5.159491E-01
228+5.916515E+01,
+1.222809E+01,
+8.100000E+01,
+1.000000E+00,
+6.754510E-01,
-6.397194E-01,
-5.528387E-01
229+3.236264E+01,
+8.635817E+00,
+3.647436E+01,
+3.600000E+01,
+1.000000E+00,
-9.078873E-01,
-1.553995E-01
230-6.905257E+01,
-2.468855E+01,
-5.757475E+01,
-5.447324E+01,
+1.000000E+02,
+1.000000E+00,
+1.096865E-01
231-4.000320E+01,
-1.444658E+01,
-3.482884E+01,
-6.526778E+00,
+7.678058E+00,
+4.900000E+01,
+1.000000E+00
239+2.000000E+00,
+1.000000E+01
242+1.000000E+00,
+6.533102E-01
243+6.533102E-01,
+1.000000E+00
245call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
247+1.000000E+00,
+0.000000E+00,
+0.000000E+00
248+6.533102E-01,
+1.000000E+00,
+0.000000E+00
249call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
251+1.000000E+00,
+4.000000E+00,
+1.306620E+01
252+6.533102E-01,
+1.000000E+00,
+1.000000E+02
255call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
257+0.000000E+00,
+1.000000E+00,
+6.533102E-01
258+0.000000E+00,
+0.000000E+00,
+1.000000E+00
259call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
261+4.000000E+00,
+1.000000E+00,
+6.533102E-01
262+1.306620E+01,
+1.000000E+02,
+1.000000E+00
270+4.000000E+00,
+2.000000E+00,
+4.000000E+00
273+1.000000E+00,
+2.043380E-01,
+4.009339E-01
274+2.043380E-01,
+1.000000E+00,
+9.075730E-01
275+4.009339E-01,
+9.075730E-01,
+1.000000E+00
277call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
279+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
280+2.043380E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
281+4.009339E-01,
+9.075730E-01,
+1.000000E+00,
+0.000000E+00
282call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
284+1.000000E+00,
+1.600000E+01,
+1.634704E+00,
+6.414942E+00
285+2.043380E-01,
+1.000000E+00,
+4.000000E+00,
+7.260584E+00
286+4.009339E-01,
+9.075730E-01,
+1.000000E+00,
+1.600000E+01
289call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
291+0.000000E+00,
+1.000000E+00,
+2.043380E-01,
+4.009339E-01
292+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+9.075730E-01
293+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
294call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
296+1.600000E+01,
+1.000000E+00,
+2.043380E-01,
+4.009339E-01
297+1.634704E+00,
+4.000000E+00,
+1.000000E+00,
+9.075730E-01
298+6.414942E+00,
+7.260584E+00,
+1.600000E+01,
+1.000000E+00
306+1.000000E+00,
+1.000000E+01,
+1.000000E+00
309+1.000000E+00,
+6.648872E-01,
-3.382753E-01
310+6.648872E-01,
+1.000000E+00,
-1.525386E-01
311-3.382753E-01,
-1.525386E-01,
+1.000000E+00
313call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
315+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
316+6.648872E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
317-3.382753E-01,
-1.525386E-01,
+1.000000E+00,
+0.000000E+00
318call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
320+1.000000E+00,
+1.000000E+00,
+6.648872E+00,
-3.382753E-01
321+6.648872E-01,
+1.000000E+00,
+1.000000E+02,
-1.525386E+00
322-3.382753E-01,
-1.525386E-01,
+1.000000E+00,
+1.000000E+00
325call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
327+0.000000E+00,
+1.000000E+00,
+6.648872E-01,
-3.382753E-01
328+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-1.525386E-01
329+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
330call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
332+1.000000E+00,
+1.000000E+00,
+6.648872E-01,
-3.382753E-01
333+6.648872E+00,
+1.000000E+02,
+1.000000E+00,
-1.525386E-01
334-3.382753E-01,
-1.525386E+00,
+1.000000E+00,
+1.000000E+00
342+8.000000E+00,
+4.000000E+00,
+1.000000E+01,
+8.000000E+00,
+5.000000E+00,
+4.000000E+00,
+9.000000E+00
345+1.000000E+00,
+9.621283E-01,
-6.292422E-01,
-7.141774E-01,
+3.926840E-01,
+2.214031E-01,
-4.519387E-01
346+9.621283E-01,
+1.000000E+00,
-7.547417E-01,
-8.737954E-01,
+4.880659E-01,
+1.197928E-01,
-5.189089E-01
347-6.292422E-01,
-7.547417E-01,
+1.000000E+00,
+8.021672E-01,
-7.937289E-01,
-2.737021E-01,
+4.582914E-01
348-7.141774E-01,
-8.737954E-01,
+8.021672E-01,
+1.000000E+00,
-5.158389E-01,
+1.213096E-01,
+5.098246E-01
349+3.926840E-01,
+4.880659E-01,
-7.937289E-01,
-5.158389E-01,
+1.000000E+00,
+9.855510E-02,
-5.040165E-01
350+2.214031E-01,
+1.197928E-01,
-2.737021E-01,
+1.213096E-01,
+9.855510E-02,
+1.000000E+00,
-1.514607E-01
351-4.519387E-01,
-5.189089E-01,
+4.582914E-01,
+5.098246E-01,
-5.040165E-01,
-1.514607E-01,
+1.000000E+00
353call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
355+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
356+9.621283E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
357-6.292422E-01,
-7.547417E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
358-7.141774E-01,
-8.737954E-01,
+8.021672E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
359+3.926840E-01,
+4.880659E-01,
-7.937289E-01,
-5.158389E-01,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00
360+2.214031E-01,
+1.197928E-01,
-2.737021E-01,
+1.213096E-01,
+9.855510E-02,
+1.000000E+00,
+0.000000E+00,
+0.000000E+00
361-4.519387E-01,
-5.189089E-01,
+4.582914E-01,
+5.098246E-01,
-5.040165E-01,
-1.514607E-01,
+1.000000E+00,
+0.000000E+00
362call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
364+1.000000E+00,
+6.400000E+01,
+3.078811E+01,
-5.033938E+01,
-4.570735E+01,
+1.570736E+01,
+7.084900E+00,
-3.253959E+01
365+9.621283E-01,
+1.000000E+00,
+1.600000E+01,
-3.018967E+01,
-2.796145E+01,
+9.761318E+00,
+1.916685E+00,
-1.868072E+01
366-6.292422E-01,
-7.547417E-01,
+1.000000E+00,
+1.000000E+02,
+6.417337E+01,
-3.968645E+01,
-1.094808E+01,
+4.124623E+01
367-7.141774E-01,
-8.737954E-01,
+8.021672E-01,
+1.000000E+00,
+6.400000E+01,
-2.063355E+01,
+3.881907E+00,
+3.670737E+01
368+3.926840E-01,
+4.880659E-01,
-7.937289E-01,
-5.158389E-01,
+1.000000E+00,
+2.500000E+01,
+1.971102E+00,
-2.268074E+01
369+2.214031E-01,
+1.197928E-01,
-2.737021E-01,
+1.213096E-01,
+9.855510E-02,
+1.000000E+00,
+1.600000E+01,
-5.452584E+00
370-4.519387E-01,
-5.189089E-01,
+4.582914E-01,
+5.098246E-01,
-5.040165E-01,
-1.514607E-01,
+1.000000E+00,
+8.100000E+01
373call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
375+0.000000E+00,
+1.000000E+00,
+9.621283E-01,
-6.292422E-01,
-7.141774E-01,
+3.926840E-01,
+2.214031E-01,
-4.519387E-01
376+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-7.547417E-01,
-8.737954E-01,
+4.880659E-01,
+1.197928E-01,
-5.189089E-01
377+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+8.021672E-01,
-7.937289E-01,
-2.737021E-01,
+4.582914E-01
378+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-5.158389E-01,
+1.213096E-01,
+5.098246E-01
379+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
+9.855510E-02,
-5.040165E-01
380+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00,
-1.514607E-01
381+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+0.000000E+00,
+1.000000E+00
382call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
384+6.400000E+01,
+1.000000E+00,
+9.621283E-01,
-6.292422E-01,
-7.141774E-01,
+3.926840E-01,
+2.214031E-01,
-4.519387E-01
385+3.078811E+01,
+1.600000E+01,
+1.000000E+00,
-7.547417E-01,
-8.737954E-01,
+4.880659E-01,
+1.197928E-01,
-5.189089E-01
386-5.033938E+01,
-3.018967E+01,
+1.000000E+02,
+1.000000E+00,
+8.021672E-01,
-7.937289E-01,
-2.737021E-01,
+4.582914E-01
387-4.570735E+01,
-2.796145E+01,
+6.417337E+01,
+6.400000E+01,
+1.000000E+00,
-5.158389E-01,
+1.213096E-01,
+5.098246E-01
388+1.570736E+01,
+9.761318E+00,
-3.968645E+01,
-2.063355E+01,
+2.500000E+01,
+1.000000E+00,
+9.855510E-02,
-5.040165E-01
389+7.084900E+00,
+1.916685E+00,
-1.094808E+01,
+3.881907E+00,
+1.971102E+00,
+1.600000E+01,
+1.000000E+00,
-1.514607E-01
390-3.253959E+01,
-1.868072E+01,
+4.124623E+01,
+3.670737E+01,
-2.268074E+01,
-5.452584E+00,
+8.100000E+01,
+1.000000E+00
403call setMatCopy(cov(
1:ndim,
1:ndim), rdpack, cor, rdpack, lowDia)
405+1.000000E+00,
+0.000000E+00
406call setCov(cov(
1:ndim,
2:ndim
+1), uppDia, cov(
1:ndim,
1:ndim), lowDia, std)
408+1.000000E+00,
+1.000000E+02
411call setMatCopy(cov(
1:ndim,
2:ndim
+1), rdpack, cor, rdpack, uppDia)
413+0.000000E+00,
+1.000000E+00
414call setCov(cov(
1:ndim,
1:ndim), lowDia, cov(
1:ndim,
2:ndim
+1), uppDia, std)
416+1.000000E+02,
+1.000000E+00
425ndim
= 2; nsam
= 10; dim
= 2
426sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
428+5.00000000,
+8.00000000,
+4.00000000,
+2.00000000,
+15.0000000,
+20.0000000,
+4.00000000,
+18.0000000,
+7.00000000,
+1.00000000
429+12.0000000,
+17.0000000,
+8.00000000,
+13.0000000,
+3.00000000,
+19.0000000,
+10.0000000,
+9.00000000,
+3.00000000,
+16.0000000
432+8.40000057,
+11.0000000
435-3.40000057,
-0.400000572,
-4.40000057,
-6.40000057,
+6.59999943,
+11.5999994,
-4.40000057,
+9.59999943,
-1.40000057,
-7.40000057
436+1.00000000,
+6.00000000,
-3.00000000,
+2.00000000,
-8.00000000,
+8.00000000,
-1.00000000,
-2.00000000,
-8.00000000,
+5.00000000
438call setCov(cov, uppDia, mean, sample, dim)
440+41.8400002,
-0.599999845
441+0.00000000,
+27.2000008
443call setCov(cov, uppDia, samShifted, dim)
445+41.8400002,
-0.599999845
446+0.00000000,
+27.2000008
449call setCov(cov, lowDia, mean, sample, dim)
451+41.8400002,
+0.00000000
452-0.599999845,
+27.2000008
454call setCov(cov, lowDia, samShifted, dim)
456+41.8400002,
+0.00000000
457-0.599999845,
+27.2000008
459'Compute the sample covariance along the first dimension.'
463call setCov(cov, uppDia, mean,
transpose(sample), dim)
465+41.8400002,
-0.599999845
466+0.00000000,
+27.2000008
468call setCov(cov, uppDia,
transpose(samShifted), dim)
470+41.8400002,
-0.599999845
471+0.00000000,
+27.2000008
474call setCov(cov, lowDia, mean,
transpose(sample), dim)
476+41.8400002,
+0.00000000
477-0.599999845,
+27.2000008
479call setCov(cov, lowDia,
transpose(samShifted), dim)
481+41.8400002,
+0.00000000
482-0.599999845,
+27.2000008
484'Compute the full sample covariance for a pair of time series.'
487call setCov(cov, mean, sample(
1,:), sample(
2,:))
489+41.8400002,
-0.599999845
490-0.599999845,
+27.2000008
492call setCov(cov, samShifted(
1,:), samShifted(
2,:))
494+41.8400002,
-0.599999845
495-0.599999845,
+27.2000008
504ndim
= 2; nsam
= 10; dim
= 2
505sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
507+18.0000000,
+11.0000000,
+5.00000000,
+18.0000000,
+17.0000000,
+13.0000000,
+13.0000000,
+9.00000000,
+9.00000000,
+17.0000000
508+16.0000000,
+11.0000000,
+12.0000000,
+19.0000000,
+5.00000000,
+17.0000000,
+15.0000000,
+20.0000000,
+11.0000000,
+16.0000000
512+6,
+3,
+7,
+1,
+7,
+7,
+1,
+1,
+3,
+2
513call setMean(mean, sample, dim, iweight, iweisum)
515+12.8157892,
+12.7894735
520+6,
+3,
+7,
+1,
+7,
+7,
+1,
+1,
+3,
+2
521call setMean(mean, sample, dim, rweight, rweisum)
523+12.8157892,
+12.7894735
533+5.18421078,
-1.81578922,
-7.81578922,
+5.18421078,
+4.18421078,
+0.184210777,
+0.184210777,
-3.81578922,
-3.81578922,
+4.18421078
534+3.21052647,
-1.78947353,
-0.789473534,
+6.21052647,
-7.78947353,
+4.21052647,
+2.21052647,
+7.21052647,
-1.78947353,
+3.21052647
536call setCov(cov, uppDia, mean, sample, dim, iweight, iweisum)
538+22.1502762,
-0.459834158
539+0.00000000,
+19.7451534
541call setCov(cov, uppDia, samShifted, dim, iweight,
sum(iweight))
543+22.1502762,
-0.459834158
544+0.00000000,
+19.7451534
547call setCov(cov, lowDia, mean, sample, dim, iweight, iweisum)
549+22.1502762,
+0.00000000
550-0.459833741,
+19.7451534
552call setCov(cov, lowDia, samShifted, dim, iweight,
sum(iweight))
554+22.1502762,
+0.00000000
555-0.459833741,
+19.7451534
557'Compute the sample covariance along the first dimension.'
561call setCov(cov, uppDia, mean,
transpose(sample), dim, iweight, iweisum)
563+22.1502762,
-0.459833443
564+0.00000000,
+19.7451534
566call setCov(cov, uppDia,
transpose(samShifted), dim, iweight,
sum(iweight))
568+22.1502762,
-0.459833443
569+0.00000000,
+19.7451534
572call setCov(cov, lowDia, mean,
transpose(sample), dim, iweight, iweisum)
574+22.1502762,
+0.00000000
575-0.459833443,
+19.7451534
577call setCov(cov, lowDia,
transpose(samShifted), dim, iweight,
sum(iweight))
579+22.1502762,
+0.00000000
580-0.459833443,
+19.7451534
582'Compute the full sample covariance for a pair of time series.'
585call setCov(cov, mean, sample(
1,:), sample(
2,:), iweight, iweisum)
587+22.1502762,
-0.459834158
588-0.459834158,
+19.7451534
590call setCov(cov, samShifted(
1,:), samShifted(
2,:), iweight,
sum(iweight))
592+22.1502762,
-0.459834158
593-0.459834158,
+19.7451534
603+5.18421078,
-1.81578922,
-7.81578922,
+5.18421078,
+4.18421078,
+0.184210777,
+0.184210777,
-3.81578922,
-3.81578922,
+4.18421078
604+3.21052647,
-1.78947353,
-0.789473534,
+6.21052647,
-7.78947353,
+4.21052647,
+2.21052647,
+7.21052647,
-1.78947353,
+3.21052647
606call setCov(cov, uppDia, mean, sample, dim, rweight, rweisum)
608+22.1502762,
-0.459834158
609+0.00000000,
+19.7451534
611call setCov(cov, uppDia, samShifted, dim, rweight,
sum(rweight))
613+22.1502762,
-0.459834158
614+0.00000000,
+19.7451534
617call setCov(cov, lowDia, mean, sample, dim, rweight, rweisum)
619+22.1502762,
+0.00000000
620-0.459833741,
+19.7451534
622call setCov(cov, lowDia, samShifted, dim, rweight,
sum(rweight))
624+22.1502762,
+0.00000000
625-0.459833741,
+19.7451534
627'Compute the sample covariance along the first dimension.'
631call setCov(cov, uppDia, mean,
transpose(sample), dim, rweight, rweisum)
633+22.1502762,
-0.459833443
634+0.00000000,
+19.7451534
636call setCov(cov, uppDia,
transpose(samShifted), dim, rweight,
sum(rweight))
638+22.1502762,
-0.459833443
639+0.00000000,
+19.7451534
642call setCov(cov, lowDia, mean,
transpose(sample), dim, rweight, rweisum)
644+22.1502762,
+0.00000000
645-0.459833443,
+19.7451534
647call setCov(cov, lowDia,
transpose(samShifted), dim, rweight,
sum(rweight))
649+22.1502762,
+0.00000000
650-0.459833443,
+19.7451534
652'Compute the full sample covariance for a pair of time series.'
655call setCov(cov, mean, sample(
1,:), sample(
2,:), rweight, rweisum)
657+22.1502762,
-0.459834158
658-0.459834158,
+19.7451534
660call setCov(cov, samShifted(
1,:), samShifted(
2,:), rweight,
sum(rweight))
662+22.1502762,
-0.459834158
663-0.459834158,
+19.7451534