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

Generate a refined version of the input array where the sequentially unweighted entries along the specified dimension of array are skipped every skip to create a refined weighted output array of size rsize.
More...

Detailed Description

Generate a refined version of the input array where the sequentially unweighted entries along the specified dimension of array are skipped every skip to create a refined weighted output array of size rsize.

Parameters
[in,out]array: The input/output contiguous array of shape (:) or (:,:) of either
  1. type character of kind any supported by the processor (e.g., SK, SKA, SKD , or SKU) or
  2. type integer of kind any supported by the processor (e.g., IK, IK8, IK16, IK32, or IK64) or
  3. type logical of kind any supported by the processor (e.g., LK) or
  4. type complex of kind any supported by the processor (e.g., CK, CK32, CK64, or CK128) or
  5. type real of kind any supported by the processor (e.g., RK, RK32, RK64, or RK128)
or scalar character of kind any supported by the processor (e.g., SK, SKA, SKD , or SKU).
On input, it contains the elements to be refined.
On output, the first rsize entries of array contain the condensed elements of the original input array, such that no two adjacent entries in array from entry 1 to entry rsize are duplicates.
[in]dim: The input scalar of type integer of default kind IK representing the axis of array(:,:) along which array must be refined.
(optional, it must be present if and only if array is of shape (:,:).)
[in,out]weight: The input/output vector of
  1. type integer of default kind IK,
of size size(array, dim) containing the weights of individual elements of the input array.
On output, the first rsize elements of weight are written with the new refined weights of the corresponding elements in the output array.
[in]skip: The input scalar of type integer of default kind IK representing the number of elements to skip in the input sequence.
[out]rsize: The output scalar of type integer of default kind IK such that,
  1. weight(1:rsize) contains the weights of the refined elements of the input array, and
  2. array(1:rsize) or array(:,1:rsize) or array(1:rsize,:) contains the refined elements of the input array.


Possible calling interfaces ⛓

use pm_arrayRefine, only: setRefined
call setRefined(array, weight(:), rsize) ! scalar character objects.
call setRefined(array(:), weight(:), rsize) ! all intrinsic array objects.
call setRefined(array(:,:), dim, weight(:), rsize) ! all intrinsic array objects.
!
pm_arrayRefine::setRefined
Generate a refined version of the input array where the sequentially unweighted entries along the spe...
Definition: pm_arrayRefine.F90:1018
pm_arrayRefine
This module contains procedures and generic interfaces for refining (thinning) (weighted) arrays of a...
Definition: pm_arrayRefine.F90:60
Warning
The condition 0 < skip must hold.
The condition dim == 1 .or. dim == 2 must hold.
The condition size(weight) == size(array) must hold when array is rank 1.
The condition size(weight) == size(array, dim) must hold when array is rank 2.
These conditions are verified only if the library is built with the preprocessor macro CHECK_ENABLED=1.
The pure procedure(s) documented herein become impure when the ParaMonte library is compiled with preprocessor macro CHECK_ENABLED=1.
By default, these procedures are pure in release build and impure in debug and testing builds.
Note
See pm_arrayCopy for refining unweighted strings and arrays.
See also
setRefined
getCompact
getVerbose
pm_arrayCopy


Example usage ⛓

