ParaMonte Fortran 2.0.0
Parallel Monte Carlo and Machine Learning Library
See the latest version documentation.
pm_distCov::setCovRand Interface Reference

Return a random positive-definite power-law-distributed (correlation) matrix.
More...

Detailed Description

Return a random positive-definite power-law-distributed (correlation) matrix.

See the documentation of pm_distCov for details.

Parameters
[in,out]rng: The input/output scalar that can be an object of,
  1. type rngf_type, implying the use of intrinsic Fortran uniform RNG.
  2. type xoshiro256ssw_type, implying the use of xoshiro256** uniform RNG.
[out]rand: The output matrix of shape (1:ndim, 1:ndim) of,
  1. type complex of kind any supported by the processor (e.g., CK, CK32, CK64, or CK128),
  2. type real of kind any supported by the processor (e.g., RK, RK32, RK64, or RK128),
containing a random (optionally power-law-distributed determinant) positive-definite matrix.
The output rand can of complex type if and only if the optional input argument method is missing.
[in,out]method: The input/output scalar constant that can be one of the following:
  1. The scalar input constant dvine implying the use of the Dvine algorithm for generating random covariance matrices whose determinants are power-law distributed with exponent eta.
    In this case, the argument method has intent(in).
  2. A scalar output variable of type onion_type such as onion implying the use of the Onion algorithm for generating random covariance matrices whose determinants are power-law distributed with exponent eta.
    In this case, the argument method has intent(out).
    If the Cholesky factorization within the Onion algorithm fails, methodinfo will be set to the order of the leading minor of the specified input subset of mat that is not positive definite, indicating the occurrence of an error and that the factorization could not be completed.
    Otherwise, the info component of the onion method is set to 0.
The resulting matrix distribution from dvine and onion are identically distributed but onion method tends to have slightly faster runtime.
The larger eta is, the more the output random matrix looks like the Identity matrix.
Setting eta = 0. corresponds to a uniform distribution of the output matrix over the space of positive-definite correlation matrices.
See the description of the output argument rand for more information on the effects of eta on the off-diagonal elements of the output positive-definite matrix.
(optional. If missing the Gram method is used for random matrix generation. It must be missing for output rand of type complex.)
[in]eta: The input non-negative scalar of type real of the same kind as the output argument rand.
The larger eta is, the more the output random matrix looks like the Identity matrix.
Setting eta = 0. corresponds to a uniform distribution of the output matrix over the space of positive-definite correlation matrices.
See the description of the output argument rand for more information on the effects of eta on the off-diagonal elements of the output positive-definite matrix.
(optional. It must be present if and only if the input argument method is also present.)
[in]scale: The input scalar or contiguous vector of size ndim of type real of the same kind as the output argument rand, representing the scale of the matrix (e.g., the standard deviation of a covariance matrix) along each dimension.
(optional. default = 1.)


Possible calling interfaces

! Default (Gram) method.
call setCovRand(rng, rand(1:ndim, 1:ndim))
call setCovRand(rng, rand(1:ndim, 1:ndim), scale)
call setCovRand(rng, rand(1:ndim, 1:ndim), scale(1:ndim))
! Other methods.
call setCovRand(rng, rand(1:ndim, 1:ndim), method, eta)
call setCovRand(rng, rand(1:ndim, 1:ndim), method, eta, scale)
call setCovRand(rng, rand(1:ndim, 1:ndim), method, eta, scale(1:ndim))
Return a random positive-definite power-law-distributed (correlation) matrix.
Definition: pm_distCov.F90:787
This module contains classes and procedures for generating random matrices distributed on the space o...
Definition: pm_distCov.F90:72
Warning
The condition 0 <= eta must hold for the corresponding input arguments.
The condition all([0 < scale]) must hold for the corresponding input arguments.
The condition size(rand, 1) == size(rand, 2) must hold for the corresponding input arguments.
The condition rank(scale) == 0 .or. all(size(scale) == shape(rand)) must hold for the corresponding input arguments.
These conditions are verified only if the library is built with the preprocessor macro CHECK_ENABLED=1.
Beware that when the input argument scale is missing, the diagonal elements of the output correlation matrix are not enforced to match 1.
As such, numerical matrix multiplication errors may lead to diagonal matrix values slightly deviating from 1. If you need such a guarantee on the diagonal elements of the output random correlation matrix, use getCovRand.


Example usage

1program example
2
3 use pm_kind, only: SK
4 use pm_kind, only: IK, LK
5 use pm_io, only: field_type
6 use pm_io, only: display_type
7 use pm_matrixChol, only: setChoLow
10 use pm_distCov, only: setCovRand, dvine, onion
12 use pm_arrayResize, only: setResized
13 use pm_matrixDet, only: getMatDet
14
15 implicit none
16
17 integer(IK) :: itry, ndim
18
19 type(display_type) :: disp
20 disp = display_type(file = SK_"main.out.F90", format = field_type(complex = SK_"math"))
21
22 call disp%skip()
23 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
24 call disp%show("!Gram method for real covariance.")
25 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
26 call disp%skip()
27
28 block
29
30 use pm_kind, only: TKG => RKS ! all kinds are supported.
31 real(TKG), allocatable :: scale(:)
32 real(TKG), allocatable :: rand(:,:)
33 do itry = 1, 5
34
35 call disp%skip()
36 call disp%show("ndim = getUnifRand(3, 9)")
37 ndim = getUnifRand(3, 9)
38 call disp%show("ndim")
39 call disp%show( ndim )
40 call disp%show("call setResized(rand, [ndim, ndim])")
41 call setResized(rand, [ndim, ndim])
42 call disp%show("scale = getUnifRand(1, 10, ndim)")
43 scale = getUnifRand(1, 10, ndim)
44 call disp%show("scale")
45 call disp%show( scale )
46 call disp%skip()
47
48 call disp%show("call setCovRand(rngf, rand)")
49 call setCovRand(rngf, rand)
50 call disp%show("rand")
51 call disp%show( rand )
52 call disp%show("isMatClass(rand, posdefmat)")
53 call disp%show( isMatClass(rand, posdefmat) )
54 call disp%show("isMatClass(rand, hermitian)")
55 call disp%show( isMatClass(rand, hermitian) )
56 call disp%show("getMatDet(rand)")
57 call disp%show( getMatDet(rand) )
58
59 call disp%show("call setCovRand(rngf, rand, scale(1))")
60 call setCovRand(rngf, rand, scale(1))
61 call disp%show("rand")
62 call disp%show( rand )
63 call disp%show("isMatClass(rand, posdefmat)")
64 call disp%show( isMatClass(rand, posdefmat) )
65 call disp%show("isMatClass(rand, hermitian)")
66 call disp%show( isMatClass(rand, hermitian) )
67 call disp%show("getMatDet(rand)")
68 call disp%show( getMatDet(rand) )
69 call disp%skip()
70
71 call disp%show("call setCovRand(rngf, rand, scale)")
72 call setCovRand(rngf, rand, scale)
73 call disp%show("rand")
74 call disp%show( rand )
75 call disp%show("isMatClass(rand, posdefmat)")
76 call disp%show( isMatClass(rand, posdefmat) )
77 call disp%show("isMatClass(rand, hermitian)")
78 call disp%show( isMatClass(rand, hermitian) )
79 call disp%show("getMatDet(rand)")
80 call disp%show( getMatDet(rand) )
81 call disp%skip()
82
83 end do
84
85 end block
86
87 call disp%skip()
88 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
89 call disp%show("!Gram method for complex covariance.")
90 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
91 call disp%skip()
92
93 block
94
95 use pm_kind, only: TKG => RKD ! all kinds are supported.
96 real(TKG), allocatable :: scale(:)
97 complex(TKG), allocatable :: rand(:,:)
98 do itry = 1, 5
99
100 call disp%skip()
101 call disp%show("ndim = getUnifRand(3, 9)")
102 ndim = getUnifRand(3, 9)
103 call disp%show("ndim")
104 call disp%show( ndim )
105 call disp%show("call setResized(rand, [ndim, ndim])")
106 call setResized(rand, [ndim, ndim])
107 call disp%show("scale = getUnifRand(1, 10, ndim)")
108 scale = getUnifRand(1, 10, ndim)
109 call disp%show("scale")
110 call disp%show( scale )
111 call disp%skip()
112
113 call disp%show("call setCovRand(rngf, rand)")
114 call setCovRand(rngf, rand)
115 call disp%show("rand")
116 call disp%show( rand )
117 call disp%show("isMatClass(rand, posdefmat)")
118 call disp%show( isMatClass(rand, posdefmat) )
119 call disp%show("isMatClass(rand, hermitian)")
120 call disp%show( isMatClass(rand, hermitian) )
121 call disp%show("getMatDet(rand)")
122 call disp%show( getMatDet(rand) )
123
124 call disp%show("call setCovRand(rngf, rand, scale(1))")
125 call setCovRand(rngf, rand, scale(1))
126 call disp%show("rand")
127 call disp%show( rand )
128 call disp%show("isMatClass(rand, posdefmat)")
129 call disp%show( isMatClass(rand, posdefmat) )
130 call disp%show("isMatClass(rand, hermitian)")
131 call disp%show( isMatClass(rand, hermitian) )
132 call disp%show("getMatDet(rand)")
133 call disp%show( getMatDet(rand) )
134 call disp%skip()
135
136 call disp%show("call setCovRand(rngf, rand, scale)")
137 call setCovRand(rngf, rand, scale)
138 call disp%show("rand")
139 call disp%show( rand )
140 call disp%show("isMatClass(rand, posdefmat)")
141 call disp%show( isMatClass(rand, posdefmat) )
142 call disp%show("isMatClass(rand, hermitian)")
143 call disp%show( isMatClass(rand, hermitian) )
144 call disp%show("getMatDet(rand)")
145 call disp%show( getMatDet(rand) )
146 call disp%skip()
147
148 end do
149
150 end block
151
152 call disp%skip()
153 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%")
154 call disp%show("!Dvine and Onion methods.")
155 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%")
156 call disp%skip()
157
158 block
159
160 use pm_kind, only: TKG => RKS ! all kinds are supported.
161 real(TKG), allocatable :: rand(:,:)
162 real(TKG) :: eta, scale
163 do itry = 1, 10
164
165 call disp%skip()
166 call disp%show("eta = getUnifRand(1, 10)")
167 eta = getUnifRand(1, 10)
168 call disp%show("ndim = getUnifRand(2, 5)")
169 ndim = getUnifRand(2, 5)
170 call disp%show("call setResized(rand, [ndim, ndim])")
171 call setResized(rand, [ndim, ndim])
172 call disp%show("scale = getUnifRand(1, 10)")
173 scale = getUnifRand(1, 10)
174 call disp%skip()
175
176 call disp%show("call setCovRand(rngf, rand, dvine, eta)")
177 call setCovRand(rngf, rand, dvine, eta)
178 call disp%show("onion%info")
179 call disp%show( onion%info )
180 call disp%show("rand")
181 call disp%show( rand )
182 call disp%show("isMatClass(rand, posdefmat)")
183 call disp%show( isMatClass(rand, posdefmat) )
184 call disp%show("isMatClass(rand, hermitian)")
185 call disp%show( isMatClass(rand, hermitian) )
186 call disp%show("getMatDet(rand)")
187 call disp%show( getMatDet(rand) )
188
189 call disp%show("call setCovRand(rngf, rand, onion, eta)")
190 call setCovRand(rngf, rand, onion, eta)
191 call disp%show("onion%info")
192 call disp%show( onion%info )
193 call disp%show("rand")
194 call disp%show( rand )
195 call disp%show("isMatClass(rand, posdefmat)")
196 call disp%show( isMatClass(rand, posdefmat) )
197 call disp%show("isMatClass(rand, hermitian)")
198 call disp%show( isMatClass(rand, hermitian) )
199 call disp%show("getMatDet(rand)")
200 call disp%show( getMatDet(rand) )
201
202 call disp%show("call setCovRand(rngf, rand, dvine, eta, scale)")
203 call setCovRand(rngf, rand, dvine, eta, scale)
204 call disp%show("rand")
205 call disp%show( rand )
206 call disp%show("isMatClass(rand, posdefmat)")
207 call disp%show( isMatClass(rand, posdefmat) )
208 call disp%show("isMatClass(rand, hermitian)")
209 call disp%show( isMatClass(rand, hermitian) )
210 call disp%show("getMatDet(rand)")
211 call disp%show( getMatDet(rand) )
212 call disp%skip()
213
214 call disp%show("call setCovRand(rngf, rand, onion, eta, scale)")
215 call setCovRand(rngf, rand, onion, eta, scale)
216 call disp%show("onion%info")
217 call disp%show( onion%info )
218 call disp%show("rand")
219 call disp%show( rand )
220 call disp%show("isMatClass(rand, posdefmat)")
221 call disp%show( isMatClass(rand, posdefmat) )
222 call disp%show("isMatClass(rand, hermitian)")
223 call disp%show( isMatClass(rand, hermitian) )
224 call disp%show("getMatDet(rand)")
225 call disp%show( getMatDet(rand) )
226 call disp%skip()
227
228 end do
229
230 call disp%skip()
231 call disp%show("ndim = getUnifRand(2, 10)")
232 ndim = getUnifRand(2, 10)
233 call disp%show("call setResized(rand, [ndim, ndim])")
234 call setResized(rand, [ndim, ndim])
235 call disp%skip()
236
237 call disp%show("call setCovRand(rngf, rand, dvine, eta = 0._TKG, scale = [(real(itry, TKG), itry = 1, ndim)])")
238 call setCovRand(rngf, rand, dvine, eta = 0._TKG, scale = [(real(itry, TKG), itry = 1, ndim)])
239 call disp%show("rand")
240 call disp%show( rand )
241 call disp%show("isMatClass(rand, posdefmat)")
242 call disp%show( isMatClass(rand, posdefmat) )
243 call disp%skip()
244
245 call disp%show("call setCovRand(rngf, rand, onion, eta = 0._TKG, scale = [(real(itry, TKG), itry = 1, ndim)])")
246 call setCovRand(rngf, rand, onion, eta = 0._TKG, scale = [(real(itry, TKG), itry = 1, ndim)])
247 call disp%show("onion%info")
248 call disp%show( onion%info )
249 call disp%show("rand")
250 call disp%show( rand )
251 call disp%show("isMatClass(rand, posdefmat)")
252 call disp%show( isMatClass(rand, posdefmat) )
253 call disp%skip()
254
255 end block
256
257end program example
Allocate or resize (shrink or expand) an input allocatable scalar string or array of rank 1....
Generate and return a (collection) of random vector(s) of size ndim from the ndim-dimensional MultiVa...
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.
Definition: pm_io.F90:11726
This is a generic method of the derived type display_type with pass attribute.
Definition: pm_io.F90:11508
[LEGACY code] Return the lower-triangle of the Cholesky factorization of the symmetric positive-def...
Generate and return .true. if and only if the input matrix is of the specified input class.
Generate and return the determinant of the input general square matrix.
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
type(onion_type) onion
The scalar module variable object of type onion_type implying the use of the Onion algorithm for gene...
Definition: pm_distCov.F90:331
type(dvine_type), parameter dvine
The scalar constant of type dvine_type implying the use of the Dvine algorithm for generating random ...
Definition: pm_distCov.F90:238
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for computing various statistical quantities related to t...
type(rngf_type) rngf
The scalar constant object of type rngf_type whose presence signified the use of the Fortran intrinsi...
This module contains classes and procedures for input/output (IO) or generic display operations on st...
Definition: pm_io.F90:252
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
Definition: pm_io.F90:11393
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
Definition: pm_kind.F90:268
integer, parameter LK
The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
Definition: pm_kind.F90:541
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
Definition: pm_kind.F90:540
integer, parameter RKD
The double precision real kind in Fortran mode. On most platforms, this is an 64-bit real kind.
Definition: pm_kind.F90:568
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
Definition: pm_kind.F90:539
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
Definition: pm_kind.F90:567
This module contains procedures and generic interfaces for computing the Cholesky factorization of po...
This module contains abstract and concrete derived types that are required for compile-time resolutio...
type(posdefmat_type), parameter posdefmat
This is a scalar parameter object of type hermitian_type that is exclusively used to signify the Herm...
type(hermitian_type), parameter hermitian
This is a scalar parameter object of type hermitian_type that is exclusively used to signify the Herm...
This module contains procedures and generic interfaces relevant to the computation of the determinant...
Generate and return an object of type display_type.
Definition: pm_io.F90:10282
The derived type that can be used for constructing containers of format or left and right delimiters ...
Definition: pm_io.F90:482

