18 type(display_type) :: disp
19 integer(IK) :: itry, ntry
= 5
20 integer(IK) :: dim, isam, ndim, nsam
24 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
25 call disp%show(
"!Transform 2D sample along the second dimension.")
26 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
31 real(RKG),
allocatable :: tlate(:), tform(:,:), sample(:,:), affinInv(:,:), affinity(:,:)
34 call disp%show(
"dim = 2; ndim = getUnifRand(1, 3); nsam = getUnifRand(4, 5)")
38 call disp%show(
"tlate = getUnifRand(-5, +5, ndim)")
42 call disp%show(
"sample = getUnifRand(-9, +9, ndim, nsam)")
46 call disp%show(
"call setResized(affinity, shape(sample, IK))")
48 call disp%show(
"call setResized(affinInv, shape(sample, IK))")
50 call disp%show(
"tform = getCovRand(1., ndim)")
54 call disp%show(
"affinity = getAffinity(sample, dim, tform) ! whole sample transformation.")
56 call disp%show(
"transpose(affinity)")
57 call disp%show(
transpose(affinity) )
59 call disp%show(
" affinity(:,isam) = getAffinity(sample(:,isam), tform) ! point-wise transformation")
62 affinity(:,isam)
= getAffinity(sample(:,isam), tform)
64 call disp%show(
"transpose(affinity)")
65 call disp%show(
transpose(affinity) )
66 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform)) ! inverse transformation.")
68 call disp%show(
"transpose(affinInv)")
69 call disp%show(
transpose(affinInv) )
70 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
73 call disp%show(
"call setMatInit(tform(2:,1:ndim-1), lowDia, vlow = 0._RKG, vdia = 0._RKG) ! make tform an upper-diagonal matrix.")
74 call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
77 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = upperDiag) ! whole sample transformation.")
78 affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
79 call disp%show(
"transpose(affinity)")
80 call disp%show(
transpose(affinity) )
82 call disp%show(
" affinity(:,isam) = getAffinity(sample(:,isam), tform, upperDiag) ! point-wise transformation")
85 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperDiag)
87 call disp%show(
"transpose(affinity)")
88 call disp%show(
transpose(affinity) )
89 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform, upperDiag), class = upperDiag) ! inverse transformation.")
91 call disp%show(
"transpose(affinInv)")
92 call disp%show(
transpose(affinInv) )
93 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
96 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = upperUnit) ! whole sample transformation.")
97 affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
98 call disp%show(
"transpose(affinity)")
99 call disp%show(
transpose(affinity) )
101 call disp%show(
" affinity(:,isam) = getAffinity(sample(:,isam), tform, upperUnit) ! point-wise transformation")
104 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperUnit)
106 call disp%show(
"transpose(affinity)")
107 call disp%show(
transpose(affinity) )
108 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform), class = upperUnit) ! inverse transformation.")
110 call disp%show(
"transpose(affinInv)")
111 call disp%show(
transpose(affinInv) )
112 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
115 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = lowerDiag) ! whole sample transformation.")
116 affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
117 call disp%show(
"transpose(affinity)")
118 call disp%show(
transpose(affinity) )
120 call disp%show(
" affinity(:,isam) = getAffinity(sample(:,isam), tform, lowerDiag) ! point-wise transformation")
123 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerDiag)
125 call disp%show(
"transpose(affinity)")
126 call disp%show(
transpose(affinity) )
127 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform), class = lowerDiag) ! inverse transformation.")
129 call disp%show(
"transpose(affinInv)")
130 call disp%show(
transpose(affinInv) )
131 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
134 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = lowerUnit) ! whole sample transformation.")
135 affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
136 call disp%show(
"transpose(affinity)")
137 call disp%show(
transpose(affinity) )
139 call disp%show(
" affinity(:,isam) = getAffinity(sample(:,isam), tform, lowerUnit) ! point-wise transformation")
142 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerUnit)
144 call disp%show(
"transpose(affinity)")
145 call disp%show(
transpose(affinity) )
146 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform), class = lowerUnit) ! inverse transformation.")
148 call disp%show(
"transpose(affinInv)")
149 call disp%show(
transpose(affinInv) )
150 call disp%show(
"transpose(sample) ! for comparison with affinInv.")
156 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
157 call disp%show(
"!Transform 2D sample along the first dimension.")
158 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
163 real(RKG),
allocatable :: tlate(:), tform(:,:), sample(:,:), affinInv(:,:), affinity(:,:)
166 call disp%show(
"dim = 1; ndim = getUnifRand(1, 3); nsam = getUnifRand(4, 5)")
170 call disp%show(
"tlate = getUnifRand(-5, +5, ndim)")
174 call disp%show(
"sample = getUnifRand(-9, +9, nsam, ndim)")
178 call disp%show(
"call setResized(affinity, shape(sample, IK))")
180 call disp%show(
"call setResized(affinInv, shape(sample, IK))")
182 call disp%show(
"tform = getCovRand(1., ndim)")
186 call disp%show(
"affinity = getAffinity(sample, dim, tform) ! whole sample transformation.")
191 call disp%show(
" affinity(isam,:) = getAffinity(sample(isam,:), tform) ! point-wise transformation")
194 affinity(isam,:)
= getAffinity(sample(isam,:), tform)
198 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform)) ! inverse transformation.")
202 call disp%show(
"sample ! for comparison with affinInv.")
205 call disp%show(
"call setMatInit(tform(2:,1:ndim-1), lowDia, vlow = 0._RKG, vdia = 0._RKG) ! make tform an upper-diagonal matrix.")
206 call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
209 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = upperDiag) ! whole sample transformation.")
210 affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
214 call disp%show(
" affinity(isam,:) = getAffinity(sample(isam,:), tform, upperDiag) ! point-wise transformation")
217 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperDiag)
221 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform, upperDiag), class = upperDiag) ! inverse transformation.")
225 call disp%show(
"sample ! for comparison with affinInv.")
228 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = upperUnit) ! whole sample transformation.")
229 affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
233 call disp%show(
" affinity(isam,:) = getAffinity(sample(isam,:), tform, upperUnit) ! point-wise transformation")
236 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperUnit)
240 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform), class = upperUnit) ! inverse transformation.")
244 call disp%show(
"sample ! for comparison with affinInv.")
247 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = lowerDiag) ! whole sample transformation.")
248 affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
252 call disp%show(
" affinity(isam,:) = getAffinity(sample(isam,:), tform, lowerDiag) ! point-wise transformation")
255 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerDiag)
259 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform), class = lowerDiag) ! inverse transformation.")
263 call disp%show(
"sample ! for comparison with affinInv.")
266 call disp%show(
"affinity = getAffinity(sample, dim, tform, class = lowerUnit) ! whole sample transformation.")
267 affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
271 call disp%show(
" affinity(isam,:) = getAffinity(sample(isam,:), tform, lowerUnit) ! point-wise transformation")
274 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerUnit)
278 call disp%show(
"affinInv = getAffinity(affinity, dim, getMatInv(tform), class = lowerUnit) ! inverse transformation.")
282 call disp%show(
"sample ! for comparison with affinInv.")
303 type(css_type),
allocatable :: shapes(:), tforms(:)
304 integer(IK) :: ishape, itform
305 integer(IK),
parameter :: nsam
= 1000, ndim
= 2, dim
= 2
306 real(RKG),
allocatable :: angle, tlate(:), cov(:,:), tform(:,:), sample(:,:), affinity(:,:), affinInv(:,:)
310 do ishape
= 1,
size(shapes)
311 if (shapes(ishape)
%val
== "circle")
then
313 elseif (shapes(ishape)
%val
== "square")
then
316 error stop "Unrecognized shape."
318 do itform
= 1,
size(tforms)
319 if (tforms(itform)
%val
== "rotation")
then
321 tform
= reshape([
cos(angle),
-sin(angle),
sin(angle),
cos(angle)], [ndim, ndim])
322 affinity
= getAffinity(sample, dim, tform, genrecmat, tlate)
323 elseif (tforms(itform)
%val
== "warp")
then
324 cov
= reshape([
1., .
5, .
5,
1.], shape
= [ndim, ndim])
326 affinity
= getAffinity(sample, dim, tform, lowerDiag, tlate)
328 error stop "Unrecognized tform."
331 if (
0 /= getErrTableWrite(SK_
"getAffinity."//tforms(itform)
%val
//SK_
"."//shapes(ishape)
%val
//SK_
".sample.txt", sample, trans))
error stop "Failed table-write."
332 if (
0 /= getErrTableWrite(SK_
"getAffinity."//tforms(itform)
%val
//SK_
"."//shapes(ishape)
%val
//SK_
".affinity.txt", affinity, trans))
error stop "Failed table-write."
333 if (
0 /= getErrTableWrite(SK_
"getAffinity."//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.
15+0.00000000,
+9.00000000,
+0.00000000,
+9.00000000,
+5.00000000
29 affinity(:,isam)
= getAffinity(sample(:,isam), tform)
51call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
54affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
62 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperDiag)
84affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
92 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperUnit)
114affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
122 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerDiag)
144affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
152 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerUnit)
179+2.00000000,
-2.00000000,
+4.00000000
182+1.00000000,
+4.00000000,
+0.00000000,
-4.00000000
183+8.00000000,
-4.00000000,
-5.00000000,
+6.00000000
184+1.00000000,
+7.00000000,
-2.00000000,
-9.00000000
189+1.00000000,
-0.683877289,
+0.683055997
190-0.683877289,
+1.00000000,
-0.219324827
191+0.683055997,
-0.219324827,
+1.00000000
194-3.78796244,
+7.09679794,
-0.715426207E-1
195+11.5169010,
-8.27078247,
+10.6095238
196+2.05327439,
-4.56135035,
-0.903375864
197-14.2507677,
+10.7094326,
-13.0481730
199 affinity(:,isam)
= getAffinity(sample(:,isam), tform)
202-3.78796244,
+7.09679794,
-0.715426207E-1
203+11.5169010,
-8.27078247,
+10.6095238
204+2.05327439,
-4.56135035,
-0.903375864
205-14.2507677,
+10.7094326,
-13.0481730
208+0.999998987,
+7.99999762,
+0.999999881
209+4.00000000,
-3.99999809,
+7.00000191
210+0.953674316E-6,
-5.00000000,
-2.00000000
211-3.99999809,
+5.99999523,
-9.00000000
213+1.00000000,
+8.00000000,
+1.00000000
214+4.00000000,
-4.00000000,
+7.00000000
215+0.00000000,
-5.00000000,
-2.00000000
216-4.00000000,
+6.00000000,
-9.00000000
218call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
220+1.00000000,
-0.683877289,
+0.683055997
221+0.00000000,
+1.00000000,
-0.219324827
222+0.00000000,
+0.00000000,
+1.00000000
223affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
225-3.78796244,
+7.78067493,
+1.00000000
226+11.5169010,
-5.53527355,
+7.00000000
227+2.05327439,
-4.56135035,
-2.00000000
228-14.2507677,
+7.97392368,
-9.00000000
230 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperDiag)
233-3.78796244,
+7.78067493,
+1.00000000
234+11.5169010,
-5.53527355,
+7.00000000
235+2.05327439,
-4.56135035,
-2.00000000
236-14.2507677,
+7.97392368,
-9.00000000
239+0.999999523,
+8.00000000,
+1.00000000
240+4.00000000,
-3.99999976,
+7.00000000
241+0.00000000,
-5.00000000,
-2.00000000
242-4.00000000,
+6.00000000,
-9.00000000
244+1.00000000,
+8.00000000,
+1.00000000
245+4.00000000,
-4.00000000,
+7.00000000
246+0.00000000,
-5.00000000,
-2.00000000
247-4.00000000,
+6.00000000,
-9.00000000
249affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
251-3.78796244,
+7.78067493,
+1.00000000
252+11.5169010,
-5.53527355,
+7.00000000
253+2.05327439,
-4.56135035,
-2.00000000
254-14.2507677,
+7.97392368,
-9.00000000
256 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperUnit)
259-3.78796244,
+7.78067493,
+1.00000000
260+11.5169010,
-5.53527355,
+7.00000000
261+2.05327439,
-4.56135035,
-2.00000000
262-14.2507677,
+7.97392368,
-9.00000000
265+0.999999523,
+8.00000000,
+1.00000000
266+4.00000000,
-3.99999976,
+7.00000000
267+0.00000000,
-5.00000000,
-2.00000000
268-4.00000000,
+6.00000000,
-9.00000000
270+1.00000000,
+8.00000000,
+1.00000000
271+4.00000000,
-4.00000000,
+7.00000000
272+0.00000000,
-5.00000000,
-2.00000000
273-4.00000000,
+6.00000000,
-9.00000000
275affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
277+1.00000000,
+8.00000000,
+1.00000000
278+4.00000000,
-4.00000000,
+7.00000000
279+0.00000000,
-5.00000000,
-2.00000000
280-4.00000000,
+6.00000000,
-9.00000000
282 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerDiag)
285+1.00000000,
+8.00000000,
+1.00000000
286+4.00000000,
-4.00000000,
+7.00000000
287+0.00000000,
-5.00000000,
-2.00000000
288-4.00000000,
+6.00000000,
-9.00000000
291+1.00000000,
+8.00000000,
+1.00000000
292+4.00000000,
-4.00000000,
+7.00000000
293+0.00000000,
-5.00000000,
-2.00000000
294-4.00000000,
+6.00000000,
-9.00000000
296+1.00000000,
+8.00000000,
+1.00000000
297+4.00000000,
-4.00000000,
+7.00000000
298+0.00000000,
-5.00000000,
-2.00000000
299-4.00000000,
+6.00000000,
-9.00000000
301affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
303+1.00000000,
+8.00000000,
+1.00000000
304+4.00000000,
-4.00000000,
+7.00000000
305+0.00000000,
-5.00000000,
-2.00000000
306-4.00000000,
+6.00000000,
-9.00000000
308 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerUnit)
311+1.00000000,
+8.00000000,
+1.00000000
312+4.00000000,
-4.00000000,
+7.00000000
313+0.00000000,
-5.00000000,
-2.00000000
314-4.00000000,
+6.00000000,
-9.00000000
317+1.00000000,
+8.00000000,
+1.00000000
318+4.00000000,
-4.00000000,
+7.00000000
319+0.00000000,
-5.00000000,
-2.00000000
320-4.00000000,
+6.00000000,
-9.00000000
322+1.00000000,
+8.00000000,
+1.00000000
323+4.00000000,
-4.00000000,
+7.00000000
324+0.00000000,
-5.00000000,
-2.00000000
325-4.00000000,
+6.00000000,
-9.00000000
335-9.00000000,
+8.00000000,
+5.00000000,
-4.00000000,
-6.00000000
349 affinity(:,isam)
= getAffinity(sample(:,isam), tform)
371call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
374affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
382 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperDiag)
404affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
412 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperUnit)
434affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
442 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerDiag)
464affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
472 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerUnit)
499+5.00000000,
-2.00000000
502-7.00000000,
+0.00000000,
+6.00000000,
+8.00000000,
-5.00000000
503-4.00000000,
+5.00000000,
-1.00000000,
+3.00000000,
+7.00000000
508+1.00000000,
+0.987879276
509+0.987879276,
+1.00000000
512-10.9515171,
-10.9151554
513+4.93939638,
+5.00000000
514+5.01212072,
+4.92727566
515+10.9636383,
+10.9030342
516+1.91515493,
+2.06060362
518 affinity(:,isam)
= getAffinity(sample(:,isam), tform)
521-10.9515171,
-10.9151554
522+4.93939638,
+5.00000000
523+5.01212072,
+4.92727566
524+10.9636383,
+10.9030342
525+1.91515493,
+2.06060362
528-6.99996948,
-4.00003052
529+0.00000000,
+5.00000000
530+6.00000000,
-1.00001526
531+8.00000000,
+2.99996948
532-5.00000000,
+7.00000000
534-7.00000000,
-4.00000000
535+0.00000000,
+5.00000000
536+6.00000000,
-1.00000000
537+8.00000000,
+3.00000000
538-5.00000000,
+7.00000000
540call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
542+1.00000000,
+0.987879276
543+0.00000000,
+1.00000000
544affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
546-10.9515171,
-4.00000000
547+4.93939638,
+5.00000000
548+5.01212072,
-1.00000000
549+10.9636383,
+3.00000000
550+1.91515493,
+7.00000000
552 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperDiag)
555-10.9515171,
-4.00000000
556+4.93939638,
+5.00000000
557+5.01212072,
-1.00000000
558+10.9636383,
+3.00000000
559+1.91515493,
+7.00000000
562-7.00000000,
-4.00000000
563+0.00000000,
+5.00000000
564+6.00000000,
-1.00000000
565+8.00000000,
+3.00000000
566-5.00000000,
+7.00000000
568-7.00000000,
-4.00000000
569+0.00000000,
+5.00000000
570+6.00000000,
-1.00000000
571+8.00000000,
+3.00000000
572-5.00000000,
+7.00000000
574affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
576-10.9515171,
-4.00000000
577+4.93939638,
+5.00000000
578+5.01212072,
-1.00000000
579+10.9636383,
+3.00000000
580+1.91515493,
+7.00000000
582 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperUnit)
585-10.9515171,
-4.00000000
586+4.93939638,
+5.00000000
587+5.01212072,
-1.00000000
588+10.9636383,
+3.00000000
589+1.91515493,
+7.00000000
592-7.00000000,
-4.00000000
593+0.00000000,
+5.00000000
594+6.00000000,
-1.00000000
595+8.00000000,
+3.00000000
596-5.00000000,
+7.00000000
598-7.00000000,
-4.00000000
599+0.00000000,
+5.00000000
600+6.00000000,
-1.00000000
601+8.00000000,
+3.00000000
602-5.00000000,
+7.00000000
604affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
606-7.00000000,
-4.00000000
607+0.00000000,
+5.00000000
608+6.00000000,
-1.00000000
609+8.00000000,
+3.00000000
610-5.00000000,
+7.00000000
612 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerDiag)
615-7.00000000,
-4.00000000
616+0.00000000,
+5.00000000
617+6.00000000,
-1.00000000
618+8.00000000,
+3.00000000
619-5.00000000,
+7.00000000
622-7.00000000,
-4.00000000
623+0.00000000,
+5.00000000
624+6.00000000,
-1.00000000
625+8.00000000,
+3.00000000
626-5.00000000,
+7.00000000
628-7.00000000,
-4.00000000
629+0.00000000,
+5.00000000
630+6.00000000,
-1.00000000
631+8.00000000,
+3.00000000
632-5.00000000,
+7.00000000
634affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
636-7.00000000,
-4.00000000
637+0.00000000,
+5.00000000
638+6.00000000,
-1.00000000
639+8.00000000,
+3.00000000
640-5.00000000,
+7.00000000
642 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerUnit)
645-7.00000000,
-4.00000000
646+0.00000000,
+5.00000000
647+6.00000000,
-1.00000000
648+8.00000000,
+3.00000000
649-5.00000000,
+7.00000000
652-7.00000000,
-4.00000000
653+0.00000000,
+5.00000000
654+6.00000000,
-1.00000000
655+8.00000000,
+3.00000000
656-5.00000000,
+7.00000000
658-7.00000000,
-4.00000000
659+0.00000000,
+5.00000000
660+6.00000000,
-1.00000000
661+8.00000000,
+3.00000000
662-5.00000000,
+7.00000000
669+3.00000000,
+0.00000000,
+5.00000000
672+2.00000000,
+7.00000000,
-9.00000000,
+4.00000000
673+4.00000000,
+9.00000000,
+5.00000000,
-6.00000000
674-8.00000000,
-1.00000000,
+2.00000000,
-9.00000000
679+1.00000000,
-0.604018271,
+0.676454782
680-0.604018271,
+1.00000000,
+0.467890501E-1
681+0.676454782,
+0.467890501E-1,
+1.00000000
684-5.82771111,
+2.41765118,
-6.45993423
685+0.887380838,
+4.72508287,
+4.15628529
686-10.6671810,
+10.5297432,
-3.85414743
687+1.53601694,
-8.83717442,
-6.57491493
689 affinity(:,isam)
= getAffinity(sample(:,isam), tform)
692-5.82771111,
+2.41765118,
-6.45993423
693+0.887380838,
+4.72508287,
+4.15628529
694-10.6671810,
+10.5297432,
-3.85414743
695+1.53601694,
-8.83717442,
-6.57491493
698+2.00000000,
+4.00000191,
-8.00000191
699+7.00000000,
+8.99999714,
-0.999998093
700-8.99999046,
+5.00000191,
+1.99999237
701+4.00000000,
-6.00000000,
-9.00000191
703+2.00000000,
+4.00000000,
-8.00000000
704+7.00000000,
+9.00000000,
-1.00000000
705-9.00000000,
+5.00000000,
+2.00000000
706+4.00000000,
-6.00000000,
-9.00000000
708call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
710+1.00000000,
-0.604018271,
+0.676454782
711+0.00000000,
+1.00000000,
+0.467890501E-1
712+0.00000000,
+0.00000000,
+1.00000000
713affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
715-5.82771111,
+3.62568760,
-8.00000000
716+0.887380600,
+8.95321083,
-1.00000000
717-10.6671820,
+5.09357834,
+2.00000000
718+1.53601694,
-6.42110157,
-9.00000000
720 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperDiag)
723-5.82771111,
+3.62568760,
-8.00000000
724+0.887380600,
+8.95321083,
-1.00000000
725-10.6671820,
+5.09357834,
+2.00000000
726+1.53601694,
-6.42110157,
-9.00000000
729+2.00000000,
+4.00000000,
-8.00000000
730+6.99999952,
+9.00000000,
-1.00000000
731-9.00000000,
+5.00000000,
+2.00000000
732+4.00000000,
-6.00000000,
-9.00000000
734+2.00000000,
+4.00000000,
-8.00000000
735+7.00000000,
+9.00000000,
-1.00000000
736-9.00000000,
+5.00000000,
+2.00000000
737+4.00000000,
-6.00000000,
-9.00000000
739affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
741-5.82771111,
+3.62568760,
-8.00000000
742+0.887380600,
+8.95321083,
-1.00000000
743-10.6671820,
+5.09357834,
+2.00000000
744+1.53601694,
-6.42110157,
-9.00000000
746 affinity(:,isam)
= getAffinity(sample(:,isam), tform, upperUnit)
749-5.82771111,
+3.62568760,
-8.00000000
750+0.887380600,
+8.95321083,
-1.00000000
751-10.6671820,
+5.09357834,
+2.00000000
752+1.53601694,
-6.42110157,
-9.00000000
755+2.00000000,
+4.00000000,
-8.00000000
756+6.99999952,
+9.00000000,
-1.00000000
757-9.00000000,
+5.00000000,
+2.00000000
758+4.00000000,
-6.00000000,
-9.00000000
760+2.00000000,
+4.00000000,
-8.00000000
761+7.00000000,
+9.00000000,
-1.00000000
762-9.00000000,
+5.00000000,
+2.00000000
763+4.00000000,
-6.00000000,
-9.00000000
765affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
767+2.00000000,
+4.00000000,
-8.00000000
768+7.00000000,
+9.00000000,
-1.00000000
769-9.00000000,
+5.00000000,
+2.00000000
770+4.00000000,
-6.00000000,
-9.00000000
772 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerDiag)
775+2.00000000,
+4.00000000,
-8.00000000
776+7.00000000,
+9.00000000,
-1.00000000
777-9.00000000,
+5.00000000,
+2.00000000
778+4.00000000,
-6.00000000,
-9.00000000
781+2.00000000,
+4.00000000,
-8.00000000
782+7.00000000,
+9.00000000,
-1.00000000
783-9.00000000,
+5.00000000,
+2.00000000
784+4.00000000,
-6.00000000,
-9.00000000
786+2.00000000,
+4.00000000,
-8.00000000
787+7.00000000,
+9.00000000,
-1.00000000
788-9.00000000,
+5.00000000,
+2.00000000
789+4.00000000,
-6.00000000,
-9.00000000
791affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
793+2.00000000,
+4.00000000,
-8.00000000
794+7.00000000,
+9.00000000,
-1.00000000
795-9.00000000,
+5.00000000,
+2.00000000
796+4.00000000,
-6.00000000,
-9.00000000
798 affinity(:,isam)
= getAffinity(sample(:,isam), tform, lowerUnit)
801+2.00000000,
+4.00000000,
-8.00000000
802+7.00000000,
+9.00000000,
-1.00000000
803-9.00000000,
+5.00000000,
+2.00000000
804+4.00000000,
-6.00000000,
-9.00000000
807+2.00000000,
+4.00000000,
-8.00000000
808+7.00000000,
+9.00000000,
-1.00000000
809-9.00000000,
+5.00000000,
+2.00000000
810+4.00000000,
-6.00000000,
-9.00000000
812+2.00000000,
+4.00000000,
-8.00000000
813+7.00000000,
+9.00000000,
-1.00000000
814-9.00000000,
+5.00000000,
+2.00000000
815+4.00000000,
-6.00000000,
-9.00000000
848 affinity(isam,:)
= getAffinity(sample(isam,:), tform)
870call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
873affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
881 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperDiag)
903affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
911 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperUnit)
933affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
941 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerDiag)
963affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
971 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerUnit)
998+3.00000000,
+0.00000000,
+2.00000000
1001-9.00000000,
-9.00000000,
-6.00000000
1002-7.00000000,
-5.00000000,
-8.00000000
1003-9.00000000,
-2.00000000,
+3.00000000
1004-8.00000000,
-3.00000000,
+5.00000000
1009+1.00000000,
+0.999079287,
-0.286785483
1010+0.999079287,
+1.00000000,
-0.324553668
1011-0.286785483,
-0.324553668,
+1.00000000
1014-16.2710018,
-16.0443916,
-0.497947693
1015-9.70111275,
-9.39712524,
-4.36973333
1016-11.8585148,
-11.9653749,
+6.23017693
1017-12.4311657,
-12.6154022,
+8.26794434
1019 affinity(isam,:)
= getAffinity(sample(isam,:), tform)
1022-16.2710018,
-16.0443916,
-0.497947693
1023-9.70111275,
-9.39712524,
-4.36973333
1024-11.8585148,
-11.9653749,
+6.23017693
1025-12.4311657,
-12.6154022,
+8.26794434
1028-9.00228119,
-8.99603271,
-5.99978352
1029-6.99859619,
-5.00042725,
-8.00002670
1030-8.99652100,
-1.99938965,
+2.99975204
1031-8.00170898,
-2.99584961,
+5.00006866
1033-9.00000000,
-9.00000000,
-6.00000000
1034-7.00000000,
-5.00000000,
-8.00000000
1035-9.00000000,
-2.00000000,
+3.00000000
1036-8.00000000,
-3.00000000,
+5.00000000
1038call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1040+1.00000000,
+0.999079287,
-0.286785483
1041+0.00000000,
+1.00000000,
-0.324553668
1042+0.00000000,
+0.00000000,
+1.00000000
1043affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
1045-16.2709999,
-7.05267811,
-6.00000000
1046-9.70111275,
-2.40357065,
-8.00000000
1047-11.8585148,
-2.97366095,
+3.00000000
1048-12.4311657,
-4.62276840,
+5.00000000
1050 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperDiag)
1053-16.2709999,
-7.05267811,
-6.00000000
1054-9.70111275,
-2.40357065,
-8.00000000
1055-11.8585148,
-2.97366095,
+3.00000000
1056-12.4311657,
-4.62276840,
+5.00000000
1059-8.99999905,
-9.00000000,
-6.00000000
1060-7.00000000,
-5.00000000,
-8.00000000
1061-9.00000000,
-2.00000000,
+3.00000000
1062-8.00000000,
-3.00000000,
+5.00000000
1064-9.00000000,
-9.00000000,
-6.00000000
1065-7.00000000,
-5.00000000,
-8.00000000
1066-9.00000000,
-2.00000000,
+3.00000000
1067-8.00000000,
-3.00000000,
+5.00000000
1069affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
1071-16.2709999,
-7.05267811,
-6.00000000
1072-9.70111275,
-2.40357065,
-8.00000000
1073-11.8585148,
-2.97366095,
+3.00000000
1074-12.4311657,
-4.62276840,
+5.00000000
1076 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperUnit)
1079-16.2709999,
-7.05267811,
-6.00000000
1080-9.70111275,
-2.40357065,
-8.00000000
1081-11.8585148,
-2.97366095,
+3.00000000
1082-12.4311657,
-4.62276840,
+5.00000000
1085-8.99999905,
-9.00000000,
-6.00000000
1086-7.00000000,
-5.00000000,
-8.00000000
1087-9.00000000,
-2.00000000,
+3.00000000
1088-8.00000000,
-3.00000000,
+5.00000000
1090-9.00000000,
-9.00000000,
-6.00000000
1091-7.00000000,
-5.00000000,
-8.00000000
1092-9.00000000,
-2.00000000,
+3.00000000
1093-8.00000000,
-3.00000000,
+5.00000000
1095affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
1097-9.00000000,
-9.00000000,
-6.00000000
1098-7.00000000,
-5.00000000,
-8.00000000
1099-9.00000000,
-2.00000000,
+3.00000000
1100-8.00000000,
-3.00000000,
+5.00000000
1102 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerDiag)
1105-9.00000000,
-9.00000000,
-6.00000000
1106-7.00000000,
-5.00000000,
-8.00000000
1107-9.00000000,
-2.00000000,
+3.00000000
1108-8.00000000,
-3.00000000,
+5.00000000
1111-9.00000000,
-9.00000000,
-6.00000000
1112-7.00000000,
-5.00000000,
-8.00000000
1113-9.00000000,
-2.00000000,
+3.00000000
1114-8.00000000,
-3.00000000,
+5.00000000
1116-9.00000000,
-9.00000000,
-6.00000000
1117-7.00000000,
-5.00000000,
-8.00000000
1118-9.00000000,
-2.00000000,
+3.00000000
1119-8.00000000,
-3.00000000,
+5.00000000
1121affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
1123-9.00000000,
-9.00000000,
-6.00000000
1124-7.00000000,
-5.00000000,
-8.00000000
1125-9.00000000,
-2.00000000,
+3.00000000
1126-8.00000000,
-3.00000000,
+5.00000000
1128 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerUnit)
1131-9.00000000,
-9.00000000,
-6.00000000
1132-7.00000000,
-5.00000000,
-8.00000000
1133-9.00000000,
-2.00000000,
+3.00000000
1134-8.00000000,
-3.00000000,
+5.00000000
1137-9.00000000,
-9.00000000,
-6.00000000
1138-7.00000000,
-5.00000000,
-8.00000000
1139-9.00000000,
-2.00000000,
+3.00000000
1140-8.00000000,
-3.00000000,
+5.00000000
1142-9.00000000,
-9.00000000,
-6.00000000
1143-7.00000000,
-5.00000000,
-8.00000000
1144-9.00000000,
-2.00000000,
+3.00000000
1145-8.00000000,
-3.00000000,
+5.00000000
1173 affinity(isam,:)
= getAffinity(sample(isam,:), tform)
1195call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1198affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
1206 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperDiag)
1228affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
1236 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperUnit)
1258affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
1266 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerDiag)
1288affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
1296 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerUnit)
1323+0.00000000,
+2.00000000,
-3.00000000
1326-7.00000000,
-1.00000000,
+0.00000000
1327+7.00000000,
-7.00000000,
-9.00000000
1328-9.00000000,
+4.00000000,
+7.00000000
1329+7.00000000,
+7.00000000,
+9.00000000
1334+1.00000000,
+0.831882730E-1,
-0.676033199
1335+0.831882730E-1,
+1.00000000,
+0.297843218
1336-0.676033199,
+0.297843218,
+1.00000000
1339-7.08318806,
-1.58231783,
+4.43438911
1340+12.5019808,
-9.09827137,
-15.8171349
1341-13.3994789,
+5.33620834,
+14.2756710
1342+1.49801922,
+10.2629070,
+6.35266972
1344 affinity(isam,:)
= getAffinity(sample(isam,:), tform)
1347-7.08318806,
-1.58231783,
+4.43438911
1348+12.5019808,
-9.09827137,
-15.8171349
1349-13.3994789,
+5.33620834,
+14.2756710
1350+1.49801922,
+10.2629070,
+6.35266972
1353-7.00000048,
-0.999999762,
+0.00000000
1354+7.00000000,
-7.00000095,
-8.99999809
1355-9.00000191,
+4.00000191,
+7.00000000
1356+7.00000000,
+7.00000048,
+8.99999905
1358-7.00000000,
-1.00000000,
+0.00000000
1359+7.00000000,
-7.00000000,
-9.00000000
1360-9.00000000,
+4.00000000,
+7.00000000
1361+7.00000000,
+7.00000000,
+9.00000000
1363call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1365+1.00000000,
+0.831882730E-1,
-0.676033199
1366+0.00000000,
+1.00000000,
+0.297843218
1367+0.00000000,
+0.00000000,
+1.00000000
1368affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
1370-7.08318806,
-1.00000000,
+0.00000000
1371+12.5019808,
-9.68058872,
-9.00000000
1372-13.3994789,
+6.08490276,
+7.00000000
1373+1.49801922,
+9.68058872,
+9.00000000
1375 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperDiag)
1378-7.08318806,
-1.00000000,
+0.00000000
1379+12.5019808,
-9.68058872,
-9.00000000
1380-13.3994789,
+6.08490276,
+7.00000000
1381+1.49801922,
+9.68058872,
+9.00000000
1384-7.00000000,
-1.00000000,
+0.00000000
1385+7.00000000,
-7.00000000,
-9.00000000
1386-9.00000000,
+4.00000000,
+7.00000000
1387+7.00000000,
+7.00000000,
+9.00000000
1389-7.00000000,
-1.00000000,
+0.00000000
1390+7.00000000,
-7.00000000,
-9.00000000
1391-9.00000000,
+4.00000000,
+7.00000000
1392+7.00000000,
+7.00000000,
+9.00000000
1394affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
1396-7.08318806,
-1.00000000,
+0.00000000
1397+12.5019808,
-9.68058872,
-9.00000000
1398-13.3994789,
+6.08490276,
+7.00000000
1399+1.49801922,
+9.68058872,
+9.00000000
1401 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperUnit)
1404-7.08318806,
-1.00000000,
+0.00000000
1405+12.5019808,
-9.68058872,
-9.00000000
1406-13.3994789,
+6.08490276,
+7.00000000
1407+1.49801922,
+9.68058872,
+9.00000000
1410-7.00000000,
-1.00000000,
+0.00000000
1411+7.00000000,
-7.00000000,
-9.00000000
1412-9.00000000,
+4.00000000,
+7.00000000
1413+7.00000000,
+7.00000000,
+9.00000000
1415-7.00000000,
-1.00000000,
+0.00000000
1416+7.00000000,
-7.00000000,
-9.00000000
1417-9.00000000,
+4.00000000,
+7.00000000
1418+7.00000000,
+7.00000000,
+9.00000000
1420affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
1422-7.00000000,
-1.00000000,
+0.00000000
1423+7.00000000,
-7.00000000,
-9.00000000
1424-9.00000000,
+4.00000000,
+7.00000000
1425+7.00000000,
+7.00000000,
+9.00000000
1427 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerDiag)
1430-7.00000000,
-1.00000000,
+0.00000000
1431+7.00000000,
-7.00000000,
-9.00000000
1432-9.00000000,
+4.00000000,
+7.00000000
1433+7.00000000,
+7.00000000,
+9.00000000
1436-7.00000000,
-1.00000000,
+0.00000000
1437+7.00000000,
-7.00000000,
-9.00000000
1438-9.00000000,
+4.00000000,
+7.00000000
1439+7.00000000,
+7.00000000,
+9.00000000
1441-7.00000000,
-1.00000000,
+0.00000000
1442+7.00000000,
-7.00000000,
-9.00000000
1443-9.00000000,
+4.00000000,
+7.00000000
1444+7.00000000,
+7.00000000,
+9.00000000
1446affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
1448-7.00000000,
-1.00000000,
+0.00000000
1449+7.00000000,
-7.00000000,
-9.00000000
1450-9.00000000,
+4.00000000,
+7.00000000
1451+7.00000000,
+7.00000000,
+9.00000000
1453 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerUnit)
1456-7.00000000,
-1.00000000,
+0.00000000
1457+7.00000000,
-7.00000000,
-9.00000000
1458-9.00000000,
+4.00000000,
+7.00000000
1459+7.00000000,
+7.00000000,
+9.00000000
1462-7.00000000,
-1.00000000,
+0.00000000
1463+7.00000000,
-7.00000000,
-9.00000000
1464-9.00000000,
+4.00000000,
+7.00000000
1465+7.00000000,
+7.00000000,
+9.00000000
1467-7.00000000,
-1.00000000,
+0.00000000
1468+7.00000000,
-7.00000000,
-9.00000000
1469-9.00000000,
+4.00000000,
+7.00000000
1470+7.00000000,
+7.00000000,
+9.00000000
1477+0.00000000,
+4.00000000
1480+0.00000000,
-8.00000000
1481-1.00000000,
-6.00000000
1482-2.00000000,
+5.00000000
1483+2.00000000,
+5.00000000
1484-2.00000000,
-1.00000000
1489+1.00000000,
-0.977806449
1490-0.977806449,
+1.00000000
1493+7.82245159,
-8.00000000
1494+4.86683846,
-5.02219343
1495-6.88903236,
+6.95561314
1496-2.88903236,
+3.04438710
1497-1.02219355,
+0.955612898
1499 affinity(isam,:)
= getAffinity(sample(isam,:), tform)
1502+7.82245159,
-8.00000000
1503+4.86683846,
-5.02219343
1504-6.88903236,
+6.95561314
1505-2.88903236,
+3.04438710
1506-1.02219355,
+0.955612898
1509+0.00000000,
-8.00000000
1510-1.00000763,
-6.00000000
1511-2.00000000,
+5.00000000
1512+2.00000000,
+4.99999237
1513-2.00000191,
-1.00000191
1515+0.00000000,
-8.00000000
1516-1.00000000,
-6.00000000
1517-2.00000000,
+5.00000000
1518+2.00000000,
+5.00000000
1519-2.00000000,
-1.00000000
1521call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1523+1.00000000,
-0.977806449
1524+0.00000000,
+1.00000000
1525affinity
= getAffinity(sample, dim, tform,
class = upperDiag)
1527+7.82245159,
-8.00000000
1528+4.86683846,
-6.00000000
1529-6.88903236,
+5.00000000
1530-2.88903236,
+5.00000000
1531-1.02219355,
-1.00000000
1533 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperDiag)
1536+7.82245159,
-8.00000000
1537+4.86683846,
-6.00000000
1538-6.88903236,
+5.00000000
1539-2.88903236,
+5.00000000
1540-1.02219355,
-1.00000000
1543+0.00000000,
-8.00000000
1544-1.00000000,
-6.00000000
1545-2.00000000,
+5.00000000
1546+2.00000000,
+5.00000000
1547-2.00000000,
-1.00000000
1549+0.00000000,
-8.00000000
1550-1.00000000,
-6.00000000
1551-2.00000000,
+5.00000000
1552+2.00000000,
+5.00000000
1553-2.00000000,
-1.00000000
1555affinity
= getAffinity(sample, dim, tform,
class = upperUnit)
1557+7.82245159,
-8.00000000
1558+4.86683846,
-6.00000000
1559-6.88903236,
+5.00000000
1560-2.88903236,
+5.00000000
1561-1.02219355,
-1.00000000
1563 affinity(isam,:)
= getAffinity(sample(isam,:), tform, upperUnit)
1566+7.82245159,
-8.00000000
1567+4.86683846,
-6.00000000
1568-6.88903236,
+5.00000000
1569-2.88903236,
+5.00000000
1570-1.02219355,
-1.00000000
1573+0.00000000,
-8.00000000
1574-1.00000000,
-6.00000000
1575-2.00000000,
+5.00000000
1576+2.00000000,
+5.00000000
1577-2.00000000,
-1.00000000
1579+0.00000000,
-8.00000000
1580-1.00000000,
-6.00000000
1581-2.00000000,
+5.00000000
1582+2.00000000,
+5.00000000
1583-2.00000000,
-1.00000000
1585affinity
= getAffinity(sample, dim, tform,
class = lowerDiag)
1587+0.00000000,
-8.00000000
1588-1.00000000,
-6.00000000
1589-2.00000000,
+5.00000000
1590+2.00000000,
+5.00000000
1591-2.00000000,
-1.00000000
1593 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerDiag)
1596+0.00000000,
-8.00000000
1597-1.00000000,
-6.00000000
1598-2.00000000,
+5.00000000
1599+2.00000000,
+5.00000000
1600-2.00000000,
-1.00000000
1603+0.00000000,
-8.00000000
1604-1.00000000,
-6.00000000
1605-2.00000000,
+5.00000000
1606+2.00000000,
+5.00000000
1607-2.00000000,
-1.00000000
1609+0.00000000,
-8.00000000
1610-1.00000000,
-6.00000000
1611-2.00000000,
+5.00000000
1612+2.00000000,
+5.00000000
1613-2.00000000,
-1.00000000
1615affinity
= getAffinity(sample, dim, tform,
class = lowerUnit)
1617+0.00000000,
-8.00000000
1618-1.00000000,
-6.00000000
1619-2.00000000,
+5.00000000
1620+2.00000000,
+5.00000000
1621-2.00000000,
-1.00000000
1623 affinity(isam,:)
= getAffinity(sample(isam,:), tform, lowerUnit)
1626+0.00000000,
-8.00000000
1627-1.00000000,
-6.00000000
1628-2.00000000,
+5.00000000
1629+2.00000000,
+5.00000000
1630-2.00000000,
-1.00000000
1633+0.00000000,
-8.00000000
1634-1.00000000,
-6.00000000
1635-2.00000000,
+5.00000000
1636+2.00000000,
+5.00000000
1637-2.00000000,
-1.00000000
1639+0.00000000,
-8.00000000
1640-1.00000000,
-6.00000000
1641-2.00000000,
+5.00000000
1642+2.00000000,
+5.00000000
1643-2.00000000,
-1.00000000
3import matplotlib.pyplot
as plt
12for kind
in [
"warp",
"rotation"]:
14 pattern =
"*."+kind+
"*.txt"
15 fileList = glob.glob(pattern)
19 df = pd.read_csv(file, delimiter =
",", header =
None)
22 left, width = 0.1, 0.65
23 bottom, height = 0.1, 0.65
27 fig = plt.figure(figsize=(8, 8))
29 plt.rcParams.update({
'font.size': fontsize - 2})
30 ax = fig.add_axes([left, bottom, width, height])
31 ax_histx = fig.add_axes([left, bottom + height + spacing, width, 0.2], sharex = ax)
32 ax_histy = fig.add_axes([left + width + spacing, bottom, 0.2, height], sharey = ax)
34 for axes
in [ax, ax_histx, ax_histy]:
35 axes.grid(visible =
True, which =
"both", axis =
"both", color =
"0.85", linestyle =
"-")
36 axes.tick_params(axis =
"y", which =
"minor")
37 axes.tick_params(axis =
"x", which =
"minor")
40 ax_histy.tick_params(axis =
"y", labelleft =
False)
41 ax_histx.tick_params(axis =
"x", labelbottom =
False)
44 ax.scatter ( df.values[:, 0]
50 ax_histx.hist(df.values[:, 0], bins = 50, zorder = 1000)
51 ax_histy.hist(df.values[:, 1], bins = 50, orientation =
"horizontal", zorder = 1000)
53 ax.set_xlabel(
"X", fontsize = 17)
54 ax.set_ylabel(
"Y", fontsize = 17)
55 ax.legend([file.split(
".")[-2]], fontsize = fontsize)
58 plt.savefig(file.replace(
".txt",
".png"))