1program example
2
3 use pm_arrayVerbose, only: getVerbose
4 use pm_arrayRefine, only: setRefined
5 use pm_arrayChoice, only: getChoice
6 use pm_distUnif, only: getUnifRand
7 use pm_distBern, only: isHead
8 use pm_kind, only: SK, IK, LK
9 use pm_io, only: display_type
10
11 implicit none
12
13 type(display_type) :: disp
14 integer(IK) :: rsize, nsam, itry, ntry = 10
15 disp = display_type(file = "main.out.F90")
16
17 call disp%skip()
18 call disp%show("!%%%%%%%%%%%%%%%%%")
19 call disp%show("! Refine 1D array.")
20 call disp%show("!%%%%%%%%%%%%%%%%%")
21 call disp%skip()
22
23 block
24
25 block
26 use pm_kind, only: TKG => SK ! all kinds are supported.
27 character(:,TKG), allocatable :: array, arref
28 integer(IK), allocatable :: weight(:), weiref(:)
29 integer(IK) :: skip
30 do itry = 1, ntry
31 call disp%show("nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
32 nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
33 call disp%show("nsam")
34 call disp%show( nsam )
35 call disp%show("skip")
36 call disp%show( skip )
37 call disp%show("weight = getUnifRand(-1, 9, nsam)")
38 weight = getUnifRand(-1, 9, nsam)
39 call disp%show("weight")
40 call disp%show( weight )
41 call disp%show("array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))")
42 array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
43 call disp%show("array")
44 call disp%show( array , deliml = TKG_"""" )
45 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0))")
46 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0)) , deliml = TKG_"""" )
47 call disp%show("arref = array")
48 arref = array
49 call disp%show("weiref = weight")
50 weiref = weight
51 call disp%show("call setRefined(arref, weiref, skip, rsize)")
52 call setRefined(arref, weiref, skip, rsize)
53 call disp%show("rsize")
54 call disp%show( rsize )
55 call disp%show("weiref(1:rsize)")
56 call disp%show( weiref(1:rsize) )
57 call disp%show("arref")
58 call disp%show( arref , deliml = TKG_"""" )
59 call disp%skip()
60 end do
61 end block
62
63 block
64 use pm_kind, only: TKG => SK ! all kinds are supported.
65 character(2,TKG), allocatable :: array(:), arref(:)
66 integer(IK), allocatable :: weight(:), weiref(:)
67 integer(IK) :: skip
68 do itry = 1, ntry
69 call disp%show("nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
70 nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
71 call disp%show("nsam")
72 call disp%show( nsam )
73 call disp%show("skip")
74 call disp%show( skip )
75 call disp%show("weight = getUnifRand(-1, 9, nsam)")
76 weight = getUnifRand(-1, 9, nsam)
77 call disp%show("weight")
78 call disp%show( weight )
79 call disp%show("array = getUnifRand('AA', 'ZZ', nsam)")
80 array = getUnifRand('AA', 'ZZ', nsam)
81 call disp%show("array")
82 call disp%show( array , deliml = TKG_"""" )
83 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0))")
84 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0)) , deliml = TKG_"""" )
85 call disp%show("arref = array")
86 arref = array
87 call disp%show("weiref = weight")
88 weiref = weight
89 call disp%show("call setRefined(arref, weiref, skip, rsize)")
90 call setRefined(arref, weiref, skip, rsize)
91 call disp%show("rsize")
92 call disp%show( rsize )
93 call disp%show("weiref(1:rsize)")
94 call disp%show( weiref(1:rsize) )
95 call disp%show("arref")
96 call disp%show( arref , deliml = TKG_"""" )
97 call disp%skip()
98 end do
99 end block
100
101 block
102 use pm_kind, only: TKG => IK ! all kinds are supported.
103 integer(TKG), allocatable :: array(:), arref(:)
104 integer(IK), allocatable :: weight(:), weiref(:)
105 integer(IK) :: skip
106 do itry = 1, ntry
107 call disp%show("nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
108 nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
109 call disp%show("nsam")
110 call disp%show( nsam )
111 call disp%show("skip")
112 call disp%show( skip )
113 call disp%show("weight = getUnifRand(-1, 9, nsam)")
114 weight = getUnifRand(-1, 9, nsam)
115 call disp%show("weight")
116 call disp%show( weight )
117 call disp%show("array = getUnifRand(0, 9, nsam)")
118 array = getUnifRand(0, 9, nsam)
119 call disp%show("array")
120 call disp%show( array )
121 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0))")
122 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0)) )
123 call disp%show("arref = array")
124 arref = array
125 call disp%show("weiref = weight")
126 weiref = weight
127 call disp%show("call setRefined(arref, weiref, skip, rsize)")
128 call setRefined(arref, weiref, skip, rsize)
129 call disp%show("rsize")
130 call disp%show( rsize )
131 call disp%show("weiref(1:rsize)")
132 call disp%show( weiref(1:rsize) )
133 call disp%show("arref(1:rsize)")
134 call disp%show( arref(1:rsize) )
135 call disp%skip()
136 end do
137 end block
138
139 block
140 use pm_kind, only: TKG => LK ! all kinds are supported.
141 logical(TKG), allocatable :: array(:), arref(:)
142 integer(IK), allocatable :: weight(:), weiref(:)
143 integer(IK) :: skip
144 do itry = 1, ntry
145 call disp%show("nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
146 nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
147 call disp%show("nsam")
148 call disp%show( nsam )
149 call disp%show("skip")
150 call disp%show( skip )
151 call disp%show("weight = getUnifRand(-1, 9, nsam)")
152 weight = getUnifRand(-1, 9, nsam)
153 call disp%show("weight")
154 call disp%show( weight )
155 call disp%show("array = getUnifRand(.false., .true., nsam)")
156 array = getUnifRand(.false., .true., nsam)
157 call disp%show("array")
158 call disp%show( array )
159 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0))")
160 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0)) )
161 call disp%show("arref = array")
162 arref = array
163 call disp%show("weiref = weight")
164 weiref = weight
165 call disp%show("call setRefined(arref, weiref, skip, rsize)")
166 call setRefined(arref, weiref, skip, rsize)
167 call disp%show("rsize")
168 call disp%show( rsize )
169 call disp%show("weiref(1:rsize)")
170 call disp%show( weiref(1:rsize) )
171 call disp%show("arref(1:rsize)")
172 call disp%show( arref(1:rsize) )
173 call disp%skip()
174 end do
175 end block
176
177 block
178 use pm_kind, only: TKG => CKS ! all kinds are supported.
179 complex(TKG), allocatable :: array(:), arref(:)
180 integer(IK), allocatable :: weight(:), weiref(:)
181 integer(IK) :: skip
182 do itry = 1, ntry
183 call disp%show("nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
184 nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
185 call disp%show("nsam")
186 call disp%show( nsam )
187 call disp%show("skip")
188 call disp%show( skip )
189 call disp%show("weight = getUnifRand(-1, 9, nsam)")
190 weight = getUnifRand(-1, 9, nsam)
191 call disp%show("weight")
192 call disp%show( weight )
193 call disp%show("array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)")
194 array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
195 call disp%show("array")
196 call disp%show( array )
197 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0))")
198 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0)) )
199 call disp%show("arref = array")
200 arref = array
201 call disp%show("weiref = weight")
202 weiref = weight
203 call disp%show("call setRefined(arref, weiref, skip, rsize)")
204 call setRefined(arref, weiref, skip, rsize)
205 call disp%show("rsize")
206 call disp%show( rsize )
207 call disp%show("weiref(1:rsize)")
208 call disp%show( weiref(1:rsize) )
209 call disp%show("arref(1:rsize)")
210 call disp%show( arref(1:rsize) )
211 call disp%skip()
212 end do
213 end block
214
215 block
216 use pm_kind, only: TKG => RKS ! all kinds are supported.
217 real(TKG), allocatable :: array(:), arref(:)
218 integer(IK), allocatable :: weight(:), weiref(:)
219 integer(IK) :: skip
220 do itry = 1, ntry
221 call disp%show("nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
222 nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
223 call disp%show("nsam")
224 call disp%show( nsam )
225 call disp%show("skip")
226 call disp%show( skip )
227 call disp%show("weight = getUnifRand(-1, 9, nsam)")
228 weight = getUnifRand(-1, 9, nsam)
229 call disp%show("weight")
230 call disp%show( weight )
231 call disp%show("array = getUnifRand(0, 9, nsam)")
232 array = getUnifRand(0, 9, nsam)
233 call disp%show("array")
234 call disp%show( array )
235 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0))")
236 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0)) )
237 call disp%show("arref = array")
238 arref = array
239 call disp%show("weiref = weight")
240 weiref = weight
241 call disp%show("call setRefined(arref, weiref, skip, rsize)")
242 call setRefined(arref, weiref, skip, rsize)
243 call disp%show("rsize")
244 call disp%show( rsize )
245 call disp%show("weiref(1:rsize)")
246 call disp%show( weiref(1:rsize) )
247 call disp%show("arref(1:rsize)")
248 call disp%show( arref(1:rsize) )
249 call disp%skip()
250 end do
251 end block
252
253 end block
254
255 call disp%skip()
256 call disp%show("!%%%%%%%%%%%%%%%%%")
257 call disp%show("! Refine 2D array.")
258 call disp%show("!%%%%%%%%%%%%%%%%%")
259 call disp%skip()
260
261 block
262
263 integer(IK) :: dim, ndim
264
265 block
266 use pm_kind, only: TKG => SK ! all kinds are supported.
267 character(2,TKG), allocatable :: array(:,:), arref(:,:)
268 integer(IK), allocatable :: weight(:), weiref(:)
269 integer(IK) :: skip
270 do itry = 1, ntry
271 call disp%show("dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
272 dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
273 call disp%show("[dim, ndim, nsam]")
274 call disp%show( [dim, ndim, nsam] )
275 call disp%show("skip")
276 call disp%show( skip )
277 call disp%show("weight = getUnifRand(-1, 9, nsam)")
278 weight = getUnifRand(-1, 9, nsam)
279 call disp%show("weight")
280 call disp%show( weight )
281 call disp%show("array = getUnifRand('AA', 'ZZ', ndim, nsam)")
282 array = getUnifRand('AA', 'ZZ', ndim, nsam)
283 call disp%show("array")
284 call disp%show( array , deliml = TKG_"""" )
285 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0), dim)")
286 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0), dim) , deliml = TKG_"""" )
287 call disp%show("arref = array")
288 arref = array
289 call disp%show("weiref = weight")
290 weiref = weight
291 call disp%show("call setRefined(arref, dim, weiref, skip, rsize)")
292 call setRefined(arref, dim, weiref, skip, rsize)
293 call disp%show("rsize")
294 call disp%show( rsize )
295 call disp%show("weiref(1:rsize)")
296 call disp%show( weiref(1:rsize) )
297 call disp%show("arref(:, 1:rsize)")
298 call disp%show( arref(:, 1:rsize) , deliml = TKG_"""" )
299 call disp%show("weiref = weight")
300 weiref = weight
301 call disp%show("arref = transpose(array)")
302 arref = transpose(array)
303 call disp%show("call setRefined(arref, 3_IK - dim, weiref, skip, rsize)")
304 call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
305 call disp%show("rsize")
306 call disp%show( rsize )
307 call disp%show("weiref(1:rsize)")
308 call disp%show( weiref(1:rsize) )
309 call disp%show("arref(1:rsize, :)")
310 call disp%show( arref(1:rsize, :) , deliml = TKG_"""" )
311 call disp%skip()
312 end do
313 end block
314
315 block
316 use pm_kind, only: TKG => IK ! all kinds are supported.
317 integer(TKG), allocatable :: array(:,:), arref(:,:)
318 integer(IK), allocatable :: weight(:), weiref(:)
319 integer(IK) :: skip
320 do itry = 1, ntry
321 call disp%show("dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
322 dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
323 call disp%show("[dim, ndim, nsam]")
324 call disp%show( [dim, ndim, nsam] )
325 call disp%show("skip")
326 call disp%show( skip )
327 call disp%show("weight = getUnifRand(-1, 9, nsam)")
328 weight = getUnifRand(-1, 9, nsam)
329 call disp%show("weight")
330 call disp%show( weight )
331 call disp%show("array = getUnifRand(0, 9, ndim, nsam)")
332 array = getUnifRand(0, 9, ndim, nsam)
333 call disp%show("array")
334 call disp%show( array )
335 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0), dim)")
336 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0), dim) )
337 call disp%show("arref = array")
338 arref = array
339 call disp%show("weiref = weight")
340 weiref = weight
341 call disp%show("call setRefined(arref, dim, weiref, skip, rsize)")
342 call setRefined(arref, dim, weiref, skip, rsize)
343 call disp%show("rsize")
344 call disp%show( rsize )
345 call disp%show("weiref(1:rsize)")
346 call disp%show( weiref(1:rsize) )
347 call disp%show("arref(:, 1:rsize)")
348 call disp%show( arref(:, 1:rsize) )
349 call disp%show("weiref = weight")
350 weiref = weight
351 call disp%show("arref = transpose(array)")
352 arref = transpose(array)
353 call disp%show("call setRefined(arref, 3_IK - dim, weiref, skip, rsize)")
354 call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
355 call disp%show("rsize")
356 call disp%show( rsize )
357 call disp%show("weiref(1:rsize)")
358 call disp%show( weiref(1:rsize) )
359 call disp%show("arref(1:rsize, :)")
360 call disp%show( arref(1:rsize, :) )
361 call disp%skip()
362 end do
363 end block
364
365 block
366 use pm_kind, only: TKG => LK ! all kinds are supported.
367 logical(TKG), allocatable :: array(:,:), arref(:,:)
368 integer(IK), allocatable :: weight(:), weiref(:)
369 integer(IK) :: skip
370 do itry = 1, ntry
371 call disp%show("dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
372 dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
373 call disp%show("[dim, ndim, nsam]")
374 call disp%show( [dim, ndim, nsam] )
375 call disp%show("skip")
376 call disp%show( skip )
377 call disp%show("weight = getUnifRand(-1, 9, nsam)")
378 weight = getUnifRand(-1, 9, nsam)
379 call disp%show("weight")
380 call disp%show( weight )
381 call disp%show("array = getUnifRand(.false., .true., ndim, nsam)")
382 array = getUnifRand(.false., .true., ndim, nsam)
383 call disp%show("array")
384 call disp%show( array )
385 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0), dim)")
386 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0), dim) )
387 call disp%show("arref = array")
388 arref = array
389 call disp%show("weiref = weight")
390 weiref = weight
391 call disp%show("call setRefined(arref, dim, weiref, skip, rsize)")
392 call setRefined(arref, dim, weiref, skip, rsize)
393 call disp%show("rsize")
394 call disp%show( rsize )
395 call disp%show("weiref(1:rsize)")
396 call disp%show( weiref(1:rsize) )
397 call disp%show("arref(:, 1:rsize)")
398 call disp%show( arref(:, 1:rsize) )
399 call disp%show("weiref = weight")
400 weiref = weight
401 call disp%show("arref = transpose(array)")
402 arref = transpose(array)
403 call disp%show("call setRefined(arref, 3_IK - dim, weiref, skip, rsize)")
404 call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
405 call disp%show("rsize")
406 call disp%show( rsize )
407 call disp%show("weiref(1:rsize)")
408 call disp%show( weiref(1:rsize) )
409 call disp%show("arref(1:rsize, :)")
410 call disp%show( arref(1:rsize, :) )
411 call disp%skip()
412 end do
413 end block
414
415 block
416 use pm_kind, only: TKG => CKS ! all kinds are supported.
417 complex(TKG), allocatable :: array(:,:), arref(:,:)
418 integer(IK), allocatable :: weight(:), weiref(:)
419 integer(IK) :: skip
420 do itry = 1, ntry
421 call disp%show("dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
422 dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
423 call disp%show("[dim, ndim, nsam]")
424 call disp%show( [dim, ndim, nsam] )
425 call disp%show("skip")
426 call disp%show( skip )
427 call disp%show("weight = getUnifRand(-1, 9, nsam)")
428 weight = getUnifRand(-1, 9, nsam)
429 call disp%show("weight")
430 call disp%show( weight )
431 call disp%show("array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)")
432 array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
433 call disp%show("array")
434 call disp%show( array )
435 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0), dim)")
436 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0), dim) )
437 call disp%show("arref = array")
438 arref = array
439 call disp%show("weiref = weight")
440 weiref = weight
441 call disp%show("call setRefined(arref, dim, weiref, skip, rsize)")
442 call setRefined(arref, dim, weiref, skip, rsize)
443 call disp%show("rsize")
444 call disp%show( rsize )
445 call disp%show("weiref(1:rsize)")
446 call disp%show( weiref(1:rsize) )
447 call disp%show("arref(:, 1:rsize)")
448 call disp%show( arref(:, 1:rsize) )
449 call disp%show("weiref = weight")
450 weiref = weight
451 call disp%show("arref = transpose(array)")
452 arref = transpose(array)
453 call disp%show("call setRefined(arref, 3_IK - dim, weiref, skip, rsize)")
454 call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
455 call disp%show("rsize")
456 call disp%show( rsize )
457 call disp%show("weiref(1:rsize)")
458 call disp%show( weiref(1:rsize) )
459 call disp%show("arref(1:rsize, :)")
460 call disp%show( arref(1:rsize, :) )
461 call disp%skip()
462 end do
463 end block
464
465 block
466 use pm_kind, only: TKG => RKS ! all kinds are supported.
467 real(TKG), allocatable :: array(:,:), arref(:,:)
468 integer(IK), allocatable :: weight(:), weiref(:)
469 integer(IK) :: skip
470 do itry = 1, ntry
471 call disp%show("dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)")
472 dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
473 call disp%show("[dim, ndim, nsam]")
474 call disp%show( [dim, ndim, nsam] )
475 call disp%show("skip")
476 call disp%show( skip )
477 call disp%show("weight = getUnifRand(-1, 9, nsam)")
478 weight = getUnifRand(-1, 9, nsam)
479 call disp%show("weight")
480 call disp%show( weight )
481 call disp%show("array = getUnifRand(0, 9, ndim, nsam)")
482 array = getUnifRand(0, 9, ndim, nsam)
483 call disp%show("array")
484 call disp%show( array )
485 call disp%show("getVerbose(array, weight, sum(weight, mask = weight > 0), dim)")
486 call disp%show( getVerbose(array, weight, sum(weight, mask = weight > 0), dim) )
487 call disp%show("arref = array")
488 arref = array
489 call disp%show("weiref = weight")
490 weiref = weight
491 call disp%show("call setRefined(arref, dim, weiref, skip, rsize)")
492 call setRefined(arref, dim, weiref, skip, rsize)
493 call disp%show("rsize")
494 call disp%show( rsize )
495 call disp%show("weiref(1:rsize)")
496 call disp%show( weiref(1:rsize) )
497 call disp%show("arref(:, 1:rsize)")
498 call disp%show( arref(:, 1:rsize) )
499 call disp%show("weiref = weight")
500 weiref = weight
501 call disp%show("arref = transpose(array)")
502 arref = transpose(array)
503 call disp%show("call setRefined(arref, 3_IK - dim, weiref, skip, rsize)")
504 call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
505 call disp%show("rsize")
506 call disp%show( rsize )
507 call disp%show("weiref(1:rsize)")
508 call disp%show( weiref(1:rsize) )
509 call disp%show("arref(1:rsize, :)")
510 call disp%show( arref(1:rsize, :) )
511 call disp%skip()
512 end do
513 end block
514
515 end block
516
517end program example
pm_arrayChoice::getChoice
Select a single (or multiple) element(s) from the input array of intrinsic type of arbitrary kind ran...
Definition: pm_arrayChoice.F90:137
pm_arrayVerbose::getVerbose
Generate an equally-weighted (verbose or flattened) array of the input weighted array of rank 1 or 2.
Definition: pm_arrayVerbose.F90:142
pm_distBern::isHead
Generate and return a scalar (or a vector of length size or an array of the same shape as the input p...
Definition: pm_distBern.F90:160
pm_distUnif::getUnifRand
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
Definition: pm_distUnif.F90:4150
pm_io::show
This is a generic method of the derived type display_type with pass attribute.
Definition: pm_io.F90:11726
pm_io::skip
This is a generic method of the derived type display_type with pass attribute.
Definition: pm_io.F90:11508
pm_arrayChoice
This module contains procedures and generic interfaces for selecting uniformly-distributed or arbitra...
Definition: pm_arrayChoice.F90:38
pm_arrayVerbose
This module contains procedures and generic interfaces for flattening (duplicating the elements of) a...
Definition: pm_arrayVerbose.F90:38
pm_distBern
This module contains classes and procedures for generating Bernoulli-distributed random numbers.
Definition: pm_distBern.F90:39
pm_distUnif
This module contains classes and procedures for computing various statistical quantities related to t...
Definition: pm_distUnif.F90:274
pm_io
This module contains classes and procedures for input/output (IO) or generic display operations on st...
Definition: pm_io.F90:252
pm_io::disp
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
Definition: pm_io.F90:11393
pm_kind
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
Definition: pm_kind.F90:268
pm_kind::LK
integer, parameter LK
The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
Definition: pm_kind.F90:541
pm_kind::CKS
integer, parameter CKS
The single-precision complex kind in Fortran mode. On most platforms, this is a 32-bit real kind.
Definition: pm_kind.F90:570
pm_kind::IK
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
pm_kind::SK
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
pm_kind::RKS
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
pm_io::display_type
Generate and return an object of type display_type.
Definition: pm_io.F90:10282

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! Refine 1D array.
4!%%%%%%%%%%%%%%%%%
5
6nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
7nsam
8+4
9skip
10+1
11weight = getUnifRand(-1, 9, nsam)
12weight
13+4, +9, -1, +6
14array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
15array
16"HIOF"
17getVerbose(array, weight, sum(weight, mask = weight > 0))
18"HHHHIIIIIIIIIFFFFFF"
19arref = array
20weiref = weight
21call setRefined(arref, weiref, skip, rsize)
22rsize
23+3
24weiref(1:rsize)
25+4, +9, +6
26arref
27"HIFF"
28
29nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
30nsam
31+8
32skip
33+2
34weight = getUnifRand(-1, 9, nsam)
35weight
36+5, +8, +3, +8, +2, +7, +7, +2
37array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
38array
39"QDMGKRMJ"
40getVerbose(array, weight, sum(weight, mask = weight > 0))
41"QQQQQDDDDDDDDMMMGGGGGGGGKKRRRRRRRMMMMMMMJJ"
42arref = array
43weiref = weight
44call setRefined(arref, weiref, skip, rsize)
45rsize
46+8
47weiref(1:rsize)
48+2, +4, +2, +4, +1, +3, +4, +1
49arref
50"QDMGKRMJ"
51
52nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
53nsam
54+7
55skip
56+4
57weight = getUnifRand(-1, 9, nsam)
58weight
59+8, +4, +6, +7, +0, +9, +1
60array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
61array
62"FVPNBNY"
63getVerbose(array, weight, sum(weight, mask = weight > 0))
64"FFFFFFFFVVVVPPPPPPNNNNNNNNNNNNNNNNY"
65arref = array
66weiref = weight
67call setRefined(arref, weiref, skip, rsize)
68rsize
69+5
70weiref(1:rsize)
71+2, +1, +1, +2, +2
72arref
73"FVPNNNY"
74
75nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
76nsam
77+7
78skip
79+3
80weight = getUnifRand(-1, 9, nsam)
81weight
82+9, +8, +0, +6, +6, +5, +9
83array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
84array
85"NVISPNO"
86getVerbose(array, weight, sum(weight, mask = weight > 0))
87"NNNNNNNNNVVVVVVVVSSSSSSPPPPPPNNNNNOOOOOOOOO"
88arref = array
89weiref = weight
90call setRefined(arref, weiref, skip, rsize)
91rsize
92+6
93weiref(1:rsize)
94+3, +2, +2, +2, +2, +3
95arref
96"NVSPNOO"
97
98nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
99nsam
100+5
101skip
102+1
103weight = getUnifRand(-1, 9, nsam)
104weight
105+1, +7, +3, +6, +8
106array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
107array
108"TWZMR"
109getVerbose(array, weight, sum(weight, mask = weight > 0))
110"TWWWWWWWZZZMMMMMMRRRRRRRR"
111arref = array
112weiref = weight
113call setRefined(arref, weiref, skip, rsize)
114rsize
115+5
116weiref(1:rsize)
117+1, +7, +3, +6, +8
118arref
119"TWZMR"
120
121nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
122nsam
123+1
124skip
125+4
126weight = getUnifRand(-1, 9, nsam)
127weight
128+5
129array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
130array
131"W"
132getVerbose(array, weight, sum(weight, mask = weight > 0))
133"WWWWW"
134arref = array
135weiref = weight
136call setRefined(arref, weiref, skip, rsize)
137rsize
138+1
139weiref(1:rsize)
140+1
141arref
142"W"
143
144nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
145nsam
146+2
147skip
148+1
149weight = getUnifRand(-1, 9, nsam)
150weight
151+4, +9
152array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
153array
154"ST"
155getVerbose(array, weight, sum(weight, mask = weight > 0))
156"SSSSTTTTTTTTT"
157arref = array
158weiref = weight
159call setRefined(arref, weiref, skip, rsize)
160rsize
161+2
162weiref(1:rsize)
163+4, +9
164arref
165"ST"
166
167nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
168nsam
169+1
170skip
171+3
172weight = getUnifRand(-1, 9, nsam)
173weight
174-1
175array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
176array
177"N"
178getVerbose(array, weight, sum(weight, mask = weight > 0))
179""
180arref = array
181weiref = weight
182call setRefined(arref, weiref, skip, rsize)
183rsize
184+0
185weiref(1:rsize)
186
187arref
188"N"
189
190nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
191nsam
192+4
193skip
194+3
195weight = getUnifRand(-1, 9, nsam)
196weight
197+3, +9, +8, +3
198array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
199array
200"LJGB"
201getVerbose(array, weight, sum(weight, mask = weight > 0))
202"LLLJJJJJJJJJGGGGGGGGBBB"
203arref = array
204weiref = weight
205call setRefined(arref, weiref, skip, rsize)
206rsize
207+4
208weiref(1:rsize)
209+1, +3, +2, +1
210arref
211"LJGB"
212
213nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
214nsam
215+2
216skip
217+3
218weight = getUnifRand(-1, 9, nsam)
219weight
220+6, +3
221array = getUnifRand(repeat('A', nsam), repeat('Z', nsam))
222array
223"AM"
224getVerbose(array, weight, sum(weight, mask = weight > 0))
225"AAAAAAMMM"
226arref = array
227weiref = weight
228call setRefined(arref, weiref, skip, rsize)
229rsize
230+2
231weiref(1:rsize)
232+2, +1
233arref
234"AM"
235
236nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
237nsam
238+5
239skip
240+3
241weight = getUnifRand(-1, 9, nsam)
242weight
243+5, +6, +0, -1, +6
244array = getUnifRand('AA', 'ZZ', nsam)
245array
246"CM", "ST", "CL", "GS", "TN"
247getVerbose(array, weight, sum(weight, mask = weight > 0))
248"CM", "CM", "CM", "CM", "CM", "ST", "ST", "ST", "ST", "ST", "ST", "TN", "TN", "TN", "TN", "TN", "TN"
249arref = array
250weiref = weight
251call setRefined(arref, weiref, skip, rsize)
252rsize
253+3
254weiref(1:rsize)
255+1, +2, +2
256arref
257"CM", "ST", "TN", "GS", "TN"
258
259nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
260nsam
261+2
262skip
263+4
264weight = getUnifRand(-1, 9, nsam)
265weight
266+4, +0
267array = getUnifRand('AA', 'ZZ', nsam)
268array
269"IY", "VD"
270getVerbose(array, weight, sum(weight, mask = weight > 0))
271"IY", "IY", "IY", "IY"
272arref = array
273weiref = weight
274call setRefined(arref, weiref, skip, rsize)
275rsize
276+1
277weiref(1:rsize)
278+1
279arref
280"IY", "VD"
281
282nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
283nsam
284+2
285skip
286+1
287weight = getUnifRand(-1, 9, nsam)
288weight
289+4, +5
290array = getUnifRand('AA', 'ZZ', nsam)
291array
292"PO", "HD"
293getVerbose(array, weight, sum(weight, mask = weight > 0))
294"PO", "PO", "PO", "PO", "HD", "HD", "HD", "HD", "HD"
295arref = array
296weiref = weight
297call setRefined(arref, weiref, skip, rsize)
298rsize
299+2
300weiref(1:rsize)
301+4, +5
302arref
303"PO", "HD"
304
305nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
306nsam
307+0
308skip
309+4
310weight = getUnifRand(-1, 9, nsam)
311weight
312
313array = getUnifRand('AA', 'ZZ', nsam)
314array
315
316getVerbose(array, weight, sum(weight, mask = weight > 0))
317
318arref = array
319weiref = weight
320call setRefined(arref, weiref, skip, rsize)
321rsize
322+0
323weiref(1:rsize)
324
325arref
326
327
328nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
329nsam
330+2
331skip
332+2
333weight = getUnifRand(-1, 9, nsam)
334weight
335+8, +0
336array = getUnifRand('AA', 'ZZ', nsam)
337array
338"WQ", "FO"
339getVerbose(array, weight, sum(weight, mask = weight > 0))
340"WQ", "WQ", "WQ", "WQ", "WQ", "WQ", "WQ", "WQ"
341arref = array
342weiref = weight
343call setRefined(arref, weiref, skip, rsize)
344rsize
345+1
346weiref(1:rsize)
347+4
348arref
349"WQ", "FO"
350
351nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
352nsam
353+3
354skip
355+4
356weight = getUnifRand(-1, 9, nsam)
357weight
358+8, +8, +6
359array = getUnifRand('AA', 'ZZ', nsam)
360array
361"HI", "SE", "WD"
362getVerbose(array, weight, sum(weight, mask = weight > 0))
363"HI", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "SE", "SE", "SE", "SE", "SE", "SE", "SE", "SE", "WD", "WD", "WD", "WD", "WD", "WD"
364arref = array
365weiref = weight
366call setRefined(arref, weiref, skip, rsize)
367rsize
368+3
369weiref(1:rsize)
370+2, +2, +1
371arref
372"HI", "SE", "WD"
373
374nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
375nsam
376+8
377skip
378+3
379weight = getUnifRand(-1, 9, nsam)
380weight
381+2, +2, +6, +4, +5, +9, +5, +7
382array = getUnifRand('AA', 'ZZ', nsam)
383array
384"XI", "OS", "QH", "RR", "TG", "HI", "AK", "OH"
385getVerbose(array, weight, sum(weight, mask = weight > 0))
386"XI", "XI", "OS", "OS", "QH", "QH", "QH", "QH", "QH", "QH", "RR", "RR", "RR", "RR", "TG", "TG", "TG", "TG", "TG", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "HI", "AK", "AK", "AK", "AK", "AK", "OH", "OH", "OH", "OH", "OH", "OH", "OH"
387arref = array
388weiref = weight
389call setRefined(arref, weiref, skip, rsize)
390rsize
391+7
392weiref(1:rsize)
393+1, +2, +1, +2, +3, +2, +2
394arref
395"OS", "QH", "RR", "TG", "HI", "AK", "OH", "OH"
396
397nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
398nsam
399+8
400skip
401+1
402weight = getUnifRand(-1, 9, nsam)
403weight
404+6, +6, +0, +3, +3, +5, +8, +7
405array = getUnifRand('AA', 'ZZ', nsam)
406array
407"PO", "VM", "UE", "KH", "OY", "RX", "PB", "IZ"
408getVerbose(array, weight, sum(weight, mask = weight > 0))
409"PO", "PO", "PO", "PO", "PO", "PO", "VM", "VM", "VM", "VM", "VM", "VM", "KH", "KH", "KH", "OY", "OY", "OY", "RX", "RX", "RX", "RX", "RX", "PB", "PB", "PB", "PB", "PB", "PB", "PB", "PB", "IZ", "IZ", "IZ", "IZ", "IZ", "IZ", "IZ"
410arref = array
411weiref = weight
412call setRefined(arref, weiref, skip, rsize)
413rsize
414+7
415weiref(1:rsize)
416+6, +6, +3, +3, +5, +8, +7
417arref
418"PO", "VM", "KH", "OY", "RX", "PB", "IZ", "IZ"
419
420nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
421nsam
422+5
423skip
424+3
425weight = getUnifRand(-1, 9, nsam)
426weight
427+8, +1, +2, +4, +0
428array = getUnifRand('AA', 'ZZ', nsam)
429array
430"OD", "JG", "TX", "WD", "CJ"
431getVerbose(array, weight, sum(weight, mask = weight > 0))
432"OD", "OD", "OD", "OD", "OD", "OD", "OD", "OD", "JG", "TX", "TX", "WD", "WD", "WD", "WD"
433arref = array
434weiref = weight
435call setRefined(arref, weiref, skip, rsize)
436rsize
437+3
438weiref(1:rsize)
439+2, +1, +2
440arref
441"OD", "JG", "WD", "WD", "CJ"
442
443nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
444nsam
445+2
446skip
447+1
448weight = getUnifRand(-1, 9, nsam)
449weight
450+8, +9
451array = getUnifRand('AA', 'ZZ', nsam)
452array
453"LC", "QP"
454getVerbose(array, weight, sum(weight, mask = weight > 0))
455"LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "QP", "QP", "QP", "QP", "QP", "QP", "QP", "QP", "QP"
456arref = array
457weiref = weight
458call setRefined(arref, weiref, skip, rsize)
459rsize
460+2
461weiref(1:rsize)
462+8, +9
463arref
464"LC", "QP"
465
466nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
467nsam
468+2
469skip
470+2
471weight = getUnifRand(-1, 9, nsam)
472weight
473+8, +5
474array = getUnifRand(0, 9, nsam)
475array
476+3, +5
477getVerbose(array, weight, sum(weight, mask = weight > 0))
478+3, +3, +3, +3, +3, +3, +3, +3, +5, +5, +5, +5, +5
479arref = array
480weiref = weight
481call setRefined(arref, weiref, skip, rsize)
482rsize
483+2
484weiref(1:rsize)
485+4, +2
486arref(1:rsize)
487+3, +5
488
489nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
490nsam
491+1
492skip
493+1
494weight = getUnifRand(-1, 9, nsam)
495weight
496+1
497array = getUnifRand(0, 9, nsam)
498array
499+6
500getVerbose(array, weight, sum(weight, mask = weight > 0))
501+6
502arref = array
503weiref = weight
504call setRefined(arref, weiref, skip, rsize)
505rsize
506+1
507weiref(1:rsize)
508+1
509arref(1:rsize)
510+6
511
512nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
513nsam
514+2
515skip
516+3
517weight = getUnifRand(-1, 9, nsam)
518weight
519+1, +9
520array = getUnifRand(0, 9, nsam)
521array
522+0, +0
523getVerbose(array, weight, sum(weight, mask = weight > 0))
524+0, +0, +0, +0, +0, +0, +0, +0, +0, +0
525arref = array
526weiref = weight
527call setRefined(arref, weiref, skip, rsize)
528rsize
529+1
530weiref(1:rsize)
531+3
532arref(1:rsize)
533+0
534
535nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
536nsam
537+2
538skip
539+2
540weight = getUnifRand(-1, 9, nsam)
541weight
542+8, +3
543array = getUnifRand(0, 9, nsam)
544array
545+8, +1
546getVerbose(array, weight, sum(weight, mask = weight > 0))
547+8, +8, +8, +8, +8, +8, +8, +8, +1, +1, +1
548arref = array
549weiref = weight
550call setRefined(arref, weiref, skip, rsize)
551rsize
552+2
553weiref(1:rsize)
554+4, +1
555arref(1:rsize)
556+8, +1
557
558nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
559nsam
560+0
561skip
562+4
563weight = getUnifRand(-1, 9, nsam)
564weight
565
566array = getUnifRand(0, 9, nsam)
567array
568
569getVerbose(array, weight, sum(weight, mask = weight > 0))
570
571arref = array
572weiref = weight
573call setRefined(arref, weiref, skip, rsize)
574rsize
575+0
576weiref(1:rsize)
577
578arref(1:rsize)
579
580
581nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
582nsam
583+6
584skip
585+4
586weight = getUnifRand(-1, 9, nsam)
587weight
588+8, +0, +6, +8, +2, +0
589array = getUnifRand(0, 9, nsam)
590array
591+8, +3, +5, +1, +9, +8
592getVerbose(array, weight, sum(weight, mask = weight > 0))
593+8, +8, +8, +8, +8, +8, +8, +8, +5, +5, +5, +5, +5, +5, +1, +1, +1, +1, +1, +1, +1, +1, +9, +9
594arref = array
595weiref = weight
596call setRefined(arref, weiref, skip, rsize)
597rsize
598+4
599weiref(1:rsize)
600+2, +1, +2, +1
601arref(1:rsize)
602+8, +5, +1, +9
603
604nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
605nsam
606+1
607skip
608+3
609weight = getUnifRand(-1, 9, nsam)
610weight
611+0
612array = getUnifRand(0, 9, nsam)
613array
614+2
615getVerbose(array, weight, sum(weight, mask = weight > 0))
616
617arref = array
618weiref = weight
619call setRefined(arref, weiref, skip, rsize)
620rsize
621+0
622weiref(1:rsize)
623
624arref(1:rsize)
625
626
627nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
628nsam
629+1
630skip
631+4
632weight = getUnifRand(-1, 9, nsam)
633weight
634+1
635array = getUnifRand(0, 9, nsam)
636array
637+5
638getVerbose(array, weight, sum(weight, mask = weight > 0))
639+5
640arref = array
641weiref = weight
642call setRefined(arref, weiref, skip, rsize)
643rsize
644+0
645weiref(1:rsize)
646
647arref(1:rsize)
648
649
650nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
651nsam
652+2
653skip
654+3
655weight = getUnifRand(-1, 9, nsam)
656weight
657+4, +2
658array = getUnifRand(0, 9, nsam)
659array
660+2, +2
661getVerbose(array, weight, sum(weight, mask = weight > 0))
662+2, +2, +2, +2, +2, +2
663arref = array
664weiref = weight
665call setRefined(arref, weiref, skip, rsize)
666rsize
667+2
668weiref(1:rsize)
669+1, +1
670arref(1:rsize)
671+2, +2
672
673nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
674nsam
675+6
676skip
677+3
678weight = getUnifRand(-1, 9, nsam)
679weight
680+5, +0, +7, +7, +6, +1
681array = getUnifRand(0, 9, nsam)
682array
683+7, +2, +5, +9, +5, +0
684getVerbose(array, weight, sum(weight, mask = weight > 0))
685+7, +7, +7, +7, +7, +5, +5, +5, +5, +5, +5, +5, +9, +9, +9, +9, +9, +9, +9, +5, +5, +5, +5, +5, +5, +0
686arref = array
687weiref = weight
688call setRefined(arref, weiref, skip, rsize)
689rsize
690+4
691weiref(1:rsize)
692+1, +3, +2, +2
693arref(1:rsize)
694+7, +5, +9, +5
695
696nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
697nsam
698+3
699skip
700+4
701weight = getUnifRand(-1, 9, nsam)
702weight
703+0, +3, +6
704array = getUnifRand(.false., .true., nsam)
705array
706F, T, F
707getVerbose(array, weight, sum(weight, mask = weight > 0))
708T, T, T, F, F, F, F, F, F
709arref = array
710weiref = weight
711call setRefined(arref, weiref, skip, rsize)
712rsize
713+1
714weiref(1:rsize)
715+2
716arref(1:rsize)
717F
718
719nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
720nsam
721+4
722skip
723+2
724weight = getUnifRand(-1, 9, nsam)
725weight
726+8, +2, +8, -1
727array = getUnifRand(.false., .true., nsam)
728array
729F, T, T, T
730getVerbose(array, weight, sum(weight, mask = weight > 0))
731F, F, F, F, F, F, F, F, T, T, T, T, T, T, T, T, T, T
732arref = array
733weiref = weight
734call setRefined(arref, weiref, skip, rsize)
735rsize
736+3
737weiref(1:rsize)
738+4, +1, +4
739arref(1:rsize)
740F, T, T
741
742nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
743nsam
744+4
745skip
746+4
747weight = getUnifRand(-1, 9, nsam)
748weight
749+6, +0, +3, +9
750array = getUnifRand(.false., .true., nsam)
751array
752T, F, F, F
753getVerbose(array, weight, sum(weight, mask = weight > 0))
754T, T, T, T, T, T, F, F, F, F, F, F, F, F, F, F, F, F
755arref = array
756weiref = weight
757call setRefined(arref, weiref, skip, rsize)
758rsize
759+3
760weiref(1:rsize)
761+1, +1, +2
762arref(1:rsize)
763T, F, F
764
765nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
766nsam
767+4
768skip
769+4
770weight = getUnifRand(-1, 9, nsam)
771weight
772+1, +0, +5, +7
773array = getUnifRand(.false., .true., nsam)
774array
775T, F, T, T
776getVerbose(array, weight, sum(weight, mask = weight > 0))
777T, T, T, T, T, T, T, T, T, T, T, T, T
778arref = array
779weiref = weight
780call setRefined(arref, weiref, skip, rsize)
781rsize
782+2
783weiref(1:rsize)
784+1, +2
785arref(1:rsize)
786T, T
787
788nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
789nsam
790+1
791skip
792+3
793weight = getUnifRand(-1, 9, nsam)
794weight
795+6
796array = getUnifRand(.false., .true., nsam)
797array
798T
799getVerbose(array, weight, sum(weight, mask = weight > 0))
800T, T, T, T, T, T
801arref = array
802weiref = weight
803call setRefined(arref, weiref, skip, rsize)
804rsize
805+1
806weiref(1:rsize)
807+2
808arref(1:rsize)
809T
810
811nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
812nsam
813+1
814skip
815+1
816weight = getUnifRand(-1, 9, nsam)
817weight
818+0
819array = getUnifRand(.false., .true., nsam)
820array
821F
822getVerbose(array, weight, sum(weight, mask = weight > 0))
823
824arref = array
825weiref = weight
826call setRefined(arref, weiref, skip, rsize)
827rsize
828+0
829weiref(1:rsize)
830
831arref(1:rsize)
832
833
834nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
835nsam
836+4
837skip
838+3
839weight = getUnifRand(-1, 9, nsam)
840weight
841+3, +3, +4, +1
842array = getUnifRand(.false., .true., nsam)
843array
844F, T, F, F
845getVerbose(array, weight, sum(weight, mask = weight > 0))
846F, F, F, T, T, T, F, F, F, F, F
847arref = array
848weiref = weight
849call setRefined(arref, weiref, skip, rsize)
850rsize
851+3
852weiref(1:rsize)
853+1, +1, +1
854arref(1:rsize)
855F, T, F
856
857nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
858nsam
859+9
860skip
861+4
862weight = getUnifRand(-1, 9, nsam)
863weight
864+5, +8, +4, -1, +5, +2, +0, +0, -1
865array = getUnifRand(.false., .true., nsam)
866array
867T, T, F, T, T, F, F, T, T
868getVerbose(array, weight, sum(weight, mask = weight > 0))
869T, T, T, T, T, T, T, T, T, T, T, T, T, F, F, F, F, T, T, T, T, T, F, F
870arref = array
871weiref = weight
872call setRefined(arref, weiref, skip, rsize)
873rsize
874+5
875weiref(1:rsize)
876+1, +2, +1, +1, +1
877arref(1:rsize)
878T, T, F, T, F
879
880nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
881nsam
882+9
883skip
884+1
885weight = getUnifRand(-1, 9, nsam)
886weight
887+2, +7, +0, +2, +8, +2, +8, +7, +5
888array = getUnifRand(.false., .true., nsam)
889array
890F, T, F, F, T, F, F, F, F
891getVerbose(array, weight, sum(weight, mask = weight > 0))
892F, F, T, T, T, T, T, T, T, F, F, T, T, T, T, T, T, T, T, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F, F
893arref = array
894weiref = weight
895call setRefined(arref, weiref, skip, rsize)
896rsize
897+8
898weiref(1:rsize)
899+2, +7, +2, +8, +2, +8, +7, +5
900arref(1:rsize)
901F, T, F, T, F, F, F, F
902
903nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
904nsam
905+0
906skip
907+1
908weight = getUnifRand(-1, 9, nsam)
909weight
910
911array = getUnifRand(.false., .true., nsam)
912array
913
914getVerbose(array, weight, sum(weight, mask = weight > 0))
915
916arref = array
917weiref = weight
918call setRefined(arref, weiref, skip, rsize)
919rsize
920+0
921weiref(1:rsize)
922
923arref(1:rsize)
924
925
926nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
927nsam
928+8
929skip
930+3
931weight = getUnifRand(-1, 9, nsam)
932weight
933-1, +5, +6, +2, +0, -1, +0, +5
934array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
935array
936(+8.00000000, +9.00000000), (+6.00000000, +0.00000000), (+2.00000000, +7.00000000), (+7.00000000, +8.00000000), (+1.00000000, +4.00000000), (+6.00000000, +1.00000000), (+8.00000000, +5.00000000), (+9.00000000, +9.00000000)
937getVerbose(array, weight, sum(weight, mask = weight > 0))
938(+6.00000000, +0.00000000), (+6.00000000, +0.00000000), (+6.00000000, +0.00000000), (+6.00000000, +0.00000000), (+6.00000000, +0.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+7.00000000, +8.00000000), (+7.00000000, +8.00000000), (+9.00000000, +9.00000000), (+9.00000000, +9.00000000), (+9.00000000, +9.00000000), (+9.00000000, +9.00000000), (+9.00000000, +9.00000000)
939arref = array
940weiref = weight
941call setRefined(arref, weiref, skip, rsize)
942rsize
943+4
944weiref(1:rsize)
945+1, +2, +1, +2
946arref(1:rsize)
947(+6.00000000, +0.00000000), (+2.00000000, +7.00000000), (+7.00000000, +8.00000000), (+9.00000000, +9.00000000)
948
949nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
950nsam
951+5
952skip
953+3
954weight = getUnifRand(-1, 9, nsam)
955weight
956-1, +4, +8, +0, +2
957array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
958array
959(+4.00000000, +9.00000000), (+2.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +7.00000000), (+6.00000000, +0.00000000)
960getVerbose(array, weight, sum(weight, mask = weight > 0))
961(+2.00000000, +2.00000000), (+2.00000000, +2.00000000), (+2.00000000, +2.00000000), (+2.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+6.00000000, +0.00000000), (+6.00000000, +0.00000000)
962arref = array
963weiref = weight
964call setRefined(arref, weiref, skip, rsize)
965rsize
966+2
967weiref(1:rsize)
968+1, +3
969arref(1:rsize)
970(+2.00000000, +2.00000000), (+3.00000000, +2.00000000)
971
972nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
973nsam
974+1
975skip
976+1
977weight = getUnifRand(-1, 9, nsam)
978weight
979-1
980array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
981array
982(+6.00000000, +8.00000000)
983getVerbose(array, weight, sum(weight, mask = weight > 0))
984
985arref = array
986weiref = weight
987call setRefined(arref, weiref, skip, rsize)
988rsize
989+0
990weiref(1:rsize)
991
992arref(1:rsize)
993
994
995nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
996nsam
997+2
998skip
999+3
1000weight = getUnifRand(-1, 9, nsam)
1001weight
1002+0, +7
1003array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1004array
1005(+8.00000000, +7.00000000), (+2.00000000, +0.00000000)
1006getVerbose(array, weight, sum(weight, mask = weight > 0))
1007(+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000)
1008arref = array
1009weiref = weight
1010call setRefined(arref, weiref, skip, rsize)
1011rsize
1012+1
1013weiref(1:rsize)
1014+2
1015arref(1:rsize)
1016(+2.00000000, +0.00000000)
1017
1018nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1019nsam
1020+1
1021skip
1022+2
1023weight = getUnifRand(-1, 9, nsam)
1024weight
1025+6
1026array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1027array
1028(+3.00000000, +2.00000000)
1029getVerbose(array, weight, sum(weight, mask = weight > 0))
1030(+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000)
1031arref = array
1032weiref = weight
1033call setRefined(arref, weiref, skip, rsize)
1034rsize
1035+1
1036weiref(1:rsize)
1037+3
1038arref(1:rsize)
1039(+3.00000000, +2.00000000)
1040
1041nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1042nsam
1043+4
1044skip
1045+4
1046weight = getUnifRand(-1, 9, nsam)
1047weight
1048+2, +0, +9, +2
1049array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1050array
1051(+3.00000000, +6.00000000), (+9.00000000, +4.00000000), (+0.00000000, +5.00000000), (+4.00000000, +5.00000000)
1052getVerbose(array, weight, sum(weight, mask = weight > 0))
1053(+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+4.00000000, +5.00000000), (+4.00000000, +5.00000000)
1054arref = array
1055weiref = weight
1056call setRefined(arref, weiref, skip, rsize)
1057rsize
1058+2
1059weiref(1:rsize)
1060+2, +1
1061arref(1:rsize)
1062(+0.00000000, +5.00000000), (+4.00000000, +5.00000000)
1063
1064nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1065nsam
1066+7
1067skip
1068+1
1069weight = getUnifRand(-1, 9, nsam)
1070weight
1071+4, +7, +1, +0, +8, +5, +8
1072array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1073array
1074(+3.00000000, +6.00000000), (+8.00000000, +7.00000000), (+0.00000000, +1.00000000), (+5.00000000, +5.00000000), (+8.00000000, +4.00000000), (+1.00000000, +4.00000000), (+0.00000000, +8.00000000)
1075getVerbose(array, weight, sum(weight, mask = weight > 0))
1076(+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+0.00000000, +1.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000), (+0.00000000, +8.00000000)
1077arref = array
1078weiref = weight
1079call setRefined(arref, weiref, skip, rsize)
1080rsize
1081+6
1082weiref(1:rsize)
1083+4, +7, +1, +8, +5, +8
1084arref(1:rsize)
1085(+3.00000000, +6.00000000), (+8.00000000, +7.00000000), (+0.00000000, +1.00000000), (+8.00000000, +4.00000000), (+1.00000000, +4.00000000), (+0.00000000, +8.00000000)
1086
1087nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1088nsam
1089+9
1090skip
1091+1
1092weight = getUnifRand(-1, 9, nsam)
1093weight
1094-1, +7, +3, +3, +3, +4, +8, +7, +5
1095array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1096array
1097(+9.00000000, +6.00000000), (+5.00000000, +0.00000000), (+4.00000000, +8.00000000), (+5.00000000, +1.00000000), (+1.00000000, +9.00000000), (+9.00000000, +1.00000000), (+1.00000000, +5.00000000), (+2.00000000, +0.00000000), (+6.00000000, +3.00000000)
1098getVerbose(array, weight, sum(weight, mask = weight > 0))
1099(+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+4.00000000, +8.00000000), (+4.00000000, +8.00000000), (+4.00000000, +8.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+1.00000000, +9.00000000), (+1.00000000, +9.00000000), (+1.00000000, +9.00000000), (+9.00000000, +1.00000000), (+9.00000000, +1.00000000), (+9.00000000, +1.00000000), (+9.00000000, +1.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000)
1100arref = array
1101weiref = weight
1102call setRefined(arref, weiref, skip, rsize)
1103rsize
1104+8
1105weiref(1:rsize)
1106+7, +3, +3, +3, +4, +8, +7, +5
1107arref(1:rsize)
1108(+5.00000000, +0.00000000), (+4.00000000, +8.00000000), (+5.00000000, +1.00000000), (+1.00000000, +9.00000000), (+9.00000000, +1.00000000), (+1.00000000, +5.00000000), (+2.00000000, +0.00000000), (+6.00000000, +3.00000000)
1109
1110nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1111nsam
1112+1
1113skip
1114+3
1115weight = getUnifRand(-1, 9, nsam)
1116weight
1117+9
1118array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1119array
1120(+1.00000000, +2.00000000)
1121getVerbose(array, weight, sum(weight, mask = weight > 0))
1122(+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000)
1123arref = array
1124weiref = weight
1125call setRefined(arref, weiref, skip, rsize)
1126rsize
1127+1
1128weiref(1:rsize)
1129+3
1130arref(1:rsize)
1131(+1.00000000, +2.00000000)
1132
1133nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1134nsam
1135+8
1136skip
1137+3
1138weight = getUnifRand(-1, 9, nsam)
1139weight
1140+7, +5, +1, -1, +6, +9, +7, +3
1141array = cmplx(getUnifRand(0, 9, nsam), getUnifRand(0, 9, nsam), TKG)
1142array
1143(+1.00000000, +5.00000000), (+8.00000000, +0.00000000), (+0.00000000, +3.00000000), (+4.00000000, +8.00000000), (+5.00000000, +6.00000000), (+5.00000000, +1.00000000), (+9.00000000, +6.00000000), (+8.00000000, +4.00000000)
1144getVerbose(array, weight, sum(weight, mask = weight > 0))
1145(+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+0.00000000, +3.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+5.00000000, +1.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000), (+8.00000000, +4.00000000)
1146arref = array
1147weiref = weight
1148call setRefined(arref, weiref, skip, rsize)
1149rsize
1150+6
1151weiref(1:rsize)
1152+2, +2, +2, +3, +2, +1
1153arref(1:rsize)
1154(+1.00000000, +5.00000000), (+8.00000000, +0.00000000), (+5.00000000, +6.00000000), (+5.00000000, +1.00000000), (+9.00000000, +6.00000000), (+8.00000000, +4.00000000)
1155
1156nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1157nsam
1158+9
1159skip
1160+1
1161weight = getUnifRand(-1, 9, nsam)
1162weight
1163-1, +5, -1, +4, +7, +5, +8, +6, +4
1164array = getUnifRand(0, 9, nsam)
1165array
1166+7.00000000, +4.00000000, +8.00000000, +7.00000000, +4.00000000, +7.00000000, +4.00000000, +5.00000000, +2.00000000
1167getVerbose(array, weight, sum(weight, mask = weight > 0))
1168+4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000
1169arref = array
1170weiref = weight
1171call setRefined(arref, weiref, skip, rsize)
1172rsize
1173+7
1174weiref(1:rsize)
1175+5, +4, +7, +5, +8, +6, +4
1176arref(1:rsize)
1177+4.00000000, +7.00000000, +4.00000000, +7.00000000, +4.00000000, +5.00000000, +2.00000000
1178
1179nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1180nsam
1181+8
1182skip
1183+2
1184weight = getUnifRand(-1, 9, nsam)
1185weight
1186+9, +0, +1, +8, +8, +8, +5, +6
1187array = getUnifRand(0, 9, nsam)
1188array
1189+1.00000000, +2.00000000, +7.00000000, +8.00000000, +5.00000000, +5.00000000, +3.00000000, +2.00000000
1190getVerbose(array, weight, sum(weight, mask = weight > 0))
1191+1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +7.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000
1192arref = array
1193weiref = weight
1194call setRefined(arref, weiref, skip, rsize)
1195rsize
1196+7
1197weiref(1:rsize)
1198+4, +1, +4, +4, +4, +2, +3
1199arref(1:rsize)
1200+1.00000000, +7.00000000, +8.00000000, +5.00000000, +5.00000000, +3.00000000, +2.00000000
1201
1202nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1203nsam
1204+0
1205skip
1206+2
1207weight = getUnifRand(-1, 9, nsam)
1208weight
1209
1210array = getUnifRand(0, 9, nsam)
1211array
1212
1213getVerbose(array, weight, sum(weight, mask = weight > 0))
1214
1215arref = array
1216weiref = weight
1217call setRefined(arref, weiref, skip, rsize)
1218rsize
1219+0
1220weiref(1:rsize)
1221
1222arref(1:rsize)
1223
1224
1225nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1226nsam
1227+6
1228skip
1229+1
1230weight = getUnifRand(-1, 9, nsam)
1231weight
1232+5, +7, -1, +4, +4, +0
1233array = getUnifRand(0, 9, nsam)
1234array
1235+4.00000000, +4.00000000, +1.00000000, +7.00000000, +8.00000000, +7.00000000
1236getVerbose(array, weight, sum(weight, mask = weight > 0))
1237+4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000
1238arref = array
1239weiref = weight
1240call setRefined(arref, weiref, skip, rsize)
1241rsize
1242+4
1243weiref(1:rsize)
1244+5, +7, +4, +4
1245arref(1:rsize)
1246+4.00000000, +4.00000000, +7.00000000, +8.00000000
1247
1248nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1249nsam
1250+6
1251skip
1252+1
1253weight = getUnifRand(-1, 9, nsam)
1254weight
1255+2, +3, +9, +7, +1, +9
1256array = getUnifRand(0, 9, nsam)
1257array
1258+2.00000000, +2.00000000, +1.00000000, +5.00000000, +0.00000000, +6.00000000
1259getVerbose(array, weight, sum(weight, mask = weight > 0))
1260+2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +0.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000
1261arref = array
1262weiref = weight
1263call setRefined(arref, weiref, skip, rsize)
1264rsize
1265+6
1266weiref(1:rsize)
1267+2, +3, +9, +7, +1, +9
1268arref(1:rsize)
1269+2.00000000, +2.00000000, +1.00000000, +5.00000000, +0.00000000, +6.00000000
1270
1271nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1272nsam
1273+4
1274skip
1275+1
1276weight = getUnifRand(-1, 9, nsam)
1277weight
1278+1, +9, +4, +5
1279array = getUnifRand(0, 9, nsam)
1280array
1281+0.00000000, +8.00000000, +4.00000000, +4.00000000
1282getVerbose(array, weight, sum(weight, mask = weight > 0))
1283+0.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000
1284arref = array
1285weiref = weight
1286call setRefined(arref, weiref, skip, rsize)
1287rsize
1288+4
1289weiref(1:rsize)
1290+1, +9, +4, +5
1291arref(1:rsize)
1292+0.00000000, +8.00000000, +4.00000000, +4.00000000
1293
1294nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1295nsam
1296+9
1297skip
1298+4
1299weight = getUnifRand(-1, 9, nsam)
1300weight
1301+9, +2, +8, +8, +4, +7, +8, +9, +3
1302array = getUnifRand(0, 9, nsam)
1303array
1304+8.00000000, +8.00000000, +2.00000000, +7.00000000, +7.00000000, +0.00000000, +3.00000000, +3.00000000, +5.00000000
1305getVerbose(array, weight, sum(weight, mask = weight > 0))
1306+8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +0.00000000, +0.00000000, +0.00000000, +0.00000000, +0.00000000, +0.00000000, +0.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +5.00000000, +5.00000000, +5.00000000
1307arref = array
1308weiref = weight
1309call setRefined(arref, weiref, skip, rsize)
1310rsize
1311+8
1312weiref(1:rsize)
1313+2, +2, +2, +1, +2, +2, +2, +1
1314arref(1:rsize)
1315+8.00000000, +2.00000000, +7.00000000, +7.00000000, +0.00000000, +3.00000000, +3.00000000, +5.00000000
1316
1317nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1318nsam
1319+4
1320skip
1321+4
1322weight = getUnifRand(-1, 9, nsam)
1323weight
1324+4, +2, +0, +6
1325array = getUnifRand(0, 9, nsam)
1326array
1327+8.00000000, +2.00000000, +5.00000000, +5.00000000
1328getVerbose(array, weight, sum(weight, mask = weight > 0))
1329+8.00000000, +8.00000000, +8.00000000, +8.00000000, +2.00000000, +2.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000
1330arref = array
1331weiref = weight
1332call setRefined(arref, weiref, skip, rsize)
1333rsize
1334+2
1335weiref(1:rsize)
1336+1, +2
1337arref(1:rsize)
1338+8.00000000, +5.00000000
1339
1340nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1341nsam
1342+0
1343skip
1344+1
1345weight = getUnifRand(-1, 9, nsam)
1346weight
1347
1348array = getUnifRand(0, 9, nsam)
1349array
1350
1351getVerbose(array, weight, sum(weight, mask = weight > 0))
1352
1353arref = array
1354weiref = weight
1355call setRefined(arref, weiref, skip, rsize)
1356rsize
1357+0
1358weiref(1:rsize)
1359
1360arref(1:rsize)
1361
1362
1363nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1364nsam
1365+5
1366skip
1367+4
1368weight = getUnifRand(-1, 9, nsam)
1369weight
1370+9, +6, +9, +1, -1
1371array = getUnifRand(0, 9, nsam)
1372array
1373+6.00000000, +3.00000000, +1.00000000, +6.00000000, +9.00000000
1374getVerbose(array, weight, sum(weight, mask = weight > 0))
1375+6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +6.00000000
1376arref = array
1377weiref = weight
1378call setRefined(arref, weiref, skip, rsize)
1379rsize
1380+3
1381weiref(1:rsize)
1382+2, +1, +3
1383arref(1:rsize)
1384+6.00000000, +3.00000000, +1.00000000
1385
1386
1387!%%%%%%%%%%%%%%%%%
1388! Refine 2D array.
1389!%%%%%%%%%%%%%%%%%
1390
1391dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1392[dim, ndim, nsam]
1393+2, +3, +6
1394skip
1395+2
1396weight = getUnifRand(-1, 9, nsam)
1397weight
1398+1, +3, +1, +2, +9, +9
1399array = getUnifRand('AA', 'ZZ', ndim, nsam)
1400array
1401"DO", "HA", "KY", "TD", "KX", "PT"
1402"CX", "TL", "MZ", "JF", "FT", "QP"
1403"AH", "XS", "DN", "LA", "ZN", "XV"
1404getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1405"DO", "HA", "HA", "HA", "KY", "TD", "TD", "KX", "KX", "KX", "KX", "KX", "KX", "KX", "KX", "KX", "PT", "PT", "PT", "PT", "PT", "PT", "PT", "PT", "PT"
1406"CX", "TL", "TL", "TL", "MZ", "JF", "JF", "FT", "FT", "FT", "FT", "FT", "FT", "FT", "FT", "FT", "QP", "QP", "QP", "QP", "QP", "QP", "QP", "QP", "QP"
1407"AH", "XS", "XS", "XS", "DN", "LA", "LA", "ZN", "ZN", "ZN", "ZN", "ZN", "ZN", "ZN", "ZN", "ZN", "XV", "XV", "XV", "XV", "XV", "XV", "XV", "XV", "XV"
1408arref = array
1409weiref = weight
1410call setRefined(arref, dim, weiref, skip, rsize)
1411rsize
1412+4
1413weiref(1:rsize)
1414+2, +1, +5, +4
1415arref(:, 1:rsize)
1416"HA", "TD", "KX", "PT"
1417"TL", "JF", "FT", "QP"
1418"XS", "LA", "ZN", "XV"
1419weiref = weight
1420arref = transpose(array)
1421call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1422rsize
1423+4
1424weiref(1:rsize)
1425+2, +1, +5, +4
1426arref(1:rsize, :)
1427"HA", "TL", "XS"
1428"TD", "JF", "LA"
1429"KX", "FT", "ZN"
1430"PT", "QP", "XV"
1431
1432dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1433[dim, ndim, nsam]
1434+2, +2, +1
1435skip
1436+2
1437weight = getUnifRand(-1, 9, nsam)
1438weight
1439-1
1440array = getUnifRand('AA', 'ZZ', ndim, nsam)
1441array
1442"ZC"
1443"OZ"
1444getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1445
1446
1447arref = array
1448weiref = weight
1449call setRefined(arref, dim, weiref, skip, rsize)
1450rsize
1451+0
1452weiref(1:rsize)
1453
1454arref(:, 1:rsize)
1455
1456
1457weiref = weight
1458arref = transpose(array)
1459call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1460rsize
1461+0
1462weiref(1:rsize)
1463
1464arref(1:rsize, :)
1465
1466dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1467[dim, ndim, nsam]
1468+2, +1, +1
1469skip
1470+2
1471weight = getUnifRand(-1, 9, nsam)
1472weight
1473+2
1474array = getUnifRand('AA', 'ZZ', ndim, nsam)
1475array
1476"US"
1477getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1478"US", "US"
1479arref = array
1480weiref = weight
1481call setRefined(arref, dim, weiref, skip, rsize)
1482rsize
1483+1
1484weiref(1:rsize)
1485+1
1486arref(:, 1:rsize)
1487"US"
1488weiref = weight
1489arref = transpose(array)
1490call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1491rsize
1492+1
1493weiref(1:rsize)
1494+1
1495arref(1:rsize, :)
1496"US"
1497
1498dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1499[dim, ndim, nsam]
1500+2, +1, +3
1501skip
1502+2
1503weight = getUnifRand(-1, 9, nsam)
1504weight
1505+2, +3, -1
1506array = getUnifRand('AA', 'ZZ', ndim, nsam)
1507array
1508"LG", "HV", "FY"
1509getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1510"LG", "LG", "HV", "HV", "HV"
1511arref = array
1512weiref = weight
1513call setRefined(arref, dim, weiref, skip, rsize)
1514rsize
1515+2
1516weiref(1:rsize)
1517+1, +1
1518arref(:, 1:rsize)
1519"LG", "HV"
1520weiref = weight
1521arref = transpose(array)
1522call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1523rsize
1524+2
1525weiref(1:rsize)
1526+1, +1
1527arref(1:rsize, :)
1528"LG"
1529"HV"
1530
1531dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1532[dim, ndim, nsam]
1533+2, +3, +4
1534skip
1535+1
1536weight = getUnifRand(-1, 9, nsam)
1537weight
1538+2, +8, +9, +8
1539array = getUnifRand('AA', 'ZZ', ndim, nsam)
1540array
1541"XD", "EP", "TL", "HG"
1542"QP", "HC", "EG", "YE"
1543"GZ", "SQ", "WO", "PW"
1544getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1545"XD", "XD", "EP", "EP", "EP", "EP", "EP", "EP", "EP", "EP", "TL", "TL", "TL", "TL", "TL", "TL", "TL", "TL", "TL", "HG", "HG", "HG", "HG", "HG", "HG", "HG", "HG"
1546"QP", "QP", "HC", "HC", "HC", "HC", "HC", "HC", "HC", "HC", "EG", "EG", "EG", "EG", "EG", "EG", "EG", "EG", "EG", "YE", "YE", "YE", "YE", "YE", "YE", "YE", "YE"
1547"GZ", "GZ", "SQ", "SQ", "SQ", "SQ", "SQ", "SQ", "SQ", "SQ", "WO", "WO", "WO", "WO", "WO", "WO", "WO", "WO", "WO", "PW", "PW", "PW", "PW", "PW", "PW", "PW", "PW"
1548arref = array
1549weiref = weight
1550call setRefined(arref, dim, weiref, skip, rsize)
1551rsize
1552+4
1553weiref(1:rsize)
1554+2, +8, +9, +8
1555arref(:, 1:rsize)
1556"XD", "EP", "TL", "HG"
1557"QP", "HC", "EG", "YE"
1558"GZ", "SQ", "WO", "PW"
1559weiref = weight
1560arref = transpose(array)
1561call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1562rsize
1563+4
1564weiref(1:rsize)
1565+2, +8, +9, +8
1566arref(1:rsize, :)
1567"XD", "QP", "GZ"
1568"EP", "HC", "SQ"
1569"TL", "EG", "WO"
1570"HG", "YE", "PW"
1571
1572dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1573[dim, ndim, nsam]
1574+2, +1, +5
1575skip
1576+2
1577weight = getUnifRand(-1, 9, nsam)
1578weight
1579+9, +7, +6, +7, +2
1580array = getUnifRand('AA', 'ZZ', ndim, nsam)
1581array
1582"CO", "KT", "LZ", "VO", "TR"
1583getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1584"CO", "CO", "CO", "CO", "CO", "CO", "CO", "CO", "CO", "KT", "KT", "KT", "KT", "KT", "KT", "KT", "LZ", "LZ", "LZ", "LZ", "LZ", "LZ", "VO", "VO", "VO", "VO", "VO", "VO", "VO", "TR", "TR"
1585arref = array
1586weiref = weight
1587call setRefined(arref, dim, weiref, skip, rsize)
1588rsize
1589+5
1590weiref(1:rsize)
1591+4, +4, +3, +3, +1
1592arref(:, 1:rsize)
1593"CO", "KT", "LZ", "VO", "TR"
1594weiref = weight
1595arref = transpose(array)
1596call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1597rsize
1598+5
1599weiref(1:rsize)
1600+4, +4, +3, +3, +1
1601arref(1:rsize, :)
1602"CO"
1603"KT"
1604"LZ"
1605"VO"
1606"TR"
1607
1608dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1609[dim, ndim, nsam]
1610+2, +2, +8
1611skip
1612+1
1613weight = getUnifRand(-1, 9, nsam)
1614weight
1615+1, +1, +3, +7, +3, +0, +8, -1
1616array = getUnifRand('AA', 'ZZ', ndim, nsam)
1617array
1618"TR", "WK", "CG", "YX", "SD", "OO", "SS", "NR"
1619"ZH", "NA", "QZ", "XY", "AF", "TL", "HO", "YA"
1620getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1621"TR", "WK", "CG", "CG", "CG", "YX", "YX", "YX", "YX", "YX", "YX", "YX", "SD", "SD", "SD", "SS", "SS", "SS", "SS", "SS", "SS", "SS", "SS"
1622"ZH", "NA", "QZ", "QZ", "QZ", "XY", "XY", "XY", "XY", "XY", "XY", "XY", "AF", "AF", "AF", "HO", "HO", "HO", "HO", "HO", "HO", "HO", "HO"
1623arref = array
1624weiref = weight
1625call setRefined(arref, dim, weiref, skip, rsize)
1626rsize
1627+6
1628weiref(1:rsize)
1629+1, +1, +3, +7, +3, +8
1630arref(:, 1:rsize)
1631"TR", "WK", "CG", "YX", "SD", "SS"
1632"ZH", "NA", "QZ", "XY", "AF", "HO"
1633weiref = weight
1634arref = transpose(array)
1635call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1636rsize
1637+6
1638weiref(1:rsize)
1639+1, +1, +3, +7, +3, +8
1640arref(1:rsize, :)
1641"TR", "ZH"
1642"WK", "NA"
1643"CG", "QZ"
1644"YX", "XY"
1645"SD", "AF"
1646"SS", "HO"
1647
1648dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1649[dim, ndim, nsam]
1650+2, +3, +3
1651skip
1652+1
1653weight = getUnifRand(-1, 9, nsam)
1654weight
1655+4, +0, +8
1656array = getUnifRand('AA', 'ZZ', ndim, nsam)
1657array
1658"UD", "WN", "LU"
1659"YS", "IL", "GH"
1660"WB", "UA", "OT"
1661getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1662"UD", "UD", "UD", "UD", "LU", "LU", "LU", "LU", "LU", "LU", "LU", "LU"
1663"YS", "YS", "YS", "YS", "GH", "GH", "GH", "GH", "GH", "GH", "GH", "GH"
1664"WB", "WB", "WB", "WB", "OT", "OT", "OT", "OT", "OT", "OT", "OT", "OT"
1665arref = array
1666weiref = weight
1667call setRefined(arref, dim, weiref, skip, rsize)
1668rsize
1669+2
1670weiref(1:rsize)
1671+4, +8
1672arref(:, 1:rsize)
1673"UD", "LU"
1674"YS", "GH"
1675"WB", "OT"
1676weiref = weight
1677arref = transpose(array)
1678call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1679rsize
1680+2
1681weiref(1:rsize)
1682+4, +8
1683arref(1:rsize, :)
1684"UD", "YS", "WB"
1685"LU", "GH", "OT"
1686
1687dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1688[dim, ndim, nsam]
1689+2, +2, +0
1690skip
1691+4
1692weight = getUnifRand(-1, 9, nsam)
1693weight
1694
1695array = getUnifRand('AA', 'ZZ', ndim, nsam)
1696array
1697
1698
1699getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1700
1701
1702arref = array
1703weiref = weight
1704call setRefined(arref, dim, weiref, skip, rsize)
1705rsize
1706+0
1707weiref(1:rsize)
1708
1709arref(:, 1:rsize)
1710
1711
1712weiref = weight
1713arref = transpose(array)
1714call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1715rsize
1716+0
1717weiref(1:rsize)
1718
1719arref(1:rsize, :)
1720
1721dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1722[dim, ndim, nsam]
1723+2, +1, +5
1724skip
1725+1
1726weight = getUnifRand(-1, 9, nsam)
1727weight
1728+6, +8, -1, -1, +7
1729array = getUnifRand('AA', 'ZZ', ndim, nsam)
1730array
1731"XH", "KY", "AZ", "JF", "PD"
1732getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1733"XH", "XH", "XH", "XH", "XH", "XH", "KY", "KY", "KY", "KY", "KY", "KY", "KY", "KY", "PD", "PD", "PD", "PD", "PD", "PD", "PD"
1734arref = array
1735weiref = weight
1736call setRefined(arref, dim, weiref, skip, rsize)
1737rsize
1738+3
1739weiref(1:rsize)
1740+6, +8, +7
1741arref(:, 1:rsize)
1742"XH", "KY", "PD"
1743weiref = weight
1744arref = transpose(array)
1745call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1746rsize
1747+3
1748weiref(1:rsize)
1749+6, +8, +7
1750arref(1:rsize, :)
1751"XH"
1752"KY"
1753"PD"
1754
1755dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1756[dim, ndim, nsam]
1757+2, +0, +1
1758skip
1759+4
1760weight = getUnifRand(-1, 9, nsam)
1761weight
1762+4
1763array = getUnifRand(0, 9, ndim, nsam)
1764array
1765getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1766arref = array
1767weiref = weight
1768call setRefined(arref, dim, weiref, skip, rsize)
1769rsize
1770+0
1771weiref(1:rsize)
1772
1773arref(:, 1:rsize)
1774weiref = weight
1775arref = transpose(array)
1776call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1777rsize
1778+1
1779weiref(1:rsize)
1780+1
1781arref(1:rsize, :)
1782
1783
1784dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1785[dim, ndim, nsam]
1786+2, +2, +8
1787skip
1788+4
1789weight = getUnifRand(-1, 9, nsam)
1790weight
1791+9, +4, +8, +9, +6, +8, +2, +0
1792array = getUnifRand(0, 9, ndim, nsam)
1793array
1794+6, +4, +3, +0, +5, +8, +7, +0
1795+8, +6, +9, +9, +8, +0, +3, +1
1796getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1797+6, +6, +6, +6, +6, +6, +6, +6, +6, +4, +4, +4, +4, +3, +3, +3, +3, +3, +3, +3, +3, +0, +0, +0, +0, +0, +0, +0, +0, +0, +5, +5, +5, +5, +5, +5, +8, +8, +8, +8, +8, +8, +8, +8, +7, +7
1798+8, +8, +8, +8, +8, +8, +8, +8, +8, +6, +6, +6, +6, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +9, +8, +8, +8, +8, +8, +8, +0, +0, +0, +0, +0, +0, +0, +0, +3, +3
1799arref = array
1800weiref = weight
1801call setRefined(arref, dim, weiref, skip, rsize)
1802rsize
1803+6
1804weiref(1:rsize)
1805+2, +1, +2, +2, +2, +2
1806arref(:, 1:rsize)
1807+6, +4, +3, +0, +5, +8
1808+8, +6, +9, +9, +8, +0
1809weiref = weight
1810arref = transpose(array)
1811call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1812rsize
1813+6
1814weiref(1:rsize)
1815+2, +1, +2, +2, +2, +2
1816arref(1:rsize, :)
1817+6, +8
1818+4, +6
1819+3, +9
1820+0, +9
1821+5, +8
1822+8, +0
1823
1824dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1825[dim, ndim, nsam]
1826+2, +1, +8
1827skip
1828+4
1829weight = getUnifRand(-1, 9, nsam)
1830weight
1831+8, +8, +9, +2, +2, +0, +2, +3
1832array = getUnifRand(0, 9, ndim, nsam)
1833array
1834+9, +5, +9, +4, +1, +8, +5, +1
1835getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1836+9, +9, +9, +9, +9, +9, +9, +9, +5, +5, +5, +5, +5, +5, +5, +5, +9, +9, +9, +9, +9, +9, +9, +9, +9, +4, +4, +1, +1, +5, +5, +1, +1, +1
1837arref = array
1838weiref = weight
1839call setRefined(arref, dim, weiref, skip, rsize)
1840rsize
1841+5
1842weiref(1:rsize)
1843+2, +2, +2, +1, +1
1844arref(:, 1:rsize)
1845+9, +5, +9, +1, +1
1846weiref = weight
1847arref = transpose(array)
1848call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1849rsize
1850+5
1851weiref(1:rsize)
1852+2, +2, +2, +1, +1
1853arref(1:rsize, :)
1854+9
1855+5
1856+9
1857+1
1858+1
1859
1860dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1861[dim, ndim, nsam]
1862+2, +1, +4
1863skip
1864+2
1865weight = getUnifRand(-1, 9, nsam)
1866weight
1867-1, +6, +0, +6
1868array = getUnifRand(0, 9, ndim, nsam)
1869array
1870+8, +9, +2, +2
1871getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1872+9, +9, +9, +9, +9, +9, +2, +2, +2, +2, +2, +2
1873arref = array
1874weiref = weight
1875call setRefined(arref, dim, weiref, skip, rsize)
1876rsize
1877+2
1878weiref(1:rsize)
1879+3, +3
1880arref(:, 1:rsize)
1881+9, +2
1882weiref = weight
1883arref = transpose(array)
1884call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1885rsize
1886+2
1887weiref(1:rsize)
1888+3, +3
1889arref(1:rsize, :)
1890+9
1891+2
1892
1893dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1894[dim, ndim, nsam]
1895+2, +3, +0
1896skip
1897+3
1898weight = getUnifRand(-1, 9, nsam)
1899weight
1900
1901array = getUnifRand(0, 9, ndim, nsam)
1902array
1903
1904
1905
1906getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1907
1908
1909
1910arref = array
1911weiref = weight
1912call setRefined(arref, dim, weiref, skip, rsize)
1913rsize
1914+0
1915weiref(1:rsize)
1916
1917arref(:, 1:rsize)
1918
1919
1920
1921weiref = weight
1922arref = transpose(array)
1923call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1924rsize
1925+0
1926weiref(1:rsize)
1927
1928arref(1:rsize, :)
1929
1930dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1931[dim, ndim, nsam]
1932+2, +1, +3
1933skip
1934+2
1935weight = getUnifRand(-1, 9, nsam)
1936weight
1937+1, +3, +8
1938array = getUnifRand(0, 9, ndim, nsam)
1939array
1940+5, +6, +7
1941getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1942+5, +6, +6, +6, +7, +7, +7, +7, +7, +7, +7, +7
1943arref = array
1944weiref = weight
1945call setRefined(arref, dim, weiref, skip, rsize)
1946rsize
1947+2
1948weiref(1:rsize)
1949+2, +4
1950arref(:, 1:rsize)
1951+6, +7
1952weiref = weight
1953arref = transpose(array)
1954call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1955rsize
1956+2
1957weiref(1:rsize)
1958+2, +4
1959arref(1:rsize, :)
1960+6
1961+7
1962
1963dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
1964[dim, ndim, nsam]
1965+2, +3, +3
1966skip
1967+1
1968weight = getUnifRand(-1, 9, nsam)
1969weight
1970+5, +8, +7
1971array = getUnifRand(0, 9, ndim, nsam)
1972array
1973+6, +7, +8
1974+0, +3, +4
1975+9, +8, +2
1976getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
1977+6, +6, +6, +6, +6, +7, +7, +7, +7, +7, +7, +7, +7, +8, +8, +8, +8, +8, +8, +8
1978+0, +0, +0, +0, +0, +3, +3, +3, +3, +3, +3, +3, +3, +4, +4, +4, +4, +4, +4, +4
1979+9, +9, +9, +9, +9, +8, +8, +8, +8, +8, +8, +8, +8, +2, +2, +2, +2, +2, +2, +2
1980arref = array
1981weiref = weight
1982call setRefined(arref, dim, weiref, skip, rsize)
1983rsize
1984+3
1985weiref(1:rsize)
1986+5, +8, +7
1987arref(:, 1:rsize)
1988+6, +7, +8
1989+0, +3, +4
1990+9, +8, +2
1991weiref = weight
1992arref = transpose(array)
1993call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
1994rsize
1995+3
1996weiref(1:rsize)
1997+5, +8, +7
1998arref(1:rsize, :)
1999+6, +0, +9
2000+7, +3, +8
2001+8, +4, +2
2002
2003dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2004[dim, ndim, nsam]
2005+2, +1, +4
2006skip
2007+3
2008weight = getUnifRand(-1, 9, nsam)
2009weight
2010+1, +0, -1, +1
2011array = getUnifRand(0, 9, ndim, nsam)
2012array
2013+6, +4, +0, +3
2014getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2015+6, +3
2016arref = array
2017weiref = weight
2018call setRefined(arref, dim, weiref, skip, rsize)
2019rsize
2020+0
2021weiref(1:rsize)
2022
2023arref(:, 1:rsize)
2024
2025weiref = weight
2026arref = transpose(array)
2027call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2028rsize
2029+0
2030weiref(1:rsize)
2031
2032arref(1:rsize, :)
2033
2034dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2035[dim, ndim, nsam]
2036+2, +1, +4
2037skip
2038+2
2039weight = getUnifRand(-1, 9, nsam)
2040weight
2041-1, +4, -1, +3
2042array = getUnifRand(0, 9, ndim, nsam)
2043array
2044+6, +9, +8, +5
2045getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2046+9, +9, +9, +9, +5, +5, +5
2047arref = array
2048weiref = weight
2049call setRefined(arref, dim, weiref, skip, rsize)
2050rsize
2051+2
2052weiref(1:rsize)
2053+2, +1
2054arref(:, 1:rsize)
2055+9, +5
2056weiref = weight
2057arref = transpose(array)
2058call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2059rsize
2060+2
2061weiref(1:rsize)
2062+2, +1
2063arref(1:rsize, :)
2064+9
2065+5
2066
2067dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2068[dim, ndim, nsam]
2069+2, +2, +4
2070skip
2071+4
2072weight = getUnifRand(-1, 9, nsam)
2073weight
2074+6, +8, +3, +9
2075array = getUnifRand(0, 9, ndim, nsam)
2076array
2077+1, +2, +3, +9
2078+8, +4, +6, +7
2079getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2080+1, +1, +1, +1, +1, +1, +2, +2, +2, +2, +2, +2, +2, +2, +3, +3, +3, +9, +9, +9, +9, +9, +9, +9, +9, +9
2081+8, +8, +8, +8, +8, +8, +4, +4, +4, +4, +4, +4, +4, +4, +6, +6, +6, +7, +7, +7, +7, +7, +7, +7, +7, +7
2082arref = array
2083weiref = weight
2084call setRefined(arref, dim, weiref, skip, rsize)
2085rsize
2086+4
2087weiref(1:rsize)
2088+1, +2, +1, +2
2089arref(:, 1:rsize)
2090+1, +2, +3, +9
2091+8, +4, +6, +7
2092weiref = weight
2093arref = transpose(array)
2094call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2095rsize
2096+4
2097weiref(1:rsize)
2098+1, +2, +1, +2
2099arref(1:rsize, :)
2100+1, +8
2101+2, +4
2102+3, +6
2103+9, +7
2104
2105dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2106[dim, ndim, nsam]
2107+2, +1, +8
2108skip
2109+2
2110weight = getUnifRand(-1, 9, nsam)
2111weight
2112+2, +5, +0, +3, -1, +2, +0, +5
2113array = getUnifRand(.false., .true., ndim, nsam)
2114array
2115F, F, T, T, F, F, T, F
2116getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2117F, F, F, F, F, F, F, T, T, T, F, F, F, F, F, F, F
2118arref = array
2119weiref = weight
2120call setRefined(arref, dim, weiref, skip, rsize)
2121rsize
2122+5
2123weiref(1:rsize)
2124+1, +2, +2, +1, +2
2125arref(:, 1:rsize)
2126F, F, T, F, F
2127weiref = weight
2128arref = transpose(array)
2129call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2130rsize
2131+5
2132weiref(1:rsize)
2133+1, +2, +2, +1, +2
2134arref(1:rsize, :)
2135F
2136F
2137T
2138F
2139F
2140
2141dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2142[dim, ndim, nsam]
2143+2, +2, +7
2144skip
2145+2
2146weight = getUnifRand(-1, 9, nsam)
2147weight
2148-1, +0, +5, -1, -1, +2, +7
2149array = getUnifRand(.false., .true., ndim, nsam)
2150array
2151T, F, F, F, T, T, T
2152T, F, F, T, F, T, T
2153getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2154F, F, F, F, F, T, T, T, T, T, T, T, T, T
2155F, F, F, F, F, T, T, T, T, T, T, T, T, T
2156arref = array
2157weiref = weight
2158call setRefined(arref, dim, weiref, skip, rsize)
2159rsize
2160+3
2161weiref(1:rsize)
2162+2, +1, +4
2163arref(:, 1:rsize)
2164F, T, T
2165F, T, T
2166weiref = weight
2167arref = transpose(array)
2168call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2169rsize
2170+3
2171weiref(1:rsize)
2172+2, +1, +4
2173arref(1:rsize, :)
2174F, F
2175T, T
2176T, T
2177
2178dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2179[dim, ndim, nsam]
2180+2, +1, +5
2181skip
2182+1
2183weight = getUnifRand(-1, 9, nsam)
2184weight
2185+1, +4, -1, -1, +3
2186array = getUnifRand(.false., .true., ndim, nsam)
2187array
2188F, T, T, F, T
2189getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2190F, T, T, T, T, T, T, T
2191arref = array
2192weiref = weight
2193call setRefined(arref, dim, weiref, skip, rsize)
2194rsize
2195+3
2196weiref(1:rsize)
2197+1, +4, +3
2198arref(:, 1:rsize)
2199F, T, T
2200weiref = weight
2201arref = transpose(array)
2202call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2203rsize
2204+3
2205weiref(1:rsize)
2206+1, +4, +3
2207arref(1:rsize, :)
2208F
2209T
2210T
2211
2212dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2213[dim, ndim, nsam]
2214+2, +2, +3
2215skip
2216+4
2217weight = getUnifRand(-1, 9, nsam)
2218weight
2219+0, -1, +4
2220array = getUnifRand(.false., .true., ndim, nsam)
2221array
2222F, T, T
2223F, T, T
2224getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2225T, T, T, T
2226T, T, T, T
2227arref = array
2228weiref = weight
2229call setRefined(arref, dim, weiref, skip, rsize)
2230rsize
2231+1
2232weiref(1:rsize)
2233+1
2234arref(:, 1:rsize)
2235T
2236T
2237weiref = weight
2238arref = transpose(array)
2239call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2240rsize
2241+1
2242weiref(1:rsize)
2243+1
2244arref(1:rsize, :)
2245T, T
2246
2247dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2248[dim, ndim, nsam]
2249+2, +1, +0
2250skip
2251+3
2252weight = getUnifRand(-1, 9, nsam)
2253weight
2254
2255array = getUnifRand(.false., .true., ndim, nsam)
2256array
2257
2258getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2259
2260arref = array
2261weiref = weight
2262call setRefined(arref, dim, weiref, skip, rsize)
2263rsize
2264+0
2265weiref(1:rsize)
2266
2267arref(:, 1:rsize)
2268
2269weiref = weight
2270arref = transpose(array)
2271call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2272rsize
2273+0
2274weiref(1:rsize)
2275
2276arref(1:rsize, :)
2277
2278dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2279[dim, ndim, nsam]
2280+2, +1, +5
2281skip
2282+4
2283weight = getUnifRand(-1, 9, nsam)
2284weight
2285+6, +8, +6, +9, +1
2286array = getUnifRand(.false., .true., ndim, nsam)
2287array
2288T, F, F, T, T
2289getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2290T, T, T, T, T, T, F, F, F, F, F, F, F, F, F, F, F, F, F, F, T, T, T, T, T, T, T, T, T, T
2291arref = array
2292weiref = weight
2293call setRefined(arref, dim, weiref, skip, rsize)
2294rsize
2295+4
2296weiref(1:rsize)
2297+1, +2, +2, +2
2298arref(:, 1:rsize)
2299T, F, F, T
2300weiref = weight
2301arref = transpose(array)
2302call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2303rsize
2304+4
2305weiref(1:rsize)
2306+1, +2, +2, +2
2307arref(1:rsize, :)
2308T
2309F
2310F
2311T
2312
2313dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2314[dim, ndim, nsam]
2315+2, +2, +7
2316skip
2317+1
2318weight = getUnifRand(-1, 9, nsam)
2319weight
2320-1, -1, -1, +1, +1, +4, +1
2321array = getUnifRand(.false., .true., ndim, nsam)
2322array
2323F, T, F, T, F, F, F
2324T, T, F, T, T, T, T
2325getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2326T, F, F, F, F, F, F
2327T, T, T, T, T, T, T
2328arref = array
2329weiref = weight
2330call setRefined(arref, dim, weiref, skip, rsize)
2331rsize
2332+4
2333weiref(1:rsize)
2334+1, +1, +4, +1
2335arref(:, 1:rsize)
2336T, F, F, F
2337T, T, T, T
2338weiref = weight
2339arref = transpose(array)
2340call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2341rsize
2342+4
2343weiref(1:rsize)
2344+1, +1, +4, +1
2345arref(1:rsize, :)
2346T, T
2347F, T
2348F, T
2349F, T
2350
2351dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2352[dim, ndim, nsam]
2353+2, +2, +1
2354skip
2355+3
2356weight = getUnifRand(-1, 9, nsam)
2357weight
2358+8
2359array = getUnifRand(.false., .true., ndim, nsam)
2360array
2361F
2362T
2363getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2364F, F, F, F, F, F, F, F
2365T, T, T, T, T, T, T, T
2366arref = array
2367weiref = weight
2368call setRefined(arref, dim, weiref, skip, rsize)
2369rsize
2370+1
2371weiref(1:rsize)
2372+2
2373arref(:, 1:rsize)
2374F
2375T
2376weiref = weight
2377arref = transpose(array)
2378call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2379rsize
2380+1
2381weiref(1:rsize)
2382+2
2383arref(1:rsize, :)
2384F, T
2385
2386dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2387[dim, ndim, nsam]
2388+2, +2, +6
2389skip
2390+1
2391weight = getUnifRand(-1, 9, nsam)
2392weight
2393+0, +2, +4, +3, +0, +5
2394array = getUnifRand(.false., .true., ndim, nsam)
2395array
2396T, T, F, T, F, F
2397T, T, T, T, F, F
2398getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2399T, T, F, F, F, F, T, T, T, F, F, F, F, F
2400T, T, T, T, T, T, T, T, T, F, F, F, F, F
2401arref = array
2402weiref = weight
2403call setRefined(arref, dim, weiref, skip, rsize)
2404rsize
2405+4
2406weiref(1:rsize)
2407+2, +4, +3, +5
2408arref(:, 1:rsize)
2409T, F, T, F
2410T, T, T, F
2411weiref = weight
2412arref = transpose(array)
2413call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2414rsize
2415+4
2416weiref(1:rsize)
2417+2, +4, +3, +5
2418arref(1:rsize, :)
2419T, T
2420F, T
2421T, T
2422F, F
2423
2424dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2425[dim, ndim, nsam]
2426+2, +3, +0
2427skip
2428+2
2429weight = getUnifRand(-1, 9, nsam)
2430weight
2431
2432array = getUnifRand(.false., .true., ndim, nsam)
2433array
2434
2435
2436
2437getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2438
2439
2440
2441arref = array
2442weiref = weight
2443call setRefined(arref, dim, weiref, skip, rsize)
2444rsize
2445+0
2446weiref(1:rsize)
2447
2448arref(:, 1:rsize)
2449
2450
2451
2452weiref = weight
2453arref = transpose(array)
2454call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2455rsize
2456+0
2457weiref(1:rsize)
2458
2459arref(1:rsize, :)
2460
2461dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2462[dim, ndim, nsam]
2463+2, +2, +8
2464skip
2465+1
2466weight = getUnifRand(-1, 9, nsam)
2467weight
2468+9, +6, +4, +3, +2, +0, +1, +3
2469array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2470array
2471(+6.00000000, +9.00000000), (+8.00000000, +0.00000000), (+0.00000000, +9.00000000), (+2.00000000, +2.00000000), (+1.00000000, +0.00000000), (+5.00000000, +3.00000000), (+0.00000000, +7.00000000), (+3.00000000, +0.00000000)
2472(+2.00000000, +8.00000000), (+7.00000000, +6.00000000), (+8.00000000, +7.00000000), (+9.00000000, +6.00000000), (+8.00000000, +6.00000000), (+5.00000000, +0.00000000), (+1.00000000, +8.00000000), (+1.00000000, +1.00000000)
2473getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2474(+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+8.00000000, +0.00000000), (+0.00000000, +9.00000000), (+0.00000000, +9.00000000), (+0.00000000, +9.00000000), (+0.00000000, +9.00000000), (+2.00000000, +2.00000000), (+2.00000000, +2.00000000), (+2.00000000, +2.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+0.00000000, +7.00000000), (+3.00000000, +0.00000000), (+3.00000000, +0.00000000), (+3.00000000, +0.00000000)
2475(+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+7.00000000, +6.00000000), (+7.00000000, +6.00000000), (+7.00000000, +6.00000000), (+7.00000000, +6.00000000), (+7.00000000, +6.00000000), (+7.00000000, +6.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+9.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+1.00000000, +8.00000000), (+1.00000000, +1.00000000), (+1.00000000, +1.00000000), (+1.00000000, +1.00000000)
2476arref = array
2477weiref = weight
2478call setRefined(arref, dim, weiref, skip, rsize)
2479rsize
2480+7
2481weiref(1:rsize)
2482+9, +6, +4, +3, +2, +1, +3
2483arref(:, 1:rsize)
2484(+6.00000000, +9.00000000), (+8.00000000, +0.00000000), (+0.00000000, +9.00000000), (+2.00000000, +2.00000000), (+1.00000000, +0.00000000), (+0.00000000, +7.00000000), (+3.00000000, +0.00000000)
2485(+2.00000000, +8.00000000), (+7.00000000, +6.00000000), (+8.00000000, +7.00000000), (+9.00000000, +6.00000000), (+8.00000000, +6.00000000), (+1.00000000, +8.00000000), (+1.00000000, +1.00000000)
2486weiref = weight
2487arref = transpose(array)
2488call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2489rsize
2490+7
2491weiref(1:rsize)
2492+9, +6, +4, +3, +2, +1, +3
2493arref(1:rsize, :)
2494(+6.00000000, +9.00000000), (+2.00000000, +8.00000000)
2495(+8.00000000, +0.00000000), (+7.00000000, +6.00000000)
2496(+0.00000000, +9.00000000), (+8.00000000, +7.00000000)
2497(+2.00000000, +2.00000000), (+9.00000000, +6.00000000)
2498(+1.00000000, +0.00000000), (+8.00000000, +6.00000000)
2499(+0.00000000, +7.00000000), (+1.00000000, +8.00000000)
2500(+3.00000000, +0.00000000), (+1.00000000, +1.00000000)
2501
2502dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2503[dim, ndim, nsam]
2504+2, +3, +8
2505skip
2506+1
2507weight = getUnifRand(-1, 9, nsam)
2508weight
2509+6, +0, +2, -1, +8, +5, +9, +4
2510array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2511array
2512(+1.00000000, +1.00000000), (+6.00000000, +3.00000000), (+2.00000000, +9.00000000), (+0.00000000, +5.00000000), (+2.00000000, +5.00000000), (+3.00000000, +3.00000000), (+1.00000000, +0.00000000), (+5.00000000, +0.00000000)
2513(+4.00000000, +0.00000000), (+6.00000000, +2.00000000), (+4.00000000, +6.00000000), (+8.00000000, +5.00000000), (+1.00000000, +4.00000000), (+0.00000000, +1.00000000), (+2.00000000, +8.00000000), (+8.00000000, +2.00000000)
2514(+8.00000000, +7.00000000), (+5.00000000, +4.00000000), (+7.00000000, +9.00000000), (+5.00000000, +6.00000000), (+6.00000000, +6.00000000), (+3.00000000, +3.00000000), (+1.00000000, +2.00000000), (+2.00000000, +0.00000000)
2515getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2516(+1.00000000, +1.00000000), (+1.00000000, +1.00000000), (+1.00000000, +1.00000000), (+1.00000000, +1.00000000), (+1.00000000, +1.00000000), (+1.00000000, +1.00000000), (+2.00000000, +9.00000000), (+2.00000000, +9.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+2.00000000, +5.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+1.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000), (+5.00000000, +0.00000000)
2517(+4.00000000, +0.00000000), (+4.00000000, +0.00000000), (+4.00000000, +0.00000000), (+4.00000000, +0.00000000), (+4.00000000, +0.00000000), (+4.00000000, +0.00000000), (+4.00000000, +6.00000000), (+4.00000000, +6.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+1.00000000, +4.00000000), (+0.00000000, +1.00000000), (+0.00000000, +1.00000000), (+0.00000000, +1.00000000), (+0.00000000, +1.00000000), (+0.00000000, +1.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+2.00000000, +8.00000000), (+8.00000000, +2.00000000), (+8.00000000, +2.00000000), (+8.00000000, +2.00000000), (+8.00000000, +2.00000000)
2518(+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+8.00000000, +7.00000000), (+7.00000000, +9.00000000), (+7.00000000, +9.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+6.00000000, +6.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+1.00000000, +2.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000), (+2.00000000, +0.00000000)
2519arref = array
2520weiref = weight
2521call setRefined(arref, dim, weiref, skip, rsize)
2522rsize
2523+6
2524weiref(1:rsize)
2525+6, +2, +8, +5, +9, +4
2526arref(:, 1:rsize)
2527(+1.00000000, +1.00000000), (+2.00000000, +9.00000000), (+2.00000000, +5.00000000), (+3.00000000, +3.00000000), (+1.00000000, +0.00000000), (+5.00000000, +0.00000000)
2528(+4.00000000, +0.00000000), (+4.00000000, +6.00000000), (+1.00000000, +4.00000000), (+0.00000000, +1.00000000), (+2.00000000, +8.00000000), (+8.00000000, +2.00000000)
2529(+8.00000000, +7.00000000), (+7.00000000, +9.00000000), (+6.00000000, +6.00000000), (+3.00000000, +3.00000000), (+1.00000000, +2.00000000), (+2.00000000, +0.00000000)
2530weiref = weight
2531arref = transpose(array)
2532call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2533rsize
2534+6
2535weiref(1:rsize)
2536+6, +2, +8, +5, +9, +4
2537arref(1:rsize, :)
2538(+1.00000000, +1.00000000), (+4.00000000, +0.00000000), (+8.00000000, +7.00000000)
2539(+2.00000000, +9.00000000), (+4.00000000, +6.00000000), (+7.00000000, +9.00000000)
2540(+2.00000000, +5.00000000), (+1.00000000, +4.00000000), (+6.00000000, +6.00000000)
2541(+3.00000000, +3.00000000), (+0.00000000, +1.00000000), (+3.00000000, +3.00000000)
2542(+1.00000000, +0.00000000), (+2.00000000, +8.00000000), (+1.00000000, +2.00000000)
2543(+5.00000000, +0.00000000), (+8.00000000, +2.00000000), (+2.00000000, +0.00000000)
2544
2545dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2546[dim, ndim, nsam]
2547+2, +1, +8
2548skip
2549+3
2550weight = getUnifRand(-1, 9, nsam)
2551weight
2552+8, +3, +4, -1, +7, +6, +4, +5
2553array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2554array
2555(+8.00000000, +6.00000000), (+1.00000000, +3.00000000), (+5.00000000, +6.00000000), (+3.00000000, +9.00000000), (+2.00000000, +3.00000000), (+6.00000000, +3.00000000), (+9.00000000, +0.00000000), (+0.00000000, +5.00000000)
2556getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2557(+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+8.00000000, +6.00000000), (+1.00000000, +3.00000000), (+1.00000000, +3.00000000), (+1.00000000, +3.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+5.00000000, +6.00000000), (+2.00000000, +3.00000000), (+2.00000000, +3.00000000), (+2.00000000, +3.00000000), (+2.00000000, +3.00000000), (+2.00000000, +3.00000000), (+2.00000000, +3.00000000), (+2.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000)
2558arref = array
2559weiref = weight
2560call setRefined(arref, dim, weiref, skip, rsize)
2561rsize
2562+7
2563weiref(1:rsize)
2564+2, +1, +2, +2, +2, +1, +2
2565arref(:, 1:rsize)
2566(+8.00000000, +6.00000000), (+1.00000000, +3.00000000), (+5.00000000, +6.00000000), (+2.00000000, +3.00000000), (+6.00000000, +3.00000000), (+9.00000000, +0.00000000), (+0.00000000, +5.00000000)
2567weiref = weight
2568arref = transpose(array)
2569call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2570rsize
2571+7
2572weiref(1:rsize)
2573+2, +1, +2, +2, +2, +1, +2
2574arref(1:rsize, :)
2575(+8.00000000, +6.00000000)
2576(+1.00000000, +3.00000000)
2577(+5.00000000, +6.00000000)
2578(+2.00000000, +3.00000000)
2579(+6.00000000, +3.00000000)
2580(+9.00000000, +0.00000000)
2581(+0.00000000, +5.00000000)
2582
2583dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2584[dim, ndim, nsam]
2585+2, +0, +5
2586skip
2587+3
2588weight = getUnifRand(-1, 9, nsam)
2589weight
2590+8, +3, +6, +1, +2
2591array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2592array
2593getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2594arref = array
2595weiref = weight
2596call setRefined(arref, dim, weiref, skip, rsize)
2597rsize
2598+0
2599weiref(1:rsize)
2600
2601arref(:, 1:rsize)
2602weiref = weight
2603arref = transpose(array)
2604call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2605rsize
2606+4
2607weiref(1:rsize)
2608+2, +1, +2, +1
2609arref(1:rsize, :)
2610
2611
2612
2613
2614
2615dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2616[dim, ndim, nsam]
2617+2, +3, +4
2618skip
2619+3
2620weight = getUnifRand(-1, 9, nsam)
2621weight
2622+3, +4, +3, +1
2623array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2624array
2625(+9.00000000, +5.00000000), (+3.00000000, +6.00000000), (+7.00000000, +3.00000000), (+4.00000000, +9.00000000)
2626(+8.00000000, +5.00000000), (+2.00000000, +9.00000000), (+1.00000000, +5.00000000), (+2.00000000, +5.00000000)
2627(+0.00000000, +0.00000000), (+7.00000000, +7.00000000), (+5.00000000, +5.00000000), (+9.00000000, +0.00000000)
2628getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2629(+9.00000000, +5.00000000), (+9.00000000, +5.00000000), (+9.00000000, +5.00000000), (+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+3.00000000, +6.00000000), (+7.00000000, +3.00000000), (+7.00000000, +3.00000000), (+7.00000000, +3.00000000), (+4.00000000, +9.00000000)
2630(+8.00000000, +5.00000000), (+8.00000000, +5.00000000), (+8.00000000, +5.00000000), (+2.00000000, +9.00000000), (+2.00000000, +9.00000000), (+2.00000000, +9.00000000), (+2.00000000, +9.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+1.00000000, +5.00000000), (+2.00000000, +5.00000000)
2631(+0.00000000, +0.00000000), (+0.00000000, +0.00000000), (+0.00000000, +0.00000000), (+7.00000000, +7.00000000), (+7.00000000, +7.00000000), (+7.00000000, +7.00000000), (+7.00000000, +7.00000000), (+5.00000000, +5.00000000), (+5.00000000, +5.00000000), (+5.00000000, +5.00000000), (+9.00000000, +0.00000000)
2632arref = array
2633weiref = weight
2634call setRefined(arref, dim, weiref, skip, rsize)
2635rsize
2636+3
2637weiref(1:rsize)
2638+1, +1, +1
2639arref(:, 1:rsize)
2640(+9.00000000, +5.00000000), (+3.00000000, +6.00000000), (+7.00000000, +3.00000000)
2641(+8.00000000, +5.00000000), (+2.00000000, +9.00000000), (+1.00000000, +5.00000000)
2642(+0.00000000, +0.00000000), (+7.00000000, +7.00000000), (+5.00000000, +5.00000000)
2643weiref = weight
2644arref = transpose(array)
2645call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2646rsize
2647+3
2648weiref(1:rsize)
2649+1, +1, +1
2650arref(1:rsize, :)
2651(+9.00000000, +5.00000000), (+8.00000000, +5.00000000), (+0.00000000, +0.00000000)
2652(+3.00000000, +6.00000000), (+2.00000000, +9.00000000), (+7.00000000, +7.00000000)
2653(+7.00000000, +3.00000000), (+1.00000000, +5.00000000), (+5.00000000, +5.00000000)
2654
2655dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2656[dim, ndim, nsam]
2657+2, +1, +8
2658skip
2659+3
2660weight = getUnifRand(-1, 9, nsam)
2661weight
2662-1, -1, +7, -1, -1, +4, +7, +8
2663array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2664array
2665(+0.00000000, +7.00000000), (+6.00000000, +3.00000000), (+9.00000000, +0.00000000), (+5.00000000, +8.00000000), (+1.00000000, +9.00000000), (+9.00000000, +3.00000000), (+2.00000000, +4.00000000), (+1.00000000, +7.00000000)
2666getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2667(+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +0.00000000), (+9.00000000, +3.00000000), (+9.00000000, +3.00000000), (+9.00000000, +3.00000000), (+9.00000000, +3.00000000), (+2.00000000, +4.00000000), (+2.00000000, +4.00000000), (+2.00000000, +4.00000000), (+2.00000000, +4.00000000), (+2.00000000, +4.00000000), (+2.00000000, +4.00000000), (+2.00000000, +4.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000), (+1.00000000, +7.00000000)
2668arref = array
2669weiref = weight
2670call setRefined(arref, dim, weiref, skip, rsize)
2671rsize
2672+4
2673weiref(1:rsize)
2674+2, +1, +3, +2
2675arref(:, 1:rsize)
2676(+9.00000000, +0.00000000), (+9.00000000, +3.00000000), (+2.00000000, +4.00000000), (+1.00000000, +7.00000000)
2677weiref = weight
2678arref = transpose(array)
2679call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2680rsize
2681+4
2682weiref(1:rsize)
2683+2, +1, +3, +2
2684arref(1:rsize, :)
2685(+9.00000000, +0.00000000)
2686(+9.00000000, +3.00000000)
2687(+2.00000000, +4.00000000)
2688(+1.00000000, +7.00000000)
2689
2690dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2691[dim, ndim, nsam]
2692+2, +1, +0
2693skip
2694+3
2695weight = getUnifRand(-1, 9, nsam)
2696weight
2697
2698array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2699array
2700
2701getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2702
2703arref = array
2704weiref = weight
2705call setRefined(arref, dim, weiref, skip, rsize)
2706rsize
2707+0
2708weiref(1:rsize)
2709
2710arref(:, 1:rsize)
2711
2712weiref = weight
2713arref = transpose(array)
2714call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2715rsize
2716+0
2717weiref(1:rsize)
2718
2719arref(1:rsize, :)
2720
2721dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2722[dim, ndim, nsam]
2723+2, +2, +7
2724skip
2725+1
2726weight = getUnifRand(-1, 9, nsam)
2727weight
2728+9, +1, +6, +0, +7, -1, +1
2729array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2730array
2731(+7.00000000, +4.00000000), (+2.00000000, +9.00000000), (+5.00000000, +3.00000000), (+1.00000000, +8.00000000), (+2.00000000, +7.00000000), (+9.00000000, +0.00000000), (+5.00000000, +3.00000000)
2732(+6.00000000, +9.00000000), (+4.00000000, +5.00000000), (+0.00000000, +5.00000000), (+6.00000000, +9.00000000), (+8.00000000, +3.00000000), (+0.00000000, +7.00000000), (+8.00000000, +4.00000000)
2733getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2734(+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+7.00000000, +4.00000000), (+2.00000000, +9.00000000), (+5.00000000, +3.00000000), (+5.00000000, +3.00000000), (+5.00000000, +3.00000000), (+5.00000000, +3.00000000), (+5.00000000, +3.00000000), (+5.00000000, +3.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+2.00000000, +7.00000000), (+5.00000000, +3.00000000)
2735(+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+6.00000000, +9.00000000), (+4.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+0.00000000, +5.00000000), (+8.00000000, +3.00000000), (+8.00000000, +3.00000000), (+8.00000000, +3.00000000), (+8.00000000, +3.00000000), (+8.00000000, +3.00000000), (+8.00000000, +3.00000000), (+8.00000000, +3.00000000), (+8.00000000, +4.00000000)
2736arref = array
2737weiref = weight
2738call setRefined(arref, dim, weiref, skip, rsize)
2739rsize
2740+5
2741weiref(1:rsize)
2742+9, +1, +6, +7, +1
2743arref(:, 1:rsize)
2744(+7.00000000, +4.00000000), (+2.00000000, +9.00000000), (+5.00000000, +3.00000000), (+2.00000000, +7.00000000), (+5.00000000, +3.00000000)
2745(+6.00000000, +9.00000000), (+4.00000000, +5.00000000), (+0.00000000, +5.00000000), (+8.00000000, +3.00000000), (+8.00000000, +4.00000000)
2746weiref = weight
2747arref = transpose(array)
2748call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2749rsize
2750+5
2751weiref(1:rsize)
2752+9, +1, +6, +7, +1
2753arref(1:rsize, :)
2754(+7.00000000, +4.00000000), (+6.00000000, +9.00000000)
2755(+2.00000000, +9.00000000), (+4.00000000, +5.00000000)
2756(+5.00000000, +3.00000000), (+0.00000000, +5.00000000)
2757(+2.00000000, +7.00000000), (+8.00000000, +3.00000000)
2758(+5.00000000, +3.00000000), (+8.00000000, +4.00000000)
2759
2760dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2761[dim, ndim, nsam]
2762+2, +1, +8
2763skip
2764+3
2765weight = getUnifRand(-1, 9, nsam)
2766weight
2767+1, -1, -1, +0, +7, +2, +2, +5
2768array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2769array
2770(+4.00000000, +3.00000000), (+9.00000000, +8.00000000), (+9.00000000, +3.00000000), (+2.00000000, +3.00000000), (+4.00000000, +7.00000000), (+1.00000000, +8.00000000), (+3.00000000, +2.00000000), (+5.00000000, +4.00000000)
2771getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2772(+4.00000000, +3.00000000), (+4.00000000, +7.00000000), (+4.00000000, +7.00000000), (+4.00000000, +7.00000000), (+4.00000000, +7.00000000), (+4.00000000, +7.00000000), (+4.00000000, +7.00000000), (+4.00000000, +7.00000000), (+1.00000000, +8.00000000), (+1.00000000, +8.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+5.00000000, +4.00000000), (+5.00000000, +4.00000000), (+5.00000000, +4.00000000), (+5.00000000, +4.00000000), (+5.00000000, +4.00000000)
2773arref = array
2774weiref = weight
2775call setRefined(arref, dim, weiref, skip, rsize)
2776rsize
2777+4
2778weiref(1:rsize)
2779+2, +1, +1, +1
2780arref(:, 1:rsize)
2781(+4.00000000, +7.00000000), (+1.00000000, +8.00000000), (+3.00000000, +2.00000000), (+5.00000000, +4.00000000)
2782weiref = weight
2783arref = transpose(array)
2784call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2785rsize
2786+4
2787weiref(1:rsize)
2788+2, +1, +1, +1
2789arref(1:rsize, :)
2790(+4.00000000, +7.00000000)
2791(+1.00000000, +8.00000000)
2792(+3.00000000, +2.00000000)
2793(+5.00000000, +4.00000000)
2794
2795dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2796[dim, ndim, nsam]
2797+2, +1, +7
2798skip
2799+4
2800weight = getUnifRand(-1, 9, nsam)
2801weight
2802+4, -1, +9, +3, +3, +9, +9
2803array = cmplx(getUnifRand(0, 9, ndim, nsam), getUnifRand(0, 9, ndim, nsam), TKG)
2804array
2805(+6.00000000, +2.00000000), (+7.00000000, +4.00000000), (+3.00000000, +2.00000000), (+9.00000000, +8.00000000), (+0.00000000, +9.00000000), (+6.00000000, +3.00000000), (+3.00000000, +3.00000000)
2806getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2807(+6.00000000, +2.00000000), (+6.00000000, +2.00000000), (+6.00000000, +2.00000000), (+6.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+3.00000000, +2.00000000), (+9.00000000, +8.00000000), (+9.00000000, +8.00000000), (+9.00000000, +8.00000000), (+0.00000000, +9.00000000), (+0.00000000, +9.00000000), (+0.00000000, +9.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+6.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000), (+3.00000000, +3.00000000)
2808arref = array
2809weiref = weight
2810call setRefined(arref, dim, weiref, skip, rsize)
2811rsize
2812+5
2813weiref(1:rsize)
2814+1, +2, +1, +3, +2
2815arref(:, 1:rsize)
2816(+6.00000000, +2.00000000), (+3.00000000, +2.00000000), (+9.00000000, +8.00000000), (+6.00000000, +3.00000000), (+3.00000000, +3.00000000)
2817weiref = weight
2818arref = transpose(array)
2819call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2820rsize
2821+5
2822weiref(1:rsize)
2823+1, +2, +1, +3, +2
2824arref(1:rsize, :)
2825(+6.00000000, +2.00000000)
2826(+3.00000000, +2.00000000)
2827(+9.00000000, +8.00000000)
2828(+6.00000000, +3.00000000)
2829(+3.00000000, +3.00000000)
2830
2831dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2832[dim, ndim, nsam]
2833+2, +2, +7
2834skip
2835+2
2836weight = getUnifRand(-1, 9, nsam)
2837weight
2838+8, -1, +2, +6, -1, +9, +8
2839array = getUnifRand(0, 9, ndim, nsam)
2840array
2841+7.00000000, +1.00000000, +8.00000000, +2.00000000, +9.00000000, +5.00000000, +6.00000000
2842+4.00000000, +4.00000000, +6.00000000, +3.00000000, +7.00000000, +2.00000000, +7.00000000
2843getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2844+7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +8.00000000, +8.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000, +6.00000000
2845+4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +6.00000000, +6.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000
2846arref = array
2847weiref = weight
2848call setRefined(arref, dim, weiref, skip, rsize)
2849rsize
2850+5
2851weiref(1:rsize)
2852+4, +1, +3, +4, +4
2853arref(:, 1:rsize)
2854+7.00000000, +8.00000000, +2.00000000, +5.00000000, +6.00000000
2855+4.00000000, +6.00000000, +3.00000000, +2.00000000, +7.00000000
2856weiref = weight
2857arref = transpose(array)
2858call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2859rsize
2860+5
2861weiref(1:rsize)
2862+4, +1, +3, +4, +4
2863arref(1:rsize, :)
2864+7.00000000, +4.00000000
2865+8.00000000, +6.00000000
2866+2.00000000, +3.00000000
2867+5.00000000, +2.00000000
2868+6.00000000, +7.00000000
2869
2870dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2871[dim, ndim, nsam]
2872+2, +0, +2
2873skip
2874+4
2875weight = getUnifRand(-1, 9, nsam)
2876weight
2877+2, +2
2878array = getUnifRand(0, 9, ndim, nsam)
2879array
2880getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2881arref = array
2882weiref = weight
2883call setRefined(arref, dim, weiref, skip, rsize)
2884rsize
2885+0
2886weiref(1:rsize)
2887
2888arref(:, 1:rsize)
2889weiref = weight
2890arref = transpose(array)
2891call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2892rsize
2893+1
2894weiref(1:rsize)
2895+1
2896arref(1:rsize, :)
2897
2898
2899dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2900[dim, ndim, nsam]
2901+2, +3, +7
2902skip
2903+1
2904weight = getUnifRand(-1, 9, nsam)
2905weight
2906+4, +0, +2, +5, +8, +9, +6
2907array = getUnifRand(0, 9, ndim, nsam)
2908array
2909+2.00000000, +9.00000000, +7.00000000, +1.00000000, +1.00000000, +4.00000000, +9.00000000
2910+4.00000000, +1.00000000, +7.00000000, +4.00000000, +4.00000000, +8.00000000, +9.00000000
2911+9.00000000, +9.00000000, +9.00000000, +9.00000000, +7.00000000, +5.00000000, +2.00000000
2912getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2913+2.00000000, +2.00000000, +2.00000000, +2.00000000, +7.00000000, +7.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000
2914+4.00000000, +4.00000000, +4.00000000, +4.00000000, +7.00000000, +7.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000
2915+9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +5.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000, +2.00000000
2916arref = array
2917weiref = weight
2918call setRefined(arref, dim, weiref, skip, rsize)
2919rsize
2920+6
2921weiref(1:rsize)
2922+4, +2, +5, +8, +9, +6
2923arref(:, 1:rsize)
2924+2.00000000, +7.00000000, +1.00000000, +1.00000000, +4.00000000, +9.00000000
2925+4.00000000, +7.00000000, +4.00000000, +4.00000000, +8.00000000, +9.00000000
2926+9.00000000, +9.00000000, +9.00000000, +7.00000000, +5.00000000, +2.00000000
2927weiref = weight
2928arref = transpose(array)
2929call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2930rsize
2931+6
2932weiref(1:rsize)
2933+4, +2, +5, +8, +9, +6
2934arref(1:rsize, :)
2935+2.00000000, +4.00000000, +9.00000000
2936+7.00000000, +7.00000000, +9.00000000
2937+1.00000000, +4.00000000, +9.00000000
2938+1.00000000, +4.00000000, +7.00000000
2939+4.00000000, +8.00000000, +5.00000000
2940+9.00000000, +9.00000000, +2.00000000
2941
2942dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2943[dim, ndim, nsam]
2944+2, +2, +2
2945skip
2946+2
2947weight = getUnifRand(-1, 9, nsam)
2948weight
2949+8, +9
2950array = getUnifRand(0, 9, ndim, nsam)
2951array
2952+9.00000000, +8.00000000
2953+3.00000000, +3.00000000
2954getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2955+9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +9.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000, +8.00000000
2956+3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000
2957arref = array
2958weiref = weight
2959call setRefined(arref, dim, weiref, skip, rsize)
2960rsize
2961+2
2962weiref(1:rsize)
2963+4, +4
2964arref(:, 1:rsize)
2965+9.00000000, +8.00000000
2966+3.00000000, +3.00000000
2967weiref = weight
2968arref = transpose(array)
2969call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
2970rsize
2971+2
2972weiref(1:rsize)
2973+4, +4
2974arref(1:rsize, :)
2975+9.00000000, +3.00000000
2976+8.00000000, +3.00000000
2977
2978dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
2979[dim, ndim, nsam]
2980+2, +1, +5
2981skip
2982+4
2983weight = getUnifRand(-1, 9, nsam)
2984weight
2985+7, +9, +4, +0, +3
2986array = getUnifRand(0, 9, ndim, nsam)
2987array
2988+4.00000000, +3.00000000, +4.00000000, +4.00000000, +6.00000000
2989getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
2990+4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +3.00000000, +4.00000000, +4.00000000, +4.00000000, +4.00000000, +6.00000000, +6.00000000, +6.00000000
2991arref = array
2992weiref = weight
2993call setRefined(arref, dim, weiref, skip, rsize)
2994rsize
2995+3
2996weiref(1:rsize)
2997+1, +3, +1
2998arref(:, 1:rsize)
2999+4.00000000, +3.00000000, +4.00000000
3000weiref = weight
3001arref = transpose(array)
3002call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
3003rsize
3004+3
3005weiref(1:rsize)
3006+1, +3, +1
3007arref(1:rsize, :)
3008+4.00000000
3009+3.00000000
3010+4.00000000
3011
3012dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
3013[dim, ndim, nsam]
3014+2, +1, +6
3015skip
3016+2
3017weight = getUnifRand(-1, 9, nsam)
3018weight
3019+6, +0, +1, +9, +2, +1
3020array = getUnifRand(0, 9, ndim, nsam)
3021array
3022+7.00000000, +1.00000000, +6.00000000, +1.00000000, +0.00000000, +1.00000000
3023getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
3024+7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +7.00000000, +6.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +0.00000000, +0.00000000, +1.00000000
3025arref = array
3026weiref = weight
3027call setRefined(arref, dim, weiref, skip, rsize)
3028rsize
3029+3
3030weiref(1:rsize)
3031+3, +5, +1
3032arref(:, 1:rsize)
3033+7.00000000, +1.00000000, +0.00000000
3034weiref = weight
3035arref = transpose(array)
3036call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
3037rsize
3038+3
3039weiref(1:rsize)
3040+3, +5, +1
3041arref(1:rsize, :)
3042+7.00000000
3043+1.00000000
3044+0.00000000
3045
3046dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
3047[dim, ndim, nsam]
3048+2, +1, +3
3049skip
3050+4
3051weight = getUnifRand(-1, 9, nsam)
3052weight
3053+3, +0, +7
3054array = getUnifRand(0, 9, ndim, nsam)
3055array
3056+9.00000000, +6.00000000, +1.00000000
3057getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
3058+9.00000000, +9.00000000, +9.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000, +1.00000000
3059arref = array
3060weiref = weight
3061call setRefined(arref, dim, weiref, skip, rsize)
3062rsize
3063+1
3064weiref(1:rsize)
3065+2
3066arref(:, 1:rsize)
3067+1.00000000
3068weiref = weight
3069arref = transpose(array)
3070call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
3071rsize
3072+1
3073weiref(1:rsize)
3074+2
3075arref(1:rsize, :)
3076+1.00000000
3077
3078dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
3079[dim, ndim, nsam]
3080+2, +0, +1
3081skip
3082+1
3083weight = getUnifRand(-1, 9, nsam)
3084weight
3085+5
3086array = getUnifRand(0, 9, ndim, nsam)
3087array
3088getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
3089arref = array
3090weiref = weight
3091call setRefined(arref, dim, weiref, skip, rsize)
3092rsize
3093+0
3094weiref(1:rsize)
3095
3096arref(:, 1:rsize)
3097weiref = weight
3098arref = transpose(array)
3099call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
3100rsize
3101+1
3102weiref(1:rsize)
3103+5
3104arref(1:rsize, :)
3105
3106
3107dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
3108[dim, ndim, nsam]
3109+2, +3, +1
3110skip
3111+2
3112weight = getUnifRand(-1, 9, nsam)
3113weight
3114+2
3115array = getUnifRand(0, 9, ndim, nsam)
3116array
3117+9.00000000
3118+6.00000000
3119+8.00000000
3120getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
3121+9.00000000, +9.00000000
3122+6.00000000, +6.00000000
3123+8.00000000, +8.00000000
3124arref = array
3125weiref = weight
3126call setRefined(arref, dim, weiref, skip, rsize)
3127rsize
3128+1
3129weiref(1:rsize)
3130+1
3131arref(:, 1:rsize)
3132+9.00000000
3133+6.00000000
3134+8.00000000
3135weiref = weight
3136arref = transpose(array)
3137call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
3138rsize
3139+1
3140weiref(1:rsize)
3141+1
3142arref(1:rsize, :)
3143+9.00000000, +6.00000000, +8.00000000
3144
3145dim = 2; ndim = merge(0, getUnifRand(1, 3), isHead(.1)); nsam = getUnifRand(0, 9); skip = getUnifRand(1, 4)
3146[dim, ndim, nsam]
3147+2, +0, +5
3148skip
3149+2
3150weight = getUnifRand(-1, 9, nsam)
3151weight
3152+8, -1, +3, +2, -1
3153array = getUnifRand(0, 9, ndim, nsam)
3154array
3155getVerbose(array, weight, sum(weight, mask = weight > 0), dim)
3156arref = array
3157weiref = weight
3158call setRefined(arref, dim, weiref, skip, rsize)
3159rsize
3160+0
3161weiref(1:rsize)
3162
3163arref(:, 1:rsize)
3164weiref = weight
3165arref = transpose(array)
3166call setRefined(arref, 3_IK - dim, weiref, skip, rsize)
3167rsize
3168+3
3169weiref(1:rsize)
3170+4, +1, +1
3171arref(1:rsize, :)
3172
3173
3174
3175
3176
Test:
test_pm_arrayRefine


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.

Copyright
Computational Data Science Lab
Author:
Amir Shahmoradi, Saturday 1:48 AM, August 20, 2016, Institute for Computational Engineering and Sciences, UT Austin, TX

Definition at line 1018 of file pm_arrayRefine.F90.

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