23 type(display_type) :: disp
24 integer(IK) :: itry, ntry
= 10
25 character(:),
allocatable :: format
29 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
30 call disp%show(
"!Compute the correlation matrix from correlation matrix.")
31 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
36 real(TKG),
allocatable :: cov(:,:), cor(:,:), std(:)
38 format = getFormat(mold
= [
0._TKG], ed
= SK_
"es", signed
= .true._LK)
42 call disp%show(
"ndim = getUnifRand(0, 4)")
46 call disp%show(
"std = getUnifRand(1, ndim, ndim)")
49 call disp%show( std ,
format = format )
50 call disp%show(
"cov = getCovRand(1._TKG, std)")
53 call disp%show( cov ,
format = format )
60 call disp%show(
"call setCor(cor, uppDia, cor, uppDia)")
61 call setCor(cor, uppDia, cor, uppDia)
63 call disp%show( cor ,
format = format )
67 call disp%show(
"call setCor(cor, uppDia, cor, lowDia)")
68 call setCor(cor, uppDia, cor, lowDia)
70 call disp%show( cor ,
format = format )
74 call disp%show(
"call setCor(cor, upp, cor, uppDia)")
75 call setCor(cor, upp, cor, uppDia)
77 call disp%show( cor ,
format = format )
81 call disp%show(
"call setCor(cor, upp, cor, lowDia)")
82 call setCor(cor, upp, cor, lowDia)
84 call disp%show( cor ,
format = format )
88 call disp%show(
"call setCor(cor, lowDia, cor, uppDia)")
89 call setCor(cor, lowDia, cor, uppDia)
91 call disp%show( cor ,
format = format )
95 call disp%show(
"call setCor(cor, lowDia, cor, lowDia)")
96 call setCor(cor, lowDia, cor, lowDia)
98 call disp%show( cor ,
format = format )
102 call disp%show(
"call setCor(cor, low, cor, uppDia)")
103 call setCor(cor, low, cor, uppDia)
105 call disp%show( cor ,
format = format )
109 call disp%show(
"call setCor(cor, low, cor, lowDia)")
110 call setCor(cor, low, cor, lowDia)
112 call disp%show( cor ,
format = format )
119 call disp%show(
"call setCor(cor, uppDia, cor, upp, stdinv = 1 / std)")
120 call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
122 call disp%show( cor ,
format = format )
126 call disp%show(
"call setCor(cor, uppDia, cor, low, stdinv = 1 / std)")
127 call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
129 call disp%show( cor ,
format = format )
133 call disp%show(
"call setCor(cor, upp, cor, upp, stdinv = 1 / std)")
134 call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
140 call disp%show(
"call setCor(cor, upp, cor, low, stdinv = 1 / std)")
141 call setCor(cor, upp, cor, low, stdinv
= 1 / std)
143 call disp%show( cor ,
format = format )
147 call disp%show(
"call setCor(cor, lowDia, cor, upp, stdinv = 1 / std)")
148 call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
150 call disp%show( cor ,
format = format )
154 call disp%show(
"call setCor(cor, lowDia, cor, low, stdinv = 1 / std)")
155 call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
157 call disp%show( cor ,
format = format )
161 call disp%show(
"call setCor(cor, low, cor, upp, stdinv = 1 / std)")
162 call setCor(cor, low, cor, upp, stdinv
= 1 / std)
164 call disp%show( cor ,
format = format )
168 call disp%show(
"call setCor(cor, low, cor, low, stdinv = 1 / std)")
169 call setCor(cor, low, cor, low, stdinv
= 1 / std)
171 call disp%show( cor ,
format = format )
178 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
179 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
180 call disp%show(
"!Compute the Pearson correlation matrix for a pair of time series.")
181 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
182 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
187 integer(IK) :: ndim, nsam, dim
188 real(TKG),
allocatable :: sample(:,:), samShifted(:,:), cor(:,:), mean(:)
189 format = getFormat(mold
= [
0._TKG], ed
= SK_
"es", signed
= .true._LK)
190 call disp%show(
"ndim = 2; nsam = 10; dim = 2")
191 ndim
= 2; nsam
= 10; dim
= 2
192 call disp%show(
"sample = reshape(getUnifRand(1, 20, ndim * nsam), shape = [ndim, nsam], order = [2, 1])")
193 sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
195 call disp%show( sample ,
format = format )
196 call disp%show(
"mean = getMean(sample, dim)")
199 call disp%show( mean ,
format = format )
200 call disp%show(
"samShifted = getShifted(sample, dim, -mean)")
203 call disp%show( samShifted ,
format = format )
204 call disp%show(
"cor = getFilled(0., ndim, ndim)")
206 call disp%show(
"call setCor(cor, uppDia, mean, sample, dim)")
207 call setCor(cor, uppDia, mean, sample, dim)
209 call disp%show( cor ,
format = format )
210 call disp%show(
"cor = getFilled(0., ndim, ndim)")
212 call disp%show(
"call setCor(cor, uppDia, samShifted, dim) ! same result as above.")
213 call setCor(cor, uppDia, samShifted, dim)
215 call disp%show( cor ,
format = format )
217 call disp%show(
"cor = getFilled(0., ndim, ndim)")
219 call disp%show(
"call setCor(cor, lowDia, mean, sample, dim)")
220 call setCor(cor, lowDia, mean, sample, dim)
222 call disp%show( cor ,
format = format )
223 call disp%show(
"cor = getFilled(0., ndim, ndim)")
225 call disp%show(
"call setCor(cor, lowDia, samShifted, dim) ! same result as above.")
226 call setCor(cor, lowDia, samShifted, dim)
228 call disp%show( cor ,
format = format )
230 call disp%show(
"Compute the sample correlation along the first dimension.", deliml
= SK_
'''')
234 call disp%show(
"cor = getFilled(0., ndim, ndim)")
236 call disp%show(
"call setCor(cor, uppDia, mean, transpose(sample), dim)")
237 call setCor(cor, uppDia, mean,
transpose(sample), dim)
239 call disp%show( cor ,
format = format )
240 call disp%show(
"cor = getFilled(0., ndim, ndim)")
242 call disp%show(
"call setCor(cor, uppDia, transpose(samShifted), dim) ! same result as above.")
243 call setCor(cor, uppDia,
transpose(samShifted), dim)
245 call disp%show( cor ,
format = format )
247 call disp%show(
"cor = getFilled(0., ndim, ndim)")
249 call disp%show(
"call setCor(cor, lowDia, mean, transpose(sample), dim)")
250 call setCor(cor, lowDia, mean,
transpose(sample), dim)
252 call disp%show( cor ,
format = format )
253 call disp%show(
"cor = getFilled(0., ndim, ndim)")
255 call disp%show(
"call setCor(cor, lowDia, transpose(samShifted), dim) ! same result as above.")
256 call setCor(cor, lowDia,
transpose(samShifted), dim)
258 call disp%show( cor ,
format = format )
260 call disp%show(
"Compute the full sample correlation for a pair of time series.", deliml
= SK_
'''')
262 call disp%show(
"call setCor(cor(1,1), mean, sample(1,:), sample(2,:))")
263 call setCor(cor(
1,
1), mean, sample(
1,:), sample(
2,:))
265 call disp%show( cor(
1,
1) ,
format = format )
266 call disp%show(
"call setCor(cor(1,1), samShifted(1,:), samShifted(2,:)) ! same result as above.")
267 call setCor(cor(
1,
1), samShifted(
1,:), samShifted(
2,:))
269 call disp%show( cor(
1,
1) ,
format = format )
274 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
275 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
276 call disp%show(
"!Compute the Pearson correlation matrix for a weighted pair of time series.")
277 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
278 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
284 integer(IK),
allocatable :: iweight(:)
285 real(TKG),
allocatable :: rweight(:)
288 integer(IK) :: iweisum
289 integer(IK) :: ndim, nsam
290 real(TKG),
allocatable :: sample(:,:), samShifted(:,:), cor(:,:), mean(:)
291 format = getFormat(mold
= [
0._TKG], ed
= SK_
"es", signed
= .true._LK)
292 call disp%show(
"ndim = 2; nsam = 10; dim = 2")
293 ndim
= 2; nsam
= 10; dim
= 2
294 call disp%show(
"sample = reshape(getUnifRand(1, 20, ndim * nsam), shape = [ndim, nsam], order = [2, 1])")
295 sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
297 call disp%show( sample ,
format = format )
298 call disp%show(
"call setResized(mean, ndim)")
300 call disp%show(
"iweight = getUnifRand(1, 10, nsam) ! integer-valued weights.")
304 call disp%show(
"call setMean(mean, sample, dim, iweight, iweisum)")
305 call setMean(mean, sample, dim, iweight, iweisum)
307 call disp%show( mean ,
format = format )
310 call disp%show(
"rweight = iweight ! or real-valued weights.")
314 call disp%show(
"call setMean(mean, sample, dim, rweight, rweisum)")
315 call setMean(mean, sample, dim, rweight, rweisum)
317 call disp%show( mean ,
format = format )
319 call disp%show( rweisum ,
format = format )
322 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
323 call disp%show(
"!Compute the correlation matrix with integer weights.")
324 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
327 call disp%show(
"samShifted = getShifted(sample, dim, -mean)")
330 call disp%show( samShifted ,
format = format )
331 call disp%show(
"cor = getFilled(0., ndim, ndim)")
333 call disp%show(
"call setCor(cor, uppDia, mean, sample, dim, iweight, iweisum)")
334 call setCor(cor, uppDia, mean, sample, dim, iweight, iweisum)
336 call disp%show( cor ,
format = format )
337 call disp%show(
"cor = getFilled(0., ndim, ndim)")
339 call disp%show(
"call setCor(cor, uppDia, samShifted, dim, iweight, iweisum) ! same result as above.")
340 call setCor(cor, uppDia, samShifted, dim, iweight, iweisum)
342 call disp%show( cor ,
format = format )
344 call disp%show(
"cor = getFilled(0., ndim, ndim)")
346 call disp%show(
"call setCor(cor, lowDia, mean, sample, dim, iweight, iweisum)")
347 call setCor(cor, lowDia, mean, sample, dim, iweight, iweisum)
349 call disp%show( cor ,
format = format )
350 call disp%show(
"cor = getFilled(0., ndim, ndim)")
352 call disp%show(
"call setCor(cor, lowDia, samShifted, dim, iweight, iweisum) ! same result as above.")
353 call setCor(cor, lowDia, samShifted, dim, iweight, iweisum)
355 call disp%show( cor ,
format = format )
356 call disp%show(
"cor = getFilled(0., ndim, ndim)")
358 call disp%show(
"call setCor(cor, lowDia, getVerbose(samShifted, iweight, sum(iweight), dim), dim) ! same result as above.")
359 call setCor(cor, lowDia,
getVerbose(samShifted, iweight,
sum(iweight), dim), dim)
361 call disp%show( cor ,
format = format )
363 call disp%show(
"Compute the sample correlation along the first dimension.", deliml
= SK_
'''')
367 call disp%show(
"cor = getFilled(0., ndim, ndim)")
369 call disp%show(
"call setCor(cor, uppDia, mean, transpose(sample), dim, iweight, iweisum)")
370 call setCor(cor, uppDia, mean,
transpose(sample), dim, iweight, iweisum)
372 call disp%show( cor ,
format = format )
373 call disp%show(
"cor = getFilled(0., ndim, ndim)")
375 call disp%show(
"call setCor(cor, uppDia, transpose(samShifted), dim, iweight, iweisum) ! same result as above.")
376 call setCor(cor, uppDia,
transpose(samShifted), dim, iweight, iweisum)
378 call disp%show( cor ,
format = format )
380 call disp%show(
"cor = getFilled(0., ndim, ndim)")
382 call disp%show(
"call setCor(cor, lowDia, mean, transpose(sample), dim, iweight, iweisum)")
383 call setCor(cor, lowDia, mean,
transpose(sample), dim, iweight, iweisum)
385 call disp%show( cor ,
format = format )
386 call disp%show(
"cor = getFilled(0., ndim, ndim)")
388 call disp%show(
"call setCor(cor, lowDia, transpose(samShifted), dim, iweight, iweisum) ! same result as above.")
389 call setCor(cor, lowDia,
transpose(samShifted), dim, iweight, iweisum)
391 call disp%show( cor ,
format = format )
393 call disp%show(
"Compute the full sample correlation for a pair of time series.", deliml
= SK_
'''')
395 call disp%show(
"call setCor(cor(1,1), mean, sample(1,:), sample(2,:), iweight, iweisum)")
396 call setCor(cor(
1,
1), mean, sample(
1,:), sample(
2,:), iweight, iweisum)
398 call disp%show( cor(
1,
1) ,
format = format )
399 call disp%show(
"call setCor(cor(1,1), samShifted(1,:), samShifted(2,:), iweight, iweisum) ! same result as above.")
400 call setCor(cor(
1,
1), samShifted(
1,:), samShifted(
2,:), iweight, iweisum)
402 call disp%show( cor(
1,
1) ,
format = format )
406 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
407 call disp%show(
"!Compute the correlation matrix with real weights.")
408 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
413 call disp%show(
"samShifted = getShifted(sample, dim, -mean)")
417 call disp%show(
"cor = getFilled(0., ndim, ndim)")
419 call disp%show(
"call setCor(cor, uppDia, mean, sample, dim, rweight, rweisum)")
420 call setCor(cor, uppDia, mean, sample, dim, rweight, rweisum)
422 call disp%show( cor ,
format = format )
423 call disp%show(
"cor = getFilled(0., ndim, ndim)")
425 call disp%show(
"call setCor(cor, uppDia, samShifted, dim, rweight, rweisum) ! same result as above.")
426 call setCor(cor, uppDia, samShifted, dim, rweight, rweisum)
428 call disp%show( cor ,
format = format )
430 call disp%show(
"cor = getFilled(0., ndim, ndim)")
432 call disp%show(
"call setCor(cor, lowDia, mean, sample, dim, rweight, rweisum)")
433 call setCor(cor, lowDia, mean, sample, dim, rweight, rweisum)
435 call disp%show( cor ,
format = format )
436 call disp%show(
"cor = getFilled(0., ndim, ndim)")
438 call disp%show(
"call setCor(cor, lowDia, samShifted, dim, rweight, rweisum) ! same result as above.")
439 call setCor(cor, lowDia, samShifted, dim, rweight, rweisum)
441 call disp%show( cor ,
format = format )
443 call disp%show(
"Compute the sample correlation along the first dimension.", deliml
= SK_
'''')
447 call disp%show(
"cor = getFilled(0., ndim, ndim)")
449 call disp%show(
"call setCor(cor, uppDia, mean, transpose(sample), dim, rweight, rweisum)")
450 call setCor(cor, uppDia, mean,
transpose(sample), dim, rweight, rweisum)
452 call disp%show( cor ,
format = format )
453 call disp%show(
"cor = getFilled(0., ndim, ndim)")
455 call disp%show(
"call setCor(cor, uppDia, transpose(samShifted), dim, rweight, rweisum) ! same result as above.")
456 call setCor(cor, uppDia,
transpose(samShifted), dim, rweight, rweisum)
458 call disp%show( cor ,
format = format )
460 call disp%show(
"cor = getFilled(0., ndim, ndim)")
462 call disp%show(
"call setCor(cor, lowDia, mean, transpose(sample), dim, rweight, rweisum)")
463 call setCor(cor, lowDia, mean,
transpose(sample), dim, rweight, rweisum)
465 call disp%show( cor ,
format = format )
466 call disp%show(
"cor = getFilled(0., ndim, ndim)")
468 call disp%show(
"call setCor(cor, lowDia, transpose(samShifted), dim, rweight, rweisum) ! same result as above.")
469 call setCor(cor, lowDia,
transpose(samShifted), dim, rweight, rweisum)
471 call disp%show( cor ,
format = format )
473 call disp%show(
"Compute the full sample correlation for a pair of time series.", deliml
= SK_
'''')
475 call disp%show(
"call setCor(cor(1,1), mean, sample(1,:), sample(2,:), rweight, rweisum)")
476 call setCor(cor(
1,
1), mean, sample(
1,:), sample(
2,:), rweight, rweisum)
478 call disp%show( cor(
1,
1) ,
format = format )
479 call disp%show(
"call setCor(cor(1,1), samShifted(1,:), samShifted(2,:), rweight, rweisum) ! same result as above.")
480 call setCor(cor(
1,
1), samShifted(
1,:), samShifted(
2,:), rweight, rweisum)
482 call disp%show( cor(
1,
1) ,
format = format )
Generate and return an array of the specified rank and shape of arbitrary intrinsic type and kind wit...
Generate minimally-spaced character, integer, real sequences or sequences at fixed intervals of size ...
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 an equally-weighted (verbose or flattened) array of the input weighted array of rank 1 or 2.
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.
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 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 generating arrays with linear or logarithm...
This module contains procedures and generic interfaces for flattening (duplicating the elements of) a...
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 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.
12+2.000000E+00,
+2.000000E+00,
+3.000000E+00
15+4.000000E+00,
+2.731849E+00,
-3.342749E+00
16+2.731849E+00,
+4.000000E+00,
-2.704294E+00
17-3.342749E+00,
-2.704294E+00,
+8.999998E+00
21call setCor(cor, uppDia, cor, uppDia)
23+1.000000E+00,
+6.829623E-01,
-5.571249E-01
24+2.731849E+00,
+1.000000E+00,
-4.507158E-01
25-3.342749E+00,
-2.704294E+00,
+1.000000E+00
28call setCor(cor, uppDia, cor, lowDia)
30+1.000000E+00,
+6.829623E-01,
-5.571249E-01
31+2.731849E+00,
+1.000000E+00,
-4.507158E-01
32-3.342749E+00,
-2.704294E+00,
+1.000000E+00
35call setCor(cor, upp, cor, uppDia)
37+4.000000E+00,
+6.829623E-01,
-5.571249E-01
38+2.731849E+00,
+4.000000E+00,
-4.507158E-01
39-3.342749E+00,
-2.704294E+00,
+8.999998E+00
42call setCor(cor, upp, cor, lowDia)
44+4.000000E+00,
+6.829623E-01,
-5.571249E-01
45+2.731849E+00,
+4.000000E+00,
-4.507158E-01
46-3.342749E+00,
-2.704294E+00,
+8.999998E+00
49call setCor(cor, lowDia, cor, uppDia)
51+1.000000E+00,
+2.731849E+00,
-3.342749E+00
52+2.731849E+00,
+1.000000E+00,
-2.704294E+00
53-3.342749E+00,
-2.704294E+00,
+1.000000E+00
56call setCor(cor, lowDia, cor, lowDia)
58+1.000000E+00,
+2.731849E+00,
-3.342749E+00
59+6.829623E-01,
+1.000000E+00,
-2.704294E+00
60-5.571249E-01,
-4.507158E-01,
+1.000000E+00
63call setCor(cor, low, cor, uppDia)
65+4.000000E+00,
+2.731849E+00,
-3.342749E+00
66+6.829623E-01,
+4.000000E+00,
-2.704294E+00
67-5.571249E-01,
-4.507158E-01,
+8.999998E+00
70call setCor(cor, low, cor, lowDia)
72+4.000000E+00,
+2.731849E+00,
-3.342749E+00
73+6.829623E-01,
+4.000000E+00,
-2.704294E+00
74-5.571249E-01,
-4.507158E-01,
+8.999998E+00
78call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
80+1.000000E+00,
+6.829623E-01,
-5.571248E-01
81+2.731849E+00,
+1.000000E+00,
-4.507157E-01
82-3.342749E+00,
-2.704294E+00,
+1.000000E+00
85call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
87+1.000000E+00,
+6.829623E-01,
-5.571248E-01
88+2.731849E+00,
+1.000000E+00,
-4.507157E-01
89-3.342749E+00,
-2.704294E+00,
+1.000000E+00
92call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
94+4.00000000,
+0.682962298,
-0.557124794
95+2.73184919,
+4.00000000,
-0.450715721
96-3.34274864,
-2.70429420,
+8.99999809
99call setCor(cor, upp, cor, low, stdinv
= 1 / std)
101+4.000000E+00,
+6.829623E-01,
-5.571248E-01
102+2.731849E+00,
+4.000000E+00,
-4.507157E-01
103-3.342749E+00,
-2.704294E+00,
+8.999998E+00
106call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
108+1.000000E+00,
+2.731849E+00,
-3.342749E+00
109+6.829623E-01,
+1.000000E+00,
-2.704294E+00
110-5.571248E-01,
-4.507157E-01,
+1.000000E+00
113call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
115+1.000000E+00,
+2.731849E+00,
-3.342749E+00
116+6.829623E-01,
+1.000000E+00,
-2.704294E+00
117-5.571248E-01,
-4.507157E-01,
+1.000000E+00
120call setCor(cor, low, cor, upp, stdinv
= 1 / std)
122+4.000000E+00,
+2.731849E+00,
-3.342749E+00
123+6.829623E-01,
+4.000000E+00,
-2.704294E+00
124-5.571248E-01,
-4.507157E-01,
+8.999998E+00
127call setCor(cor, low, cor, low, stdinv
= 1 / std)
129+4.000000E+00,
+2.731849E+00,
-3.342749E+00
130+6.829623E-01,
+4.000000E+00,
-2.704294E+00
131-5.571248E-01,
-4.507157E-01,
+8.999998E+00
145call setCor(cor, uppDia, cor, uppDia)
149call setCor(cor, uppDia, cor, lowDia)
153call setCor(cor, upp, cor, uppDia)
157call setCor(cor, upp, cor, lowDia)
161call setCor(cor, lowDia, cor, uppDia)
165call setCor(cor, lowDia, cor, lowDia)
169call setCor(cor, low, cor, uppDia)
173call setCor(cor, low, cor, lowDia)
178call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
182call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
186call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
190call setCor(cor, upp, cor, low, stdinv
= 1 / std)
194call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
198call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
202call setCor(cor, low, cor, upp, stdinv
= 1 / std)
206call setCor(cor, low, cor, low, stdinv
= 1 / std)
215+2.000000E+00,
+4.000000E+00,
+4.000000E+00,
+4.000000E+00
218+4.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
219+6.281831E+00,
+1.600000E+01,
+1.212170E+01,
-1.554213E+01
220+6.033614E+00,
+1.212170E+01,
+1.600000E+01,
-1.003023E+01
221-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.600000E+01
225call setCor(cor, uppDia, cor, uppDia)
227+1.000000E+00,
+7.852287E-01,
+7.542017E-01,
-6.940523E-01
228+6.281831E+00,
+1.000000E+00,
+7.576063E-01,
-9.713830E-01
229+6.033614E+00,
+1.212170E+01,
+1.000000E+00,
-6.268895E-01
230-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.000000E+00
233call setCor(cor, uppDia, cor, lowDia)
235+1.000000E+00,
+7.852287E-01,
+7.542017E-01,
-6.940523E-01
236+6.281831E+00,
+1.000000E+00,
+7.576063E-01,
-9.713830E-01
237+6.033614E+00,
+1.212170E+01,
+1.000000E+00,
-6.268895E-01
238-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.000000E+00
241call setCor(cor, upp, cor, uppDia)
243+4.000000E+00,
+7.852287E-01,
+7.542017E-01,
-6.940523E-01
244+6.281831E+00,
+1.600000E+01,
+7.576063E-01,
-9.713830E-01
245+6.033614E+00,
+1.212170E+01,
+1.600000E+01,
-6.268895E-01
246-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.600000E+01
249call setCor(cor, upp, cor, lowDia)
251+4.000000E+00,
+7.852287E-01,
+7.542017E-01,
-6.940523E-01
252+6.281831E+00,
+1.600000E+01,
+7.576063E-01,
-9.713830E-01
253+6.033614E+00,
+1.212170E+01,
+1.600000E+01,
-6.268895E-01
254-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.600000E+01
257call setCor(cor, lowDia, cor, uppDia)
259+1.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
260+6.281831E+00,
+1.000000E+00,
+1.212170E+01,
-1.554213E+01
261+6.033614E+00,
+1.212170E+01,
+1.000000E+00,
-1.003023E+01
262-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.000000E+00
265call setCor(cor, lowDia, cor, lowDia)
267+1.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
268+7.852287E-01,
+1.000000E+00,
+1.212170E+01,
-1.554213E+01
269+7.542017E-01,
+7.576063E-01,
+1.000000E+00,
-1.003023E+01
270-6.940523E-01,
-9.713830E-01,
-6.268895E-01,
+1.000000E+00
273call setCor(cor, low, cor, uppDia)
275+4.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
276+7.852287E-01,
+1.600000E+01,
+1.212170E+01,
-1.554213E+01
277+7.542017E-01,
+7.576063E-01,
+1.600000E+01,
-1.003023E+01
278-6.940523E-01,
-9.713830E-01,
-6.268895E-01,
+1.600000E+01
281call setCor(cor, low, cor, lowDia)
283+4.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
284+7.852287E-01,
+1.600000E+01,
+1.212170E+01,
-1.554213E+01
285+7.542017E-01,
+7.576063E-01,
+1.600000E+01,
-1.003023E+01
286-6.940523E-01,
-9.713830E-01,
-6.268895E-01,
+1.600000E+01
290call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
292+1.000000E+00,
+7.852288E-01,
+7.542017E-01,
-6.940523E-01
293+6.281831E+00,
+1.000000E+00,
+7.576064E-01,
-9.713831E-01
294+6.033614E+00,
+1.212170E+01,
+1.000000E+00,
-6.268895E-01
295-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.000000E+00
298call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
300+1.000000E+00,
+7.852288E-01,
+7.542017E-01,
-6.940523E-01
301+6.281831E+00,
+1.000000E+00,
+7.576064E-01,
-9.713831E-01
302+6.033614E+00,
+1.212170E+01,
+1.000000E+00,
-6.268895E-01
303-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.000000E+00
306call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
308+4.00000000,
+0.785228848,
+0.754201710,
-0.694052339
309+6.28183079,
+16.0000038,
+0.757606387,
-0.971383095
310+6.03361368,
+12.1217022,
+16.0000019,
-0.626889467
311-5.55241871,
-15.5421295,
-10.0302315,
+16.0000019
314call setCor(cor, upp, cor, low, stdinv
= 1 / std)
316+4.000000E+00,
+7.852288E-01,
+7.542017E-01,
-6.940523E-01
317+6.281831E+00,
+1.600000E+01,
+7.576064E-01,
-9.713831E-01
318+6.033614E+00,
+1.212170E+01,
+1.600000E+01,
-6.268895E-01
319-5.552419E+00,
-1.554213E+01,
-1.003023E+01,
+1.600000E+01
322call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
324+1.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
325+7.852288E-01,
+1.000000E+00,
+1.212170E+01,
-1.554213E+01
326+7.542017E-01,
+7.576064E-01,
+1.000000E+00,
-1.003023E+01
327-6.940523E-01,
-9.713831E-01,
-6.268895E-01,
+1.000000E+00
330call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
332+1.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
333+7.852288E-01,
+1.000000E+00,
+1.212170E+01,
-1.554213E+01
334+7.542017E-01,
+7.576064E-01,
+1.000000E+00,
-1.003023E+01
335-6.940523E-01,
-9.713831E-01,
-6.268895E-01,
+1.000000E+00
338call setCor(cor, low, cor, upp, stdinv
= 1 / std)
340+4.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
341+7.852288E-01,
+1.600000E+01,
+1.212170E+01,
-1.554213E+01
342+7.542017E-01,
+7.576064E-01,
+1.600000E+01,
-1.003023E+01
343-6.940523E-01,
-9.713831E-01,
-6.268895E-01,
+1.600000E+01
346call setCor(cor, low, cor, low, stdinv
= 1 / std)
348+4.000000E+00,
+6.281831E+00,
+6.033614E+00,
-5.552419E+00
349+7.852288E-01,
+1.600000E+01,
+1.212170E+01,
-1.554213E+01
350+7.542017E-01,
+7.576064E-01,
+1.600000E+01,
-1.003023E+01
351-6.940523E-01,
-9.713831E-01,
-6.268895E-01,
+1.600000E+01
359+3.000000E+00,
+4.000000E+00,
+2.000000E+00,
+4.000000E+00
362+9.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
363+8.984529E+00,
+1.600000E+01,
-3.727107E+00,
+1.758899E+00
364+6.489483E-01,
-3.727107E+00,
+4.000000E+00,
+4.784607E+00
365+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.600000E+01
369call setCor(cor, uppDia, cor, uppDia)
371+1.000000E+00,
+7.487109E-01,
+1.081581E-01,
+3.671595E-01
372+8.984529E+00,
+1.000000E+00,
-4.658885E-01,
+1.099312E-01
373+6.489483E-01,
-3.727107E+00,
+1.000000E+00,
+5.980760E-01
374+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.000000E+00
377call setCor(cor, uppDia, cor, lowDia)
379+1.000000E+00,
+7.487109E-01,
+1.081581E-01,
+3.671595E-01
380+8.984529E+00,
+1.000000E+00,
-4.658885E-01,
+1.099312E-01
381+6.489483E-01,
-3.727107E+00,
+1.000000E+00,
+5.980760E-01
382+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.000000E+00
385call setCor(cor, upp, cor, uppDia)
387+9.000000E+00,
+7.487109E-01,
+1.081581E-01,
+3.671595E-01
388+8.984529E+00,
+1.600000E+01,
-4.658885E-01,
+1.099312E-01
389+6.489483E-01,
-3.727107E+00,
+4.000000E+00,
+5.980760E-01
390+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.600000E+01
393call setCor(cor, upp, cor, lowDia)
395+9.000000E+00,
+7.487109E-01,
+1.081581E-01,
+3.671595E-01
396+8.984529E+00,
+1.600000E+01,
-4.658885E-01,
+1.099312E-01
397+6.489483E-01,
-3.727107E+00,
+4.000000E+00,
+5.980760E-01
398+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.600000E+01
401call setCor(cor, lowDia, cor, uppDia)
403+1.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
404+8.984529E+00,
+1.000000E+00,
-3.727107E+00,
+1.758899E+00
405+6.489483E-01,
-3.727107E+00,
+1.000000E+00,
+4.784607E+00
406+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.000000E+00
409call setCor(cor, lowDia, cor, lowDia)
411+1.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
412+7.487109E-01,
+1.000000E+00,
-3.727107E+00,
+1.758899E+00
413+1.081581E-01,
-4.658885E-01,
+1.000000E+00,
+4.784607E+00
414+3.671595E-01,
+1.099312E-01,
+5.980760E-01,
+1.000000E+00
417call setCor(cor, low, cor, uppDia)
419+9.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
420+7.487109E-01,
+1.600000E+01,
-3.727107E+00,
+1.758899E+00
421+1.081581E-01,
-4.658885E-01,
+4.000000E+00,
+4.784607E+00
422+3.671595E-01,
+1.099312E-01,
+5.980760E-01,
+1.600000E+01
425call setCor(cor, low, cor, lowDia)
427+9.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
428+7.487109E-01,
+1.600000E+01,
-3.727107E+00,
+1.758899E+00
429+1.081581E-01,
-4.658885E-01,
+4.000000E+00,
+4.784607E+00
430+3.671595E-01,
+1.099312E-01,
+5.980760E-01,
+1.600000E+01
434call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
436+1.000000E+00,
+7.487108E-01,
+1.081581E-01,
+3.671595E-01
437+8.984529E+00,
+1.000000E+00,
-4.658884E-01,
+1.099312E-01
438+6.489483E-01,
-3.727107E+00,
+1.000000E+00,
+5.980759E-01
439+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.000000E+00
442call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
444+1.000000E+00,
+7.487108E-01,
+1.081581E-01,
+3.671595E-01
445+8.984529E+00,
+1.000000E+00,
-4.658884E-01,
+1.099312E-01
446+6.489483E-01,
-3.727107E+00,
+1.000000E+00,
+5.980759E-01
447+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.000000E+00
450call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
452+9.00000000,
+0.748710811,
+0.108158052,
+0.367159456
453+8.98452950,
+15.9999981,
-0.465888381,
+0.109931186
454+0.648948312,
-3.72710705,
+3.99999952,
+0.598075867
455+4.40591335,
+1.75889897,
+4.78460693,
+15.9999990
458call setCor(cor, upp, cor, low, stdinv
= 1 / std)
460+9.000000E+00,
+7.487108E-01,
+1.081581E-01,
+3.671595E-01
461+8.984529E+00,
+1.600000E+01,
-4.658884E-01,
+1.099312E-01
462+6.489483E-01,
-3.727107E+00,
+4.000000E+00,
+5.980759E-01
463+4.405913E+00,
+1.758899E+00,
+4.784607E+00,
+1.600000E+01
466call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
468+1.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
469+7.487108E-01,
+1.000000E+00,
-3.727107E+00,
+1.758899E+00
470+1.081581E-01,
-4.658884E-01,
+1.000000E+00,
+4.784607E+00
471+3.671595E-01,
+1.099312E-01,
+5.980759E-01,
+1.000000E+00
474call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
476+1.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
477+7.487108E-01,
+1.000000E+00,
-3.727107E+00,
+1.758899E+00
478+1.081581E-01,
-4.658884E-01,
+1.000000E+00,
+4.784607E+00
479+3.671595E-01,
+1.099312E-01,
+5.980759E-01,
+1.000000E+00
482call setCor(cor, low, cor, upp, stdinv
= 1 / std)
484+9.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
485+7.487108E-01,
+1.600000E+01,
-3.727107E+00,
+1.758899E+00
486+1.081581E-01,
-4.658884E-01,
+4.000000E+00,
+4.784607E+00
487+3.671595E-01,
+1.099312E-01,
+5.980759E-01,
+1.600000E+01
490call setCor(cor, low, cor, low, stdinv
= 1 / std)
492+9.000000E+00,
+8.984529E+00,
+6.489483E-01,
+4.405913E+00
493+7.487108E-01,
+1.600000E+01,
-3.727107E+00,
+1.758899E+00
494+1.081581E-01,
-4.658884E-01,
+4.000000E+00,
+4.784607E+00
495+3.671595E-01,
+1.099312E-01,
+5.980759E-01,
+1.600000E+01
503+2.000000E+00,
+2.000000E+00,
+1.000000E+00
506+4.000000E+00,
-2.356627E+00,
+6.051202E-01
507-2.356627E+00,
+4.000000E+00,
+1.020138E+00
508+6.051202E-01,
+1.020138E+00,
+1.000000E+00
512call setCor(cor, uppDia, cor, uppDia)
514+1.000000E+00,
-5.891567E-01,
+3.025601E-01
515-2.356627E+00,
+1.000000E+00,
+5.100691E-01
516+6.051202E-01,
+1.020138E+00,
+1.000000E+00
519call setCor(cor, uppDia, cor, lowDia)
521+1.000000E+00,
-5.891567E-01,
+3.025601E-01
522-2.356627E+00,
+1.000000E+00,
+5.100691E-01
523+6.051202E-01,
+1.020138E+00,
+1.000000E+00
526call setCor(cor, upp, cor, uppDia)
528+4.000000E+00,
-5.891567E-01,
+3.025601E-01
529-2.356627E+00,
+4.000000E+00,
+5.100691E-01
530+6.051202E-01,
+1.020138E+00,
+1.000000E+00
533call setCor(cor, upp, cor, lowDia)
535+4.000000E+00,
-5.891567E-01,
+3.025601E-01
536-2.356627E+00,
+4.000000E+00,
+5.100691E-01
537+6.051202E-01,
+1.020138E+00,
+1.000000E+00
540call setCor(cor, lowDia, cor, uppDia)
542+1.000000E+00,
-2.356627E+00,
+6.051202E-01
543-2.356627E+00,
+1.000000E+00,
+1.020138E+00
544+6.051202E-01,
+1.020138E+00,
+1.000000E+00
547call setCor(cor, lowDia, cor, lowDia)
549+1.000000E+00,
-2.356627E+00,
+6.051202E-01
550-5.891567E-01,
+1.000000E+00,
+1.020138E+00
551+3.025601E-01,
+5.100691E-01,
+1.000000E+00
554call setCor(cor, low, cor, uppDia)
556+4.000000E+00,
-2.356627E+00,
+6.051202E-01
557-5.891567E-01,
+4.000000E+00,
+1.020138E+00
558+3.025601E-01,
+5.100691E-01,
+1.000000E+00
561call setCor(cor, low, cor, lowDia)
563+4.000000E+00,
-2.356627E+00,
+6.051202E-01
564-5.891567E-01,
+4.000000E+00,
+1.020138E+00
565+3.025601E-01,
+5.100691E-01,
+1.000000E+00
569call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
571+1.000000E+00,
-5.891566E-01,
+3.025601E-01
572-2.356627E+00,
+1.000000E+00,
+5.100690E-01
573+6.051202E-01,
+1.020138E+00,
+1.000000E+00
576call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
578+1.000000E+00,
-5.891566E-01,
+3.025601E-01
579-2.356627E+00,
+1.000000E+00,
+5.100690E-01
580+6.051202E-01,
+1.020138E+00,
+1.000000E+00
583call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
585+4.00000000,
-0.589156628,
+0.302560091
586-2.35662651,
+3.99999952,
+0.510069013
587+0.605120182,
+1.02013803,
+1.00000000
590call setCor(cor, upp, cor, low, stdinv
= 1 / std)
592+4.000000E+00,
-5.891566E-01,
+3.025601E-01
593-2.356627E+00,
+4.000000E+00,
+5.100690E-01
594+6.051202E-01,
+1.020138E+00,
+1.000000E+00
597call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
599+1.000000E+00,
-2.356627E+00,
+6.051202E-01
600-5.891566E-01,
+1.000000E+00,
+1.020138E+00
601+3.025601E-01,
+5.100690E-01,
+1.000000E+00
604call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
606+1.000000E+00,
-2.356627E+00,
+6.051202E-01
607-5.891566E-01,
+1.000000E+00,
+1.020138E+00
608+3.025601E-01,
+5.100690E-01,
+1.000000E+00
611call setCor(cor, low, cor, upp, stdinv
= 1 / std)
613+4.000000E+00,
-2.356627E+00,
+6.051202E-01
614-5.891566E-01,
+4.000000E+00,
+1.020138E+00
615+3.025601E-01,
+5.100690E-01,
+1.000000E+00
618call setCor(cor, low, cor, low, stdinv
= 1 / std)
620+4.000000E+00,
-2.356627E+00,
+6.051202E-01
621-5.891566E-01,
+4.000000E+00,
+1.020138E+00
622+3.025601E-01,
+5.100690E-01,
+1.000000E+00
630+3.000000E+00,
+2.000000E+00,
+2.000000E+00
633+9.000000E+00,
+5.503739E+00,
+3.622772E+00
634+5.503739E+00,
+3.999999E+00,
+1.997328E+00
635+3.622772E+00,
+1.997328E+00,
+3.999999E+00
639call setCor(cor, uppDia, cor, uppDia)
641+1.000000E+00,
+9.172900E-01,
+6.037955E-01
642+5.503739E+00,
+1.000000E+00,
+4.993322E-01
643+3.622772E+00,
+1.997328E+00,
+1.000000E+00
646call setCor(cor, uppDia, cor, lowDia)
648+1.000000E+00,
+9.172900E-01,
+6.037955E-01
649+5.503739E+00,
+1.000000E+00,
+4.993322E-01
650+3.622772E+00,
+1.997328E+00,
+1.000000E+00
653call setCor(cor, upp, cor, uppDia)
655+9.000000E+00,
+9.172900E-01,
+6.037955E-01
656+5.503739E+00,
+3.999999E+00,
+4.993322E-01
657+3.622772E+00,
+1.997328E+00,
+3.999999E+00
660call setCor(cor, upp, cor, lowDia)
662+9.000000E+00,
+9.172900E-01,
+6.037955E-01
663+5.503739E+00,
+3.999999E+00,
+4.993322E-01
664+3.622772E+00,
+1.997328E+00,
+3.999999E+00
667call setCor(cor, lowDia, cor, uppDia)
669+1.000000E+00,
+5.503739E+00,
+3.622772E+00
670+5.503739E+00,
+1.000000E+00,
+1.997328E+00
671+3.622772E+00,
+1.997328E+00,
+1.000000E+00
674call setCor(cor, lowDia, cor, lowDia)
676+1.000000E+00,
+5.503739E+00,
+3.622772E+00
677+9.172900E-01,
+1.000000E+00,
+1.997328E+00
678+6.037955E-01,
+4.993322E-01,
+1.000000E+00
681call setCor(cor, low, cor, uppDia)
683+9.000000E+00,
+5.503739E+00,
+3.622772E+00
684+9.172900E-01,
+3.999999E+00,
+1.997328E+00
685+6.037955E-01,
+4.993322E-01,
+3.999999E+00
688call setCor(cor, low, cor, lowDia)
690+9.000000E+00,
+5.503739E+00,
+3.622772E+00
691+9.172900E-01,
+3.999999E+00,
+1.997328E+00
692+6.037955E-01,
+4.993322E-01,
+3.999999E+00
696call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
698+1.000000E+00,
+9.172899E-01,
+6.037954E-01
699+5.503739E+00,
+1.000000E+00,
+4.993321E-01
700+3.622772E+00,
+1.997328E+00,
+1.000000E+00
703call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
705+1.000000E+00,
+9.172899E-01,
+6.037954E-01
706+5.503739E+00,
+1.000000E+00,
+4.993321E-01
707+3.622772E+00,
+1.997328E+00,
+1.000000E+00
710call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
712+9.00000000,
+0.917289913,
+0.603795409
713+5.50373936,
+3.99999928,
+0.499332070
714+3.62277222,
+1.99732828,
+3.99999905
717call setCor(cor, upp, cor, low, stdinv
= 1 / std)
719+9.000000E+00,
+9.172899E-01,
+6.037954E-01
720+5.503739E+00,
+3.999999E+00,
+4.993321E-01
721+3.622772E+00,
+1.997328E+00,
+3.999999E+00
724call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
726+1.000000E+00,
+5.503739E+00,
+3.622772E+00
727+9.172899E-01,
+1.000000E+00,
+1.997328E+00
728+6.037954E-01,
+4.993321E-01,
+1.000000E+00
731call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
733+1.000000E+00,
+5.503739E+00,
+3.622772E+00
734+9.172899E-01,
+1.000000E+00,
+1.997328E+00
735+6.037954E-01,
+4.993321E-01,
+1.000000E+00
738call setCor(cor, low, cor, upp, stdinv
= 1 / std)
740+9.000000E+00,
+5.503739E+00,
+3.622772E+00
741+9.172899E-01,
+3.999999E+00,
+1.997328E+00
742+6.037954E-01,
+4.993321E-01,
+3.999999E+00
745call setCor(cor, low, cor, low, stdinv
= 1 / std)
747+9.000000E+00,
+5.503739E+00,
+3.622772E+00
748+9.172899E-01,
+3.999999E+00,
+1.997328E+00
749+6.037954E-01,
+4.993321E-01,
+3.999999E+00
757+2.000000E+00,
+1.000000E+00,
+2.000000E+00
760+4.000000E+00,
+7.658734E-01,
+2.311712E+00
761+7.658734E-01,
+1.000000E+00,
-9.103572E-02
762+2.311712E+00,
-9.103572E-02,
+4.000000E+00
766call setCor(cor, uppDia, cor, uppDia)
768+1.000000E+00,
+3.829367E-01,
+5.779281E-01
769+7.658734E-01,
+1.000000E+00,
-4.551786E-02
770+2.311712E+00,
-9.103572E-02,
+1.000000E+00
773call setCor(cor, uppDia, cor, lowDia)
775+1.000000E+00,
+3.829367E-01,
+5.779281E-01
776+7.658734E-01,
+1.000000E+00,
-4.551786E-02
777+2.311712E+00,
-9.103572E-02,
+1.000000E+00
780call setCor(cor, upp, cor, uppDia)
782+4.000000E+00,
+3.829367E-01,
+5.779281E-01
783+7.658734E-01,
+1.000000E+00,
-4.551786E-02
784+2.311712E+00,
-9.103572E-02,
+4.000000E+00
787call setCor(cor, upp, cor, lowDia)
789+4.000000E+00,
+3.829367E-01,
+5.779281E-01
790+7.658734E-01,
+1.000000E+00,
-4.551786E-02
791+2.311712E+00,
-9.103572E-02,
+4.000000E+00
794call setCor(cor, lowDia, cor, uppDia)
796+1.000000E+00,
+7.658734E-01,
+2.311712E+00
797+7.658734E-01,
+1.000000E+00,
-9.103572E-02
798+2.311712E+00,
-9.103572E-02,
+1.000000E+00
801call setCor(cor, lowDia, cor, lowDia)
803+1.000000E+00,
+7.658734E-01,
+2.311712E+00
804+3.829367E-01,
+1.000000E+00,
-9.103572E-02
805+5.779281E-01,
-4.551786E-02,
+1.000000E+00
808call setCor(cor, low, cor, uppDia)
810+4.000000E+00,
+7.658734E-01,
+2.311712E+00
811+3.829367E-01,
+1.000000E+00,
-9.103572E-02
812+5.779281E-01,
-4.551786E-02,
+4.000000E+00
815call setCor(cor, low, cor, lowDia)
817+4.000000E+00,
+7.658734E-01,
+2.311712E+00
818+3.829367E-01,
+1.000000E+00,
-9.103572E-02
819+5.779281E-01,
-4.551786E-02,
+4.000000E+00
823call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
825+1.000000E+00,
+3.829367E-01,
+5.779281E-01
826+7.658734E-01,
+1.000000E+00,
-4.551786E-02
827+2.311712E+00,
-9.103572E-02,
+1.000000E+00
830call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
832+1.000000E+00,
+3.829367E-01,
+5.779281E-01
833+7.658734E-01,
+1.000000E+00,
-4.551786E-02
834+2.311712E+00,
-9.103572E-02,
+1.000000E+00
837call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
839+4.00000000,
+0.382936716,
+0.577928066
840+0.765873432,
+1.00000000,
-0.455178618E-1
841+2.31171227,
-0.910357237E-1,
+4.00000000
844call setCor(cor, upp, cor, low, stdinv
= 1 / std)
846+4.000000E+00,
+3.829367E-01,
+5.779281E-01
847+7.658734E-01,
+1.000000E+00,
-4.551786E-02
848+2.311712E+00,
-9.103572E-02,
+4.000000E+00
851call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
853+1.000000E+00,
+7.658734E-01,
+2.311712E+00
854+3.829367E-01,
+1.000000E+00,
-9.103572E-02
855+5.779281E-01,
-4.551786E-02,
+1.000000E+00
858call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
860+1.000000E+00,
+7.658734E-01,
+2.311712E+00
861+3.829367E-01,
+1.000000E+00,
-9.103572E-02
862+5.779281E-01,
-4.551786E-02,
+1.000000E+00
865call setCor(cor, low, cor, upp, stdinv
= 1 / std)
867+4.000000E+00,
+7.658734E-01,
+2.311712E+00
868+3.829367E-01,
+1.000000E+00,
-9.103572E-02
869+5.779281E-01,
-4.551786E-02,
+4.000000E+00
872call setCor(cor, low, cor, low, stdinv
= 1 / std)
874+4.000000E+00,
+7.658734E-01,
+2.311712E+00
875+3.829367E-01,
+1.000000E+00,
-9.103572E-02
876+5.779281E-01,
-4.551786E-02,
+4.000000E+00
891call setCor(cor, uppDia, cor, uppDia)
896call setCor(cor, uppDia, cor, lowDia)
901call setCor(cor, upp, cor, uppDia)
906call setCor(cor, upp, cor, lowDia)
911call setCor(cor, lowDia, cor, uppDia)
916call setCor(cor, lowDia, cor, lowDia)
921call setCor(cor, low, cor, uppDia)
926call setCor(cor, low, cor, lowDia)
932call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
937call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
942call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
947call setCor(cor, upp, cor, low, stdinv
= 1 / std)
952call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
957call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
962call setCor(cor, low, cor, upp, stdinv
= 1 / std)
967call setCor(cor, low, cor, low, stdinv
= 1 / std)
977+2.000000E+00,
+2.000000E+00
980+4.000000E+00,
-1.403646E-01
981-1.403646E-01,
+3.999999E+00
985call setCor(cor, uppDia, cor, uppDia)
987+1.000000E+00,
-3.509115E-02
988-1.403646E-01,
+1.000000E+00
991call setCor(cor, uppDia, cor, lowDia)
993+1.000000E+00,
-3.509115E-02
994-1.403646E-01,
+1.000000E+00
997call setCor(cor, upp, cor, uppDia)
999+4.000000E+00,
-3.509115E-02
1000-1.403646E-01,
+3.999999E+00
1003call setCor(cor, upp, cor, lowDia)
1005+4.000000E+00,
-3.509115E-02
1006-1.403646E-01,
+3.999999E+00
1009call setCor(cor, lowDia, cor, uppDia)
1011+1.000000E+00,
-1.403646E-01
1012-1.403646E-01,
+1.000000E+00
1015call setCor(cor, lowDia, cor, lowDia)
1017+1.000000E+00,
-1.403646E-01
1018-3.509115E-02,
+1.000000E+00
1021call setCor(cor, low, cor, uppDia)
1023+4.000000E+00,
-1.403646E-01
1024-3.509115E-02,
+3.999999E+00
1027call setCor(cor, low, cor, lowDia)
1029+4.000000E+00,
-1.403646E-01
1030-3.509115E-02,
+3.999999E+00
1034call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
1036+1.000000E+00,
-3.509115E-02
1037-1.403646E-01,
+1.000000E+00
1040call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
1042+1.000000E+00,
-3.509115E-02
1043-1.403646E-01,
+1.000000E+00
1046call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
1048+4.00000000,
-0.350911506E-1
1049-0.140364602,
+3.99999928
1052call setCor(cor, upp, cor, low, stdinv
= 1 / std)
1054+4.000000E+00,
-3.509115E-02
1055-1.403646E-01,
+3.999999E+00
1058call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
1060+1.000000E+00,
-1.403646E-01
1061-3.509115E-02,
+1.000000E+00
1064call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
1066+1.000000E+00,
-1.403646E-01
1067-3.509115E-02,
+1.000000E+00
1070call setCor(cor, low, cor, upp, stdinv
= 1 / std)
1072+4.000000E+00,
-1.403646E-01
1073-3.509115E-02,
+3.999999E+00
1076call setCor(cor, low, cor, low, stdinv
= 1 / std)
1078+4.000000E+00,
-1.403646E-01
1079-3.509115E-02,
+3.999999E+00
1093call setCor(cor, uppDia, cor, uppDia)
1097call setCor(cor, uppDia, cor, lowDia)
1101call setCor(cor, upp, cor, uppDia)
1105call setCor(cor, upp, cor, lowDia)
1109call setCor(cor, lowDia, cor, uppDia)
1113call setCor(cor, lowDia, cor, lowDia)
1117call setCor(cor, low, cor, uppDia)
1121call setCor(cor, low, cor, lowDia)
1126call setCor(cor, uppDia, cor, upp, stdinv
= 1 / std)
1130call setCor(cor, uppDia, cor, low, stdinv
= 1 / std)
1134call setCor(cor, upp, cor, upp, stdinv
= 1 / std)
1138call setCor(cor, upp, cor, low, stdinv
= 1 / std)
1142call setCor(cor, lowDia, cor, upp, stdinv
= 1 / std)
1146call setCor(cor, lowDia, cor, low, stdinv
= 1 / std)
1150call setCor(cor, low, cor, upp, stdinv
= 1 / std)
1154call setCor(cor, low, cor, low, stdinv
= 1 / std)
1164ndim
= 2; nsam
= 10; dim
= 2
1165sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
1167+1.300000E+01,
+3.000000E+00,
+1.700000E+01,
+1.400000E+01,
+8.000000E+00,
+2.000000E+00,
+2.000000E+00,
+8.000000E+00,
+1.800000E+01,
+2.000000E+01
1168+1.000000E+01,
+9.000000E+00,
+1.500000E+01,
+1.700000E+01,
+1.700000E+01,
+9.000000E+00,
+6.000000E+00,
+5.000000E+00,
+1.000000E+01,
+1.000000E+01
1171+1.050000E+01,
+1.080000E+01
1174+2.500000E+00,
-7.500000E+00,
+6.500000E+00,
+3.500000E+00,
-2.500000E+00,
-8.500000E+00,
-8.500000E+00,
-2.500000E+00,
+7.500000E+00,
+9.500000E+00
1175-8.000002E-01,
-1.800000E+00,
+4.200000E+00,
+6.200000E+00,
+6.200000E+00,
-1.800000E+00,
-4.800000E+00,
-5.800000E+00,
-8.000002E-01,
-8.000002E-01
1177call setCor(cor, uppDia, mean, sample, dim)
1179+1.000000E+00,
+3.937320E-01
1180+0.000000E+00,
+1.000000E+00
1182call setCor(cor, uppDia, samShifted, dim)
1184+1.000000E+00,
+3.937320E-01
1185+0.000000E+00,
+1.000000E+00
1188call setCor(cor, lowDia, mean, sample, dim)
1190+1.000000E+00,
+0.000000E+00
1191+3.937320E-01,
+1.000000E+00
1193call setCor(cor, lowDia, samShifted, dim)
1195+1.000000E+00,
+0.000000E+00
1196+3.937320E-01,
+1.000000E+00
1198'Compute the sample correlation along the first dimension.'
1202call setCor(cor, uppDia, mean,
transpose(sample), dim)
1204+1.000000E+00,
+3.937320E-01
1205+0.000000E+00,
+1.000000E+00
1207call setCor(cor, uppDia,
transpose(samShifted), dim)
1209+1.000000E+00,
+3.937320E-01
1210+0.000000E+00,
+1.000000E+00
1213call setCor(cor, lowDia, mean,
transpose(sample), dim)
1215+1.000000E+00,
+0.000000E+00
1216+3.937320E-01,
+1.000000E+00
1218call setCor(cor, lowDia,
transpose(samShifted), dim)
1220+1.000000E+00,
+0.000000E+00
1221+3.937320E-01,
+1.000000E+00
1223'Compute the full sample correlation for a pair of time series.'
1225call setCor(cor(
1,
1), mean, sample(
1,:), sample(
2,:))
1228call setCor(cor(
1,
1), samShifted(
1,:), samShifted(
2,:))
1239ndim
= 2; nsam
= 10; dim
= 2
1240sample
= reshape(
getUnifRand(
1,
20, ndim
* nsam), shape
= [ndim, nsam], order
= [
2,
1])
1242+1.300000E+01,
+5.000000E+00,
+2.000000E+01,
+1.400000E+01,
+4.000000E+00,
+9.000000E+00,
+1.700000E+01,
+1.000000E+00,
+6.000000E+00,
+8.000000E+00
1243+2.000000E+00,
+6.000000E+00,
+1.000000E+01,
+1.700000E+01,
+2.000000E+00,
+1.200000E+01,
+1.300000E+01,
+4.000000E+00,
+1.700000E+01,
+1.200000E+01
1247+7,
+3,
+7,
+6,
+1,
+1,
+8,
+1,
+4,
+10
1248call setMean(mean, sample, dim, iweight, iweisum)
1250+1.216667E+01,
+1.070833E+01
1255+7,
+3,
+7,
+6,
+1,
+1,
+8,
+1,
+4,
+10
1256call setMean(mean, sample, dim, rweight, rweisum)
1258+1.216667E+01,
+1.070833E+01
1268+8.333330E-01,
-7.166667E+00,
+7.833333E+00,
+1.833333E+00,
-8.166667E+00,
-3.166667E+00,
+4.833333E+00,
-1.116667E+01,
-6.166667E+00,
-4.166667E+00
1269-8.708334E+00,
-4.708334E+00,
-7.083340E-01,
+6.291666E+00,
-8.708334E+00,
+1.291666E+00,
+2.291666E+00,
-6.708334E+00,
+6.291666E+00,
+1.291666E+00
1271call setCor(cor, uppDia, mean, sample, dim, iweight, iweisum)
1273+1.000000E+00,
+8.198504E-02
1274+0.000000E+00,
+1.000000E+00
1276call setCor(cor, uppDia, samShifted, dim, iweight, iweisum)
1278+1.000000E+00,
+8.198504E-02
1279+0.000000E+00,
+1.000000E+00
1282call setCor(cor, lowDia, mean, sample, dim, iweight, iweisum)
1284+1.000000E+00,
+0.000000E+00
1285+8.198504E-02,
+1.000000E+00
1287call setCor(cor, lowDia, samShifted, dim, iweight, iweisum)
1289+1.000000E+00,
+0.000000E+00
1290+8.198504E-02,
+1.000000E+00
1292call setCor(cor, lowDia,
getVerbose(samShifted, iweight,
sum(iweight), dim), dim)
1294+1.000000E+00,
+0.000000E+00
1295+8.198504E-02,
+1.000000E+00
1297'Compute the sample correlation along the first dimension.'
1301call setCor(cor, uppDia, mean,
transpose(sample), dim, iweight, iweisum)
1303+1.000000E+00,
+8.198504E-02
1304+0.000000E+00,
+1.000000E+00
1306call setCor(cor, uppDia,
transpose(samShifted), dim, iweight, iweisum)
1308+1.000000E+00,
+8.198504E-02
1309+0.000000E+00,
+1.000000E+00
1312call setCor(cor, lowDia, mean,
transpose(sample), dim, iweight, iweisum)
1314+1.000000E+00,
+0.000000E+00
1315+8.198504E-02,
+1.000000E+00
1317call setCor(cor, lowDia,
transpose(samShifted), dim, iweight, iweisum)
1319+1.000000E+00,
+0.000000E+00
1320+8.198504E-02,
+1.000000E+00
1322'Compute the full sample correlation for a pair of time series.'
1324call setCor(cor(
1,
1), mean, sample(
1,:), sample(
2,:), iweight, iweisum)
1327call setCor(cor(
1,
1), samShifted(
1,:), samShifted(
2,:), iweight, iweisum)
1339+0.833333015,
-7.16666698,
+7.83333302,
+1.83333302,
-8.16666698,
-3.16666698,
+4.83333302,
-11.1666670,
-6.16666698,
-4.16666698
1340-8.70833397,
-4.70833397,
-0.708333969,
+6.29166603,
-8.70833397,
+1.29166603,
+2.29166603,
-6.70833397,
+6.29166603,
+1.29166603
1342call setCor(cor, uppDia, mean, sample, dim, rweight, rweisum)
1344+1.000000E+00,
+8.198504E-02
1345+0.000000E+00,
+1.000000E+00
1347call setCor(cor, uppDia, samShifted, dim, rweight, rweisum)
1349+1.000000E+00,
+8.198504E-02
1350+0.000000E+00,
+1.000000E+00
1353call setCor(cor, lowDia, mean, sample, dim, rweight, rweisum)
1355+1.000000E+00,
+0.000000E+00
1356+8.198504E-02,
+1.000000E+00
1358call setCor(cor, lowDia, samShifted, dim, rweight, rweisum)
1360+1.000000E+00,
+0.000000E+00
1361+8.198504E-02,
+1.000000E+00
1363'Compute the sample correlation along the first dimension.'
1367call setCor(cor, uppDia, mean,
transpose(sample), dim, rweight, rweisum)
1369+1.000000E+00,
+8.198504E-02
1370+0.000000E+00,
+1.000000E+00
1372call setCor(cor, uppDia,
transpose(samShifted), dim, rweight, rweisum)
1374+1.000000E+00,
+8.198504E-02
1375+0.000000E+00,
+1.000000E+00
1378call setCor(cor, lowDia, mean,
transpose(sample), dim, rweight, rweisum)
1380+1.000000E+00,
+0.000000E+00
1381+8.198504E-02,
+1.000000E+00
1383call setCor(cor, lowDia,
transpose(samShifted), dim, rweight, rweisum)
1385+1.000000E+00,
+0.000000E+00
1386+8.198504E-02,
+1.000000E+00
1388'Compute the full sample correlation for a pair of time series.'
1390call setCor(cor(
1,
1), mean, sample(
1,:), sample(
2,:), rweight, rweisum)
1393call setCor(cor(
1,
1), samShifted(
1,:), samShifted(
2,:), rweight, rweisum)