19 type(display_type) :: disp
20 integer(IK) :: itry, ntry
= 5
21 integer(IK) :: dim, isam, ndim, nsam
25 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
26 call disp%show(
"!Transform 2D sample along the second dimension.")
27 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
32 real(RKG),
allocatable :: tlate(:), tform(:,:), sample(:,:), affinInv(:,:), affinity(:,:)
35 call disp%show(
"dim = 2; ndim = getUnifRand(1, 3); nsam = getUnifRand(4, 5)")
39 call disp%show(
"tlate = getUnifRand(-5, +5, ndim)")
43 call disp%show(
"sample = getUnifRand(-9, +9, ndim, nsam)")
47 call disp%show(
"call setResized(affinity, shape(sample, IK))")
49 call disp%show(
"call setResized(affinInv, shape(sample, IK))")
51 call disp%show(
"tform = getCovRand(1., ndim)")
55 call disp%show(
"call setAffinity(affinity, sample, dim, tform, genrecmat) ! whole sample transformation.")
56 call setAffinity(affinity, sample, dim, tform, genrecmat)
57 call disp%show(
"transpose(affinity)")
58 call disp%show(
transpose(affinity) )
60 call disp%show(
" call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat) ! point-wise transformation")
63 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
65 call disp%show(
"transpose(affinity)")
66 call disp%show(
transpose(affinity) )
67 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), genrecmat) ! inverse transformation.")
69 call disp%show(
"transpose(affinInv)")
70 call disp%show(
transpose(affinInv) )
71 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
74 call disp%show(
"call setMatInit(tform(2:,1:ndim-1), lowDia, vlow = 0._RKG, vdia = 0._RKG) ! make tform an upper-diagonal matrix.")
75 call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
78 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = upperDiag) ! whole sample transformation.")
79 call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
80 call disp%show(
"transpose(affinity)")
81 call disp%show(
transpose(affinity) )
83 call disp%show(
" call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag) ! point-wise transformation")
86 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
88 call disp%show(
"transpose(affinity)")
89 call disp%show(
transpose(affinity) )
90 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform, upperDiag), class = upperDiag) ! inverse transformation.")
92 call disp%show(
"transpose(affinInv)")
93 call disp%show(
transpose(affinInv) )
94 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
97 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = upperUnit) ! whole sample transformation.")
98 call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
99 call disp%show(
"transpose(affinity)")
100 call disp%show(
transpose(affinity) )
102 call disp%show(
" call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit) ! point-wise transformation")
105 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
107 call disp%show(
"transpose(affinity)")
108 call disp%show(
transpose(affinity) )
109 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), class = upperUnit) ! inverse transformation.")
111 call disp%show(
"transpose(affinInv)")
112 call disp%show(
transpose(affinInv) )
113 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
116 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = lowerDiag) ! whole sample transformation.")
117 call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
118 call disp%show(
"transpose(affinity)")
119 call disp%show(
transpose(affinity) )
121 call disp%show(
" call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag) ! point-wise transformation")
124 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
126 call disp%show(
"transpose(affinity)")
127 call disp%show(
transpose(affinity) )
128 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), class = lowerDiag) ! inverse transformation.")
130 call disp%show(
"transpose(affinInv)")
131 call disp%show(
transpose(affinInv) )
132 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
135 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = lowerUnit) ! whole sample transformation.")
136 call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
137 call disp%show(
"transpose(affinity)")
138 call disp%show(
transpose(affinity) )
140 call disp%show(
" call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit) ! point-wise transformation")
143 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
145 call disp%show(
"transpose(affinity)")
146 call disp%show(
transpose(affinity) )
147 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), class = lowerUnit) ! inverse transformation.")
149 call disp%show(
"transpose(affinInv)")
150 call disp%show(
transpose(affinInv) )
151 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
157 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
158 call disp%show(
"!Transform 2D sample along the first dimension.")
159 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
164 real(RKG),
allocatable :: tlate(:), tform(:,:), sample(:,:), affinInv(:,:), affinity(:,:)
167 call disp%show(
"dim = 1; ndim = getUnifRand(1, 3); nsam = getUnifRand(4, 5)")
171 call disp%show(
"tlate = getUnifRand(-5, +5, ndim)")
175 call disp%show(
"sample = getUnifRand(-9, +9, nsam, ndim)")
179 call disp%show(
"call setResized(affinity, shape(sample, IK))")
181 call disp%show(
"call setResized(affinInv, shape(sample, IK))")
183 call disp%show(
"tform = getCovRand(1., ndim)")
187 call disp%show(
"call setAffinity(affinity, sample, dim, tform) ! whole sample transformation.")
188 call setAffinity(affinity, sample, dim, tform, genrecmat)
192 call disp%show(
" call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat) ! point-wise transformation")
195 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
199 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), genrecmat) ! inverse transformation.")
203 call disp%show(
"sample ! for comparison with affinInv.")
206 call disp%show(
"call setMatInit(tform(2:,1:ndim-1), lowDia, vlow = 0._RKG, vdia = 0._RKG) ! make tform an upper-diagonal matrix.")
207 call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
210 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = upperDiag) ! whole sample transformation.")
211 call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
215 call disp%show(
" call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag) ! point-wise transformation")
218 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
222 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform, upperDiag), class = upperDiag) ! inverse transformation.")
226 call disp%show(
"sample ! for comparison with affinInv.")
229 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = upperUnit) ! whole sample transformation.")
230 call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
234 call disp%show(
" call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit) ! point-wise transformation")
237 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
241 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), class = upperUnit) ! inverse transformation.")
245 call disp%show(
"sample ! for comparison with affinInv.")
248 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = lowerDiag) ! whole sample transformation.")
249 call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
253 call disp%show(
" call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag) ! point-wise transformation")
256 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
260 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), class = lowerDiag) ! inverse transformation.")
264 call disp%show(
"sample ! for comparison with affinInv.")
267 call disp%show(
"call setAffinity(affinity, sample, dim, tform, class = lowerUnit) ! whole sample transformation.")
268 call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
272 call disp%show(
" call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit) ! point-wise transformation")
275 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)
279 call disp%show(
"call setAffinity(affinInv, affinity, dim, getMatInv(tform), class = lowerUnit) ! inverse transformation.")
283 call disp%show(
"sample ! for comparison with affinInv.")
304 type(css_type),
allocatable :: shapes(:), tforms(:)
305 integer(IK) :: ishape, itform
306 integer(IK),
parameter :: nsam
= 1000, ndim
= 2, dim
= 2
307 real(RKG),
allocatable :: angle, tlate(:), cov(:,:), tform(:,:), sample(:,:), affinity(:,:), affinInv(:,:)
311 do ishape
= 1,
size(shapes)
312 if (shapes(ishape)
%val
== "circle")
then
314 elseif (shapes(ishape)
%val
== "square")
then
317 error stop "Unrecognized shape."
321 do itform
= 1,
size(tforms)
322 if (tforms(itform)
%val
== "rotation")
then
324 tform
= reshape([
cos(angle),
-sin(angle),
sin(angle),
cos(angle)], [ndim, ndim])
325 call setAffinity(affinity, sample, dim, tform, genrecmat, tlate)
326 elseif (tforms(itform)
%val
== "warp")
then
327 cov
= reshape([
1., .
5, .
5,
1.], shape
= [ndim, ndim])
329 call setAffinity(affinity, sample, dim, tform, lowerDiag, tlate)
331 error stop "Unrecognized tform."
334 if (
0 /= getErrTableWrite(SK_
"setAffinity."//tforms(itform)
%val
//SK_
"."//shapes(ishape)
%val
//SK_
".sample.txt", sample, trans))
error stop "Failed table-write."
335 if (
0 /= getErrTableWrite(SK_
"setAffinity."//tforms(itform)
%val
//SK_
"."//shapes(ishape)
%val
//SK_
".affinity.txt", affinity, trans))
error stop "Failed table-write."
336 if (
0 /= getErrTableWrite(SK_
"setAffinity."//tforms(itform)
%val
//SK_
"."//shapes(ishape)
%val
//SK_
".affinInv.txt", affinInv, trans))
error stop "Failed table-write."
Generate and return an array of the specified rank and shape of arbitrary intrinsic type and kind wit...
Allocate or resize (shrink or expand) an input allocatable scalar string or array of rank 1....
Generate count evenly spaced points over the interval [x1, x2] if x1 < x2, or [x2,...
Generate and return a random positive-definite (correlation or covariance) matrix using the Gram meth...
Generate and return a (collection) of random vector(s) of size ndim from the ndim-dimensional MultiVa...
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
Generate and return the iostat code resulting from writing the input table of rank 1 or 2 to the spec...
This is a generic method of the derived type display_type with pass attribute.
This is a generic method of the derived type display_type with pass attribute.
Generate and return the upper or the lower Cholesky factorization of the input symmetric positive-def...
Generate and return a matrix of shape (shape(1), shape(2)) with the upper/lower triangle and the diag...
Set the upper/lower triangle and the diagonal elements of the input matrix of arbitrary shape (:,...
Generate and return the full inverse of an input matrix of general or triangular form directly or thr...
Generate a sample of shape (nsam), or (ndim, nsam) or (nsam, ndim) that is shifted by the specified i...
This module contains procedures and generic interfaces for convenient allocation and filling of array...
This module contains procedures and generic interfaces for resizing allocatable arrays of various typ...
This module contains procedures and generic interfaces for generating arrays with linear or logarithm...
This module contains the derived types for generating allocatable containers of scalar,...
This module contains classes and procedures for generating random matrices distributed on the space o...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for computing various statistical quantities related to t...
This module contains classes and procedures for input/output (IO) or generic display operations on st...
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
integer, parameter LK
The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
integer, parameter RKS
The single-precision real kind in Fortran mode. On most platforms, this is an 32-bit real kind.
This module contains relevant mathematical constants.
real(RKB), parameter PI
The scalar real constant of kind with highest available precision RKB representing the irrational num...
This module contains procedures and generic interfaces for computing the Cholesky factorization of po...
This module contains procedures and generic interfaces relevant to generating and initializing matric...
This module contains abstract and concrete derived types and procedures related to the inversion of s...
This module contains classes and procedures for shifting univariate or multivariate samples by arbitr...
This is the css_type type for generating instances of container of scalar of string objects.
Generate and return an object of type display_type.
12+3.00000000,
+2.00000000,
+2.00000000
15+1.00000000,
+1.00000000,
-3.00000000,
+1.00000000,
-9.00000000
16+3.00000000,
+9.00000000,
-3.00000000,
+0.00000000,
+1.00000000
17-7.00000000,
-4.00000000,
-3.00000000,
+0.00000000,
+9.00000000
22+1.00000000,
-0.941470921,
-0.876164198
23-0.941470921,
+1.00000000,
+0.983183086
24-0.876164198,
+0.983183086,
+1.00000000
25call setAffinity(affinity, sample, dim, tform, genrecmat)
27+4.30873632,
-4.82375240,
-4.92661476
28-3.96858120,
+4.12579632,
+3.97248363
29+2.45290542,
-3.12513638,
-3.32105660
30+1.00000000,
-0.941470921,
-0.876164198
31-17.8269501,
+18.3218861,
+17.8686600
33 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
36+4.30873632,
-4.82375240,
-4.92661476
37-3.96858120,
+4.12579632,
+3.97248363
38+2.45290542,
-3.12513638,
-3.32105660
39+1.00000000,
-0.941470921,
-0.876164198
40-17.8269501,
+18.3218861,
+17.8686600
43+1.00000000,
+3.00000000,
-6.99996948
44+0.999984741,
+9.00000000,
-4.00003052
45-2.99997711,
-2.99996948,
-2.99996948
46+1.00000000,
-0.762939453E-5,
+0.762939453E-5
47-9.00006104,
+1.00000000,
+8.99987793
49+1.00000000,
+3.00000000,
-7.00000000
50+1.00000000,
+9.00000000,
-4.00000000
51-3.00000000,
-3.00000000,
-3.00000000
52+1.00000000,
+0.00000000,
+0.00000000
53-9.00000000,
+1.00000000,
+9.00000000
55call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
57+1.00000000,
-0.941470921,
-0.876164198
58+0.00000000,
+1.00000000,
+0.983183086
59+0.00000000,
+0.00000000,
+1.00000000
60call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
62+4.30873632,
-3.88228178,
-7.00000000
63-3.96858120,
+5.06726742,
-4.00000000
64+2.45290565,
-5.94954920,
-3.00000000
65+1.00000000,
+0.00000000,
+0.00000000
66-17.8269501,
+9.84864807,
+9.00000000
68 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
71+4.30873632,
-3.88228178,
-7.00000000
72-3.96858120,
+5.06726742,
-4.00000000
73+2.45290565,
-5.94954920,
-3.00000000
74+1.00000000,
+0.00000000,
+0.00000000
75-17.8269501,
+9.84864807,
+9.00000000
78+0.999999523,
+3.00000000,
-7.00000000
79+1.00000000,
+9.00000000,
-4.00000000
80-2.99999952,
-3.00000000,
-3.00000000
81+1.00000000,
+0.00000000,
+0.00000000
82-9.00000095,
+1.00000000,
+9.00000000
84+1.00000000,
+3.00000000,
-7.00000000
85+1.00000000,
+9.00000000,
-4.00000000
86-3.00000000,
-3.00000000,
-3.00000000
87+1.00000000,
+0.00000000,
+0.00000000
88-9.00000000,
+1.00000000,
+9.00000000
90call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
92+4.30873632,
-3.88228178,
-7.00000000
93-3.96858120,
+5.06726742,
-4.00000000
94+2.45290565,
-5.94954920,
-3.00000000
95+1.00000000,
+0.00000000,
+0.00000000
96-17.8269501,
+9.84864807,
+9.00000000
98 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
101+4.30873632,
-3.88228178,
-7.00000000
102-3.96858120,
+5.06726742,
-4.00000000
103+2.45290565,
-5.94954920,
-3.00000000
104+1.00000000,
+0.00000000,
+0.00000000
105-17.8269501,
+9.84864807,
+9.00000000
108+0.999999523,
+3.00000000,
-7.00000000
109+1.00000000,
+9.00000000,
-4.00000000
110-2.99999952,
-3.00000000,
-3.00000000
111+1.00000000,
+0.00000000,
+0.00000000
112-9.00000095,
+1.00000000,
+9.00000000
114+1.00000000,
+3.00000000,
-7.00000000
115+1.00000000,
+9.00000000,
-4.00000000
116-3.00000000,
-3.00000000,
-3.00000000
117+1.00000000,
+0.00000000,
+0.00000000
118-9.00000000,
+1.00000000,
+9.00000000
120call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
122+1.00000000,
+3.00000000,
-7.00000000
123+1.00000000,
+9.00000000,
-4.00000000
124-3.00000000,
-3.00000000,
-3.00000000
125+1.00000000,
+0.00000000,
+0.00000000
126-9.00000000,
+1.00000000,
+9.00000000
128 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
131+1.00000000,
+3.00000000,
-7.00000000
132+1.00000000,
+9.00000000,
-4.00000000
133-3.00000000,
-3.00000000,
-3.00000000
134+1.00000000,
+0.00000000,
+0.00000000
135-9.00000000,
+1.00000000,
+9.00000000
138+1.00000000,
+3.00000000,
-7.00000000
139+1.00000000,
+9.00000000,
-4.00000000
140-3.00000000,
-3.00000000,
-3.00000000
141+1.00000000,
+0.00000000,
+0.00000000
142-9.00000000,
+1.00000000,
+9.00000000
144+1.00000000,
+3.00000000,
-7.00000000
145+1.00000000,
+9.00000000,
-4.00000000
146-3.00000000,
-3.00000000,
-3.00000000
147+1.00000000,
+0.00000000,
+0.00000000
148-9.00000000,
+1.00000000,
+9.00000000
150call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
152+1.00000000,
+3.00000000,
-7.00000000
153+1.00000000,
+9.00000000,
-4.00000000
154-3.00000000,
-3.00000000,
-3.00000000
155+1.00000000,
+0.00000000,
+0.00000000
156-9.00000000,
+1.00000000,
+9.00000000
158 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
161+1.00000000,
+3.00000000,
-7.00000000
162+1.00000000,
+9.00000000,
-4.00000000
163-3.00000000,
-3.00000000,
-3.00000000
164+1.00000000,
+0.00000000,
+0.00000000
165-9.00000000,
+1.00000000,
+9.00000000
168+1.00000000,
+3.00000000,
-7.00000000
169+1.00000000,
+9.00000000,
-4.00000000
170-3.00000000,
-3.00000000,
-3.00000000
171+1.00000000,
+0.00000000,
+0.00000000
172-9.00000000,
+1.00000000,
+9.00000000
174+1.00000000,
+3.00000000,
-7.00000000
175+1.00000000,
+9.00000000,
-4.00000000
176-3.00000000,
-3.00000000,
-3.00000000
177+1.00000000,
+0.00000000,
+0.00000000
178-9.00000000,
+1.00000000,
+9.00000000
185+2.00000000,
-3.00000000
188+7.00000000,
-4.00000000,
+0.00000000,
+1.00000000,
-8.00000000
189-9.00000000,
+9.00000000,
-2.00000000,
+4.00000000,
-6.00000000
194+1.00000000,
+0.319688737
195+0.319688737,
+1.00000000
196call setAffinity(affinity, sample, dim, tform, genrecmat)
198+4.12280130,
-6.76217890
199-1.12280130,
+7.72124481
200-0.639377475,
-2.00000000
201+2.27875495,
+4.31968880
202-9.91813278,
-8.55751038
204 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
207+4.12280130,
-6.76217890
208-1.12280130,
+7.72124481
209-0.639377475,
-2.00000000
210+2.27875495,
+4.31968880
211-9.91813278,
-8.55751038
214+6.99999952,
-9.00000000
215-3.99999952,
+9.00000000
216+0.00000000,
-1.99999988
217+0.999999881,
+3.99999976
218-8.00000000,
-6.00000000
220+7.00000000,
-9.00000000
221-4.00000000,
+9.00000000
222+0.00000000,
-2.00000000
223+1.00000000,
+4.00000000
224-8.00000000,
-6.00000000
226call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
228+1.00000000,
+0.319688737
229+0.00000000,
+1.00000000
230call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
232+4.12280130,
-9.00000000
233-1.12280130,
+9.00000000
234-0.639377475,
-2.00000000
235+2.27875495,
+4.00000000
236-9.91813278,
-6.00000000
238 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
241+4.12280130,
-9.00000000
242-1.12280130,
+9.00000000
243-0.639377475,
-2.00000000
244+2.27875495,
+4.00000000
245-9.91813278,
-6.00000000
248+7.00000000,
-9.00000000
249-4.00000000,
+9.00000000
250+0.00000000,
-2.00000000
251+1.00000000,
+4.00000000
252-8.00000000,
-6.00000000
254+7.00000000,
-9.00000000
255-4.00000000,
+9.00000000
256+0.00000000,
-2.00000000
257+1.00000000,
+4.00000000
258-8.00000000,
-6.00000000
260call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
262+4.12280130,
-9.00000000
263-1.12280130,
+9.00000000
264-0.639377475,
-2.00000000
265+2.27875495,
+4.00000000
266-9.91813278,
-6.00000000
268 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
271+4.12280130,
-9.00000000
272-1.12280130,
+9.00000000
273-0.639377475,
-2.00000000
274+2.27875495,
+4.00000000
275-9.91813278,
-6.00000000
278+7.00000000,
-9.00000000
279-4.00000000,
+9.00000000
280+0.00000000,
-2.00000000
281+1.00000000,
+4.00000000
282-8.00000000,
-6.00000000
284+7.00000000,
-9.00000000
285-4.00000000,
+9.00000000
286+0.00000000,
-2.00000000
287+1.00000000,
+4.00000000
288-8.00000000,
-6.00000000
290call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
292+7.00000000,
-9.00000000
293-4.00000000,
+9.00000000
294+0.00000000,
-2.00000000
295+1.00000000,
+4.00000000
296-8.00000000,
-6.00000000
298 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
301+7.00000000,
-9.00000000
302-4.00000000,
+9.00000000
303+0.00000000,
-2.00000000
304+1.00000000,
+4.00000000
305-8.00000000,
-6.00000000
308+7.00000000,
-9.00000000
309-4.00000000,
+9.00000000
310+0.00000000,
-2.00000000
311+1.00000000,
+4.00000000
312-8.00000000,
-6.00000000
314+7.00000000,
-9.00000000
315-4.00000000,
+9.00000000
316+0.00000000,
-2.00000000
317+1.00000000,
+4.00000000
318-8.00000000,
-6.00000000
320call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
322+7.00000000,
-9.00000000
323-4.00000000,
+9.00000000
324+0.00000000,
-2.00000000
325+1.00000000,
+4.00000000
326-8.00000000,
-6.00000000
328 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
331+7.00000000,
-9.00000000
332-4.00000000,
+9.00000000
333+0.00000000,
-2.00000000
334+1.00000000,
+4.00000000
335-8.00000000,
-6.00000000
338+7.00000000,
-9.00000000
339-4.00000000,
+9.00000000
340+0.00000000,
-2.00000000
341+1.00000000,
+4.00000000
342-8.00000000,
-6.00000000
344+7.00000000,
-9.00000000
345-4.00000000,
+9.00000000
346+0.00000000,
-2.00000000
347+1.00000000,
+4.00000000
348-8.00000000,
-6.00000000
355-5.00000000,
-2.00000000,
-1.00000000
358+7.00000000,
-5.00000000,
+9.00000000,
+7.00000000
359-2.00000000,
-3.00000000,
+0.00000000,
+1.00000000
360+9.00000000,
+1.00000000,
+2.00000000,
+1.00000000
365+1.00000000,
+0.737021267,
-0.490478843
366+0.737021267,
+1.00000000,
+0.224671006
367-0.490478843,
+0.224671006,
+1.00000000
368call setAffinity(affinity, sample, dim, tform, genrecmat)
370+1.11164808,
+5.18118763,
+5.11730576
371-7.70154285,
-6.46043539,
+2.77838135
372+8.01904202,
+7.08253384,
-2.41430950
373+7.24654245,
+6.38381958,
-2.20868111
375 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
378+1.11164808,
+5.18118763,
+5.11730576
379-7.70154285,
-6.46043539,
+2.77838135
380+8.01904202,
+7.08253384,
-2.41430950
381+7.24654245,
+6.38381958,
-2.20868111
384+7.00006104,
-2.00000000,
+9.00000000
385-4.99987793,
-3.00024414,
+0.999969482
386+8.99954224,
+0.335693359E-3,
+1.99993896
387+6.99981689,
+1.00012207,
+1.00006104
389+7.00000000,
-2.00000000,
+9.00000000
390-5.00000000,
-3.00000000,
+1.00000000
391+9.00000000,
+0.00000000,
+2.00000000
392+7.00000000,
+1.00000000,
+1.00000000
394call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
396+1.00000000,
+0.737021267,
-0.490478843
397+0.00000000,
+1.00000000,
+0.224671006
398+0.00000000,
+0.00000000,
+1.00000000
399call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
401+1.11164808,
+0.220389366E-1,
+9.00000000
402-7.70154285,
-2.77532911,
+1.00000000
403+8.01904202,
+0.449342012,
+2.00000000
404+7.24654245,
+1.22467101,
+1.00000000
406 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
409+1.11164808,
+0.220389366E-1,
+9.00000000
410-7.70154285,
-2.77532911,
+1.00000000
411+8.01904202,
+0.449342012,
+2.00000000
412+7.24654245,
+1.22467101,
+1.00000000
415+7.00000048,
-2.00000000,
+9.00000000
416-5.00000000,
-3.00000000,
+1.00000000
417+9.00000000,
+0.00000000,
+2.00000000
418+7.00000000,
+1.00000000,
+1.00000000
420+7.00000000,
-2.00000000,
+9.00000000
421-5.00000000,
-3.00000000,
+1.00000000
422+9.00000000,
+0.00000000,
+2.00000000
423+7.00000000,
+1.00000000,
+1.00000000
425call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
427+1.11164808,
+0.220389366E-1,
+9.00000000
428-7.70154285,
-2.77532911,
+1.00000000
429+8.01904202,
+0.449342012,
+2.00000000
430+7.24654245,
+1.22467101,
+1.00000000
432 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
435+1.11164808,
+0.220389366E-1,
+9.00000000
436-7.70154285,
-2.77532911,
+1.00000000
437+8.01904202,
+0.449342012,
+2.00000000
438+7.24654245,
+1.22467101,
+1.00000000
441+7.00000048,
-2.00000000,
+9.00000000
442-5.00000000,
-3.00000000,
+1.00000000
443+9.00000000,
+0.00000000,
+2.00000000
444+7.00000000,
+1.00000000,
+1.00000000
446+7.00000000,
-2.00000000,
+9.00000000
447-5.00000000,
-3.00000000,
+1.00000000
448+9.00000000,
+0.00000000,
+2.00000000
449+7.00000000,
+1.00000000,
+1.00000000
451call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
453+7.00000000,
-2.00000000,
+9.00000000
454-5.00000000,
-3.00000000,
+1.00000000
455+9.00000000,
+0.00000000,
+2.00000000
456+7.00000000,
+1.00000000,
+1.00000000
458 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
461+7.00000000,
-2.00000000,
+9.00000000
462-5.00000000,
-3.00000000,
+1.00000000
463+9.00000000,
+0.00000000,
+2.00000000
464+7.00000000,
+1.00000000,
+1.00000000
467+7.00000000,
-2.00000000,
+9.00000000
468-5.00000000,
-3.00000000,
+1.00000000
469+9.00000000,
+0.00000000,
+2.00000000
470+7.00000000,
+1.00000000,
+1.00000000
472+7.00000000,
-2.00000000,
+9.00000000
473-5.00000000,
-3.00000000,
+1.00000000
474+9.00000000,
+0.00000000,
+2.00000000
475+7.00000000,
+1.00000000,
+1.00000000
477call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
479+7.00000000,
-2.00000000,
+9.00000000
480-5.00000000,
-3.00000000,
+1.00000000
481+9.00000000,
+0.00000000,
+2.00000000
482+7.00000000,
+1.00000000,
+1.00000000
484 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
487+7.00000000,
-2.00000000,
+9.00000000
488-5.00000000,
-3.00000000,
+1.00000000
489+9.00000000,
+0.00000000,
+2.00000000
490+7.00000000,
+1.00000000,
+1.00000000
493+7.00000000,
-2.00000000,
+9.00000000
494-5.00000000,
-3.00000000,
+1.00000000
495+9.00000000,
+0.00000000,
+2.00000000
496+7.00000000,
+1.00000000,
+1.00000000
498+7.00000000,
-2.00000000,
+9.00000000
499-5.00000000,
-3.00000000,
+1.00000000
500+9.00000000,
+0.00000000,
+2.00000000
501+7.00000000,
+1.00000000,
+1.00000000
508-2.00000000,
-4.00000000,
+3.00000000
511-3.00000000,
-8.00000000,
-3.00000000,
+5.00000000
512+4.00000000,
+5.00000000,
+4.00000000,
-7.00000000
513-1.00000000,
-9.00000000,
+5.00000000,
+3.00000000
518+1.00000000,
+0.447494864,
-0.525756478
519+0.447494864,
+1.00000000,
-0.206653520
520-0.525756478,
-0.206653520,
+1.00000000
521call setAffinity(affinity, sample, dim, tform, genrecmat)
523-0.684264064,
+2.86416912,
-0.249344647
524-1.03071737,
+3.27992272,
-5.82721567
525-3.83880281,
+1.62424791,
+5.75065517
526+0.290266633,
-5.38248634,
+1.81779242
528 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
531-0.684264064,
+2.86416912,
-0.249344647
532-1.03071737,
+3.27992272,
-5.82721567
533-3.83880281,
+1.62424791,
+5.75065517
534+0.290266633,
-5.38248634,
+1.81779242
537-3.00000000,
+3.99999976,
-1.00000000
538-7.99999952,
+4.99999905,
-9.00000000
539-3.00000000,
+3.99999976,
+4.99999952
540+5.00000000,
-6.99999952,
+3.00000048
542-3.00000000,
+4.00000000,
-1.00000000
543-8.00000000,
+5.00000000,
-9.00000000
544-3.00000000,
+4.00000000,
+5.00000000
545+5.00000000,
-7.00000000,
+3.00000000
547call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
549+1.00000000,
+0.447494864,
-0.525756478
550+0.00000000,
+1.00000000,
-0.206653520
551+0.00000000,
+0.00000000,
+1.00000000
552call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
554-0.684264183,
+4.20665359,
-1.00000000
555-1.03071737,
+6.85988140,
-9.00000000
556-3.83880281,
+2.96673250,
+5.00000000
557+0.290266514,
-7.61996078,
+3.00000000
559 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
562-0.684264183,
+4.20665359,
-1.00000000
563-1.03071737,
+6.85988140,
-9.00000000
564-3.83880281,
+2.96673250,
+5.00000000
565+0.290266514,
-7.61996078,
+3.00000000
568-3.00000000,
+4.00000000,
-1.00000000
569-7.99999952,
+5.00000000,
-9.00000000
570-3.00000000,
+4.00000000,
+5.00000000
571+5.00000000,
-7.00000000,
+3.00000000
573-3.00000000,
+4.00000000,
-1.00000000
574-8.00000000,
+5.00000000,
-9.00000000
575-3.00000000,
+4.00000000,
+5.00000000
576+5.00000000,
-7.00000000,
+3.00000000
578call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
580-0.684264183,
+4.20665359,
-1.00000000
581-1.03071737,
+6.85988140,
-9.00000000
582-3.83880281,
+2.96673250,
+5.00000000
583+0.290266514,
-7.61996078,
+3.00000000
585 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
588-0.684264183,
+4.20665359,
-1.00000000
589-1.03071737,
+6.85988140,
-9.00000000
590-3.83880281,
+2.96673250,
+5.00000000
591+0.290266514,
-7.61996078,
+3.00000000
594-3.00000000,
+4.00000000,
-1.00000000
595-7.99999952,
+5.00000000,
-9.00000000
596-3.00000000,
+4.00000000,
+5.00000000
597+5.00000000,
-7.00000000,
+3.00000000
599-3.00000000,
+4.00000000,
-1.00000000
600-8.00000000,
+5.00000000,
-9.00000000
601-3.00000000,
+4.00000000,
+5.00000000
602+5.00000000,
-7.00000000,
+3.00000000
604call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
606-3.00000000,
+4.00000000,
-1.00000000
607-8.00000000,
+5.00000000,
-9.00000000
608-3.00000000,
+4.00000000,
+5.00000000
609+5.00000000,
-7.00000000,
+3.00000000
611 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
614-3.00000000,
+4.00000000,
-1.00000000
615-8.00000000,
+5.00000000,
-9.00000000
616-3.00000000,
+4.00000000,
+5.00000000
617+5.00000000,
-7.00000000,
+3.00000000
620-3.00000000,
+4.00000000,
-1.00000000
621-8.00000000,
+5.00000000,
-9.00000000
622-3.00000000,
+4.00000000,
+5.00000000
623+5.00000000,
-7.00000000,
+3.00000000
625-3.00000000,
+4.00000000,
-1.00000000
626-8.00000000,
+5.00000000,
-9.00000000
627-3.00000000,
+4.00000000,
+5.00000000
628+5.00000000,
-7.00000000,
+3.00000000
630call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
632-3.00000000,
+4.00000000,
-1.00000000
633-8.00000000,
+5.00000000,
-9.00000000
634-3.00000000,
+4.00000000,
+5.00000000
635+5.00000000,
-7.00000000,
+3.00000000
637 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
640-3.00000000,
+4.00000000,
-1.00000000
641-8.00000000,
+5.00000000,
-9.00000000
642-3.00000000,
+4.00000000,
+5.00000000
643+5.00000000,
-7.00000000,
+3.00000000
646-3.00000000,
+4.00000000,
-1.00000000
647-8.00000000,
+5.00000000,
-9.00000000
648-3.00000000,
+4.00000000,
+5.00000000
649+5.00000000,
-7.00000000,
+3.00000000
651-3.00000000,
+4.00000000,
-1.00000000
652-8.00000000,
+5.00000000,
-9.00000000
653-3.00000000,
+4.00000000,
+5.00000000
654+5.00000000,
-7.00000000,
+3.00000000
661-2.00000000,
-2.00000000
664+5.00000000,
-4.00000000,
+9.00000000,
+6.00000000
665+9.00000000,
-2.00000000,
+4.00000000,
+6.00000000
670+1.00000000,
-0.981980860
671-0.981980860,
+1.00000000
672call setAffinity(affinity, sample, dim, tform, genrecmat)
674-3.83782768,
+4.09009552
675-2.03603840,
+1.92792344
676+5.07207680,
-4.83782768
677+0.108114719,
+0.108114719
679 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
682-3.83782768,
+4.09009552
683-2.03603840,
+1.92792344
684+5.07207680,
-4.83782768
685+0.108114719,
+0.108114719
688+5.00000000,
+8.99999237
689-3.99999619,
-1.99999619
690+8.99998474,
+4.00000000
691+5.99998951,
+5.99998951
693+5.00000000,
+9.00000000
694-4.00000000,
-2.00000000
695+9.00000000,
+4.00000000
696+6.00000000,
+6.00000000
698call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
700+1.00000000,
-0.981980860
701+0.00000000,
+1.00000000
702call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
704-3.83782768,
+9.00000000
705-2.03603840,
-2.00000000
706+5.07207680,
+4.00000000
707+0.108114719,
+6.00000000
709 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
712-3.83782768,
+9.00000000
713-2.03603840,
-2.00000000
714+5.07207680,
+4.00000000
715+0.108114719,
+6.00000000
718+5.00000000,
+9.00000000
719-4.00000000,
-2.00000000
720+9.00000000,
+4.00000000
721+6.00000000,
+6.00000000
723+5.00000000,
+9.00000000
724-4.00000000,
-2.00000000
725+9.00000000,
+4.00000000
726+6.00000000,
+6.00000000
728call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
730-3.83782768,
+9.00000000
731-2.03603840,
-2.00000000
732+5.07207680,
+4.00000000
733+0.108114719,
+6.00000000
735 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
738-3.83782768,
+9.00000000
739-2.03603840,
-2.00000000
740+5.07207680,
+4.00000000
741+0.108114719,
+6.00000000
744+5.00000000,
+9.00000000
745-4.00000000,
-2.00000000
746+9.00000000,
+4.00000000
747+6.00000000,
+6.00000000
749+5.00000000,
+9.00000000
750-4.00000000,
-2.00000000
751+9.00000000,
+4.00000000
752+6.00000000,
+6.00000000
754call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
756+5.00000000,
+9.00000000
757-4.00000000,
-2.00000000
758+9.00000000,
+4.00000000
759+6.00000000,
+6.00000000
761 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
764+5.00000000,
+9.00000000
765-4.00000000,
-2.00000000
766+9.00000000,
+4.00000000
767+6.00000000,
+6.00000000
770+5.00000000,
+9.00000000
771-4.00000000,
-2.00000000
772+9.00000000,
+4.00000000
773+6.00000000,
+6.00000000
775+5.00000000,
+9.00000000
776-4.00000000,
-2.00000000
777+9.00000000,
+4.00000000
778+6.00000000,
+6.00000000
780call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
782+5.00000000,
+9.00000000
783-4.00000000,
-2.00000000
784+9.00000000,
+4.00000000
785+6.00000000,
+6.00000000
787 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
790+5.00000000,
+9.00000000
791-4.00000000,
-2.00000000
792+9.00000000,
+4.00000000
793+6.00000000,
+6.00000000
796+5.00000000,
+9.00000000
797-4.00000000,
-2.00000000
798+9.00000000,
+4.00000000
799+6.00000000,
+6.00000000
801+5.00000000,
+9.00000000
802-4.00000000,
-2.00000000
803+9.00000000,
+4.00000000
804+6.00000000,
+6.00000000
816+5.00000000,
-5.00000000
819-4.00000000,
+7.00000000
820+2.00000000,
+8.00000000
821-4.00000000,
+4.00000000
822-1.00000000,
+5.00000000
823-8.00000000,
-6.00000000
828+1.00000000,
+0.165688619E-1
829+0.165688619E-1,
+1.00000000
832-3.88401794,
+6.93372440
833+2.13255095,
+8.03313732
834-3.93372464,
+3.93372464
835-0.917155683,
+4.98343134
836-8.09941292,
-6.13255072
838 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
841-3.88401794,
+6.93372440
842+2.13255095,
+8.03313732
843-3.93372464,
+3.93372464
844-0.917155683,
+4.98343134
845-8.09941292,
-6.13255072
848-4.00000000,
+7.00000048
849+2.00000024,
+8.00000000
850-4.00000048,
+4.00000048
851-1.00000000,
+5.00000048
852-8.00000095,
-6.00000000
854-4.00000000,
+7.00000000
855+2.00000000,
+8.00000000
856-4.00000000,
+4.00000000
857-1.00000000,
+5.00000000
858-8.00000000,
-6.00000000
860call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
862+1.00000000,
+0.165688619E-1
863+0.00000000,
+1.00000000
864call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
866-3.88401794,
+7.00000000
867+2.13255095,
+8.00000000
868-3.93372464,
+4.00000000
869-0.917155683,
+5.00000000
870-8.09941292,
-6.00000000
872 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
875-3.88401794,
+7.00000000
876+2.13255095,
+8.00000000
877-3.93372464,
+4.00000000
878-0.917155683,
+5.00000000
879-8.09941292,
-6.00000000
882-4.00000000,
+7.00000000
883+2.00000000,
+8.00000000
884-4.00000000,
+4.00000000
885-1.00000000,
+5.00000000
886-7.99999952,
-6.00000000
888-4.00000000,
+7.00000000
889+2.00000000,
+8.00000000
890-4.00000000,
+4.00000000
891-1.00000000,
+5.00000000
892-8.00000000,
-6.00000000
894call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
896-3.88401794,
+7.00000000
897+2.13255095,
+8.00000000
898-3.93372464,
+4.00000000
899-0.917155683,
+5.00000000
900-8.09941292,
-6.00000000
902 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
905-3.88401794,
+7.00000000
906+2.13255095,
+8.00000000
907-3.93372464,
+4.00000000
908-0.917155683,
+5.00000000
909-8.09941292,
-6.00000000
912-4.00000000,
+7.00000000
913+2.00000000,
+8.00000000
914-4.00000000,
+4.00000000
915-1.00000000,
+5.00000000
916-7.99999952,
-6.00000000
918-4.00000000,
+7.00000000
919+2.00000000,
+8.00000000
920-4.00000000,
+4.00000000
921-1.00000000,
+5.00000000
922-8.00000000,
-6.00000000
924call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
926-4.00000000,
+7.00000000
927+2.00000000,
+8.00000000
928-4.00000000,
+4.00000000
929-1.00000000,
+5.00000000
930-8.00000000,
-6.00000000
932 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
935-4.00000000,
+7.00000000
936+2.00000000,
+8.00000000
937-4.00000000,
+4.00000000
938-1.00000000,
+5.00000000
939-8.00000000,
-6.00000000
942-4.00000000,
+7.00000000
943+2.00000000,
+8.00000000
944-4.00000000,
+4.00000000
945-1.00000000,
+5.00000000
946-8.00000000,
-6.00000000
948-4.00000000,
+7.00000000
949+2.00000000,
+8.00000000
950-4.00000000,
+4.00000000
951-1.00000000,
+5.00000000
952-8.00000000,
-6.00000000
954call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
956-4.00000000,
+7.00000000
957+2.00000000,
+8.00000000
958-4.00000000,
+4.00000000
959-1.00000000,
+5.00000000
960-8.00000000,
-6.00000000
962 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)
965-4.00000000,
+7.00000000
966+2.00000000,
+8.00000000
967-4.00000000,
+4.00000000
968-1.00000000,
+5.00000000
969-8.00000000,
-6.00000000
972-4.00000000,
+7.00000000
973+2.00000000,
+8.00000000
974-4.00000000,
+4.00000000
975-1.00000000,
+5.00000000
976-8.00000000,
-6.00000000
978-4.00000000,
+7.00000000
979+2.00000000,
+8.00000000
980-4.00000000,
+4.00000000
981-1.00000000,
+5.00000000
982-8.00000000,
-6.00000000
989-1.00000000,
+1.00000000
992+2.00000000,
+1.00000000
993-1.00000000,
-1.00000000
994-1.00000000,
+6.00000000
995+5.00000000,
-4.00000000
1000+1.00000000,
-0.720322132
1001-0.720322132,
+1.00000000
1004+1.27967787,
-0.440644264
1005-0.279677868,
-0.279677868
1006-5.32193279,
+6.72032213
1007+7.88128853,
-7.60161066
1009 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
1012+1.27967787,
-0.440644264
1013-0.279677868,
-0.279677868
1014-5.32193279,
+6.72032213
1015+7.88128853,
-7.60161066
1018+2.00000000,
+1.00000024
1019-1.00000000,
-1.00000000
1020-1.00000095,
+5.99999905
1021+5.00000191,
-3.99999809
1023+2.00000000,
+1.00000000
1024-1.00000000,
-1.00000000