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.
15+8.00000000,
+5.00000000,
-2.00000000,
-2.00000000,
+8.00000000
21call setAffinity(affinity, sample, dim, tform, genrecmat)
29 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
51call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
54call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
62 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
84call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
92 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
114call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
122 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
144call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
152 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
179-2.00000000,
+3.00000000
182-6.00000000,
+0.00000000,
-8.00000000,
-2.00000000,
-6.00000000
183-1.00000000,
+9.00000000,
+3.00000000,
-1.00000000,
+8.00000000
188+1.00000000,
-0.881028533
189-0.881028533,
+1.00000000
190call setAffinity(affinity, sample, dim, tform, genrecmat)
192-5.11897135,
+4.28617096
193-7.92925692,
+9.00000000
194-10.6430855,
+10.0482283
195-1.11897147,
+0.762057066
196-13.0482283,
+13.2861710
198 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
201-5.11897135,
+4.28617096
202-7.92925692,
+9.00000000
203-10.6430855,
+10.0482283
204-1.11897147,
+0.762057066
205-13.0482283,
+13.2861710
208-6.00000000,
-1.00000000
209+0.00000000,
+8.99999809
210-8.00000000,
+3.00000000
211-1.99999976,
-0.999999762
212-6.00000381,
+7.99999619
214-6.00000000,
-1.00000000
215+0.00000000,
+9.00000000
216-8.00000000,
+3.00000000
217-2.00000000,
-1.00000000
218-6.00000000,
+8.00000000
220call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
222+1.00000000,
-0.881028533
223+0.00000000,
+1.00000000
224call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
226-5.11897135,
-1.00000000
227-7.92925692,
+9.00000000
228-10.6430855,
+3.00000000
229-1.11897147,
-1.00000000
230-13.0482283,
+8.00000000
232 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
235-5.11897135,
-1.00000000
236-7.92925692,
+9.00000000
237-10.6430855,
+3.00000000
238-1.11897147,
-1.00000000
239-13.0482283,
+8.00000000
242-6.00000000,
-1.00000000
243+0.00000000,
+9.00000000
244-8.00000000,
+3.00000000
245-2.00000000,
-1.00000000
246-6.00000000,
+8.00000000
248-6.00000000,
-1.00000000
249+0.00000000,
+9.00000000
250-8.00000000,
+3.00000000
251-2.00000000,
-1.00000000
252-6.00000000,
+8.00000000
254call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
256-5.11897135,
-1.00000000
257-7.92925692,
+9.00000000
258-10.6430855,
+3.00000000
259-1.11897147,
-1.00000000
260-13.0482283,
+8.00000000
262 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
265-5.11897135,
-1.00000000
266-7.92925692,
+9.00000000
267-10.6430855,
+3.00000000
268-1.11897147,
-1.00000000
269-13.0482283,
+8.00000000
272-6.00000000,
-1.00000000
273+0.00000000,
+9.00000000
274-8.00000000,
+3.00000000
275-2.00000000,
-1.00000000
276-6.00000000,
+8.00000000
278-6.00000000,
-1.00000000
279+0.00000000,
+9.00000000
280-8.00000000,
+3.00000000
281-2.00000000,
-1.00000000
282-6.00000000,
+8.00000000
284call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
286-6.00000000,
-1.00000000
287+0.00000000,
+9.00000000
288-8.00000000,
+3.00000000
289-2.00000000,
-1.00000000
290-6.00000000,
+8.00000000
292 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
295-6.00000000,
-1.00000000
296+0.00000000,
+9.00000000
297-8.00000000,
+3.00000000
298-2.00000000,
-1.00000000
299-6.00000000,
+8.00000000
302-6.00000000,
-1.00000000
303+0.00000000,
+9.00000000
304-8.00000000,
+3.00000000
305-2.00000000,
-1.00000000
306-6.00000000,
+8.00000000
308-6.00000000,
-1.00000000
309+0.00000000,
+9.00000000
310-8.00000000,
+3.00000000
311-2.00000000,
-1.00000000
312-6.00000000,
+8.00000000
314call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
316-6.00000000,
-1.00000000
317+0.00000000,
+9.00000000
318-8.00000000,
+3.00000000
319-2.00000000,
-1.00000000
320-6.00000000,
+8.00000000
322 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
325-6.00000000,
-1.00000000
326+0.00000000,
+9.00000000
327-8.00000000,
+3.00000000
328-2.00000000,
-1.00000000
329-6.00000000,
+8.00000000
332-6.00000000,
-1.00000000
333+0.00000000,
+9.00000000
334-8.00000000,
+3.00000000
335-2.00000000,
-1.00000000
336-6.00000000,
+8.00000000
338-6.00000000,
-1.00000000
339+0.00000000,
+9.00000000
340-8.00000000,
+3.00000000
341-2.00000000,
-1.00000000
342-6.00000000,
+8.00000000
349-4.00000000,
+1.00000000
352-2.00000000,
+7.00000000,
-8.00000000,
+1.00000000
353-1.00000000,
+3.00000000,
-1.00000000,
-2.00000000
358+1.00000000,
+0.995931149
359+0.995931149,
+1.00000000
360call setAffinity(affinity, sample, dim, tform, genrecmat)
362-2.99593115,
-2.99186230
363+9.98779297,
+9.97151756
364-8.99593163,
-8.96744919
365-0.991862297,
-1.00406885
367 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
370-2.99593115,
-2.99186230
371+9.98779297,
+9.97151756
372-8.99593163,
-8.96744919
373-0.991862297,
-1.00406885
376-2.00003052,
-0.999969482
377+7.00000000,
+3.00000000
378-8.00000000,
-0.999877930
379+1.00000000,
-2.00000000
381-2.00000000,
-1.00000000
382+7.00000000,
+3.00000000
383-8.00000000,
-1.00000000
384+1.00000000,
-2.00000000
386call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
388+1.00000000,
+0.995931149
389+0.00000000,
+1.00000000
390call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
392-2.99593115,
-1.00000000
393+9.98779297,
+3.00000000
394-8.99593163,
-1.00000000
395-0.991862297,
-2.00000000
397 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
400-2.99593115,
-1.00000000
401+9.98779297,
+3.00000000
402-8.99593163,
-1.00000000
403-0.991862297,
-2.00000000
406-2.00000000,
-1.00000000
407+6.99999952,
+3.00000000
408-8.00000000,
-1.00000000
409+1.00000000,
-2.00000000
411-2.00000000,
-1.00000000
412+7.00000000,
+3.00000000
413-8.00000000,
-1.00000000
414+1.00000000,
-2.00000000
416call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
418-2.99593115,
-1.00000000
419+9.98779297,
+3.00000000
420-8.99593163,
-1.00000000
421-0.991862297,
-2.00000000
423 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
426-2.99593115,
-1.00000000
427+9.98779297,
+3.00000000
428-8.99593163,
-1.00000000
429-0.991862297,
-2.00000000
432-2.00000000,
-1.00000000
433+6.99999952,
+3.00000000
434-8.00000000,
-1.00000000
435+1.00000000,
-2.00000000
437-2.00000000,
-1.00000000
438+7.00000000,
+3.00000000
439-8.00000000,
-1.00000000
440+1.00000000,
-2.00000000
442call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
444-2.00000000,
-1.00000000
445+7.00000000,
+3.00000000
446-8.00000000,
-1.00000000
447+1.00000000,
-2.00000000
449 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
452-2.00000000,
-1.00000000
453+7.00000000,
+3.00000000
454-8.00000000,
-1.00000000
455+1.00000000,
-2.00000000
458-2.00000000,
-1.00000000
459+7.00000000,
+3.00000000
460-8.00000000,
-1.00000000
461+1.00000000,
-2.00000000
463-2.00000000,
-1.00000000
464+7.00000000,
+3.00000000
465-8.00000000,
-1.00000000
466+1.00000000,
-2.00000000
468call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
470-2.00000000,
-1.00000000
471+7.00000000,
+3.00000000
472-8.00000000,
-1.00000000
473+1.00000000,
-2.00000000
475 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
478-2.00000000,
-1.00000000
479+7.00000000,
+3.00000000
480-8.00000000,
-1.00000000
481+1.00000000,
-2.00000000
484-2.00000000,
-1.00000000
485+7.00000000,
+3.00000000
486-8.00000000,
-1.00000000
487+1.00000000,
-2.00000000
489-2.00000000,
-1.00000000
490+7.00000000,
+3.00000000
491-8.00000000,
-1.00000000
492+1.00000000,
-2.00000000
499-1.00000000,
-5.00000000,
+4.00000000
502+6.00000000,
+9.00000000,
+5.00000000,
+6.00000000
503+1.00000000,
+7.00000000,
-4.00000000,
+8.00000000
504+0.00000000,
+9.00000000,
+0.00000000,
-5.00000000
509+1.00000000,
-0.774100065,
+0.233371958
510-0.774100065,
+1.00000000,
-0.511827469
511+0.233371958,
-0.511827469,
+1.00000000
512call setAffinity(affinity, sample, dim, tform, genrecmat)
514+5.22589970,
-3.64460039,
+0.888404250
515+5.68164730,
-4.57334805,
+7.51755524
516+8.09640026,
-7.87050056,
+3.21416950
517-1.35966027,
+5.91453695,
-7.69438791
519 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
522+5.22589970,
-3.64460039,
+0.888404250
523+5.68164730,
-4.57334805,
+7.51755524
524+8.09640026,
-7.87050056,
+3.21416950
525-1.35966027,
+5.91453695,
-7.69438791
528+6.00000095,
+1.00000024,
-0.476837158E-6
529+9.00000095,
+6.99999905,
+9.00000000
530+5.00000334,
-3.99999905,
-0.953674316E-6
531+5.99999809,
+7.99999809,
-5.00000000
533+6.00000000,
+1.00000000,
+0.00000000
534+9.00000000,
+7.00000000,
+9.00000000
535+5.00000000,
-4.00000000,
+0.00000000
536+6.00000000,
+8.00000000,
-5.00000000
538call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
540+1.00000000,
-0.774100065,
+0.233371958
541+0.00000000,
+1.00000000,
-0.511827469
542+0.00000000,
+0.00000000,
+1.00000000
543call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
545+5.22589970,
+1.00000000,
+0.00000000
546+5.68164730,
+2.39355278,
+9.00000000
547+8.09640026,
-4.00000000,
+0.00000000
548-1.35966015,
+10.5591373,
-5.00000000
550 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
553+5.22589970,
+1.00000000,
+0.00000000
554+5.68164730,
+2.39355278,
+9.00000000
555+8.09640026,
-4.00000000,
+0.00000000
556-1.35966015,
+10.5591373,
-5.00000000
559+6.00000000,
+1.00000000,
+0.00000000
560+9.00000000,
+7.00000000,
+9.00000000
561+5.00000000,
-4.00000000,
+0.00000000
562+6.00000048,
+8.00000000,
-5.00000000
564+6.00000000,
+1.00000000,
+0.00000000
565+9.00000000,
+7.00000000,
+9.00000000
566+5.00000000,
-4.00000000,
+0.00000000
567+6.00000000,
+8.00000000,
-5.00000000
569call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
571+5.22589970,
+1.00000000,
+0.00000000
572+5.68164730,
+2.39355278,
+9.00000000
573+8.09640026,
-4.00000000,
+0.00000000
574-1.35966015,
+10.5591373,
-5.00000000
576 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
579+5.22589970,
+1.00000000,
+0.00000000
580+5.68164730,
+2.39355278,
+9.00000000
581+8.09640026,
-4.00000000,
+0.00000000
582-1.35966015,
+10.5591373,
-5.00000000
585+6.00000000,
+1.00000000,
+0.00000000
586+9.00000000,
+7.00000000,
+9.00000000
587+5.00000000,
-4.00000000,
+0.00000000
588+6.00000048,
+8.00000000,
-5.00000000
590+6.00000000,
+1.00000000,
+0.00000000
591+9.00000000,
+7.00000000,
+9.00000000
592+5.00000000,
-4.00000000,
+0.00000000
593+6.00000000,
+8.00000000,
-5.00000000
595call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
597+6.00000000,
+1.00000000,
+0.00000000
598+9.00000000,
+7.00000000,
+9.00000000
599+5.00000000,
-4.00000000,
+0.00000000
600+6.00000000,
+8.00000000,
-5.00000000
602 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
605+6.00000000,
+1.00000000,
+0.00000000
606+9.00000000,
+7.00000000,
+9.00000000
607+5.00000000,
-4.00000000,
+0.00000000
608+6.00000000,
+8.00000000,
-5.00000000
611+6.00000000,
+1.00000000,
+0.00000000
612+9.00000000,
+7.00000000,
+9.00000000
613+5.00000000,
-4.00000000,
+0.00000000
614+6.00000000,
+8.00000000,
-5.00000000
616+6.00000000,
+1.00000000,
+0.00000000
617+9.00000000,
+7.00000000,
+9.00000000
618+5.00000000,
-4.00000000,
+0.00000000
619+6.00000000,
+8.00000000,
-5.00000000
621call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
623+6.00000000,
+1.00000000,
+0.00000000
624+9.00000000,
+7.00000000,
+9.00000000
625+5.00000000,
-4.00000000,
+0.00000000
626+6.00000000,
+8.00000000,
-5.00000000
628 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
631+6.00000000,
+1.00000000,
+0.00000000
632+9.00000000,
+7.00000000,
+9.00000000
633+5.00000000,
-4.00000000,
+0.00000000
634+6.00000000,
+8.00000000,
-5.00000000
637+6.00000000,
+1.00000000,
+0.00000000
638+9.00000000,
+7.00000000,
+9.00000000
639+5.00000000,
-4.00000000,
+0.00000000
640+6.00000000,
+8.00000000,
-5.00000000
642+6.00000000,
+1.00000000,
+0.00000000
643+9.00000000,
+7.00000000,
+9.00000000
644+5.00000000,
-4.00000000,
+0.00000000
645+6.00000000,
+8.00000000,
-5.00000000
655-2.00000000,
-2.00000000,
+6.00000000,
-9.00000000
661call setAffinity(affinity, sample, dim, tform, genrecmat)
668 call setAffinity(affinity(:,isam), sample(:,isam), tform, genrecmat)
687call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
690call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
697 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperDiag)
716call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
723 call setAffinity(affinity(:,isam), sample(:,isam), tform, upperUnit)
742call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
749 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerDiag)
768call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
775 call setAffinity(affinity(:,isam), sample(:,isam), tform, lowerUnit)
804-2.00000000,
-2.00000000
807+6.00000000,
+6.00000000
808+7.00000000,
+4.00000000
809+4.00000000,
-1.00000000
810+1.00000000,
-4.00000000
811-1.00000000,
-1.00000000
816+1.00000000,
-0.973061204
817-0.973061204,
+1.00000000
820+0.161632538,
+0.161632538
821+3.10775518,
-2.81142855
822+4.97306108,
-4.89224482
823+4.89224482,
-4.97306108
824-0.269387960E-1,
-0.269387960E-1
826 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
829+0.161632538,
+0.161632538
830+3.10775518,
-2.81142855
831+4.97306108,
-4.89224482
832+4.89224482,
-4.97306108
833-0.269387960E-1,
-0.269387960E-1
836+5.99998999,
+5.99998999
837+6.99999619,
+4.00000000
838+4.00000000,
-1.00000000
839+1.00000000,
-4.00000000
840-0.999999762,
-0.999999762
842+6.00000000,
+6.00000000
843+7.00000000,
+4.00000000
844+4.00000000,
-1.00000000
845+1.00000000,
-4.00000000
846-1.00000000,
-1.00000000
848call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
850+1.00000000,
-0.973061204
851+0.00000000,
+1.00000000
852call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
854+0.161632538,
+6.00000000
855+3.10775518,
+4.00000000
856+4.97306108,
-1.00000000
857+4.89224482,
-4.00000000
858-0.269387960E-1,
-1.00000000
860 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
863+0.161632538,
+6.00000000
864+3.10775518,
+4.00000000
865+4.97306108,
-1.00000000
866+4.89224482,
-4.00000000
867-0.269387960E-1,
-1.00000000
870+6.00000000,
+6.00000000
871+7.00000000,
+4.00000000
872+4.00000000,
-1.00000000
873+1.00000000,
-4.00000000
874-1.00000000,
-1.00000000
876+6.00000000,
+6.00000000
877+7.00000000,
+4.00000000
878+4.00000000,
-1.00000000
879+1.00000000,
-4.00000000
880-1.00000000,
-1.00000000
882call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
884+0.161632538,
+6.00000000
885+3.10775518,
+4.00000000
886+4.97306108,
-1.00000000
887+4.89224482,
-4.00000000
888-0.269387960E-1,
-1.00000000
890 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
893+0.161632538,
+6.00000000
894+3.10775518,
+4.00000000
895+4.97306108,
-1.00000000
896+4.89224482,
-4.00000000
897-0.269387960E-1,
-1.00000000
900+6.00000000,
+6.00000000
901+7.00000000,
+4.00000000
902+4.00000000,
-1.00000000
903+1.00000000,
-4.00000000
904-1.00000000,
-1.00000000
906+6.00000000,
+6.00000000
907+7.00000000,
+4.00000000
908+4.00000000,
-1.00000000
909+1.00000000,
-4.00000000
910-1.00000000,
-1.00000000
912call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
914+6.00000000,
+6.00000000
915+7.00000000,
+4.00000000
916+4.00000000,
-1.00000000
917+1.00000000,
-4.00000000
918-1.00000000,
-1.00000000
920 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
923+6.00000000,
+6.00000000
924+7.00000000,
+4.00000000
925+4.00000000,
-1.00000000
926+1.00000000,
-4.00000000
927-1.00000000,
-1.00000000
930+6.00000000,
+6.00000000
931+7.00000000,
+4.00000000
932+4.00000000,
-1.00000000
933+1.00000000,
-4.00000000
934-1.00000000,
-1.00000000
936+6.00000000,
+6.00000000
937+7.00000000,
+4.00000000
938+4.00000000,
-1.00000000
939+1.00000000,
-4.00000000
940-1.00000000,
-1.00000000
942call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
944+6.00000000,
+6.00000000
945+7.00000000,
+4.00000000
946+4.00000000,
-1.00000000
947+1.00000000,
-4.00000000
948-1.00000000,
-1.00000000
950 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)
953+6.00000000,
+6.00000000
954+7.00000000,
+4.00000000
955+4.00000000,
-1.00000000
956+1.00000000,
-4.00000000
957-1.00000000,
-1.00000000
960+6.00000000,
+6.00000000
961+7.00000000,
+4.00000000
962+4.00000000,
-1.00000000
963+1.00000000,
-4.00000000
964-1.00000000,
-1.00000000
966+6.00000000,
+6.00000000
967+7.00000000,
+4.00000000
968+4.00000000,
-1.00000000
969+1.00000000,
-4.00000000
970-1.00000000,
-1.00000000
977+4.00000000,
+4.00000000
980-5.00000000,
+5.00000000
981+7.00000000,
-8.00000000
982+7.00000000,
-3.00000000
983+1.00000000,
-1.00000000
988+1.00000000,
+0.168912858
989+0.168912858,
+1.00000000
992-4.15543556,
+4.15543556
993+5.64869690,
-6.81760979
994+6.49326134,
-1.81761003
995+0.831087112,
-0.831087112
997 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
1000-4.15543556,
+4.15543556
1001+5.64869690,
-6.81760979
1002+6.49326134,
-1.81761003
1003+0.831087112,
-0.831087112
1006-5.00000000,
+5.00000000
1007+6.99999952,
-7.99999952
1008+7.00000000,
-3.00000000
1009+0.999999881,
-0.999999940
1011-5.00000000,
+5.00000000
1012+7.00000000,
-8.00000000
1013+7.00000000,
-3.00000000
1014+1.00000000,
-1.00000000
1016call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1018+1.00000000,
+0.168912858
1019+0.00000000,
+1.00000000
1020call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
1022-4.15543556,
+5.00000000
1023+5.64869690,
-8.00000000
1024+6.49326134,
-3.00000000
1025+0.831087112,
-1.00000000
1027 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
1030-4.15543556,
+5.00000000
1031+5.64869690,
-8.00000000
1032+6.49326134,
-3.00000000
1033+0.831087112,
-1.00000000
1036-5.00000000,
+5.00000000
1037+7.00000000,
-8.00000000
1038+7.00000000,
-3.00000000
1039+1.00000000,
-1.00000000
1041-5.00000000,
+5.00000000
1042+7.00000000,
-8.00000000
1043+7.00000000,
-3.00000000
1044+1.00000000,
-1.00000000
1046call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
1048-4.15543556,
+5.00000000
1049+5.64869690,
-8.00000000
1050+6.49326134,
-3.00000000
1051+0.831087112,
-1.00000000
1053 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
1056-4.15543556,
+5.00000000
1057+5.64869690,
-8.00000000
1058+6.49326134,
-3.00000000
1059+0.831087112,
-1.00000000
1062-5.00000000,
+5.00000000
1063+7.00000000,
-8.00000000
1064+7.00000000,
-3.00000000
1065+1.00000000,
-1.00000000
1067-5.00000000,
+5.00000000
1068+7.00000000,
-8.00000000
1069+7.00000000,
-3.00000000
1070+1.00000000,
-1.00000000
1072call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
1074-5.00000000,
+5.00000000
1075+7.00000000,
-8.00000000
1076+7.00000000,
-3.00000000
1077+1.00000000,
-1.00000000
1079 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
1082-5.00000000,
+5.00000000
1083+7.00000000,
-8.00000000
1084+7.00000000,
-3.00000000
1085+1.00000000,
-1.00000000
1088-5.00000000,
+5.00000000
1089+7.00000000,
-8.00000000
1090+7.00000000,
-3.00000000
1091+1.00000000,
-1.00000000
1093-5.00000000,
+5.00000000
1094+7.00000000,
-8.00000000
1095+7.00000000,
-3.00000000
1096+1.00000000,
-1.00000000
1098call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
1100-5.00000000,
+5.00000000
1101+7.00000000,
-8.00000000
1102+7.00000000,
-3.00000000
1103+1.00000000,
-1.00000000
1105 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)
1108-5.00000000,
+5.00000000
1109+7.00000000,
-8.00000000
1110+7.00000000,
-3.00000000
1111+1.00000000,
-1.00000000
1114-5.00000000,
+5.00000000
1115+7.00000000,
-8.00000000
1116+7.00000000,
-3.00000000
1117+1.00000000,
-1.00000000
1119-5.00000000,
+5.00000000
1120+7.00000000,
-8.00000000
1121+7.00000000,
-3.00000000
1122+1.00000000,
-1.00000000
1129+0.00000000,
+4.00000000,
-3.00000000
1132+4.00000000,
-3.00000000,
-3.00000000
1133-2.00000000,
+0.00000000,
+6.00000000
1134+0.00000000,
-5.00000000,
+5.00000000
1135+1.00000000,
+6.00000000,
+0.00000000
1140+1.00000000,
-0.634198010,
-0.598320663
1141-0.634198010,
+1.00000000,
+0.953590631
1142-0.598320663,
+0.953590631,
+1.00000000
1145+7.69755602,
-8.39756393,
-8.25405502
1146-5.58992386,
+6.98993969,
+7.19664145
1147+0.179386616,
-0.232047081,
+0.232047081
1148-2.80518818,
+5.36580181,
+5.12322330
1150 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
1153+7.69755602,
-8.39756393,
-8.25405502
1154-5.58992386,
+6.98993969,
+7.19664145
1155+0.179386616,
-0.232047081,
+0.232047081
1156-2.80518818,
+5.36580181,
+5.12322330
1159+3.99999380,
-3.00000000,
-3.00000000
1160-1.99999571,
+0.762939453E-5,
+6.00000763
1161-0.413507223E-6,
-5.00000954,
+5.00000906
1162+1.00000334,
+6.00000381,
+0.00000000
1164+4.00000000,
-3.00000000,
-3.00000000
1165-2.00000000,
+0.00000000,
+6.00000000
1166+0.00000000,
-5.00000000,
+5.00000000
1167+1.00000000,
+6.00000000,
+0.00000000
1169call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1171+1.00000000,
-0.634198010,
-0.598320663
1172+0.00000000,
+1.00000000,
+0.953590631
1173+0.00000000,
+0.00000000,
+1.00000000
1174call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
1176+7.69755602,
-5.86077213,
-3.00000000
1177-5.58992386,
+5.72154379,
+6.00000000
1178+0.179386616,
-0.232047081,
+5.00000000
1179-2.80518818,
+6.00000000,
+0.00000000
1181 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
1184+7.69755602,
-5.86077213,
-3.00000000
1185-5.58992386,
+5.72154379,
+6.00000000
1186+0.179386616,
-0.232047081,
+5.00000000
1187-2.80518818,
+6.00000000,
+0.00000000
1190+3.99999976,
-3.00000024,
-3.00000000
1191-2.00000000,
+0.00000000,
+6.00000000
1192-0.342726707E-6,
-5.00000000,
+5.00000000
1193+1.00000000,
+6.00000000,
+0.00000000
1195+4.00000000,
-3.00000000,
-3.00000000
1196-2.00000000,
+0.00000000,
+6.00000000
1197+0.00000000,
-5.00000000,
+5.00000000
1198+1.00000000,
+6.00000000,
+0.00000000
1200call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
1202+7.69755602,
-5.86077213,
-3.00000000
1203-5.58992386,
+5.72154379,
+6.00000000
1204+0.179386616,
-0.232047081,
+5.00000000
1205-2.80518818,
+6.00000000,
+0.00000000
1207 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
1210+7.69755602,
-5.86077213,
-3.00000000
1211-5.58992386,
+5.72154379,
+6.00000000
1212+0.179386616,
-0.232047081,
+5.00000000
1213-2.80518818,
+6.00000000,
+0.00000000
1216+3.99999976,
-3.00000024,
-3.00000000
1217-2.00000000,
+0.00000000,
+6.00000000
1218-0.342726707E-6,
-5.00000000,
+5.00000000
1219+1.00000000,
+6.00000000,
+0.00000000
1221+4.00000000,
-3.00000000,
-3.00000000
1222-2.00000000,
+0.00000000,
+6.00000000
1223+0.00000000,
-5.00000000,
+5.00000000
1224+1.00000000,
+6.00000000,
+0.00000000
1226call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
1228+4.00000000,
-3.00000000,
-3.00000000
1229-2.00000000,
+0.00000000,
+6.00000000
1230+0.00000000,
-5.00000000,
+5.00000000
1231+1.00000000,
+6.00000000,
+0.00000000
1233 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
1236+4.00000000,
-3.00000000,
-3.00000000
1237-2.00000000,
+0.00000000,
+6.00000000
1238+0.00000000,
-5.00000000,
+5.00000000
1239+1.00000000,
+6.00000000,
+0.00000000
1242+4.00000000,
-3.00000000,
-3.00000000
1243-2.00000000,
+0.00000000,
+6.00000000
1244+0.00000000,
-5.00000000,
+5.00000000
1245+1.00000000,
+6.00000000,
+0.00000000
1247+4.00000000,
-3.00000000,
-3.00000000
1248-2.00000000,
+0.00000000,
+6.00000000
1249+0.00000000,
-5.00000000,
+5.00000000
1250+1.00000000,
+6.00000000,
+0.00000000
1252call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
1254+4.00000000,
-3.00000000,
-3.00000000
1255-2.00000000,
+0.00000000,
+6.00000000
1256+0.00000000,
-5.00000000,
+5.00000000
1257+1.00000000,
+6.00000000,
+0.00000000
1259 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)
1262+4.00000000,
-3.00000000,
-3.00000000
1263-2.00000000,
+0.00000000,
+6.00000000
1264+0.00000000,
-5.00000000,
+5.00000000
1265+1.00000000,
+6.00000000,
+0.00000000
1268+4.00000000,
-3.00000000,
-3.00000000
1269-2.00000000,
+0.00000000,
+6.00000000
1270+0.00000000,
-5.00000000,
+5.00000000
1271+1.00000000,
+6.00000000,
+0.00000000
1273+4.00000000,
-3.00000000,
-3.00000000
1274-2.00000000,
+0.00000000,
+6.00000000
1275+0.00000000,
-5.00000000,
+5.00000000
1276+1.00000000,
+6.00000000,
+0.00000000
1283-5.00000000,
-1.00000000
1286-2.00000000,
+9.00000000
1287-5.00000000,
+1.00000000
1288+0.00000000,
+4.00000000
1289+1.00000000,
-7.00000000
1294+1.00000000,
+0.873878956
1295+0.873878956,
+1.00000000
1298+5.86491060,
+7.25224209
1299-4.12612104,
-3.36939478
1300+3.49551582,
+4.00000000
1301-5.11715269,
-6.12612104
1303 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
1306+5.86491060,
+7.25224209
1307-4.12612104,
-3.36939478
1308+3.49551582,
+4.00000000
1309-5.11715269,
-6.12612104
1312-2.00000000,
+9.00000000
1313-5.00000095,
+1.00000000
1314+0.00000000,
+4.00000095
1315+1.00000000,
-7.00000191
1317-2.00000000,
+9.00000000
1318-5.00000000,
+1.00000000
1319+0.00000000,
+4.00000000
1320+1.00000000,
-7.00000000
1322call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1324+1.00000000,
+0.873878956
1325+0.00000000,
+1.00000000
1326call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
1328+5.86491060,
+9.00000000
1329-4.12612104,
+1.00000000
1330+3.49551582,
+4.00000000
1331-5.11715269,
-7.00000000
1333 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
1336+5.86491060,
+9.00000000
1337-4.12612104,
+1.00000000
1338+3.49551582,
+4.00000000
1339-5.11715269,
-7.00000000
1342-2.00000000,
+9.00000000
1343-5.00000000,
+1.00000000
1344+0.00000000,
+4.00000000
1345+1.00000000,
-7.00000000
1347-2.00000000,
+9.00000000
1348-5.00000000,
+1.00000000
1349+0.00000000,
+4.00000000
1350+1.00000000,
-7.00000000
1352call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
1354+5.86491060,
+9.00000000
1355-4.12612104,
+1.00000000
1356+3.49551582,
+4.00000000
1357-5.11715269,
-7.00000000
1359 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
1362+5.86491060,
+9.00000000
1363-4.12612104,
+1.00000000
1364+3.49551582,
+4.00000000
1365-5.11715269,
-7.00000000
1368-2.00000000,
+9.00000000
1369-5.00000000,
+1.00000000
1370+0.00000000,
+4.00000000
1371+1.00000000,
-7.00000000
1373-2.00000000,
+9.00000000
1374-5.00000000,
+1.00000000
1375+0.00000000,
+4.00000000
1376+1.00000000,
-7.00000000
1378call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
1380-2.00000000,
+9.00000000
1381-5.00000000,
+1.00000000
1382+0.00000000,
+4.00000000
1383+1.00000000,
-7.00000000
1385 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
1388-2.00000000,
+9.00000000
1389-5.00000000,
+1.00000000
1390+0.00000000,
+4.00000000
1391+1.00000000,
-7.00000000
1394-2.00000000,
+9.00000000
1395-5.00000000,
+1.00000000
1396+0.00000000,
+4.00000000
1397+1.00000000,
-7.00000000
1399-2.00000000,
+9.00000000
1400-5.00000000,
+1.00000000
1401+0.00000000,
+4.00000000
1402+1.00000000,
-7.00000000
1404call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
1406-2.00000000,
+9.00000000
1407-5.00000000,
+1.00000000
1408+0.00000000,
+4.00000000
1409+1.00000000,
-7.00000000
1411 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)
1414-2.00000000,
+9.00000000
1415-5.00000000,
+1.00000000
1416+0.00000000,
+4.00000000
1417+1.00000000,
-7.00000000
1420-2.00000000,
+9.00000000
1421-5.00000000,
+1.00000000
1422+0.00000000,
+4.00000000
1423+1.00000000,
-7.00000000
1425-2.00000000,
+9.00000000
1426-5.00000000,
+1.00000000
1427+0.00000000,
+4.00000000
1428+1.00000000,
-7.00000000
1454 call setAffinity(affinity(isam,:), sample(isam,:), tform, genrecmat)
1473call setMatInit(tform(
2:,
1:ndim
-1), lowDia, vlow
= 0._RKG, vdia
= 0._RKG)
1476call setAffinity(affinity, sample, dim, tform,
class = upperDiag)
1483 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperDiag)
1502call setAffinity(affinity, sample, dim, tform,
class = upperUnit)
1509 call setAffinity(affinity(isam,:), sample(isam,:), tform, upperUnit)
1528call setAffinity(affinity, sample, dim, tform,
class = lowerDiag)
1535 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerDiag)
1554call setAffinity(affinity, sample, dim, tform,
class = lowerUnit)
1561 call setAffinity(affinity(isam,:), sample(isam,:), tform, lowerUnit)