Example Unix compile command via Intel ifort compiler
1#!/usr/bin/env sh
2rm main.exe
3ifort -fpp -standard-semantics -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte* -o main.exe
4./main.exe

Example Windows Batch compile command via Intel ifort compiler
1del main.exe
2set PATH=..\..\..\lib;%PATH%
3ifort /fpp /standard-semantics /O3 /I:..\..\..\include main.F90 ..\..\..\lib\libparamonte*.lib /exe:main.exe
4main.exe

Example Unix / MinGW compile command via GNU gfortran compiler
1#!/usr/bin/env sh
2rm main.exe
3gfortran -cpp -ffree-line-length-none -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte* -o main.exe
4./main.exe

Example output
1
2!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3!Gram method for real covariance.
4!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5
6
7ndim = getUnifRand(3, 9)
8ndim
9+9
10call setResized(rand, [ndim, ndim])
11scale = getUnifRand(1, 10, ndim)
12scale
13+5.00000000, +8.00000000, +2.00000000, +1.00000000, +10.0000000, +1.00000000, +5.00000000, +9.00000000, +10.0000000
14
15call setCovRand(rngf, rand)
16rand
17+0.999999881, -0.739881396, -0.804915726, +0.538025916, +0.398764789, +0.309528589, +0.442448050, +0.493463755, +0.441832125
18-0.739881396, +0.999999881, +0.296902210, -0.724600196, -0.263106436, +0.152529180, -0.675511360E-1, -0.668727100, -0.623100400
19-0.804915726, +0.296902210, +0.999999940, -0.341024101, -0.350007147, -0.690403700, -0.535106003, -0.102386631, -0.482003391E-2
20+0.538025916, -0.724600196, -0.341024101, +1.00000000, -0.101771802, -0.273033798, -0.159997135, +0.950107276E-1, +0.472058654
21+0.398764789, -0.263106436, -0.350007147, -0.101771802, +1.00000012, +0.497617096, +0.173769295, +0.666463375, +0.261322558
22+0.309528589, +0.152529180, -0.690403700, -0.273033798, +0.497617096, +1.00000012, +0.652712524, +0.649232343E-1, -0.516765356
23+0.442448050, -0.675511360E-1, -0.535106003, -0.159997135, +0.173769295, +0.652712524, +1.00000000, -0.101931170E-1, -0.264938653
24+0.493463755, -0.668727100, -0.102386631, +0.950107276E-1, +0.666463375, +0.649232343E-1, -0.101931170E-1, +1.00000000, +0.405425191
25+0.441832125, -0.623100400, -0.482003391E-2, +0.472058654, +0.261322558, -0.516765356, -0.264938653, +0.405425191, +1.00000000
27T
29T
30getMatDet(rand)
31+0.234202190E-6
32call setCovRand(rngf, rand, scale(1))
33rand
34+61.0346375, -10.0135965, +16.9959717, +8.01034546, -5.05790472, -27.6150284, -17.2427444, +30.9637642, +1.74493396
35-10.0135965, +96.9907684, +2.80410504, +26.5298996, +4.88708305, -0.655987859, +25.8009071, -11.9765596, -2.40745330
36+16.9959717, +2.80410504, +75.4460678, +24.5599308, -42.5905991, +13.8779697, +11.9145842, +20.8507996, -21.9599953
37+8.01034546, +26.5298996, +24.5599308, +56.4027023, -18.4717484, +3.13681293, +8.51517487, +2.92348003, +3.38907242
38-5.05790472, +4.88708305, -42.5905991, -18.4717484, +39.8450775, -7.24597883, -8.25260830, +0.948182464, +7.83119440
39-27.6150284, -0.655987859, +13.8779697, +3.13681293, -7.24597883, +37.6607590, +8.24943066, -6.02743912, -17.7582092
40-17.2427444, +25.8009071, +11.9145842, +8.51517487, -8.25260830, +8.24943066, +34.5391998, -15.7915573, +7.22107363
41+30.9637642, -11.9765596, +20.8507996, +2.92348003, +0.948182464, -6.02743912, -15.7915573, +31.5244446, -10.7031727
42+1.74493396, -2.40745330, -21.9599953, +3.38907242, +7.83119440, -17.7582092, +7.22107363, -10.7031727, +25.0000038
44T
46T
47getMatDet(rand)
48+0.423529808E+12
49
50call setCovRand(rngf, rand, scale)
51rand
52+40.5061951, -23.0232124, +4.13454056, +2.28239012, +2.36792278, +3.78461409, -1.11944640, +13.0608749, +20.6375237
53-23.0232124, +74.0647736, -20.3389320, -5.65863943, -54.4920540, +1.56228399, +6.09199953, -5.95409727, -12.9539738
54+4.13454056, -20.3389320, +19.7489681, +2.94602442, +103.552551, -6.32460117, -14.0377922, +43.6238136, +7.50013828
55+2.28239012, -5.65863943, +2.94602442, +1.25752735, +14.5796700, -0.347866327, +0.419407785, +1.32313454, +1.36859727
56+2.36792278, -54.4920540, +103.552551, +14.5796700, +690.292480, -37.4460602, -51.5100708, +227.427185, -25.1074066
57+3.78461409, +1.56228399, -6.32460117, -0.347866327, -37.4460602, +4.87174797, +1.33646905, -22.8349628, -2.51463223
58-1.11944640, +6.09199953, -14.0377922, +0.419407785, -51.5100708, +1.33646905, +45.4131699, -30.8680153, -4.68578482
59+13.0608749, -5.95409727, +43.6238136, +1.32313454, +227.427185, -22.8349628, -30.8680153, +247.720123, +70.3257294
60+20.6375237, -12.9539738, +7.50013828, +1.36859727, -25.1074066, -2.51463223, -4.68578482, +70.3257294, +99.9999924
62T
64T
65getMatDet(rand)
66+0.616349082E+11
67
68
69ndim = getUnifRand(3, 9)
70ndim
71+7
72call setResized(rand, [ndim, ndim])
73scale = getUnifRand(1, 10, ndim)
74scale
75+5.00000000, +3.00000000, +1.00000000, +1.00000000, +6.00000000, +4.00000000, +10.0000000
76
77call setCovRand(rngf, rand)
78rand
79+1.00000000, +0.869409561, -0.749419808, -0.195486188, -0.572602272, -0.416085541, +0.354812711
80+0.869409561, +1.00000000, -0.871373296, +0.513090491E-1, -0.349975944, -0.624143958, +0.109967470
81-0.749419808, -0.871373296, +0.999999881, -0.312165260, +0.166473687E-1, +0.312029839, -0.193620026E-1
82-0.195486188, +0.513090491E-1, -0.312165260, +1.00000000, +0.434048206, -0.865594447E-1, -0.234513819
83-0.572602272, -0.349975944, +0.166473687E-1, +0.434048206, +0.999999940, +0.630527973, -0.154537350
84-0.416085541, -0.624143958, +0.312029839, -0.865594447E-1, +0.630527973, +1.00000000, +0.266300499
85+0.354812711, +0.109967470, -0.193620026E-1, -0.234513819, -0.154537350, +0.266300499, +1.00000012
87T
89T
90getMatDet(rand)
91+0.990256103E-5
92call setCovRand(rngf, rand, scale(1))
93rand
94+61.0848732, -14.3712673, +37.0940132, +44.5056076, +8.52080345, +4.14040899, -3.01028705
95-14.3712673, +104.208389, -24.4708366, -21.8200912, +36.9947319, +14.5860672, -7.01430130
96+37.0940132, -24.4708366, +81.0135498, +44.8614311, +3.33841825, +13.2448359, -1.96349764
97+44.5056076, -21.8200912, +44.8614311, +71.7445984, -8.89900780, +7.54883194, +21.2461147
98+8.52080345, +36.9947319, +3.33841825, -8.89900780, +30.5617924, +17.5111866, -12.2559528
99+4.14040899, +14.5860672, +13.2448359, +7.54883194, +17.5111866, +27.0288830, +2.81372905
100-3.01028705, -7.01430130, -1.96349764, +21.2461147, -12.2559528, +2.81372905, +25.0000057
102T
104T
105getMatDet(rand)
106+0.721505485E+10
107
108call setCovRand(rngf, rand, scale)
109rand
110+39.7695770, -26.5509968, +2.51350141, +7.48580217, +0.290671550E-1, -24.8027153, -1.29713869
111-26.5509968, +73.1710815, -3.43008327, +8.98910904, -7.58647871, +13.5712662, -18.0960865
112+2.51350141, -3.43008327, +1.16550708, -0.200181246, -4.36577606, -0.174233735, -0.916450322
113+7.48580217, +8.98910904, -0.200181246, +5.61936760, +0.878764629, -7.40337849, -5.93482780
114+0.290671550E-1, -7.58647871, -4.36577606, +0.878764629, +37.6346054, -10.0206528, -19.8253746
115-24.8027153, +13.5712662, -0.174233735, -7.40337849, -10.0206528, +25.6715145, +4.66086578
116-1.29713869, -18.0960865, -0.916450322, -5.93482780, -19.8253746, +4.66086578, +99.9999924
118T
120T
121getMatDet(rand)
122+940080.875
123
124
125ndim = getUnifRand(3, 9)
126ndim
127+7
128call setResized(rand, [ndim, ndim])
129scale = getUnifRand(1, 10, ndim)
130scale
131+3.00000000, +5.00000000, +9.00000000, +9.00000000, +4.00000000, +1.00000000, +3.00000000
132
133call setCovRand(rngf, rand)
134rand
135+617.858093, -312.881348, +4.46586895, +190.990952, +274.088501, +2.21625376, +9.06907654
136-312.881348, +662.771851, +46.3033295, -2.50942993, +215.569397, -96.7078705, -13.8033905
137+4.46586895, +46.3033295, +20.9703140, +3.51977253, +19.4590225, -7.24011040, +0.804553688
138+190.990952, -2.50942993, +3.51977253, +91.8044205, +192.545731, -25.2380543, -1.69498229
139+274.088501, +215.569397, +19.4590225, +192.545731, +494.458954, -89.2690506, -8.78002357
140+2.21625376, -96.7078705, -7.24011040, -25.2380543, -89.2690506, +22.7236691, +3.14544582
141+9.06907654, -13.8033905, +0.804553688, -1.69498229, -8.78002357, +3.14544582, +1.00000012
143T
145T
146getMatDet(rand)
147+63426.3984
148call setCovRand(rngf, rand, scale(1))
149rand
150+13.9390383, -20.0662899, +9.04753399, -188.819778, +5.28256273, +8.12878323, -2.65992475
151-20.0662899, +1527.81421, -736.960205, +2720.85205, -33.1264610, +106.516724, +44.7382660
152+9.04753399, -736.960205, +386.258636, -1263.07959, +27.6776295, -72.1289368, -18.3278656
153-188.819778, +2720.85205, -1263.07959, +10366.5498, -46.0284576, -74.5005417, +63.5371208
154+5.28256273, -33.1264610, +27.6776295, -46.0284576, +9.82776070, -6.76881933, -2.03864050
155+8.12878323, +106.516724, -72.1289368, -74.5005417, -6.76881933, +34.1964340, +1.63812590
156-2.65992475, +44.7382660, -18.3278656, +63.5371208, -2.03864050, +1.63812590, +8.99999905
158T
160T
161getMatDet(rand)
162+0.115144755E+11
163
164call setCovRand(rngf, rand, scale)
165rand
166+227.928268, -607.185181, -3012.11621, -30.8781872, -65.4889221, -5.31984520, +2.29313517
167-607.185181, +4858.24170, +34753.7070, -175.558517, +343.589111, +37.0453720, +102.521461
168-3012.11621, +34753.7070, +280521.500, -1993.22046, +2047.24426, +319.004730, +750.816101
169-30.8781872, -175.558517, -1993.22046, +86.8605042, -12.6130610, -3.58376765, -1.93249154
170-65.4889221, +343.589111, +2047.24426, -12.6130610, +38.4054451, +4.03139162, +0.849420309
171-5.31984520, +37.0453720, +319.004730, -3.58376765, +4.03139162, +1.11486793, -0.966748714
172+2.29313517, +102.521461, +750.816101, -1.93249154, +0.849420309, -0.966748714, +8.99999809
174T
176T
177getMatDet(rand)
178+0.132806672E+13
179
180
181ndim = getUnifRand(3, 9)
182ndim
183+5
184call setResized(rand, [ndim, ndim])
185scale = getUnifRand(1, 10, ndim)
186scale
187+6.00000000, +6.00000000, +8.00000000, +6.00000000, +1.00000000
188
189call setCovRand(rngf, rand)
190rand
191+1.00000000, +0.236704782, -0.556507409, +0.524604499, +0.730599225
192+0.236704782, +1.00000000, +0.491989881, -0.670927286, +0.893746614E-1
193-0.556507409, +0.491989881, +0.999999881, -0.857951105, -0.462085903
194+0.524604499, -0.670927286, -0.857951105, +1.00000000, +0.310104609
195+0.730599225, +0.893746614E-1, -0.462085903, +0.310104609, +0.999999940
197T
199T
200getMatDet(rand)
201+0.505568867E-3
202call setCovRand(rngf, rand, scale(1))
203rand
204+45.9075546, +8.71592712, -22.0295620, +14.4144192, -13.0143461
205+8.71592712, +77.5395203, -36.2339630, +15.3402424, +30.7150993
206-22.0295620, -36.2339630, +36.7687035, +5.96755362, -10.1159983
207+14.4144192, +15.3402424, +5.96755362, +40.2152901, -6.58344221
208-13.0143461, +30.7150993, -10.1159983, -6.58344221, +36.0000038
210T
212T
213getMatDet(rand)
214+4724999.50
215
216call setCovRand(rngf, rand, scale)
217rand
218+38.8931961, +30.9693489, +44.3616180, -20.3788643, -3.21896076
219+30.9693489, +50.3003426, +23.7925568, +3.73630953, -2.25873351
220+44.3616180, +23.7925568, +96.0005798, -10.9213295, -8.27246094
221-20.3788643, +3.73630953, -10.9213295, +54.4054871, +0.557608604
222-3.21896076, -2.25873351, -8.27246094, +0.557608604, +1.00000000
224T
226T
227getMatDet(rand)
228+118958.555
229
230
231ndim = getUnifRand(3, 9)
232ndim
233+5
234call setResized(rand, [ndim, ndim])
235scale = getUnifRand(1, 10, ndim)
236scale
237+1.00000000, +5.00000000, +5.00000000, +9.00000000, +1.00000000
238
239call setCovRand(rngf, rand)
240rand
241+11.3617086, +6.54946613, +27.1131077, -1.26753402, -1.87211275
242+6.54946613, +586.147644, -43.1725922, +12.0227509, -14.7685471
243+27.1131077, -43.1725922, +188.709045, -11.2402802, -3.92123723
244-1.26753402, +12.0227509, -11.2402802, +1.31092739, -0.450334996
245-1.87211275, -14.7685471, -3.92123723, -0.450334996, +1.00000000
247T
249T
250getMatDet(rand)
251+7156.83740
252call setCovRand(rngf, rand, scale(1))
253rand
254+6.69470596, -10.1037893, -25.5022507, -0.843474448, +1.06889188
255-10.1037893, +2108.61084, -1698.81311, -40.7184143, +31.1774330
256-25.5022507, -1698.81311, +1571.93213, +36.8347931, -31.5178909
257-0.843474448, -40.7184143, +36.8347931, +1.03592753, -0.897042453
258+1.06889188, +31.1774330, -31.5178909, -0.897042453, +0.999999821
260T
262T
263getMatDet(rand)
264+117.119003
265
266call setCovRand(rngf, rand, scale)
267rand
268+3.27774715, -480.930725, +69.0719833, +7.50692463, +1.14822602
269-480.930725, +172131.859, +2378.35693, -3777.72461, -169.189835
270+69.0719833, +2378.35693, +6977.16113, +32.7830734, +32.5407333
271+7.50692463, -3777.72461, +32.7830734, +138.289062, +5.43580532
272+1.14822602, -169.189835, +32.5407333, +5.43580532, +1.00000000
274T
276T
277getMatDet(rand)
278+0.234586092E+11
279
280
281!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
282!Gram method for complex covariance.
283!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
284
285
286ndim = getUnifRand(3, 9)
287ndim
288+9
289call setResized(rand, [ndim, ndim])
290scale = getUnifRand(1, 10, ndim)
291scale
292+3.0000000000000000, +9.0000000000000000, +7.0000000000000000, +6.0000000000000000, +6.0000000000000000, +9.0000000000000000, +9.0000000000000000, +8.0000000000000000, +10.000000000000000
293
294call setCovRand(rngf, rand)
295rand
296+0.99999999999999989+0.0000000000000000i, +0.63217179887472008-0.33655488205228445i, -0.31886290753877666E-1+0.73468987116192752i, +0.28951683812391221+0.10414744687830052i, +0.38397383660360435+0.15989257442504839i, +0.51065747326686750-0.28997634513941206i, +0.22471899717558727+0.40960050406135695i, +0.14911362472802780+0.48945924559874388i, -0.40014393149821220E-1+0.22938147563270625i
297+0.63217179887472008+0.33655488205228445i, +1.0000000000000000+0.0000000000000000i, -0.42292861377502750+0.80704320842942678i, +0.49297932606648009+0.10799358266187534i, +0.36149434592047142-0.29808226195065679E-1i, +0.54874134650883633-0.13183855712112208i, -0.12425548507018996+0.16953776346992652i, +0.10936545944214343+0.35943312209129719i, +0.15315875861449668+0.24662973406056476i
298-0.31886290753877666E-1-0.73468987116192752i, -0.42292861377502750-0.80704320842942678i, +1.0000000000000000+0.0000000000000000i, +0.60858033402506839E-1-0.45736040699502212i, -0.31404708098929213-0.29026600664223090i, -0.40557056272212905-0.41474869164753686i, +0.10722898980115499-0.12262876026140533i, +0.18158469477473618-0.16246107187850400i, +0.14518611628562791-0.22469084197327802i
299+0.28951683812391221-0.10414744687830052i, +0.49297932606648009-0.10799358266187534i, +0.60858033402506839E-1+0.45736040699502212i, +1.0000000000000000+0.0000000000000000i, +0.39703584449336282-0.46895172512744623E-1i, +0.47659387633061712E-1-0.25722635375941050i, -0.24834159518425231-0.22852535317586947i, +0.27285449235406878+0.21726037472438064E-1i, +0.84250520470181639E-1+0.40623324807940303i
300+0.38397383660360435-0.15989257442504839i, +0.36149434592047142+0.29808226195065679E-1i, -0.31404708098929213+0.29026600664223090i, +0.39703584449336282+0.46895172512744623E-1i, +1.0000000000000000+0.0000000000000000i, +0.27674583625914728-0.68700668815060845E-1i, -0.25646642946846376E-1+0.89911635321156050E-1i, +0.47257151373924011+0.13315702704404092i, -0.11461667318785559+0.36589898683335209i
301+0.51065747326686750+0.28997634513941206i, +0.54874134650883633+0.13183855712112208i, -0.40557056272212905+0.41474869164753686i, +0.47659387633061712E-1+0.25722635375941050i, +0.27674583625914728+0.68700668815060845E-1i, +0.99999999999999989+0.0000000000000000i, -0.12334227331484443+0.62561515506547094E-1i, +0.14418866258292126+0.19088928587837942i, -0.91343747422469729E-1+0.31995844572819460i
302+0.22471899717558727-0.40960050406135695i, -0.12425548507018996-0.16953776346992652i, +0.10722898980115499+0.12262876026140533i, -0.24834159518425231+0.22852535317586947i, -0.25646642946846376E-1-0.89911635321156050E-1i, -0.12334227331484443-0.62561515506547094E-1i, +0.99999999999999978+0.0000000000000000i, +0.12683800720222232+0.50866724385910153E-1i, +0.13445439147865806+0.46502279685910863E-1i
303+0.14911362472802780-0.48945924559874388i, +0.10936545944214343-0.35943312209129719i, +0.18158469477473618+0.16246107187850400i, +0.27285449235406878-0.21726037472438064E-1i, +0.47257151373924011-0.13315702704404092i, +0.14418866258292126-0.19088928587837942i, +0.12683800720222232-0.50866724385910153E-1i, +0.99999999999999989+0.0000000000000000i, +0.13060253119857190+0.37623279273856725i
304-0.40014393149821220E-1-0.22938147563270625i, +0.15315875861449668-0.24662973406056476i, +0.14518611628562791+0.22469084197327802i, +0.84250520470181639E-1-0.40623324807940303i, -0.11461667318785559-0.36589898683335209i, -0.91343747422469729E-1-0.31995844572819460i, +0.13445439147865806-0.46502279685910863E-1i, +0.13060253119857190-0.37623279273856725i, +0.99999999999999989+0.0000000000000000i
306T
308T
309getMatDet(rand)
310+0.13277313303183156E-4-0.61833405149563925E-19i
311call setCovRand(rngf, rand, scale(1))
312rand
313+24.188444556243386+0.0000000000000000i, +10.729497493282803+4.8351992810785145i, -3.0680262530919524+3.1973337986868842i, -4.8157109776571740+6.1468466731640827i, -1.5846469587996326+3.0414784156301535i, -8.2913609317042010+0.66949107400570140i, +2.8826110579721274-5.4371720961617909i, -1.5820427626140303+2.1735404853463738i, +1.2326056931307467-1.5309451128209579i
314+10.729497493282803-4.8351992810785145i, +30.022406641402601+0.0000000000000000i, +1.5281626091344487-4.6455097816736206i, -2.9978168624431363+1.1659510336942716i, -2.5897780401906418+1.6512169812963331i, -3.0094356080880571-0.35705631154192630i, -2.8976058974161982-2.8420234970030833i, +2.1511984202530101+1.4859576014260103i, +1.5186217267564404+2.6666794857209002i
315-3.0680262530919524-3.1973337986868842i, +1.5281626091344487+4.6455097816736206i, +33.984829165845859+0.0000000000000000i, -4.5797220436831356+4.1011313455231555i, +3.6009894361075530+19.891611007748011i, -3.9065613196536244+0.67744469532300566i, -6.5257971711599856-4.8187531964625565i, +1.5082094936886277-2.2420039719935767i, -7.1611737561207178-6.3705951212391181i
316-4.8157109776571740-6.1468466731640827i, -2.9978168624431363-1.1659510336942716i, -4.5797220436831356-4.1011313455231555i, +14.983965387587626+0.0000000000000000i, +4.0404933549897510-9.3271119221148808i, +10.442943429803833-1.9114085702166399i, -6.1924565010868520+1.3043359088481616i, +4.1914143252602774-4.5678223370122328i, -2.3170547681399012+6.1887721679998542i
317-1.5846469587996326-3.0414784156301535i, -2.5897780401906418-1.6512169812963331i, +3.6009894361075530-19.891611007748011i, +4.0404933549897510+9.3271119221148808i, +30.114352209990802+0.0000000000000000i, +0.78615396853304298+11.054496091138613i, -3.2022607998766821+2.1234678573834032i, +4.7058006226404103+1.3565266607826358i, -9.0395438184283563+2.6248917739553281i
318-8.2913609317042010-0.66949107400570140i, -3.0094356080880571+0.35705631154192630i, -3.9065613196536244-0.67744469532300566i, +10.442943429803833+1.9114085702166399i, +0.78615396853304298-11.054496091138613i, +17.188294531210811+0.0000000000000000i, -5.0991381441067114+0.69023612784138177i, +6.1185026389273727-1.3190080295481599i, -2.3820526225320329+3.8463243344782363i
319+2.8826110579721274+5.4371720961617909i, -2.8976058974161982+2.8420234970030833i, -6.5257971711599856+4.8187531964625565i, -6.1924565010868520-1.3043359088481616i, -3.2022607998766821-2.1234678573834032i, -5.0991381441067114-0.69023612784138177i, +13.146477546202149+0.0000000000000000i, -0.79004801037656142E-1+4.6952969860195761i, +3.3557152822420395-1.2719590361880588i
320-1.5820427626140303-2.1735404853463738i, +2.1511984202530101-1.4859576014260103i, +1.5082094936886277+2.2420039719935767i, +4.1914143252602774+4.5678223370122328i, +4.7058006226404103-1.3565266607826358i, +6.1185026389273727+1.3190080295481599i, -0.79004801037656142E-1-4.6952969860195761i, +11.207652263677558+0.0000000000000000i, -2.4950542202274164+0.54609857384073179i
321+1.2326056931307467+1.5309451128209579i, +1.5186217267564404-2.6666794857209002i, -7.1611737561207178+6.3705951212391181i, -2.3170547681399012-6.1887721679998542i, -9.0395438184283563-2.6248917739553281i, -2.3820526225320329-3.8463243344782363i, +3.3557152822420395+1.2719590361880588i, -2.4950542202274164-0.54609857384073179i, +9.0000000000000000+0.0000000000000000i
323T
325T
326getMatDet(rand)
327+516031462.14512599+0.34272670745849609E-6i
328
329call setCovRand(rngf, rand, scale)
330rand
331+12.238469969341139+0.0000000000000000i, -17.122709299538496+12.433117500028706i, -5.1678870130208718-11.943161369012730i, -10.967238824193885-2.1581772015360765i, +1.4342817862335977+4.7710036434539633i, +4.2367627719394330+8.2944957428049033i, -19.753637041863154-6.4403854719300702i, -12.418826704187948+1.8740539564889522i, -5.8974144920238922+6.1947007117793635i
332-17.122709299538496-12.433117500028706i, +95.027263825388175+0.0000000000000000i, -35.940374010560731+46.343722926120037i, -9.8014060601206587+4.5052960992039477i, +5.2696899872174550-3.2217954465403880i, +15.009884630939451-28.978777091247764i, +38.480751130168812+18.496394454451039i, -2.1791264972714561+0.65394149534533863i, -19.833924418280343-18.640910791844082i
333-5.1678870130208718+11.943161369012730i, -35.940374010560731-46.343722926120037i, +946.48061236621140+0.0000000000000000i, +323.52511367485067+101.62235617396979i, -7.4733556015496827+38.623627731727716i, -13.754794029196834-30.784422864954401i, -354.76771125851980+106.31266250490719i, +49.643874834233600+31.779608761161040i, +12.528249874277602+23.538128676198124i
334-10.967238824193885+2.1581772015360765i, -9.8014060601206587-4.5052960992039477i, +323.52511367485067-101.62235617396979i, +278.80153207352566+0.0000000000000000i, -58.826247737649140+1.3745237829863437i, -11.572769916390680+63.236674375132061i, -100.21276103467625+109.71224824472785i, +39.883159029416163+21.638052784092871i, +63.460634453716089-39.603922518872523i
335+1.4342817862335977-4.7710036434539633i, +5.2696899872174550+3.2217954465403880i, -7.4733556015496827-38.623627731727716i, -58.826247737649140-1.3745237829863437i, +55.169796603573879+0.0000000000000000i, +7.8641031552955907-65.880946877034887i, -7.1074512323468708-15.242700531288479i, -7.9953892168440426+6.2003569155513372i, -12.900246861120564+17.810898180574554i
336+4.2367627719394330-8.2944957428049033i, +15.009884630939451+28.978777091247764i, -13.754794029196834+30.784422864954401i, -11.572769916390680-63.236674375132061i, +7.8641031552955907+65.880946877034887i, +135.77149748367887+0.0000000000000000i, -7.8625842420657488+4.8568864036183150i, -4.7903251449892146-25.749762090128407i, -40.829341329977552-25.511198823989751i
337-19.753637041863154+6.4403854719300702i, +38.480751130168812-18.496394454451039i, -354.76771125851980-106.31266250490719i, -100.21276103467625-109.71224824472785i, -7.1074512323468708+15.242700531288479i, -7.8625842420657488-4.8568864036183150i, +369.26089418097018+0.0000000000000000i, -25.731195236814060-35.798636596645352i, +21.216757956139261-36.070484486310718i
338-12.418826704187948-1.8740539564889522i, -2.1791264972714561-0.65394149534533863i, +49.643874834233600-31.779608761161040i, +39.883159029416163-21.638052784092871i, -7.9953892168440426-6.2003569155513372i, -4.7903251449892146+25.749762090128407i, -25.731195236814060+35.798636596645352i, +65.159459821234890+0.0000000000000000i, +14.938901595582838-10.295846708869442i
339-5.8974144920238922-6.1947007117793635i, -19.833924418280343+18.640910791844082i, +12.528249874277602-23.538128676198124i, +63.460634453716089+39.603922518872523i, -12.900246861120564-17.810898180574554i, -40.829341329977552+25.511198823989751i, +21.216757956139261+36.070484486310718i, +14.938901595582838+10.295846708869442i, +100.00000000000000+0.0000000000000000i
341T
343T
344getMatDet(rand)
345+453719258309644.62-4.4375000000000000i
346
347
348ndim = getUnifRand(3, 9)
349ndim
350+9
351call setResized(rand, [ndim, ndim])
352scale = getUnifRand(1, 10, ndim)
353scale
354+3.0000000000000000, +4.0000000000000000, +3.0000000000000000, +1.0000000000000000, +6.0000000000000000, +6.0000000000000000, +1.0000000000000000, +3.0000000000000000, +6.0000000000000000
355
356call setCovRand(rngf, rand)
357rand
358+1069.4305652954893+0.0000000000000000i, -1304.2893423475321+904.47243133511643i, +1921.9806013591067+5426.5167750171995i, +598.03211968956202+2669.1882162116863i, +15.808782278554816+177.08521436717854i, +224.77868721046630-594.38763626992431i, -91.967776574443690-576.34791902523477i, -156.23363222270297+41.190040518616129i, -10.214153798115436-0.50547490205611112E-1i
359-1304.2893423475321-904.47243133511643i, +7228.9613092745385+0.0000000000000000i, -14651.459873291289-7140.3009902515969i, -4739.7844923404045-8419.9793930120795i, +1439.9420178456810+1811.8563745686351i, +1045.4329265136541-132.27908008600957i, +252.58400243981617-1209.7654844759024i, -108.21049602338907-482.68154346084657i, +9.2742861654706470+20.200738124941598i
360+1921.9806013591067-5426.5167750171995i, -14651.459873291289+7140.3009902515969i, +259029.08271424440+0.0000000000000000i, +48471.155555161764+29476.127048200695i, +2998.1848883363109-8282.8207824112178i, +999.86704199230689+1389.9958210621826i, -2857.7364676516650+2091.1398454658033i, +239.47518624290856+425.34852547284152i, +44.557667935171160-37.127953530388631i
361+598.03211968956202-2669.1882162116863i, -4739.7844923404045+8419.9793930120795i, +48471.155555161764-29476.127048200695i, +37330.192123793306+0.0000000000000000i, -6904.9557870920353-3609.7565351441353i, -998.91123745515085+3808.1138688080873i, +1027.0806763775574+2291.5573894960226i, +1381.7886001299105+5.7714873110471032i, -25.071309633812369+33.928339285914731i
362+15.808782278554816-177.08521436717854i, +1439.9420178456810-1811.8563745686351i, +2998.1848883363109+8282.8207824112178i, -6904.9557870920353+3609.7565351441353i, +5272.0142000442502+0.0000000000000000i, -53.733121417640703-1140.4316343993014i, -915.41392696598871-500.11211170243803i, -362.72879051458165+163.93539010864365i, +24.115319150557511-6.2875781355768643i
363+224.77868721046630+594.38763626992431i, +1045.4329265136541+132.27908008600957i, +999.86704199230689-1389.9958210621826i, -998.91123745515085-3808.1138688080873i, -53.733121417640703+1140.4316343993014i, +3090.2634180711125+0.0000000000000000i, +1103.4618582469495-1522.1584628540236i, -353.22910208844837-810.58131355856290i, +15.644889971595592-4.9197291236728971i
364-91.967776574443690+576.34791902523477i, +252.58400243981617+1209.7654844759024i, -2857.7364676516650-2091.1398454658033i, +1027.0806763775574-2291.5573894960226i, -915.41392696598871+500.11211170243803i, +1103.4618582469495+1522.1584628540236i, +3695.8674597402996+0.0000000000000000i, +671.04923610718322-326.43325684506112i, -11.037991279931191+20.419530398029107i
365-156.23363222270297-41.190040518616129i, -108.21049602338907+482.68154346084657i, +239.47518624290856-425.34852547284152i, +1381.7886001299105-5.7714873110471032i, -362.72879051458165-163.93539010864365i, -353.22910208844837+810.58131355856290i, +671.04923610718322+326.43325684506112i, +330.17524033504515+0.0000000000000000i, -2.2420246313264114+8.2975418908786729i
366-10.214153798115436+0.50547490205611112E-1i, +9.2742861654706470-20.200738124941598i, +44.557667935171160+37.127953530388631i, -25.071309633812369-33.928339285914731i, +24.115319150557511+6.2875781355768643i, +15.644889971595592+4.9197291236728971i, -11.037991279931191-20.419530398029107i, -2.2420246313264114-8.2975418908786729i, +1.0000000000000002+0.0000000000000000i
368T
370T
371getMatDet(rand)
372+0.25799527933889667E+24-650133897216.00000i
373call setCovRand(rngf, rand, scale(1))
374rand
375+352.13537145559866+0.0000000000000000i, -454.05699832015750+57.144441712489680i, +104355.54735688977+29893.375187356087i, +474.82308743044416-312.23967238423302i, -503.21386059641600+2256.0930910019670i, +61.758687972671808+53.388482214176022i, -131.81040890000929-680.98103660165816i, -19.334248771053471+17.616026903362929i, +22.194092455787686+8.1625068322459686i
376-454.05699832015750-57.144441712489680i, +797.27148706693993+0.0000000000000000i, -162016.22294854699-49546.235164533602i, -764.19120410950165+447.71200190342654i, +873.24994060147196-3517.5484921287261i, -86.696801103274794-89.945115739765200i, +225.02064538840207+1019.1953296296520i, +8.2237208984393853-11.218571720660130i, -21.383473785706425-15.919836202070332i
377+104355.54735688977-29893.375187356087i, -162016.22294854699+49546.235164533602i, +41931750.319800243+0.0000000000000000i, +141151.67575177565-152410.83961411682i, +41516.110069339498+886166.94427728700i, +26585.407629455087+13991.799237143965i, -136247.97575548169-236243.32156854140i, -455.75883116900297+7727.0962665817697i, +8450.7313942786259-196.98282883504919i
378+474.82308743044416+312.23967238423302i, -764.19120410950165-447.71200190342654i, +141151.67575177565+152410.83961411682i, +1088.1273156826442+0.0000000000000000i, -3080.8346573194335+3184.4643833505270i, +36.139885851266996+152.39699521384094i, +392.10026183285800-1247.5594748682977i, -28.119676710149452+19.357903994570826i, +23.505500023641236+22.297048756962518i
379-503.21386059641600-2256.0930910019670i, +873.24994060147196+3517.5484921287261i, +41516.110069339498-886166.94427728700i, -3080.8346573194335-3184.4643833505270i, +18861.339730768312+0.0000000000000000i, +319.61298965853075-549.07310706603187i, -5100.7877712313702+2663.4867589833598i, +159.98931044324124+18.852098640821794i, +3.6033791245610338-173.56938276316674i
380+61.758687972671808-53.388482214176022i, -86.696801103274794+89.945115739765200i, +26585.407629455087-13991.799237143965i, +36.139885851266996-152.39699521384094i, +319.61298965853075+549.07310706603187i, +28.386722574337327+0.0000000000000000i, -157.39135835836100-107.57184992303607i, +0.78312030772047070+3.5224994129637470i, +2.7511621821195327+0.81218095358348197i
381-131.81040890000929+680.98103660165816i, +225.02064538840207-1019.1953296296520i, -136247.97575548169+236243.32156854140i, +392.10026183285800+1247.5594748682977i, -5100.7877712313702-2663.4867589833598i, -157.39135835836100+107.57184992303607i, +1828.6755229741248+0.0000000000000000i, -49.061052032403147-34.084037001562045i, -26.803564492970317+53.845564193654567i
382-19.334248771053471-17.616026903362929i, +8.2237208984393853+11.218571720660130i, -455.75883116900297-7727.0962665817697i, -28.119676710149452-19.357903994570826i, +159.98931044324124-18.852098640821794i, +0.78312030772047070-3.5224994129637470i, -49.061052032403147+34.084037001562045i, +9.0000640097880815+0.0000000000000000i, -1.6137193692771696-3.1811494964323117i
383+22.194092455787686-8.1625068322459686i, -21.383473785706425+15.919836202070332i, +8450.7313942786259+196.98282883504919i, +23.505500023641236-22.297048756962518i, +3.6033791245610338+173.56938276316674i, +2.7511621821195327-0.81218095358348197i, -26.803564492970317-53.845564193654567i, -1.6137193692771696+3.1811494964323117i, +8.9999999999999964+0.0000000000000000i
385T
387T
388getMatDet(rand)
389+1143227842179555.0+55.187500000000000i
390
391call setCovRand(rngf, rand, scale)
392rand
393+155100030.77876353+0.0000000000000000i, +7266240.2750910055+6737848.3721818365i, -25586702.958410259-25432336.587381463i, +192572.06328152466-5356364.1018221797i, +529998.88172831538+1641823.0488274698i, +2047071.0510740583-1985334.4042129130i, -369549.02992404590+1105052.8900277587i, +30619.972224723802-904146.87123442790i, -3026.3617838228010+6557.0295855333834i
394+7266240.2750910055-6737848.3721818365i, +633135.76230546902+0.0000000000000000i, -2303545.6748278551-79933.289977486362i, -223670.70255219718-259304.78068852684i, +96152.211043776595+53887.097197861476i, +9664.5160516162141-181933.26816567389i, +30693.281066086576+67826.717952106716i, -37846.423650215838-43689.645839472534i, +144.09907740852191+439.35746039042294i
395-25586702.958410259+25432336.587381463i, -2303545.6748278551+79933.289977486362i, +8392120.7056018692+0.0000000000000000i, +846583.75740124227+915267.42735936888i, -356773.51640780683-184059.51780810877i, -12330.459690614760+663340.96224512020i, -120252.94951722957-242908.59840369486i, +143226.71006185163+154165.29738442414i, -647.87836672144761-1554.3422616321736i
396+192572.06328152466+5356364.1018221797i, -223670.70255219718+259304.78068852684i, +846583.75740124227-915267.42735936888i, +185234.23661834604+0.0000000000000000i, -56033.999868330502+20328.802716384820i, +71087.816024238171+68242.234749524781i, -38624.715467285518-11390.609480434174i, +31267.682195298632-68.374314668388052i, -235.38572548026477-90.225859137721628i
397+529998.88172831538-1641823.0488274698i, +96152.211043776595-53887.097197861476i, -356773.51640780683+184059.51780810877i, -56033.999868330502-20328.802716384820i, +19377.826212137428+0.0000000000000000i, -14077.028345410534-28583.026425028642i, +10432.854276063659+7682.0799844302637i, -9456.7972032149464-3407.7440481170470i, +61.567503766448652+20.124738250852221i
398+2047071.0510740583+1985334.4042129130i, +9664.5160516162141+181933.26816567389i, -12330.459690614760-663340.96224512020i, +71087.816024238171-68242.234749524781i, -14077.028345410534+28583.026425028642i, +52612.544097875274+0.0000000000000000i, -19016.145477565060+9858.0809101611612i, +11965.609847278023-11526.309123835414i, -106.60995664553425+60.524205419122765i
399-369549.02992404590-1105052.8900277587i, +30693.281066086576-67826.717952106716i, -120252.94951722957+242908.59840369486i, -38624.715467285518+11390.609480434174i, +10432.854276063659-7682.0799844302637i, -19016.145477565060-9858.0809101611612i, +8755.6278243789511+0.0000000000000000i, -6515.7638713618908+1937.6740935903613i, +55.592047703524031+3.7465548175883789i
400+30619.972224723802+904146.87123442790i, -37846.423650215838+43689.645839472534i, +143226.71006185163-154165.29738442414i, +31267.682195298632+68.374314668388052i, -9456.7972032149464+3407.7440481170470i, +11965.609847278023+11526.309123835414i, -6515.7638713618908-1937.6740935903613i, +5287.3227903407542+0.0000000000000000i, -41.779855200171959-6.5931965490545492i
401-3026.3617838228010-6557.0295855333834i, +144.09907740852191-439.35746039042294i, -647.87836672144761+1554.3422616321736i, -235.38572548026477+90.225859137721628i, +61.567503766448652-20.124738250852221i, -106.60995664553425-60.524205419122765i, +55.592047703524031-3.7465548175883789i, -41.779855200171959+6.5931965490545492i, +36.000000000000014+0.0000000000000000i
403T
405T
406getMatDet(rand)
407+71235295923925992.+556652.00000000000i
408
409
410ndim = getUnifRand(3, 9)
411ndim
412+9
413call setResized(rand, [ndim, ndim])
414scale = getUnifRand(1, 10, ndim)
415scale
416+3.0000000000000000, +9.0000000000000000, +3.0000000000000000, +5.0000000000000000, +2.0000000000000000, +3.0000000000000000, +5.0000000000000000, +4.0000000000000000, +9.0000000000000000
417
418call setCovRand(rngf, rand)
419rand
420+1343495546965861.8+0.0000000000000000i, +62941258211740.117+58364405795764.641i, -6931972462720.0430-6889360917707.6533i, +16446534628.162292-458763734515.16888i, +35359206316.112068+109471550634.52771i, +28729744551.082150-27837064273.828739i, -1946899294.5011809+5823710276.6244688i, -17258006.765676238-2276762.9773725546i, -15314374.644660067-5217201.0336584318i
421+62941258211740.117-58364405795764.641i, +5484205635109.8047+0.0000000000000000i, -624038733582.89014-21621988311.131931i, -19158255440.695168-22206899176.625473i, +6411959429.1008234+3592718029.9040461i, +136648819.05585250-2552310164.3987365i, +161800793.64528432+357425043.59001803i, -907472.73331286036+643041.74764465610i, -944110.16338271392+420870.19704606902i
422-6931972462720.0430+6889360917707.6533i, -624038733582.89014+21621988311.131931i, +2273506066580.9067+0.0000000000000000i, +72521031564.188400+78406346532.948364i, -23787223893.432125-12271082759.071785i, -171254068.88009977+9300491364.5677319i, -634544434.92613733-1281557035.0038218i, +154484.99498740624-3595.8280210322264i, +805787.34428771422-94620.676250700446i
423+16446534628.162292+458763734515.16888i, -19158255440.695168+22206899176.625473i, +72521031564.188400-78406346532.948364i, +15862900354.153978+0.0000000000000000i, -3737197626.8435054+1356350311.2241161i, +997196634.64228415+956995181.12911737i, -203890752.21665555-60161439.865317933i, -2935.1228004332934+2356.7735582905807i, +912.22323729643017-46712.307196102469i
424+35359206316.112068-109471550634.52771i, +6411959429.1008234-3592718029.9040461i, -23787223893.432125+12271082759.071785i, -3737197626.8435054-1356350311.2241161i, +1284058857.2869618+0.0000000000000000i, -196561753.36367717-398720627.30818295i, +55025369.478341900+40529824.480642825i, +7150.1520692091090+5285.6095413207931i, -5164.1939488204125+16640.869881813029i
425+28729744551.082150+27837064273.828739i, +136648819.05585250+2552310164.3987365i, -171254068.88009977-9300491364.5677319i, +997196634.64228415-956995181.12911737i, -196561753.36367717+398720627.30818295i, +734842198.95298791+0.0000000000000000i, -100299365.07042469+51926507.260917030i, -4896.0195919282341+10527.510725671065i, -11814.965868272926-8166.8003550483709i
426-1946899294.5011809-5823710276.6244688i, +161800793.64528432-357425043.59001803i, -634544434.92613733+1281557035.0038218i, -203890752.21665555+60161439.865317933i, +55025369.478341900-40529824.480642825i, -100299365.07042469-51926507.260917030i, +46212865.232756995+0.0000000000000000i, +69.688650099512415-1809.1685839202623i, +609.63835637448130+2214.7650109107603i
427-17258006.765676238+2276762.9773725546i, -907472.73331286036-643041.74764465610i, +154484.99498740624+3595.8280210322264i, -2935.1228004332934-2356.7735582905807i, +7150.1520692091090-5285.6095413207931i, -4896.0195919282341-10527.510725671065i, +69.688650099512415+1809.1685839202623i, +1790.0265412818007+0.0000000000000000i, -4.2203466171684303-11.658045811746822i
428-15314374.644660067+5217201.0336584318i, -944110.16338271392-420870.19704606902i, +805787.34428771422+94620.676250700446i, +912.22323729643017+46712.307196102469i, -5164.1939488204125-16640.869881813029i, -11814.965868272926+8166.8003550483709i, +609.63835637448130-2214.7650109107603i, -4.2203466171684303+11.658045811746822i, +1.0000000000000002+0.0000000000000000i
430T
432T
433getMatDet(rand)
434+0.14594538903704364E+65-0.22886761546550783E+62i
435call setCovRand(rngf, rand, scale(1))
436rand
437+62319169300259712.+0.0000000000000000i, +5291747653765.9395+4238548030292.4390i, +929860377763.26416+257712204909.74655i, -24213238325.597000+9421923491.6046467i, +104038734336.63643+257737607635.79474i, -1083250038303447.9-365741004135863.94i, +1055093869868.1617-434289380330.82611i, -33963807.378569633-307948748.97563136i, -280314649.69378889+165275403.61698794i
438+5291747653765.9395-4238548030292.4390i, +738529468.47354007+0.0000000000000000i, +96391369.219072416-41193076.412463069i, -1414286.5989651291+2447083.0994010386i, +26358781.288845666+14814168.130165059i, -116850662392.13374+42588644248.787796i, +60069029.909681939-108610295.93742841i, -22743.694823822269-23106.097634172489i, -13089.185108289408+32537.185351567940i
439+929860377763.26416-257712204909.74655i, +96391369.219072416+41193076.412463069i, +14980922.099127794+0.0000000000000000i, -322397.91054642003+240524.87273952196i, +2619625.7972072540+3415911.4312568670i, -17682036844.290199-975714125.06156015i, +13950609.928086052-10848818.187349189i, -1747.6907294862933-4724.2282275594534i, -3567.4727586757608+3798.0763975727814i
440-24213238325.597000-9421923491.6046467i, -1414286.5989651291-2447083.0994010386i, -322397.91054642003-240524.87273952196i, +10839.762557678147+0.0000000000000000i, -1460.4643884795978-115865.83648643459i, +365589019.60281754+305844816.79328948i, -475582.23297354649+9247.4159575370813i, -33.758465535164589+125.90978280719079i, +133.99714213857220-20.602711105885664i
441+104038734336.63643-257737607635.79474i, +26358781.288845666-14814168.130165059i, +2619625.7972072540-3415911.4312568670i, -1460.4643884795978+115865.83648643459i, +1239689.7911372704+0.0000000000000000i, -3321257284.0576429+3869621899.6208539i, -34632.932559507368-5088887.0711300094i, -1333.4122204860798-384.57560438523672i, +215.28408083875254+1443.5335835124104i
442-1083250038303447.9+365741004135863.94i, -116850662392.13374-42588644248.787796i, -17682036844.290199+975714125.06156015i, +365589019.60281754-305844816.79328948i, -3321257284.0576429-3869621899.6208539i, +20976963499394.426+0.0000000000000000i, -15791995753.789490+13741847737.332785i, +2379889.3084241510+5198336.3694421342i, +3923113.0468479469-4540228.2140845638i
443+1055093869868.1617+434289380330.82611i, +60069029.909681939+108610295.93742841i, +13950609.928086052+10848818.187349189i, -475582.23297354649-9247.4159575370813i, -34632.932559507368+5088887.0711300094i, -15791995753.789490-13741847737.332785i, +20890820.518268209+0.0000000000000000i, +1614.1255525549884-5464.4880697867939i, -5931.0728000616682+848.61638959850393i
444-33963807.378569633+307948748.97563136i, -22743.694823822269+23106.097634172489i, -1747.6907294862933+4724.2282275594534i, -33.758465535164589-125.90978280719079i, -1333.4122204860798+384.57560438523672i, +2379889.3084241510-5198336.3694421342i, +1614.1255525549884+5464.4880697867939i, +9.0367383717977745+0.0000000000000000i, -3.9175602106874097-1.2981218058136337i
445-280314649.69378889-165275403.61698794i, -13089.185108289408-32537.185351567940i, -3567.4727586757608-3798.0763975727814i, +133.99714213857220+20.602711105885664i, +215.28408083875254-1443.5335835124104i, +3923113.0468479469+4540228.2140845638i, -5931.0728000616682-848.61638959850393i, -3.9175602106874097+1.2981218058136337i, +9.0000000000000000+0.0000000000000000i
447T
449T
450getMatDet(rand)
451+0.20348577898473199E+34+0.10300898960243798E+26i
452
453call setCovRand(rngf, rand, scale)
454rand
455+0.33019291825899974E+18+0.0000000000000000i, +0.67970196302899994E+18+0.26460804969820512E+18i, -5270788900759.4199-14252296738249.086i, -68563961991896.719-36412473633396.547i, -216819600695105.94-67928858672056.625i, +106042976361812.81-162520992438087.06i, -1031103280.5721830+149769486.44465709i, +834343811.26794457+299443921.28693718i, -966143094.80625308-2246628367.2408776i
456+0.67970196302899994E+18-0.26460804969820512E+18i, +0.16112162718696796E+19+0.0000000000000000i, -22271337628238.383-25114473501479.172i, -170318833321591.75-20009552231795.938i, -500759791379022.88+33921998296789.586i, +88049823787643.484-419530238933174.88i, -2002498421.3045659+1134597679.4090116i, +1957713572.3387547-52008695.128760509i, -3789677208.3991385-3850236963.7607374i
457-5270788900759.4199+14252296738249.086i, -22271337628238.383+25114473501479.172i, +699315997.87038493+0.0000000000000000i, +2666159270.8870912-2378218982.2697535i, +6393089491.9060278-8274376349.2001705i, +5322257724.0163002+7171480791.3806496i, +9998.2149866345826-46894.351848533472i, -26250.700363024636+31233.427691073186i, +112388.74998831221-5862.7745003185728i
458-68563961991896.719+36412473633396.547i, -170318833321591.75+20009552231795.938i, +2666159270.8870912+2378218982.2697535i, +18252881726.372681+0.0000000000000000i, +52512721723.626221-9804304571.0223255i, -4098255762.9880791+45439695332.715134i, +197598.27493009280-144810.27760669758i, -205705.09624892124+29586.307728589712i, +448320.26482212852+360970.69658912154i
459-216819600695105.94+67928858672056.625i, -500759791379022.88-33921998296789.586i, +6393089491.9060278+8274376349.2001705i, +52512721723.626221+9804304571.0223255i, +156349714884.34955+0.0000000000000000i, -36199370146.299194+128538260571.24850i, +646181.83801426529-310489.19043506583i, -610723.75463700842-25608.162878323412i, +1098662.8584059924+1274937.0166720184i
460+106042976361812.81+162520992438087.06i, +88049823787643.484+419530238933174.88i, +5322257724.0163002-7171480791.3806496i, -4098255762.9880791-45439695332.715134i, -36199370146.299194-128538260571.24850i, +114060628202.32614+0.0000000000000000i, -404890.50809690414-459247.68186258565i, +120127.99028365943+510727.69058938592i, +789786.12557069119-1200920.0972746499i
461-1031103280.5721830-149769486.44465709i, -2002498421.3045659-1134597679.4090116i, +9998.2149866345826+46894.351848533472i, +197598.27493009280+144810.27760669758i, +646181.83801426529+310489.19043506583i, -404890.50809690414+459247.68186258565i, +25.003632912258922+0.0000000000000000i, +0.52160968853186873+4.4097718201342566i, +3.5326334624192879+2.6500022088032873i
462+834343811.26794457-299443921.28693718i, +1957713572.3387547+52008695.128760509i, -26250.700363024636-31233.427691073186i, -205705.09624892124-29586.307728589712i, -610723.75463700842+25608.162878323412i, +120127.99028365943-510727.69058938592i, +0.52160968853186873-4.4097718201342566i, +16.000006586061719+0.0000000000000000i, -0.73448648795992277+4.8927346387868687i
463-966143094.80625308+2246628367.2408776i, -3789677208.3991385+3850236963.7607374i, +112388.74998831221+5862.7745003185728i, +448320.26482212852-360970.69658912154i, +1098662.8584059924-1274937.0166720184i, +789786.12557069119+1200920.0972746499i, +3.5326334624192879-2.6500022088032873i, -0.73448648795992277-4.8927346387868687i, +81.000000000000000+0.0000000000000000i
465T
467T
468getMatDet(rand)
469+0.65575783712460318E+38+0.96099466828368688E+34i
470
471
472ndim = getUnifRand(3, 9)
473ndim
474+5
475call setResized(rand, [ndim, ndim])
476scale = getUnifRand(1, 10, ndim)
477scale
478+1.0000000000000000, +5.0000000000000000, +2.0000000000000000, +3.0000000000000000, +1.0000000000000000
479
480call setCovRand(rngf, rand)
481rand
482+874791074915.72778+0.0000000000000000i, +1654863531.8330736+3456174867.6806126i, +100205245782.11807-19153093036.027203i, +539237016443.10834-120987236876.35603i, +446047.05885539198+215054.88667077123i
483+1654863531.8330736-3456174867.6806126i, +18020017.609215338+0.0000000000000000i, +121174119.69610503-480027086.35118186i, +581999420.15579212-2532978925.9001756i, +1409.2497263603198-1373.0816775984924i
484+100205245782.11807+19153093036.027203i, +121174119.69610503+480027086.35118186i, +33897971589062.527+0.0000000000000000i, +69106266407.221359-2392483333.9750562i, -1734274.9275035011-192012.64822718716i
485+539237016443.10834+120987236876.35603i, +581999420.15579212+2532978925.9001756i, +69106266407.221359+2392483333.9750562i, +374844756844.12793+0.0000000000000000i, +238681.10418605231+153717.24889554008i
486+446047.05885539198-215054.88667077123i, +1409.2497263603198+1373.0816775984924i, -1734274.9275035011+192012.64822718716i, +238681.10418605231-153717.24889554008i, +1.0000000000000000+0.0000000000000000i
488T
490T
491getMatDet(rand)
492+0.21144876388100787E+35-0.58411921786521521E+26i
493call setCovRand(rngf, rand, scale(1))
494rand
495+3.5397684569324062+0.0000000000000000i, +0.89924732755146686+1.4556013750471977i, -3448449.5197038865-536710.62355260528i, +0.15549008314429316+0.17229448832761005i, -0.74544948763553620+0.40927538701499976E-1i
496+0.89924732755146686-1.4556013750471977i, +1.0032702081018967+0.0000000000000000i, +998939.97049478150-413430.74702006503i, +0.15580355521096620+0.30617546720045868i, -0.26946105186388952+0.92018674075671852E-1i
497-3448449.5197038865+536710.62355260528i, +998939.97049478150+413430.74702006503i, +357410501630686.62+0.0000000000000000i, -6540798.3869468458+9261193.9913189206i, +5192212.4592452655-9535294.0528666209i
498+0.15549008314429316-0.17229448832761005i, +0.15580355521096620-0.30617546720045868i, -6540798.3869468458-9261193.9913189206i, +1.0000000000002194+0.0000000000000000i, -0.88142117977167334+0.20531418873610907i
499-0.74544948763553620-0.40927538701499976E-1i, -0.26946105186388952-0.92018674075671852E-1i, +5192212.4592452655+9535294.0528666209i, -0.88142117977167334-0.20531418873610907i, +0.99999999999999978+0.0000000000000000i
501T
503T
504getMatDet(rand)
505+1967920239645.0513+0.94909667968750000E-2i
506
507call setCovRand(rngf, rand, scale)
508rand
509+4160611639961.8486+0.0000000000000000i, +32581648738717.344-78824782939362.000i, +126193.37428494007+1973588.9401069998i, +916309.74610039848-2463238.1665149205i, +1107544.2427382448-19879.273728671957i
510+32581648738717.344+78824782939362.000i, +1748519417219247.5+0.0000000000000000i, -36402388.228523031+17845954.720489126i, +53842836.715489067-1929524.3886041278i, +9049805.9193919953+20827273.952775516i
511+126193.37428494007-1973588.9401069998i, -36402388.228523031-17845954.720489126i, +4.0000000000000044+0.0000000000000000i, -0.24782584983110861+1.4558889862780409i, -0.75490045887146440-0.98516431398038240i
512+916309.74610039848+2463238.1665149205i, +53842836.715489067+1929524.3886041278i, -0.24782584983110861-1.4558889862780409i, +15.870086499414956+0.0000000000000000i, +0.79014269897867673+0.40484907578558721i
513+1107544.2427382448+19879.273728671957i, +9049805.9193919953-20827273.952775516i, -0.75490045887146440+0.98516431398038240i, +0.79014269897867673-0.40484907578558721i, +1.0000000000000000+0.0000000000000000i
515T
517T
518getMatDet(rand)
519+13367343050791844.-7310620.0000000000i
520
521
522ndim = getUnifRand(3, 9)
523ndim
524+6
525call setResized(rand, [ndim, ndim])
526scale = getUnifRand(1, 10, ndim)
527scale
528+7.0000000000000000, +8.0000000000000000, +8.0000000000000000, +5.0000000000000000, +10.000000000000000, +2.0000000000000000
529
530call setCovRand(rngf, rand)
531rand
532+0.22141657866709462E+30+0.0000000000000000i, +0.10742048844852886E+21-0.11607529079002026E+20i, +0.23164442052151706E+21+0.26051205725235059E+21i, +0.23400593315388293E+24+0.98409299523304603E+24i, +1366588227906098.2+6687994319106210.0i, +118360593800977.78-240086954570624.62i
533+0.10742048844852886E+21+0.11607529079002026E+20i, +0.30842305736230348E+23+0.0000000000000000i, +1140514135867766.5-5348560453786729.0i, +61499457235268.969+489190309802430.62i, -21021486270.468246+10568671546.176046i, +35943662483.353699-102430489119.38293i
534+0.23164442052151706E+21-0.26051205725235059E+21i, +1140514135867766.5+5348560453786729.0i, +576268644451.60718+0.0000000000000000i, +1402689432545201.2+754232989207227.88i, +10266178.571503295+4982015.5660802592i, -112319.24461889527-376235.83404748584i
535+0.23400593315388293E+24-0.98409299523304603E+24i, +61499457235268.969-489190309802430.62i, +1402689432545201.2-754232989207227.88i, +0.46211435930851820E+19+0.0000000000000000i, +31169965008.791569+993825859.23412669i, -941958668.73280084-779784797.35779166i
536+1366588227906098.2-6687994319106210.0i, -21021486270.468246-10568671546.176046i, +10266178.571503295-4982015.5660802592i, +31169965008.791569-993825859.23412669i, +257.00021075401838+0.0000000000000000i, -6.0671984106243650-3.8650922676576807i
537+118360593800977.78+240086954570624.62i, +35943662483.353699+102430489119.38293i, -112319.24461889527+376235.83404748584i, -941958668.73280084+779784797.35779166i, -6.0671984106243650+3.8650922676576807i, +1.0000000000000002+0.0000000000000000i
539T
541T
542getMatDet(rand)
543+0.36530188947842479E+70+0.17526200881160303E+65i
544call setCovRand(rngf, rand, scale(1))
545rand
546+0.33526965031969882E+32+0.0000000000000000i, -0.39072310038743802E+20-0.11134656723865398E+21i, -0.40654808846107268E+30+0.37291669761398095E+30i, +6544409104033978.0+14813818982119444.i, -7809777315259260.0+2918205428920772.0i, -24023883390915136.+6412312850641423.0i
547-0.39072310038743802E+20+0.11134656723865398E+21i, +415328529.65191370+0.0000000000000000i, -0.76470450465255398E+18-0.17847857937133164E+19i, -56810.106229374178+4482.9342033948406i, -611.08661303962879-29341.674907332661i, +6706.6409785075557-87263.045224404763i
548-0.40654808846107268E+30-0.37291669761398095E+30i, -0.76470450465255398E+18+0.17847857937133164E+19i, +0.90777143502388622E+28+0.0000000000000000i, +85414932210627.734-252425151184458.09i, +127160259238908.88+51481114088113.305i, +362637130463179.38+189459252237015.06i
549+6544409104033978.0-14813818982119444.i, -56810.106229374178-4482.9342033948406i, +85414932210627.734+252425151184458.09i, +49.000000001022499+0.0000000000000000i, -28.628400864323567-0.99313159031529330i, -5.7030155561205911+28.189458850181897i
550-7809777315259260.0-2918205428920772.0i, -611.08661303962879+29341.674907332661i, +127160259238908.88-51481114088113.305i, -28.628400864323567+0.99313159031529330i, +55.553196984231498+0.0000000000000000i, +4.6035830367324895-19.508741371244433i
551-24023883390915136.-6412312850641423.0i, +6706.6409785075557+87263.045224404763i, +362637130463179.38-189459252237015.06i, -5.7030155561205911-28.189458850181897i, +4.6035830367324895+19.508741371244433i, +49.000000000000000+0.0000000000000000i
553T
555T
556getMatDet(rand)
557+0.98487591527881179E+45+0.15604037845907035E+50i
558
559call setCovRand(rngf, rand, scale)
560rand
561+0.79261801806871948E+33+0.0000000000000000i, +0.28541225886348336E+33+0.41655692836993947E+31i, +0.75396562123884413E+32-0.20943945402712286E+33i, +12498650602882804.+69097419466622000.i, -0.41785264033382054E+22+0.35473359063505431E+22i, -15658129675411800.+6151994233370259.0i
562+0.28541225886348336E+33-0.41655692836993947E+31i, +0.10279542934885601E+33+0.0000000000000000i, +0.26048699473342772E+32-0.75812883258333625E+32i, +4863752906876198.0+24815467385529364.i, -0.14859944629173265E+22+0.12993132531307630E+22i, -5605973496387436.0+2297550084862569.5i
563+0.75396562123884413E+32+0.20943945402712286E+33i, +0.26048699473342772E+32+0.75812883258333625E+32i, +0.62513752341905719E+32+0.0000000000000000i, -17069218984921004.+9875398059225988.0i, -0.13348153265161324E+22-0.76668871767484151E+21i, -3115043324749349.5-3552267113049426.5i
564+12498650602882804.-69097419466622000.i, +4863752906876198.0-24815467385529364.i, -17069218984921004.-9875398059225988.0i, +25.006768964601676+0.0000000000000000i, +474712.80528594519-850003.17377562390i, -0.34860366770042917+3.5681986790034874i
565-0.41785264033382054E+22-0.35473359063505431E+22i, -0.14859944629173265E+22-0.12993132531307630E+22i, -0.13348153265161324E+22+0.76668871767484151E+21i, +474712.80528594519+850003.17377562390i, +140029122415312.53+0.0000000000000000i, +1892728.1196094402+11509953.921628647i
566-15658129675411800.-6151994233370259.0i, -5605973496387436.0-2297550084862569.5i, -3115043324749349.5+3552267113049426.5i, -0.34860366770042917-3.5681986790034874i, +1892728.1196094402-11509953.921628647i, +4.0000000000000009+0.0000000000000000i
568F
570T
571getMatDet(rand)
572-0.56111642697928132E+81-0.99853355198952856E+80i
573
574
575!%%%%%%%%%%%%%%%%%%%%%%%%
576!Dvine and Onion methods.
577!%%%%%%%%%%%%%%%%%%%%%%%%
578
579
580eta = getUnifRand(1, 10)
581ndim = getUnifRand(2, 5)
582call setResized(rand, [ndim, ndim])
583scale = getUnifRand(1, 10)
584
585call setCovRand(rngf, rand, dvine, eta)
586onion%info
587+0
588rand
589+1.00000000, -0.182436705E-1, -0.344355106
590-0.182436705E-1, +1.00000000, +0.374031126
591-0.344355106, +0.374031126, +1.00000000
593T
595T
596getMatDet(rand)
597+0.745887041
598call setCovRand(rngf, rand, onion, eta)
599onion%info
600+0
601rand
602+1.00000000, +0.161169767E-1, +0.482837111E-1
603+0.161169767E-1, +1.00000000, -0.135480821
604+0.482837111E-1, -0.135480821, +1.00000000
606T
608T
609getMatDet(rand)
610+0.978843033
611call setCovRand(rngf, rand, dvine, eta, scale)
612rand
613+9.00000000, +3.29877591, -1.48765397
614+3.29877591, +9.00000000, +3.56839824
615-1.48765397, +3.56839824, +9.00000000
617T
619T
620getMatDet(rand)
621+461.520020
622
623call setCovRand(rngf, rand, onion, eta, scale)
624onion%info
625+0
626rand
627+9.00000000, +0.283311009, +2.90053558
628+0.283311009, +9.00000000, +11.0620708
629+2.90053558, +11.0620708, +9.00000000
631F
633T
634getMatDet(rand)
635-430.584534
636
637
638eta = getUnifRand(1, 10)
639ndim = getUnifRand(2, 5)
640call setResized(rand, [ndim, ndim])
641scale = getUnifRand(1, 10)
642
643call setCovRand(rngf, rand, dvine, eta)
644onion%info
645+0
646rand
647+1.00000000, -0.232560039E-1, +0.390677571
648-0.232560039E-1, +1.00000000, +0.120144054
649+0.390677571, +0.120144054, +1.00000000
651T
653T
654getMatDet(rand)
655+0.830212474
656call setCovRand(rngf, rand, onion, eta)
657onion%info
658+0
659rand
660+1.00000000, +0.421261787E-2, -0.119264908
661+0.421261787E-2, +1.00000000, +0.560555197E-1
662-0.119264908, +0.560555197E-1, +1.00000000
664T
666T
667getMatDet(rand)
668+0.982559562
669call setCovRand(rngf, rand, dvine, eta, scale)
670rand
671+64.0000000, +15.8274841, +5.32799530
672+15.8274841, +64.0000000, +21.2583618
673+5.32799530, +21.2583618, +64.0000000
675T
677T
678getMatDet(rand)
679+218957.234
680
681call setCovRand(rngf, rand, onion, eta, scale)
682onion%info
683+0
684rand
685+64.0000000, +37.5149689, -15.8536568
686+37.5149689, +64.0000000, -239.270523
687-15.8536568, -239.270523, +64.0000000
689F
691T
692getMatDet(rand)
693-3223426.25
694
695
696eta = getUnifRand(1, 10)
697ndim = getUnifRand(2, 5)
698call setResized(rand, [ndim, ndim])
699scale = getUnifRand(1, 10)
700
701call setCovRand(rngf, rand, dvine, eta)
702onion%info
703+0
704rand
705+1.00000000, -0.520092070, +0.461997986
706-0.520092070, +1.00000000, -0.379322737
707+0.461997986, -0.379322737, +1.00000000
709T
711T
712getMatDet(rand)
713+0.554464817
714call setCovRand(rngf, rand, onion, eta)
715onion%info
716+0
717rand
718+1.00000000, +0.327518344, +0.583361201E-1
719+0.327518344, +1.00000000, +0.845330536
720+0.583361201E-1, +0.845330536, +1.00000000
722T
724T
725getMatDet(rand)
726+0.207046866
727call setCovRand(rngf, rand, dvine, eta, scale)
728rand
729+64.0000000, -8.00318909, -15.8723068
730-8.00318909, +64.0000000, +26.2841110
731-15.8723068, +26.2841110, +64.0000000
733T
735T
736getMatDet(rand)
737+204384.203
738
739call setCovRand(rngf, rand, onion, eta, scale)
740onion%info
741+0
742rand
743+64.0000000, +25.5941772, +9.11215687
744+25.5941772, +64.0000000, -85.0748978
745+9.11215687, -85.0748978, +64.0000000
747F
749T
750getMatDet(rand)
751-287991.281
752
753
754eta = getUnifRand(1, 10)
755ndim = getUnifRand(2, 5)
756call setResized(rand, [ndim, ndim])
757scale = getUnifRand(1, 10)
758
759call setCovRand(rngf, rand, dvine, eta)
760onion%info
761+0
762rand
763+1.00000000, -0.432123780
764-0.432123780, +1.00000000
766T
768T
769getMatDet(rand)
770+0.813269019
771call setCovRand(rngf, rand, onion, eta)
772onion%info
773+0
774rand
775+1.00000000, +0.155968428
776+0.155968428, +1.00000000
778T
780T
781getMatDet(rand)
782+0.975673854
783call setCovRand(rngf, rand, dvine, eta, scale)
784rand
785+49.0000000, -3.63959265
786-3.63959265, +49.0000000
788T
790T
791getMatDet(rand)
792+2387.75342
793
794call setCovRand(rngf, rand, onion, eta, scale)
795onion%info
796+0
797rand
798+49.0000000, -4.23254156
799-4.23254156, +49.0000000
801T
803T
804getMatDet(rand)
805+2383.08545
806
807
808eta = getUnifRand(1, 10)
809ndim = getUnifRand(2, 5)
810call setResized(rand, [ndim, ndim])
811scale = getUnifRand(1, 10)
812
813call setCovRand(rngf, rand, dvine, eta)
814onion%info
815+0
816rand
817+1.00000000, -0.174169600, +0.262220860
818-0.174169600, +1.00000000, -0.186970353
819+0.262220860, -0.186970353, +1.00000000
821T
823T
824getMatDet(rand)
825+0.883025467
826call setCovRand(rngf, rand, onion, eta)
827onion%info
828+0
829rand
830+1.00000000, +0.522613525, +0.670238212E-1
831+0.522613525, +1.00000000, -0.122490957
832+0.670238212E-1, -0.122490957, +1.00000000
834T
836T
837getMatDet(rand)
838+0.698797762
839call setCovRand(rngf, rand, dvine, eta, scale)
840rand
841+25.0000000, -5.14816952, +1.68904364
842-5.14816952, +25.0000000, -0.894779325
843+1.68904364, -0.894779325, +25.0000000
845T
847T
848getMatDet(rand)
849+14886.6318
850
851call setCovRand(rngf, rand, onion, eta, scale)
852onion%info
853+0
854rand
855+25.0000000, +6.03189182, +4.14547825
856+6.03189182, +25.0000000, -17.3841953
857+4.14547825, -17.3841953, +25.0000000
859T
861T
862getMatDet(rand)
863+5861.14062
864
865
866eta = getUnifRand(1, 10)
867ndim = getUnifRand(2, 5)
868call setResized(rand, [ndim, ndim])
869scale = getUnifRand(1, 10)
870
871call setCovRand(rngf, rand, dvine, eta)
872onion%info
873+0
874rand
875+1.00000000, -0.327304602E-1
876-0.327304602E-1, +1.00000000
878T
880T
881getMatDet(rand)
882+0.998928726
883call setCovRand(rngf, rand, onion, eta)
884onion%info
885+0
886rand
887+1.00000000, +0.193414688E-1
888+0.193414688E-1, +1.00000000
890T
892T
893getMatDet(rand)
894+0.999625921
895call setCovRand(rngf, rand, dvine, eta, scale)
896rand
897+36.0000000, -2.42128372
898-2.42128372, +36.0000000
900T
902T
903getMatDet(rand)
904+1290.13745
905
906call setCovRand(rngf, rand, onion, eta, scale)
907onion%info
908+0
909rand
910+36.0000000, +16.2747421
911+16.2747421, +36.0000000
913T
915T
916getMatDet(rand)
917+1031.13281
918
919
920eta = getUnifRand(1, 10)
921ndim = getUnifRand(2, 5)
922call setResized(rand, [ndim, ndim])
923scale = getUnifRand(1, 10)
924
925call setCovRand(rngf, rand, dvine, eta)
926onion%info
927+0
928rand
929+1.00000000, -0.201402664
930-0.201402664, +1.00000000
932T
934T
935getMatDet(rand)
936+0.959436953
937call setCovRand(rngf, rand, onion, eta)
938onion%info
939+0
940rand
941+1.00000000, -0.232133389
942-0.232133389, +1.00000000
944T
946T
947getMatDet(rand)
948+0.946114063
949call setCovRand(rngf, rand, dvine, eta, scale)
950rand
951+36.0000000, +0.454919815
952+0.454919815, +36.0000000
954T
956T
957getMatDet(rand)
958+1295.79309
959
960call setCovRand(rngf, rand, onion, eta, scale)
961onion%info
962+0
963rand
964+36.0000000, +20.6934013
965+20.6934013, +36.0000000
967T
969T
970getMatDet(rand)
971+867.783142
972
973
974eta = getUnifRand(1, 10)
975ndim = getUnifRand(2, 5)
976call setResized(rand, [ndim, ndim])
977scale = getUnifRand(1, 10)
978
979call setCovRand(rngf, rand, dvine, eta)
980onion%info
981+0
982rand
983+1.00000000, -0.372969985, -0.374505222, +0.109540582
984-0.372969985, +1.00000000, +0.165067613E-1, +0.218678236
985-0.374505222, +0.165067613E-1, +1.00000000, -0.583541036
986+0.109540582, +0.218678236, -0.583541036, +1.00000000
988T
990T
991getMatDet(rand)
992+0.439671546
993call setCovRand(rngf, rand, onion, eta)
994onion%info
995+0
996rand
997+1.00000000, -0.114963531, -0.298758388, -0.104634173
998-0.114963531, +1.00000000, -0.330263108, -0.511745036
999-0.298758388, -0.330263108, +1.00000000, -0.540939808
1000-0.104634173, -0.511745036, -0.540939808, +1.00000000
1001isMatClass(rand, posdefmat)
1002F
1003isMatClass(rand, hermitian)
1004T
1005getMatDet(rand)
1006-0.341053568E-1
1007call setCovRand(rngf, rand, dvine, eta, scale)
1008rand
1009+81.0000000, +26.6339645, +10.7979622, +31.2080002
1010+26.6339645, +81.0000000, +2.52577424, +32.4428215
1011+10.7979622, +2.52577424, +81.0000000, +14.3352013
1012+31.2080002, +32.4428215, +14.3352013, +81.0000000
1013isMatClass(rand, posdefmat)
1014T
1015isMatClass(rand, hermitian)
1016T
1017getMatDet(rand)
1018+28292602.0
1019
1020call setCovRand(rngf, rand, onion, eta, scale)
1021onion%info
1022+3
1023rand
1024+1.00000000, -0.240146160, -0.407497138E-1, +31.2080002
1025+26.6339645, +1.00000000, -1.17803288, +32.4428215
1026+10.7979622, +2.52577424, +1.00000000, +14.3352013
1027+31.2080002, +32.4428215, +14.3352013, +81.0000000
1028isMatClass(rand, posdefmat)
1029F
1030isMatClass(rand, hermitian)
1031F
1032getMatDet(rand)
1033+10628.1533
1034
1035
1036eta = getUnifRand(1, 10)
1037ndim = getUnifRand(2, 5)
1038call setResized(rand, [ndim, ndim])
1039scale = getUnifRand(1, 10)
1040
1041call setCovRand(rngf, rand, dvine, eta)
1042onion%info
1043+3
1044rand
1045+1.00000000, +0.127711296, +0.515108705, -0.743637085E-1
1046+0.127711296, +1.00000000, +0.865820795E-2, -0.664231703E-1
1047+0.515108705, +0.865820795E-2, +1.00000000, +0.640287697E-1
1048-0.743637085E-1, -0.664231703E-1, +0.640287697E-1, +1.00000000
1049isMatClass(rand, posdefmat)
1050T
1051isMatClass(rand, hermitian)
1052T
1053getMatDet(rand)
1054+0.703422129
1055call setCovRand(rngf, rand, onion, eta)
1056onion%info
1057+0
1058rand
1059+1.00000000, -0.393271923, -0.177916437, +0.103122868
1060-0.393271923, +1.00000000, -0.244356006, -0.391173720
1061-0.177916437, -0.244356006, +1.00000000, +0.610120595E-3
1062+0.103122868, -0.391173720, +0.610120595E-3, +1.00000000
1063isMatClass(rand, posdefmat)
1064T
1065isMatClass(rand, hermitian)
1066T
1067getMatDet(rand)
1068+0.596956670
1069call setCovRand(rngf, rand, dvine, eta, scale)
1070rand
1071+81.0000000, +31.7667694, -19.5949402, +9.43964005
1072+31.7667694, +81.0000000, -48.8774033, +18.2937279
1073-19.5949402, -48.8774033, +81.0000000, +20.7111053
1074+9.43964005, +18.2937279, +20.7111053, +81.0000000
1075isMatClass(rand, posdefmat)
1076T
1077isMatClass(rand, hermitian)
1078T
1079getMatDet(rand)
1080+16356823.0
1081
1082call setCovRand(rngf, rand, onion, eta, scale)
1083onion%info
1084+3
1085rand
1086+1.00000000, +0.143210888E-1, +0.346610427, +9.43964005
1087+31.7667694, +1.00000000, +11.1167889, +18.2937279
1088-19.5949402, -48.8774033, +1.00000000, +20.7111053
1089+9.43964005, +18.2937279, +20.7111053, +81.0000000
1090isMatClass(rand, posdefmat)
1091F
1092isMatClass(rand, hermitian)
1093F
1094getMatDet(rand)
1095+276834.656
1096
1097
1098eta = getUnifRand(1, 10)
1099ndim = getUnifRand(2, 5)
1100call setResized(rand, [ndim, ndim])
1101scale = getUnifRand(1, 10)
1102
1103call setCovRand(rngf, rand, dvine, eta)
1104onion%info
1105+3
1106rand
1107+1.00000000, -0.801032782E-1, +0.163490295
1108-0.801032782E-1, +1.00000000, +0.158034563
1109+0.163490295, +0.158034563, +1.00000000
1110isMatClass(rand, posdefmat)
1111T
1112isMatClass(rand, hermitian)
1113T
1114getMatDet(rand)
1115+0.937740088
1116call setCovRand(rngf, rand, onion, eta)
1117onion%info
1118+0
1119rand
1120+1.00000000, -0.379363418, +0.492337257
1121-0.379363418, +1.00000000, -0.303344548
1122+0.492337257, -0.303344548, +1.00000000
1123isMatClass(rand, posdefmat)
1124T
1125isMatClass(rand, hermitian)
1126T
1127getMatDet(rand)
1128+0.634983659
1129call setCovRand(rngf, rand, dvine, eta, scale)
1130rand
1131+36.0000000, -13.0127354, -2.69560933
1132-13.0127354, +36.0000000, +2.05713034
1133-2.69560933, +2.05713034, +36.0000000
1134isMatClass(rand, posdefmat)
1135T
1136isMatClass(rand, hermitian)
1137T
1138getMatDet(rand)
1139+40290.4570
1140
1141call setCovRand(rngf, rand, onion, eta, scale)
1142onion%info
1143+0
1144rand
1145+36.0000000, +1.91980505, +3.56779456
1146+1.91980505, +36.0000000, -34.6922646
1147+3.56779456, -34.6922646, +36.0000000
1148isMatClass(rand, posdefmat)
1149T
1150isMatClass(rand, hermitian)
1151T
1152getMatDet(rand)
1153+2261.90503
1154
1155
1156ndim = getUnifRand(2, 10)
1157call setResized(rand, [ndim, ndim])
1158
1159call setCovRand(rngf, rand, dvine, eta = 0._TKG, scale = [(real(itry, TKG), itry = 1, ndim)])
1160rand
1161+1.00000000, +0.197524786, -1.08652520, -0.899771214, +1.60182416, +0.350251436, +1.98801780, -4.89459229, -5.81478357
1162+0.197524786, +4.00000000, -2.07108426, +1.63782883, +2.10013771, -5.61528635, -2.64949918, -4.22243834, -6.41157246
1163-1.08652520, -2.07108426, +9.00000000, -3.22147989, -0.655206084, +6.93056011, -7.32249546, +3.19860888, +11.3008451
1164-0.899771214, +1.63782883, -3.22147989, +16.0000000, +2.31580687, -9.37684059, -7.61445665, -1.33025646, +2.76267791
1165+1.60182416, +2.10013771, -0.655206084, +2.31580687, +25.0000000, -5.97565365, +7.52172184, -21.8173752, +4.94239283
1166+0.350251436, -5.61528635, +6.93056011, -9.37684059, -5.97565365, +36.0000000, +15.2414961, +13.5586033, +20.8919735
1167+1.98801780, -2.64949918, -7.32249546, -7.61445665, +7.52172184, +15.2414961, +49.0000000, -11.2109098, +13.1339788
1168-4.89459229, -4.22243834, +3.19860888, -1.33025646, -21.8173752, +13.5586033, -11.2109098, +64.0000000, +16.8706017
1169-5.81478357, -6.41157246, +11.3008451, +2.76267791, +4.94239283, +20.8919735, +13.1339788, +16.8706017, +81.0000000
1170isMatClass(rand, posdefmat)
1171T
1172
1173call setCovRand(rngf, rand, onion, eta = 0._TKG, scale = [(real(itry, TKG), itry = 1, ndim)])
1174onion%info
1175+5
1176rand
1177+1.00000000, +0.860533714E-1, +0.167458639, -0.874029770E-1, -0.472449660, +0.350251436, +1.98801780, -4.89459229, -5.81478357
1178+0.197524786, +1.00000000, -0.565134287E-1, +0.320119888, -0.114935264, -5.61528635, -2.64949918, -4.22243834, -6.41157246
1179-1.08652520, -2.07108426, +1.00000000, -0.381791294, +0.627015650, +6.93056011, -7.32249546, +3.19860888, +11.3008451
1180-0.899771214, +1.63782883, -3.22147989, +1.00000000, -1.98168969, -9.37684059, -7.61445665, -1.33025646, +2.76267791
1181+1.60182416, +2.10013771, -0.655206084, +2.31580687, +1.00000000, -5.97565365, +7.52172184, -21.8173752, +4.94239283
1182+0.350251436, -5.61528635, +6.93056011, -9.37684059, -5.97565365, +36.0000000, +15.2414961, +13.5586033, +20.8919735
1183+1.98801780, -2.64949918, -7.32249546, -7.61445665, +7.52172184, +15.2414961, +49.0000000, -11.2109098, +13.1339788
1184-4.89459229, -4.22243834, +3.19860888, -1.33025646, -21.8173752, +13.5586033, -11.2109098, +64.0000000, +16.8706017
1185-5.81478357, -6.41157246, +11.3008451, +2.76267791, +4.94239283, +20.8919735, +13.1339788, +16.8706017, +81.0000000
1186isMatClass(rand, posdefmat)
1187F
1188
1189
Test:
test_pm_distCov
Todo:
High Priority: The current implementation of this generic interface uses a naive method of computing the Cholesky factorization with a default matrix packing for the Onion method.
The RFP packing format must be also implemented for this generic interface.
Todo:
High Priority: The current implementation of the Gram method can be significantly improved, both computationally and functionally.


Final Remarks


If you believe this algorithm or its documentation can be improved, we appreciate your contribution and help to edit this page's documentation and source file on GitHub.
For details on the naming abbreviations, see this page.
For details on the naming conventions, see this page.
This software is distributed under the MIT license with additional terms outlined below.

  1. If you use any parts or concepts from this library to any extent, please acknowledge the usage by citing the relevant publications of the ParaMonte library.
  2. If you regenerate any parts/ideas from this library in a programming environment other than those currently supported by this ParaMonte library (i.e., other than C, C++, Fortran, MATLAB, Python, R), please also ask the end users to cite this original ParaMonte library.

This software is available to the public under a highly permissive license.
Help us justify its continued development and maintenance by acknowledging its benefit to society, distributing it, and contributing to it.

Author:
Amir Shahmoradi, Monday March 6, 2017, 3:22 pm, Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin.

Definition at line 787 of file pm_distCov.F90.


The documentation for this interface was generated from the following file: