26 integer(IK):: itry, info, ndim, ntry
= 10
27 character(:, SK),
allocatable :: cform, rform
29 type(display_type) :: disp
32 cform
= getFormat(
mold = [(
0._TKG,
0._TKG)], ed
= SK_
"f", ndigit
= 2_IK, signed
= .true._LK)
33 rform
= getFormat(mold
= [
0._TKG], ed
= SK_
"f", ndigit
= 2_IK, signed
= .true._LK)
36 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
37 call disp%show(
"! Compute inverse of an upperDiag triangular real matrix.")
38 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
43 real(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
46 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
50 call disp%show(
"mat = getUnifRand(1_IK, 9_IK, ndim, ndim)")
52 call disp%show(
"call setMatInit(mat(2 : ndim, 1 : ndim - 1), lowDia, 0._TKG, 0._TKG)")
53 call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
55 call disp%show( mat ,
format = rform )
56 call disp%show(
"call setResized(inv, shape(mat, IK))")
58 call disp%show(
"call setMatInv(inv, mat, upperDiag)")
61 call disp%show( inv ,
format = rform )
62 call disp%show(
"mul = matmul(mat, inv)")
63 mul
= matmul(mat, inv)
65 call disp%show( mul ,
format = rform )
71 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
72 call disp%show(
"! Compute inverse of a upperDiag triangular complex matrix.")
73 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
78 complex(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
81 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
85 call disp%show(
"mat = getUnifRand((1., 1.), (2., 2.), ndim, ndim)")
87 call disp%show(
"call setMatInit(mat(2 : ndim, 1 : ndim - 1), lowDia, (0._TKG, 0._TKG), (0._TKG, 0._TKG))")
88 call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
90 call disp%show( mat ,
format = cform )
91 call disp%show(
"call setResized(inv, shape(mat, IK))")
93 call disp%show(
"call setMatInv(inv, mat, upperDiag)")
96 call disp%show( inv ,
format = cform )
97 call disp%show(
"mul = matmul(mat, inv)")
98 mul
= matmul(mat, inv)
100 call disp%show( mul ,
format = cform )
106 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
107 call disp%show(
"! Compute inverse of a lowerDiag triangular real matrix.")
108 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
113 real(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
116 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
118 call disp%show(
"ndim ! matrix rank")
120 call disp%show(
"mat = getUnifRand(1_IK, 9_IK, ndim, ndim)")
122 call disp%show(
"call setMatInit(mat(1 : ndim - 1, 2 : ndim), uppDia, 0._TKG, 0._TKG)")
123 call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
125 call disp%show( mat ,
format = rform )
126 call disp%show(
"call setResized(inv, shape(mat, IK))")
128 call disp%show(
"call setMatInv(inv, mat, lowerDiag)")
131 call disp%show( inv ,
format = rform )
132 call disp%show(
"mul = matmul(mat, inv)")
133 mul
= matmul(mat, inv)
135 call disp%show( mul ,
format = rform )
141 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
142 call disp%show(
"! Compute inverse of a lowerDiag triangular complex matrix.")
143 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
148 complex(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
151 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
153 call disp%show(
"ndim ! matrix rank")
155 call disp%show(
"mat = getUnifRand((1., 1.), (2., 2.), ndim, ndim)")
157 call disp%show(
"call setMatInit(mat(1 : ndim - 1, 2 : ndim), uppDia, (0._TKG, 0._TKG), (0._TKG, 0._TKG))")
158 call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
160 call disp%show( mat ,
format = cform )
161 call disp%show(
"call setResized(inv, shape(mat, IK))")
163 call disp%show(
"call setMatInv(inv, mat, lowerDiag)")
166 call disp%show( inv ,
format = cform )
167 call disp%show(
"mul = matmul(mat, inv)")
168 mul
= matmul(mat, inv)
170 call disp%show( mul ,
format = cform )
176 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
177 call disp%show(
"! Compute inverse of a upperUnit triangular real matrix.")
178 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
183 real(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
186 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
188 call disp%show(
"ndim ! matrix rank")
190 call disp%show(
"mat = getUnifRand(1_IK, 9_IK, ndim, ndim)")
192 call disp%show(
"call setMatInit(mat, lowDia, 0._TKG, 1._TKG)")
195 call disp%show( mat ,
format = rform )
196 call disp%show(
"call setResized(inv, shape(mat, IK))")
198 call disp%show(
"call setMatInv(inv, mat, upperUnit)")
201 call disp%show( inv ,
format = rform )
202 call disp%show(
"mul = matmul(mat, inv)")
203 mul
= matmul(mat, inv)
205 call disp%show( mul ,
format = rform )
211 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
212 call disp%show(
"! Compute inverse of a upperUnit triangular complex matrix.")
213 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
218 complex(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
221 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
223 call disp%show(
"ndim ! matrix rank")
225 call disp%show(
"mat = getUnifRand((-1., -1.), (+1., +1.), ndim, ndim)")
226 mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
227 call disp%show(
"call setMatInit(mat, lowDia, (0._TKG, 0._TKG), (1._TKG, 0._TKG))")
228 call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
230 call disp%show( mat ,
format = cform )
231 call disp%show(
"call setResized(inv, shape(mat, IK))")
233 call disp%show(
"call setMatInv(inv, mat, upperUnit)")
236 call disp%show( inv ,
format = cform )
237 call disp%show(
"mul = matmul(mat, inv)")
238 mul
= matmul(mat, inv)
240 call disp%show( mul ,
format = cform )
246 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
247 call disp%show(
"! Compute inverse of a lowerUnit triangular real matrix.")
248 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
253 real(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
256 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
258 call disp%show(
"ndim ! matrix rank")
260 call disp%show(
"mat = getUnifRand(1_IK, 9_IK, ndim, ndim)")
262 call disp%show(
"call setMatInit(mat, uppDia, 0._TKG, 1._TKG)")
265 call disp%show( mat ,
format = rform )
266 call disp%show(
"call setResized(inv, shape(mat, IK))")
268 call disp%show(
"call setMatInv(inv, mat, lowerUnit)")
271 call disp%show( inv ,
format = rform )
272 call disp%show(
"mul = matmul(mat, inv)")
273 mul
= matmul(mat, inv)
275 call disp%show( mul ,
format = rform )
281 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
282 call disp%show(
"! Compute inverse of a lowerUnit triangular complex matrix.")
283 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
288 complex(TKG),
allocatable :: mat(:,:), inv(:,:), mul(:,:)
291 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
293 call disp%show(
"ndim ! matrix rank")
295 call disp%show(
"mat = getUnifRand((-1., -1.), (+1., +1.), ndim, ndim)")
296 mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
297 call disp%show(
"call setMatInit(mat, uppDia, (0._TKG, 0._TKG), (1._TKG, 0._TKG))")
298 call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
300 call disp%show( mat ,
format = cform )
301 call disp%show(
"call setResized(inv, shape(mat, IK))")
303 call disp%show(
"call setMatInv(inv, mat, lowerUnit)")
306 call disp%show( inv ,
format = cform )
307 call disp%show(
"mul = matmul(mat, inv)")
308 mul
= matmul(mat, inv)
310 call disp%show( mul ,
format = cform )
316 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
317 call disp%show(
"! Compute inverse of a general real matrix by passing its LUP factorization.")
318 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
323 integer(IK),
allocatable :: rperm(:)
324 real(TKG),
allocatable :: mat(:,:), lup(:,:), inv(:,:), mul(:,:)
327 call disp%show(
"mat = reshape([1, 0, 2, -1, 5, 0, 0, 3, -9], shape = [3,3], order = [2, 1])")
328 mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
330 call disp%show( mat ,
format = rform )
331 call disp%show(
"call setResized(rperm, size(mat, 1, IK))")
336 call disp%show( lup ,
format = rform )
337 call disp%show(
"call setMatLUP(lup, rperm, info) ! compute the LUP factorization of the matrix.")
339 call disp%show(
"if (info /= 0) error stop 'LUP factorization failed.'")
340 if (info
/= 0)
error stop 'LUP factorization failed.'
342 call disp%show( lup ,
format = rform )
343 call disp%show(
"call setResized(inv, shape(mat, IK))")
345 call disp%show(
"call setMatInv(inv, lup, rperm) ! compute the inverse of the matrix by passing its LUP factorization.")
348 call disp%show( inv ,
format = rform )
349 call disp%show(
"mul = matmul(mat, inv)")
350 mul
= matmul(mat, inv)
352 call disp%show( mul ,
format = rform )
358 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
359 call disp%show(
"! Compute inverse of a general complex matrix by passing its LUP factorization.")
360 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
365 integer(IK),
allocatable :: rperm(:)
366 complex(TKG),
parameter :: mat(
*,
*)
= reshape(
&
367 [ (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0), (
3.2,
-1.0), (
3.6,
-1.0), (
4.0,
-1.0), (
4.4,
-1.0), (
4.8,
-1.0), (
5.2,
-1.0)
&
368 , (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0), (
3.2,
-1.0), (
3.6,
-1.0), (
4.0,
-1.0), (
4.4,
-1.0), (
4.8,
-1.0)
&
369 , (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0), (
3.2,
-1.0), (
3.6,
-1.0), (
4.0,
-1.0), (
4.4,
-1.0)
&
370 , (
3.2,
1.0), (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0), (
3.2,
-1.0), (
3.6,
-1.0), (
4.0,
-1.0)
&
371 , (
3.6,
1.0), (
3.2,
1.0), (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0), (
3.2,
-1.0), (
3.6,
-1.0)
&
372 , (
4.0,
1.0), (
3.6,
1.0), (
3.2,
1.0), (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0), (
3.2,
-1.0)
&
373 , (
4.4,
1.0), (
4.0,
1.0), (
3.6,
1.0), (
3.2,
1.0), (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0), (
2.8,
-1.0)
&
374 , (
4.8,
1.0), (
4.4,
1.0), (
4.0,
1.0), (
3.6,
1.0), (
3.2,
1.0), (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0), (
2.4,
-1.0)
&
375 , (
5.2,
1.0), (
4.8,
1.0), (
4.4,
1.0), (
4.0,
1.0), (
3.6,
1.0), (
3.2,
1.0), (
2.8,
1.0), (
2.4,
1.0), (
2.0,
1.0)
&
376 ], shape
= [
9,
9], order
= [
2,
1])
377 complex(TKG), dimension(
size(mat,
1),
size(mat,
2)) :: inv, lup, mul
379 call disp%show(
"call setResized(rperm, size(mat, 1, IK))")
384 call disp%show( lup ,
format = cform )
385 call disp%show(
"call setMatLUP(lup, rperm, info) ! compute the LUP factorization of the matrix.")
387 call disp%show(
"if (info /= 0) error stop 'LUP factorization failed.'")
388 if (info
/= 0)
error stop 'LUP factorization failed.'
390 call disp%show( lup ,
format = cform )
391 call disp%show(
"call setMatInv(inv, lup, rperm) ! compute the inverse of the matrix by passing its LUP factorization.")
394 call disp%show( inv ,
format = cform )
395 call disp%show(
"mul = matmul(mat, inv)")
396 mul
= matmul(mat, inv)
398 call disp%show( mul ,
format = cform )
403 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
404 call disp%show(
"! Compute the inverse of a positive-definite real matrix by passing its Cholesky factorization.")
405 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
410 real(TKG),
allocatable :: mat(:,:), chol(:,:), inv(:,:), mul(:,:)
413 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
415 call disp%show(
"ndim ! matrix rank")
417 call disp%show(
"mat = getCovRand(mold = 1._TKG, ndim = ndim)")
420 call disp%show( mat ,
format = rform )
421 call disp%show(
"chol = getMatChol(mat, subset = lowDia)")
424 call disp%show( chol ,
format = rform )
425 call disp%show(
"call setResized(inv, shape(mat, IK))")
427 call disp%show(
"call setMatInv(inv, chol, auxil = choLow)")
428 call setMatInv(inv, chol, auxil
= choLow)
430 call disp%show( inv ,
format = rform )
431 call disp%show(
"mul = matmul(mat, inv)")
432 mul
= matmul(mat, inv)
434 call disp%show( mul ,
format = rform )
440 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
441 call disp%show(
"! Compute the inverse of a positive-definite complex matrix by passing its Cholesky factorization.")
442 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
447 integer(IK),
parameter :: ndim
= 3_IK
448 complex(TKG),
allocatable :: chol(:,:), inv(:,:), mul(:,:)
449 complex(TKG),
parameter :: mat(
*,
*)
= reshape( [ (
9.0,
0.0), (
3.0,
3.0), (
3.0,
-3.0)
&
450 , (
3.0,
-3.0),(
18.0,
0.0), (
8.0,
-6.0)
&
451 , (
3.0,
3.0), (
8.0,
6.0),(
43.0,
0.0)
&
452 ], shape
= [ndim, ndim], order
= [
2,
1])
455 call disp%show( mat ,
format = cform )
456 call disp%show(
"chol = getMatChol(mat, subset = lowDia)")
459 call disp%show( chol ,
format = cform )
460 call disp%show(
"call setResized(inv, shape(mat, IK))")
462 call disp%show(
"call setMatInv(inv, chol, auxil = choLow)")
463 call setMatInv(inv, chol, auxil
= choLow)
465 call disp%show( inv ,
format = cform )
466 call disp%show(
"mul = matmul(mat, inv)")
467 mul
= matmul(mat, inv)
469 call disp%show( mul ,
format = cform )
475 integer(IK),
parameter :: ndim
= 3_IK
476 complex(TKG),
allocatable :: chol(:,:), inv(:,:), mul(:,:)
477 complex(TKG),
parameter :: mat(
*,
*)
= reshape( [ (
9.0,
0.0), (
3.0,
3.0), (
3.0,
-3.0)
&
478 , (
3.0,
-3.0),(
18.0,
0.0), (
8.0,
-6.0)
&
479 , (
3.0,
3.0), (
8.0,
6.0),(
43.0,
0.0)
&
480 ], shape
= [ndim, ndim], order
= [
2,
1])
483 call disp%show( mat ,
format = cform )
484 call disp%show(
"chol = getMatChol(mat, subset = uppDia)")
487 call disp%show( chol ,
format = cform )
488 call disp%show(
"call setResized(inv, shape(mat, IK))")
490 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp)")
491 call setMatInv(inv, chol, auxil
= choUpp)
493 call disp%show( inv ,
format = cform )
494 call disp%show(
"mul = matmul(mat, inv)")
495 mul
= matmul(mat, inv)
497 call disp%show( mul ,
format = cform )
502 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
503 call disp%show(
"! Compute a subset of the inverse of a positive-definite real matrix by passing its Cholesky factorization.")
504 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
509 real(TKG),
allocatable :: mat(:,:), chol(:,:), inv(:,:), mul(:,:)
512 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
514 call disp%show(
"ndim ! matrix rank")
516 call disp%show(
"mat = getCovRand(mold = 1._TKG, ndim = ndim)")
519 call disp%show( mat ,
format = rform )
520 call disp%show(
"chol = getMatChol(mat, subset = lowDia)")
523 call disp%show( chol ,
format = rform )
524 call disp%show(
"call setResized(inv, shape(mat, IK))")
529 call disp%show(
"call setMatInv(inv, chol, auxil = choLow, subset = uppDia)")
530 call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
532 call disp%show( inv ,
format = rform )
533 call disp%show(
"call setMatCopy(inv(2:,1:), rdpack, inv(1:,2:), rdpack, uppDia, transHerm) ! symmetrize `inv` for multiplication.")
534 call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
536 call disp%show( inv ,
format = rform )
537 call disp%show(
"mul = matmul(mat, inv)")
538 mul
= matmul(mat, inv)
540 call disp%show( mul ,
format = rform )
544 call disp%show(
"call setMatInv(inv, chol, auxil = choLow, subset = lowDia)")
545 call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
547 call disp%show( inv ,
format = rform )
548 call disp%show(
"call setMatCopy(inv(1:,2:), rdpack, inv(2:,1:), rdpack, lowDia, transHerm) ! symmetrize `inv` for multiplication.")
549 call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
551 call disp%show( inv ,
format = rform )
552 call disp%show(
"mul = matmul(mat, inv)")
553 mul
= matmul(mat, inv)
555 call disp%show( mul ,
format = rform )
562 real(TKG),
allocatable :: mat(:,:), chol(:,:), inv(:,:), mul(:,:)
565 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
567 call disp%show(
"ndim ! matrix rank")
569 call disp%show(
"mat = getCovRand(mold = 1._TKG, ndim = ndim)")
572 call disp%show( mat ,
format = rform )
573 call disp%show(
"chol = getMatChol(mat, subset = uppDia)")
576 call disp%show( chol ,
format = rform )
577 call disp%show(
"call setResized(inv, shape(mat, IK))")
582 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp, subset = uppDia)")
583 call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
585 call disp%show( inv ,
format = rform )
586 call disp%show(
"call setMatCopy(inv(2:,1:), rdpack, inv(1:,2:), rdpack, uppDia, transHerm) ! symmetrize `inv` for multiplication.")
587 call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
589 call disp%show( inv ,
format = rform )
590 call disp%show(
"mul = matmul(mat, inv)")
591 mul
= matmul(mat, inv)
593 call disp%show( mul ,
format = rform )
597 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp, subset = lowDia)")
598 call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
600 call disp%show( inv ,
format = rform )
601 call disp%show(
"call setMatCopy(inv(1:,2:), rdpack, inv(2:,1:), rdpack, lowDia, transHerm) ! symmetrize `inv` for multiplication.")
602 call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
604 call disp%show( inv ,
format = rform )
605 call disp%show(
"mul = matmul(mat, inv)")
606 mul
= matmul(mat, inv)
608 call disp%show( mul ,
format = rform )
614 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
615 call disp%show(
"! Compute a subset of the inverse of a positive-definite complex matrix by passing its Cholesky factorization.")
616 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
621 integer(IK),
parameter :: ndim
= 3_IK
622 complex(TKG),
allocatable :: chol(:,:), inv(:,:), mul(:,:)
623 complex(TKG),
parameter :: mat(
*,
*)
= reshape( [ (
9.0,
0.0), (
3.0,
3.0), (
3.0,
-3.0)
&
624 , (
3.0,
-3.0),(
18.0,
0.0), (
8.0,
-6.0)
&
625 , (
3.0,
3.0), (
8.0,
6.0),(
43.0,
0.0)
&
626 ], shape
= [ndim, ndim], order
= [
2,
1])
629 call disp%show( mat ,
format = cform )
630 call disp%show(
"chol = getMatChol(mat, subset = lowDia)")
633 call disp%show( chol ,
format = cform )
634 call disp%show(
"call setResized(inv, shape(mat, IK))")
638 call disp%show(
"call setMatInv(inv, chol, auxil = choLow, subset = uppDia)")
639 call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
641 call disp%show( inv ,
format = cform )
642 call disp%show(
"call setMatCopy(inv(2:,1:), rdpack, inv(1:,2:), rdpack, uppDia, transHerm) ! symmetrize `inv` for multiplication.")
643 call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
645 call disp%show( inv ,
format = cform )
646 call disp%show(
"mul = matmul(mat, inv)")
647 mul
= matmul(mat, inv)
649 call disp%show( mul ,
format = cform )
653 call disp%show(
"call setMatInv(inv, chol, auxil = choLow, subset = lowDia)")
654 call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
656 call disp%show( inv ,
format = cform )
657 call disp%show(
"call setMatCopy(inv(1:,2:), rdpack, inv(2:,1:), rdpack, lowDia, transHerm) ! symmetrize `inv` for multiplication.")
658 call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
660 call disp%show( inv ,
format = cform )
661 call disp%show(
"mul = matmul(mat, inv)")
662 mul
= matmul(mat, inv)
664 call disp%show( mul ,
format = cform )
670 integer(IK),
parameter :: ndim
= 3_IK
671 complex(TKG),
allocatable :: chol(:,:), inv(:,:), mul(:,:)
672 complex(TKG),
parameter :: mat(
*,
*)
= reshape( [ (
9.0,
0.0), (
3.0,
3.0), (
3.0,
-3.0)
&
673 , (
3.0,
-3.0),(
18.0,
0.0), (
8.0,
-6.0)
&
674 , (
3.0,
3.0), (
8.0,
6.0),(
43.0,
0.0)
&
675 ], shape
= [ndim, ndim], order
= [
2,
1])
678 call disp%show( mat ,
format = cform )
679 call disp%show(
"chol = getMatChol(mat, subset = uppDia)")
682 call disp%show( chol ,
format = cform )
683 call disp%show(
"call setResized(inv, shape(mat, IK))")
687 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp, subset = lowDia)")
688 call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
690 call disp%show( inv ,
format = cform )
691 call disp%show(
"call setMatCopy(inv(1:,2:), rdpack, inv(2:,1:), rdpack, lowDia, transHerm) ! symmetrize `inv` for multiplication.")
692 call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
694 call disp%show( inv ,
format = cform )
695 call disp%show(
"mul = matmul(mat, inv)")
696 mul
= matmul(mat, inv)
698 call disp%show( mul ,
format = cform )
702 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp, subset = uppDia)")
703 call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
705 call disp%show( inv ,
format = cform )
706 call disp%show(
"call setMatCopy(inv(2:,1:), rdpack, inv(1:,2:), rdpack, uppDia, transHerm) ! symmetrize `inv` for multiplication.")
707 call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
709 call disp%show( inv ,
format = cform )
710 call disp%show(
"mul = matmul(mat, inv)")
711 mul
= matmul(mat, inv)
713 call disp%show( mul ,
format = cform )
719 complex(TKG),
allocatable :: low(:,:), upp(:,:), mat(:,:), chol(:,:), inv(:,:), mul(:,:)
722 call disp%show(
"ndim = getUnifRand(1_IK, 6_IK)")
724 call disp%show(
"ndim ! matrix rank")
726 low
= getUnifRand((
1._TKG,
-1._TKG), (
+2._TKG,
+1._TKG), ndim, ndim)
727 call setMatInit(low, uppDia, (
0._TKG,
0._TKG),
cmplx(
getUnifRand(
1._TKG,
2._TKG, ndim),
0._TKG, TKG))
728 upp
= transpose(
conjg(low))
729 mat
= matmul(low, upp)
731 call disp%show( mat ,
format = cform )
732 call disp%show(
"chol = getMatChol(mat, subset = uppDia)")
735 call disp%show( chol ,
format = cform )
736 call disp%show(
"call setResized(inv, shape(mat, IK))")
740 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp, subset = lowDia)")
741 call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
743 call disp%show( inv ,
format = cform )
744 call disp%show(
"call setMatCopy(inv(1:,2:), rdpack, inv(2:,1:), rdpack, lowDia, transHerm) ! symmetrize `inv` for multiplication.")
745 call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
747 call disp%show( inv ,
format = cform )
748 call disp%show(
"mul = matmul(mat, inv)")
749 mul
= matmul(mat, inv)
751 call disp%show( mul ,
format = cform )
755 call disp%show(
"call setMatInv(inv, chol, auxil = choUpp, subset = uppDia)")
756 call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
758 call disp%show( inv ,
format = cform )
759 call disp%show(
"call setMatCopy(inv(2:,1:), rdpack, inv(1:,2:), rdpack, uppDia, transHerm) ! symmetrize `inv` for multiplication.")
760 call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
762 call disp%show( inv ,
format = cform )
763 call disp%show(
"mul = matmul(mat, inv)")
764 mul
= matmul(mat, inv)
766 call disp%show( mul ,
format = cform )
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 and return a random positive-definite (correlation or covariance) matrix using the Gram meth...
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
Verify the input assertion holds and if it does not, print the (optional) input message on stdout and...
Generate and return an object of type stop_type with the user-specified input attributes.
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...
Copy a desired subset of the input source matrix of arbitrary shape (:) or (:,:) to the target subset...
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 (:,...
Return the LU-Pivoted decomposition of the input square matrix mat(ndim,ndim).
Generate and return the conversion of the input value to an output Fortran string,...
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 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 reporting and handling errors.
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 procedures and generic interfaces for computing the Cholesky factorization of po...
This module contains procedures and generic interfaces relevant to copying (diagonal or upper/lower t...
This module contains procedures and generic interfaces relevant to generating and initializing matric...
This module contains procedures and generic interfaces relevant to the partially LU Pivoted decomposi...
This module contains the generic procedures for converting values of different types and kinds to For...
Generate and return an object of type display_type.
11call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
30call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
46call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
48 +7.00,
+9.00,
+2.00,
+4.00,
+3.00
49 +0.00,
+6.00,
+1.00,
+8.00,
+1.00
50 +0.00,
+0.00,
+5.00,
+6.00,
+1.00
51 +0.00,
+0.00,
+0.00,
+8.00,
+7.00
52 +0.00,
+0.00,
+0.00,
+0.00,
+4.00
56 +0.14,
-0.21,
-0.01,
+0.15,
-0.32
57 +0.00,
+0.17,
-0.03,
-0.14,
+0.21
58 +0.00,
+0.00,
+0.20,
-0.15,
+0.21
59 +0.00,
+0.00,
+0.00,
+0.12,
-0.22
60 +0.00,
+0.00,
+0.00,
+0.00,
+0.25
63 +1.00,
-0.00,
+0.00,
+0.00,
-0.00
64 +0.00,
+1.00,
-0.00,
+0.00,
+0.00
65 +0.00,
+0.00,
+1.00,
+0.00,
+0.00
66 +0.00,
+0.00,
+0.00,
+1.00,
+0.00
67 +0.00,
+0.00,
+0.00,
+0.00,
+1.00
74call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
76 +6.00,
+9.00,
+3.00,
+3.00
77 +0.00,
+7.00,
+4.00,
+7.00
78 +0.00,
+0.00,
+5.00,
+4.00
79 +0.00,
+0.00,
+0.00,
+3.00
83 +0.17,
-0.21,
+0.07,
+0.24
84 +0.00,
+0.14,
-0.11,
-0.18
85 +0.00,
+0.00,
+0.20,
-0.27
86 +0.00,
+0.00,
+0.00,
+0.33
89 +1.00,
+0.00,
+0.00,
+0.00
90 +0.00,
+1.00,
-0.00,
+0.00
91 +0.00,
+0.00,
+1.00,
+0.00
92 +0.00,
+0.00,
+0.00,
+1.00
99call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
101 +3.00,
+7.00,
+2.00,
+6.00,
+3.00,
+7.00
102 +0.00,
+9.00,
+7.00,
+6.00,
+8.00,
+1.00
103 +0.00,
+0.00,
+5.00,
+7.00,
+8.00,
+3.00
104 +0.00,
+0.00,
+0.00,
+1.00,
+5.00,
+4.00
105 +0.00,
+0.00,
+0.00,
+0.00,
+1.00,
+2.00
106 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+8.00
110 +0.33,
-0.26,
+0.23,
-2.05,
+9.50,
-1.69
111 +0.00,
+0.11,
-0.16,
+0.42,
-1.76,
+0.27
112 +0.00,
+0.00,
+0.20,
-1.40,
+5.40,
-0.73
113 +0.00,
+0.00,
+0.00,
+1.00,
-5.00,
+0.75
114 +0.00,
+0.00,
+0.00,
+0.00,
+1.00,
-0.25
115 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+0.12
116mul
= matmul(mat, inv)
118 +1.00,
-0.00,
+0.00,
+0.00,
+0.00,
-0.00
119 +0.00,
+1.00,
-0.00,
-0.00,
+0.00,
-0.00
120 +0.00,
+0.00,
+1.00,
+0.00,
+0.00,
+0.00
121 +0.00,
+0.00,
+0.00,
+1.00,
+0.00,
+0.00
122 +0.00,
+0.00,
+0.00,
+0.00,
+1.00,
+0.00
123 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+1.00
130call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
132 +5.00,
+7.00,
+3.00,
+5.00,
+8.00
133 +0.00,
+6.00,
+3.00,
+1.00,
+4.00
134 +0.00,
+0.00,
+8.00,
+7.00,
+7.00
135 +0.00,
+0.00,
+0.00,
+9.00,
+6.00
136 +0.00,
+0.00,
+0.00,
+0.00,
+7.00
140 +0.20,
-0.23,
+0.01,
-0.09,
-0.03
141 +0.00,
+0.17,
-0.06,
+0.03,
-0.06
142 +0.00,
+0.00,
+0.12,
-0.10,
-0.04
143 +0.00,
+0.00,
+0.00,
+0.11,
-0.10
144 +0.00,
+0.00,
+0.00,
+0.00,
+0.14
145mul
= matmul(mat, inv)
147 +1.00,
+0.00,
-0.00,
+0.00,
+0.00
148 +0.00,
+1.00,
+0.00,
-0.00,
+0.00
149 +0.00,
+0.00,
+1.00,
+0.00,
+0.00
150 +0.00,
+0.00,
+0.00,
+1.00,
+0.00
151 +0.00,
+0.00,
+0.00,
+0.00,
+1.00
158call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
169mul
= matmul(mat, inv)
180call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
191mul
= matmul(mat, inv)
202call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
213mul
= matmul(mat, inv)
224call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia,
0._TKG,
0._TKG)
226 +2.00,
+9.00,
+9.00,
+1.00
227 +0.00,
+7.00,
+4.00,
+4.00
228 +0.00,
+0.00,
+2.00,
+1.00
229 +0.00,
+0.00,
+0.00,
+7.00
233 +0.50,
-0.64,
-0.96,
+0.43
234 +0.00,
+0.14,
-0.29,
-0.04
235 +0.00,
+0.00,
+0.50,
-0.07
236 +0.00,
+0.00,
+0.00,
+0.14
237mul
= matmul(mat, inv)
239 +1.00,
+0.00,
+0.00,
-0.00
240 +0.00,
+1.00,
+0.00,
+0.00
241 +0.00,
+0.00,
+1.00,
+0.00
242 +0.00,
+0.00,
+0.00,
+1.00
254call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
261mul
= matmul(mat, inv)
270call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
277mul
= matmul(mat, inv)
286call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
293mul
= matmul(mat, inv)
302call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
304(
+1.35,
+1.60), (
+1.87,
+1.77)
305(
+0.00,
+0.00), (
+1.42,
+1.20)
309(
+0.31,
-0.36), (
-0.45,
+0.48)
310(
+0.00,
+0.00), (
+0.41,
-0.35)
311mul
= matmul(mat, inv)
313(
+1.00,
+0.00), (
+0.00,
-0.00)
314(
+0.00,
+0.00), (
+1.00,
+0.00)
321call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
323(
+1.52,
+1.72), (
+1.59,
+1.66), (
+1.27,
+1.72), (
+1.84,
+1.17), (
+1.82,
+1.09)
324(
+0.00,
+0.00), (
+1.70,
+1.45), (
+1.08,
+1.77), (
+1.16,
+1.75), (
+1.11,
+1.45)
325(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.41,
+1.96), (
+1.22,
+1.04), (
+1.12,
+1.97)
326(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.34,
+1.20), (
+1.71,
+1.51)
327(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.38,
+1.20)
331(
+0.29,
-0.33), (
-0.33,
+0.30), (
+0.05,
+0.05), (
+0.13,
+0.15), (
-0.08,
-0.08)
332(
+0.00,
+0.00), (
+0.34,
-0.29), (
-0.31,
+0.23), (
-0.20,
+0.02), (
+0.31,
+0.06)
333(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.24,
-0.34), (
-0.21,
+0.30), (
-0.16,
-0.08)
334(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.41,
-0.37), (
-0.52,
+0.46)
335(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.41,
-0.36)
336mul
= matmul(mat, inv)
338(
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
-0.00)
339(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
340(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00)
341(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
342(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
349call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
356mul
= matmul(mat, inv)
365call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
372mul
= matmul(mat, inv)
381call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
383(
+1.21,
+1.81), (
+1.24,
+1.08), (
+1.74,
+1.67), (
+1.68,
+1.36), (
+1.68,
+1.82), (
+1.48,
+1.24)
384(
+0.00,
+0.00), (
+1.46,
+1.53), (
+1.52,
+1.20), (
+1.09,
+1.63), (
+1.70,
+1.38), (
+1.22,
+1.45)
385(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.31,
+1.66), (
+1.55,
+1.03), (
+1.61,
+1.21), (
+1.44,
+1.64)
386(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.91,
+1.66), (
+1.60,
+1.93), (
+1.85,
+1.55)
387(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.43,
+1.08), (
+1.74,
+1.62)
388(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.82,
+1.40)
392(
+0.26,
-0.38), (
-0.17,
+0.31), (
-0.15,
+0.16), (
+0.08,
+0.01), (
-0.16,
-0.13), (
+0.20,
+0.11)
393(
+0.00,
+0.00), (
+0.33,
-0.34), (
-0.22,
+0.37), (
-0.19,
-0.11), (
+0.06,
+0.19), (
+0.11,
-0.17)
394(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.29,
-0.37), (
-0.18,
+0.30), (
+0.01,
+0.07), (
-0.14,
-0.09)
395(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.30,
-0.26), (
-0.49,
+0.26), (
+0.20,
+0.05)
396(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.45,
-0.34), (
-0.49,
+0.30)
397(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.35,
-0.27)
398mul
= matmul(mat, inv)
400(
+1.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
401(
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
402(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00)
403(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
404(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
405(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
412call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
414(
+1.85,
+1.11), (
+1.95,
+1.50), (
+1.24,
+1.55), (
+1.43,
+1.41)
415(
+0.00,
+0.00), (
+1.62,
+1.98), (
+1.00,
+1.65), (
+1.79,
+1.38)
416(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.64,
+1.13), (
+1.14,
+1.87)
417(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.07,
+1.06)
421(
+0.40,
-0.24), (
-0.32,
+0.31), (
-0.04,
-0.03), (
-0.06,
-0.14)
422(
+0.00,
+0.00), (
+0.25,
-0.30), (
-0.34,
+0.17), (
+0.22,
+0.38)
423(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.41,
-0.28), (
-0.68,
+0.26)
424(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.47,
-0.47)
425mul
= matmul(mat, inv)
427(
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
-0.00)
428(
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
429(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
-0.00)
430(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
437call setMatInit(mat(
2 : ndim,
1 : ndim
- 1), lowDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
439(
+1.24,
+1.36), (
+1.39,
+1.85), (
+1.38,
+1.12), (
+1.90,
+1.62), (
+1.92,
+1.42), (
+1.14,
+1.85)
440(
+0.00,
+0.00), (
+1.52,
+1.43), (
+1.04,
+1.28), (
+1.44,
+1.90), (
+1.01,
+1.47), (
+1.57,
+1.81)
441(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.22,
+1.44), (
+1.82,
+1.80), (
+1.90,
+1.10), (
+1.17,
+1.05)
442(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+2.00,
+1.37), (
+1.62,
+1.09), (
+1.75,
+1.60)
443(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+1.52), (
+2.00,
+1.69)
444(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.83,
+1.81)
448(
+0.37,
-0.40), (
-0.48,
+0.37), (
+0.15,
+0.12), (
+0.00,
+0.02), (
-0.05,
+0.22), (
+0.02,
-0.37)
449(
+0.00,
+0.00), (
+0.35,
-0.33), (
-0.31,
+0.28), (
-0.05,
-0.04), (
-0.08,
-0.13), (
-0.02,
+0.25)
450(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.34,
-0.41), (
-0.43,
+0.36), (
+0.13,
+0.09), (
+0.05,
-0.17)
451(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.34,
-0.23), (
-0.24,
+0.37), (
-0.09,
-0.17)
452(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.30,
-0.46), (
-0.27,
+0.49)
453(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.28,
-0.27)
454mul
= matmul(mat, inv)
456(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
+0.00)
457(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
458(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00)
459(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00)
460(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
461(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
473call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
482mul
= matmul(mat, inv)
492call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
494 +4.00,
+0.00,
+0.00,
+0.00
495 +6.00,
+5.00,
+0.00,
+0.00
496 +5.00,
+2.00,
+6.00,
+0.00
497 +4.00,
+1.00,
+6.00,
+7.00
501 +0.25,
+0.00,
+0.00,
+0.00
502 -0.30,
+0.20,
+0.00,
+0.00
503 -0.11,
-0.07,
+0.17,
+0.00
504 -0.01,
+0.03,
-0.14,
+0.14
505mul
= matmul(mat, inv)
507 +1.00,
+0.00,
+0.00,
+0.00
508 +0.00,
+1.00,
+0.00,
+0.00
509 +0.00,
-0.00,
+1.00,
+0.00
510 +0.00,
+0.00,
+0.00,
+1.00
517call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
524mul
= matmul(mat, inv)
533call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
544mul
= matmul(mat, inv)
555call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
562mul
= matmul(mat, inv)
571call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
580mul
= matmul(mat, inv)
590call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
597mul
= matmul(mat, inv)
606call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
615mul
= matmul(mat, inv)
625call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
627 +2.00,
+0.00,
+0.00,
+0.00
628 +5.00,
+3.00,
+0.00,
+0.00
629 +8.00,
+6.00,
+1.00,
+0.00
630 +4.00,
+5.00,
+3.00,
+7.00
634 +0.50,
+0.00,
+0.00,
+0.00
635 -0.83,
+0.33,
+0.00,
+0.00
636 +1.00,
-2.00,
+1.00,
+0.00
637 -0.12,
+0.62,
-0.43,
+0.14
638mul
= matmul(mat, inv)
640 +1.00,
+0.00,
+0.00,
+0.00
641 +0.00,
+1.00,
+0.00,
+0.00
642 +0.00,
+0.00,
+1.00,
+0.00
643 -0.00,
+0.00,
-0.00,
+1.00
650call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia,
0._TKG,
0._TKG)
659mul
= matmul(mat, inv)
674call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
681mul
= matmul(mat, inv)
690call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
697mul
= matmul(mat, inv)
706call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
708(
+1.33,
+1.66), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
709(
+1.73,
+1.11), (
+1.15,
+1.61), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
710(
+1.79,
+1.84), (
+1.41,
+1.63), (
+1.22,
+1.68), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
711(
+1.93,
+1.62), (
+1.68,
+1.30), (
+1.44,
+1.40), (
+1.25,
+1.33), (
+0.00,
+0.00), (
+0.00,
+0.00)
712(
+1.38,
+1.76), (
+1.55,
+1.03), (
+1.14,
+1.10), (
+1.26,
+1.15), (
+1.62,
+1.94), (
+0.00,
+0.00)
713(
+1.92,
+1.77), (
+1.70,
+1.73), (
+1.38,
+1.78), (
+1.18,
+1.76), (
+1.94,
+1.53), (
+1.81,
+1.76)
717(
+0.29,
-0.37), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
718(
-0.14,
+0.47), (
+0.29,
-0.41), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
719(
-0.19,
+0.01), (
-0.27,
+0.45), (
+0.28,
-0.39), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
720(
-0.06,
-0.03), (
+0.01,
+0.02), (
-0.29,
+0.44), (
+0.38,
-0.40), (
+0.00,
+0.00), (
+0.00,
+0.00)
721(
-0.11,
-0.05), (
+0.01,
+0.04), (
+0.01,
-0.06), (
-0.22,
+0.30), (
+0.25,
-0.30), (
+0.00,
+0.00)
722(
+0.18,
+0.03), (
-0.01,
-0.03), (
+0.02,
+0.06), (
-0.20,
-0.05), (
-0.22,
+0.32), (
+0.28,
-0.28)
723mul
= matmul(mat, inv)
725(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
726(
-0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
727(
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
728(
+0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
729(
-0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
730(
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
-0.00)
737call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
739(
+1.19,
+1.09), (
+0.00,
+0.00)
740(
+1.50,
+1.06), (
+1.97,
+1.03)
744(
+0.46,
-0.42), (
+0.00,
+0.00)
745(
-0.42,
+0.29), (
+0.40,
-0.21)
746mul
= matmul(mat, inv)
748(
+1.00,
+0.00), (
+0.00,
+0.00)
749(
+0.00,
+0.00), (
+1.00,
-0.00)
756call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
758(
+1.70,
+1.45), (
+0.00,
+0.00)
759(
+1.86,
+1.93), (
+1.41,
+1.57)
763(
+0.34,
-0.29), (
+0.00,
+0.00)
764(
-0.42,
+0.38), (
+0.32,
-0.35)
765mul
= matmul(mat, inv)
767(
+1.00,
+0.00), (
+0.00,
+0.00)
768(
+0.00,
+0.00), (
+1.00,
+0.00)
775call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
777(
+1.70,
+1.35), (
+0.00,
+0.00)
778(
+1.78,
+1.51), (
+1.64,
+1.63)
782(
+0.36,
-0.29), (
+0.00,
+0.00)
783(
-0.34,
+0.32), (
+0.31,
-0.30)
784mul
= matmul(mat, inv)
786(
+1.00,
+0.00), (
+0.00,
+0.00)
787(
+0.00,
+0.00), (
+1.00,
-0.00)
794call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
796(
+1.96,
+1.20), (
+0.00,
+0.00), (
+0.00,
+0.00)
797(
+1.87,
+1.57), (
+1.81,
+1.64), (
+0.00,
+0.00)
798(
+1.16,
+1.13), (
+1.57,
+1.32), (
+1.66,
+1.55)
802(
+0.37,
-0.23), (
+0.00,
+0.00), (
+0.00,
+0.00)
803(
-0.36,
+0.24), (
+0.30,
-0.27), (
+0.00,
+0.00)
804(
+0.05,
-0.08), (
-0.26,
+0.26), (
+0.32,
-0.30)
805mul
= matmul(mat, inv)
807(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
808(
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
+0.00)
809(
+0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
+0.00)
816call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
823mul
= matmul(mat, inv)
832call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
839mul
= matmul(mat, inv)
848call setMatInit(mat(
1 : ndim
- 1,
2 : ndim), uppDia, (
0._TKG,
0._TKG), (
0._TKG,
0._TKG))
850(
+1.01,
+1.91), (
+0.00,
+0.00)
851(
+1.13,
+1.04), (
+1.07,
+1.38)
855(
+0.22,
-0.41), (
+0.00,
+0.00)
856(
-0.13,
+0.39), (
+0.35,
-0.45)
857mul
= matmul(mat, inv)
859(
+1.00,
-0.00), (
+0.00,
+0.00)
860(
+0.00,
+0.00), (
+1.00,
+0.00)
879mul
= matmul(mat, inv)
897mul
= matmul(mat, inv)
916mul
= matmul(mat, inv)
933mul
= matmul(mat, inv)
953mul
= matmul(mat, inv)
973mul
= matmul(mat, inv)
994mul
= matmul(mat, inv)
1007 +1.00,
+5.00,
+8.00,
+5.00,
+7.00
1008 +0.00,
+1.00,
+1.00,
+7.00,
+7.00
1009 +0.00,
+0.00,
+1.00,
+7.00,
+7.00
1010 +0.00,
+0.00,
+0.00,
+1.00,
+9.00
1011 +0.00,
+0.00,
+0.00,
+0.00,
+1.00
1015 +1.00,
-5.00,
-3.00,
+51.00,
-410.00
1016 +0.00,
+1.00,
-1.00,
-0.00,
-0.00
1017 +0.00,
+0.00,
+1.00,
-7.00,
+56.00
1018 +0.00,
+0.00,
+0.00,
+1.00,
-9.00
1019 +0.00,
+0.00,
+0.00,
+0.00,
+1.00
1020mul
= matmul(mat, inv)
1022 +1.00,
+0.00,
+0.00,
+0.00,
+0.00
1023 +0.00,
+1.00,
+0.00,
+0.00,
+0.00
1024 +0.00,
+0.00,
+1.00,
+0.00,
+0.00
1025 +0.00,
+0.00,
+0.00,
+1.00,
+0.00
1026 +0.00,
+0.00,
+0.00,
+0.00,
+1.00
1035 +1.00,
+7.00,
+9.00,
+6.00,
+9.00,
+6.00
1036 +0.00,
+1.00,
+6.00,
+4.00,
+2.00,
+7.00
1037 +0.00,
+0.00,
+1.00,
+9.00,
+9.00,
+4.00
1038 +0.00,
+0.00,
+0.00,
+1.00,
+7.00,
+7.00
1039 +0.00,
+0.00,
+0.00,
+0.00,
+1.00,
+1.00
1040 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+1.00
1044 +1.00,
-7.00,
+33.00,
-275.00,
+1633.00,
+203.00
1045 +0.00,
+1.00,
-6.00,
+50.00,
-298.00,
-35.00
1046 +0.00,
+0.00,
+1.00,
-9.00,
+54.00,
+5.00
1047 +0.00,
+0.00,
+0.00,
+1.00,
-7.00,
-0.00
1048 +0.00,
+0.00,
+0.00,
+0.00,
+1.00,
-1.00
1049 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+1.00
1050mul
= matmul(mat, inv)
1052 +1.00,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
1053 +0.00,
+1.00,
+0.00,
+0.00,
+0.00,
+0.00
1054 +0.00,
+0.00,
+1.00,
+0.00,
+0.00,
+0.00
1055 +0.00,
+0.00,
+0.00,
+1.00,
+0.00,
+0.00
1056 +0.00,
+0.00,
+0.00,
+0.00,
+1.00,
+0.00
1057 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+1.00
1075mul
= matmul(mat, inv)
1090mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1091call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1093(
+1.00,
+0.00), (
-0.90,
+0.64), (
+0.91,
-0.30), (
-0.46,
-0.88), (
-0.27,
+0.89), (
-0.53,
-0.47)
1094(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.86,
+0.77), (
-0.61,
+0.95), (
-0.04,
-0.57), (
+0.83,
+0.86)
1095(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.39,
-0.06), (
+0.44,
-0.99), (
+0.86,
+0.23)
1096(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.41,
-0.48), (
-0.73,
+0.04)
1097(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.15,
-0.76)
1098(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1102(
+1.00,
+0.00), (
+0.90,
-0.64), (
-2.18,
+0.16), (
-0.46,
-0.43), (
+1.86,
-2.69), (
+3.13,
+1.29)
1103(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.86,
-0.77), (
+0.32,
-1.30), (
+1.68,
+0.73), (
-0.38,
+0.43)
1104(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.39,
+0.06), (
-0.62,
+1.16), (
-1.55,
-0.50)
1105(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.41,
+0.48), (
+0.30,
-0.27)
1106(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.15,
+0.76)
1107(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1108mul
= matmul(mat, inv)
1110(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
-0.00)
1111(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
1112(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1113(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1114(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1115(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1121mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1122call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1124(
+1.00,
+0.00), (
+0.89,
+0.08), (
+0.06,
+0.02), (
-0.88,
-0.59), (
-0.35,
+0.52)
1125(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.80,
+0.06), (
-0.18,
+0.87), (
-0.15,
+0.37)
1126(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.84,
-0.13), (
-0.50,
-0.46)
1127(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.91,
+0.68)
1128(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1132(
+1.00,
+0.00), (
-0.89,
-0.08), (
-0.78,
-0.02), (
+1.32,
+1.26), (
-0.54,
-2.62)
1133(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.80,
-0.06), (
-0.49,
-0.71), (
+0.55,
+0.95)
1134(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.84,
+0.13), (
+1.36,
+0.91)
1135(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.91,
-0.68)
1136(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1137mul
= matmul(mat, inv)
1139(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
1140(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
-0.00)
1141(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1142(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1143(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1149mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1150call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1152(
+1.00,
+0.00), (
-0.72,
+0.40), (
+0.07,
+0.51), (
+0.12,
+0.93), (
+0.68,
+0.66), (
-0.38,
-0.11)
1153(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.80,
-0.80), (
+0.92,
-0.78), (
+0.15,
+0.61), (
-0.17,
-0.14)
1154(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.53,
-0.32), (
-0.65,
-0.67), (
-0.22,
+0.33)
1155(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.76,
-0.07), (
-0.57,
+0.62)
1156(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.44,
+0.85)
1157(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1161(
+1.00,
+0.00), (
+0.72,
-0.40), (
+0.83,
-0.26), (
-0.84,
+0.41), (
+0.28,
-1.02), (
-0.31,
-0.14)
1162(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.80,
+0.80), (
-1.60,
+0.61), (
+1.01,
-0.13), (
+0.42,
+0.47)
1163(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.53,
+0.32), (
+1.03,
+0.39), (
+0.91,
-0.53)
1164(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.76,
+0.07), (
+0.29,
+0.06)
1165(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.44,
-0.85)
1166(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1167mul
= matmul(mat, inv)
1169(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
-0.00,
-0.00)
1170(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00)
1171(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
1172(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1173(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1174(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1180mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1181call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1183(
+1.00,
+0.00), (
+0.00,
+0.51), (
-0.34,
+0.44), (
+0.93,
-0.98), (
+0.87,
-0.36)
1184(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.30,
+0.87), (
+0.39,
-0.70), (
+0.79,
+0.38)
1185(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.22,
-0.66), (
-0.40,
-0.32)
1186(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.04,
+0.69)
1187(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1191(
+1.00,
+0.00), (
-0.00,
-0.51), (
-0.10,
-0.28), (
-0.36,
+1.17), (
-0.22,
+0.91)
1192(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.30,
-0.87), (
+0.26,
+0.70), (
-0.14,
-0.98)
1193(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.22,
+0.66), (
+0.85,
+0.50)
1194(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.04,
-0.69)
1195(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1196mul
= matmul(mat, inv)
1198(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00)
1199(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1200(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
1201(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1202(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1208mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1209call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1216mul
= matmul(mat, inv)
1224mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1225call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1227(
+1.00,
+0.00), (
+0.08,
+0.57), (
+0.93,
-0.60), (
+0.26,
-0.06)
1228(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.38,
+0.30), (
-0.36,
-0.87)
1229(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.32,
+0.21)
1230(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1234(
+1.00,
+0.00), (
-0.08,
-0.57), (
-1.13,
+0.41), (
-0.06,
+0.15)
1235(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.38,
-0.30), (
+0.42,
+0.70)
1236(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.32,
-0.21)
1237(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1238mul
= matmul(mat, inv)
1240(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00)
1241(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1242(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1243(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1249mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1250call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1252(
+1.00,
+0.00), (
+0.72,
-0.80)
1253(
+0.00,
+0.00), (
+1.00,
+0.00)
1257(
+1.00,
+0.00), (
-0.72,
+0.80)
1258(
+0.00,
+0.00), (
+1.00,
+0.00)
1259mul
= matmul(mat, inv)
1261(
+1.00,
+0.00), (
+0.00,
+0.00)
1262(
+0.00,
+0.00), (
+1.00,
+0.00)
1268mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1269call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1271(
+1.00,
+0.00), (
+0.56,
+0.68), (
+0.30,
-0.20), (
-0.24,
+0.58), (
+0.08,
+0.24)
1272(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.47,
+0.53), (
-0.52,
+0.41), (
-0.19,
-0.81)
1273(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.17,
+0.05), (
+0.85,
+0.39)
1274(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.82,
+0.91)
1275(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1279(
+1.00,
+0.00), (
-0.56,
-0.68), (
-0.93,
+0.18), (
-0.48,
-0.63), (
+0.26,
-0.69)
1280(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.47,
-0.53), (
+0.57,
-0.52), (
-0.43,
+0.13)
1281(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.17,
-0.05), (
-0.76,
-0.58)
1282(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.82,
-0.91)
1283(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1284mul
= matmul(mat, inv)
1286(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
-0.00)
1287(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1288(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1289(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1290(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1296mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1297call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1299(
+1.00,
+0.00), (
+0.29,
-0.20), (
+0.78,
+0.90), (
+0.12,
-0.05)
1300(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.76,
+0.42), (
-0.80,
-0.38)
1301(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.81,
+0.15)
1302(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1306(
+1.00,
+0.00), (
-0.29,
+0.20), (
-0.48,
-0.93), (
-0.18,
+0.93)
1307(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.76,
-0.42), (
+1.35,
+0.83)
1308(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.81,
-0.15)
1309(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1310mul
= matmul(mat, inv)
1312(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1313(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1314(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1315(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1321mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1322call setMatInit(mat, lowDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1324(
+1.00,
+0.00), (
-0.43,
+0.31), (
+0.78,
+0.34)
1325(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.07,
+0.21)
1326(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1330(
+1.00,
+0.00), (
+0.43,
-0.31), (
-0.88,
-0.41)
1331(
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.07,
-0.21)
1332(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1333mul
= matmul(mat, inv)
1335(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1336(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1337(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1360mul
= matmul(mat, inv)
1373 +1.00,
+0.00,
+0.00,
+0.00
1374 +9.00,
+1.00,
+0.00,
+0.00
1375 +9.00,
+6.00,
+1.00,
+0.00
1376 +6.00,
+4.00,
+4.00,
+1.00
1380 +1.00,
+0.00,
+0.00,
+0.00
1381 -9.00,
+1.00,
+0.00,
+0.00
1382 +45.00,
-6.00,
+1.00,
+0.00
1383 -150.00,
+20.00,
-4.00,
+1.00
1384mul
= matmul(mat, inv)
1386 +1.00,
+0.00,
+0.00,
+0.00
1387 +0.00,
+1.00,
+0.00,
+0.00
1388 +0.00,
+0.00,
+1.00,
+0.00
1389 +0.00,
+0.00,
+0.00,
+1.00
1406 +22.00,
-4.00,
+1.00
1407mul
= matmul(mat, inv)
1427mul
= matmul(mat, inv)
1444mul
= matmul(mat, inv)
1455 +1.00,
+0.00,
+0.00,
+0.00
1456 +6.00,
+1.00,
+0.00,
+0.00
1457 +9.00,
+8.00,
+1.00,
+0.00
1458 +6.00,
+9.00,
+4.00,
+1.00
1462 +1.00,
+0.00,
+0.00,
+0.00
1463 -6.00,
+1.00,
+0.00,
+0.00
1464 +39.00,
-8.00,
+1.00,
+0.00
1465 -108.00,
+23.00,
-4.00,
+1.00
1466mul
= matmul(mat, inv)
1468 +1.00,
+0.00,
+0.00,
+0.00
1469 +0.00,
+1.00,
+0.00,
+0.00
1470 +0.00,
+0.00,
+1.00,
+0.00
1471 +0.00,
+0.00,
+0.00,
+1.00
1487mul
= matmul(mat, inv)
1499 +1.00,
+0.00,
+0.00,
+0.00
1500 +4.00,
+1.00,
+0.00,
+0.00
1501 +3.00,
+2.00,
+1.00,
+0.00
1502 +2.00,
+8.00,
+1.00,
+1.00
1506 +1.00,
+0.00,
+0.00,
+0.00
1507 -4.00,
+1.00,
+0.00,
+0.00
1508 +5.00,
-2.00,
+1.00,
+0.00
1509 +25.00,
-6.00,
-1.00,
+1.00
1510mul
= matmul(mat, inv)
1512 +1.00,
+0.00,
+0.00,
+0.00
1513 +0.00,
+1.00,
+0.00,
+0.00
1514 +0.00,
+0.00,
+1.00,
+0.00
1515 +0.00,
+0.00,
+0.00,
+1.00
1531mul
= matmul(mat, inv)
1550mul
= matmul(mat, inv)
1564mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1565call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1567(
+1.00,
+0.00), (
+0.00,
+0.00)
1568(
-0.25,
+0.08), (
+1.00,
+0.00)
1572(
+1.00,
+0.00), (
+0.00,
+0.00)
1573(
+0.25,
-0.08), (
+1.00,
+0.00)
1574mul
= matmul(mat, inv)
1576(
+1.00,
+0.00), (
+0.00,
+0.00)
1577(
+0.00,
+0.00), (
+1.00,
+0.00)
1583mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1584call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1586(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1587(
+0.89,
-0.85), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1588(
-0.71,
-0.07), (
+0.86,
+0.67), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1589(
+0.91,
-0.44), (
+0.81,
-0.25), (
-0.17,
-0.61), (
+1.00,
+0.00), (
+0.00,
+0.00)
1590(
-0.93,
+0.52), (
-0.37,
+0.70), (
-0.33,
+0.09), (
+0.75,
+0.93), (
+1.00,
+0.00)
1594(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1595(
-0.89,
+0.85), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1596(
+2.05,
-0.07), (
-0.86,
-0.67), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1597(
-0.00,
+0.77), (
-0.55,
-0.39), (
+0.17,
+0.61), (
+1.00,
+0.00), (
+0.00,
+0.00)
1598(
+2.60,
-0.37), (
+0.07,
-0.04), (
+0.77,
-0.71), (
-0.75,
-0.93), (
+1.00,
+0.00)
1599mul
= matmul(mat, inv)
1601(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1602(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1603(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1604(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1605(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1611mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1612call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1614(
+1.00,
+0.00), (
+0.00,
+0.00)
1615(
-0.73,
-0.67), (
+1.00,
+0.00)
1619(
+1.00,
+0.00), (
+0.00,
+0.00)
1620(
+0.73,
+0.67), (
+1.00,
+0.00)
1621mul
= matmul(mat, inv)
1623(
+1.00,
+0.00), (
+0.00,
+0.00)
1624(
+0.00,
+0.00), (
+1.00,
+0.00)
1630mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1631call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1633(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1634(
+0.83,
-0.65), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1635(
+0.68,
+0.59), (
+0.26,
+0.03), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1636(
-0.23,
+0.46), (
-0.73,
-0.67), (
-0.12,
+0.06), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1637(
-0.52,
+0.86), (
-0.28,
+0.53), (
-0.97,
+0.25), (
-0.69,
-0.96), (
+1.00,
+0.00), (
+0.00,
+0.00)
1638(
+0.79,
-0.97), (
-0.57,
+0.03), (
+0.08,
-0.50), (
+0.51,
-0.20), (
+0.42,
-0.80), (
+1.00,
+0.00)
1642(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1643(
-0.83,
+0.65), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1644(
-0.45,
-0.74), (
-0.26,
-0.03), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1645(
-0.91,
-0.59), (
+0.70,
+0.68), (
+0.12,
-0.06), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1646(
-0.05,
-2.12), (
-0.15,
+0.64), (
+1.11,
-0.18), (
+0.69,
+0.96), (
+1.00,
+0.00), (
+0.00,
+0.00)
1647(
+1.46,
+2.16), (
-0.34,
-0.75), (
-0.44,
+1.52), (
-1.56,
+0.35), (
-0.42,
+0.80), (
+1.00,
+0.00)
1648mul
= matmul(mat, inv)
1650(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1651(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1652(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1653(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1654(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1655(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1661mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1662call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1664(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1665(
-0.69,
-0.99), (
+1.00,
+0.00), (
+0.00,
+0.00)
1666(
-0.70,
-0.75), (
-0.61,
-0.42), (
+1.00,
+0.00)
1670(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1671(
+0.69,
+0.99), (
+1.00,
+0.00), (
+0.00,
+0.00)
1672(
+0.71,
+1.64), (
+0.61,
+0.42), (
+1.00,
+0.00)
1673mul
= matmul(mat, inv)
1675(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1676(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1677(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1683mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1684call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1686(
+1.00,
+0.00), (
+0.00,
+0.00)
1687(
-0.15,
+0.92), (
+1.00,
+0.00)
1691(
+1.00,
+0.00), (
+0.00,
+0.00)
1692(
+0.15,
-0.92), (
+1.00,
+0.00)
1693mul
= matmul(mat, inv)
1695(
+1.00,
+0.00), (
+0.00,
+0.00)
1696(
+0.00,
+0.00), (
+1.00,
+0.00)
1702mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1703call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1705(
+1.00,
+0.00), (
+0.00,
+0.00)
1706(
+0.72,
+0.26), (
+1.00,
+0.00)
1710(
+1.00,
+0.00), (
+0.00,
+0.00)
1711(
-0.72,
-0.26), (
+1.00,
+0.00)
1712mul
= matmul(mat, inv)
1714(
+1.00,
+0.00), (
+0.00,
+0.00)
1715(
+0.00,
+0.00), (
+1.00,
+0.00)
1721mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1722call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1729mul
= matmul(mat, inv)
1737mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1738call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1740(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1741(
+0.86,
+0.50), (
+1.00,
+0.00), (
+0.00,
+0.00)
1742(
+0.56,
+0.92), (
+0.16,
-0.14), (
+1.00,
+0.00)
1746(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1747(
-0.86,
-0.50), (
+1.00,
+0.00), (
+0.00,
+0.00)
1748(
-0.35,
-0.96), (
-0.16,
+0.14), (
+1.00,
+0.00)
1749mul
= matmul(mat, inv)
1751(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1752(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1753(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1759mat
= getUnifRand((
-1.,
-1.), (
+1.,
+1.), ndim, ndim)
1760call setMatInit(mat, uppDia, (
0._TKG,
0._TKG), (
1._TKG,
0._TKG))
1762(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1763(
-0.12,
-0.09), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1764(
-0.82,
+0.26), (
+0.41,
-0.99), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1765(
-0.57,
+0.69), (
+0.77,
+0.43), (
-0.45,
+0.44), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1766(
+0.77,
-0.42), (
+0.37,
-0.81), (
+0.33,
+0.69), (
-0.38,
-0.34), (
+1.00,
+0.00), (
+0.00,
+0.00)
1767(
+0.09,
+0.68), (
+0.83,
+0.06), (
+0.65,
+0.12), (
+0.08,
-0.06), (
+0.94,
-0.59), (
+1.00,
+0.00)
1771(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1772(
+0.12,
+0.09), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1773(
+0.68,
-0.17), (
-0.41,
+0.99), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1774(
+0.74,
-1.19), (
-0.52,
+0.19), (
+0.45,
-0.44), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1775(
-0.54,
-0.13), (
+0.19,
+0.66), (
-0.01,
-0.71), (
+0.38,
+0.34), (
+1.00,
+0.00), (
+0.00,
+0.00)
1776(
-0.04,
-0.79), (
-0.99,
-1.21), (
-0.23,
+0.61), (
-0.64,
-0.03), (
-0.94,
+0.59), (
+1.00,
+0.00)
1777mul
= matmul(mat, inv)
1779(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1780(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1781(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1782(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
1783(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
1784(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
1792mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1804if (info
/= 0)
error stop 'LUP factorization failed.'
1808 +0.00,
+0.60,
-10.20
1815mul
= matmul(mat, inv)
1822mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1834if (info
/= 0)
error stop 'LUP factorization failed.'
1838 +0.00,
+0.60,
-10.20
1845mul
= matmul(mat, inv)
1852mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1864if (info
/= 0)
error stop 'LUP factorization failed.'
1868 +0.00,
+0.60,
-10.20
1875mul
= matmul(mat, inv)
1882mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1894if (info
/= 0)
error stop 'LUP factorization failed.'
1898 +0.00,
+0.60,
-10.20
1905mul
= matmul(mat, inv)
1912mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1924if (info
/= 0)
error stop 'LUP factorization failed.'
1928 +0.00,
+0.60,
-10.20
1935mul
= matmul(mat, inv)
1942mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1954if (info
/= 0)
error stop 'LUP factorization failed.'
1958 +0.00,
+0.60,
-10.20
1965mul
= matmul(mat, inv)
1972mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
1984if (info
/= 0)
error stop 'LUP factorization failed.'
1988 +0.00,
+0.60,
-10.20
1995mul
= matmul(mat, inv)
2002mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
2014if (info
/= 0)
error stop 'LUP factorization failed.'
2018 +0.00,
+0.60,
-10.20
2025mul
= matmul(mat, inv)
2032mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
2044if (info
/= 0)
error stop 'LUP factorization failed.'
2048 +0.00,
+0.60,
-10.20
2055mul
= matmul(mat, inv)
2062mat
= reshape([
1,
0,
2,
-1,
5,
0,
0,
3,
-9], shape
= [
3,
3], order
= [
2,
1])
2074if (info
/= 0)
error stop 'LUP factorization failed.'
2078 +0.00,
+0.60,
-10.20
2085mul
= matmul(mat, inv)
2100(
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00), (
+3.20,
-1.00), (
+3.60,
-1.00), (
+4.00,
-1.00), (
+4.40,
-1.00), (
+4.80,
-1.00), (
+5.20,
-1.00)
2101(
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00), (
+3.20,
-1.00), (
+3.60,
-1.00), (
+4.00,
-1.00), (
+4.40,
-1.00), (
+4.80,
-1.00)
2102(
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00), (
+3.20,
-1.00), (
+3.60,
-1.00), (
+4.00,
-1.00), (
+4.40,
-1.00)
2103(
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00), (
+3.20,
-1.00), (
+3.60,
-1.00), (
+4.00,
-1.00)
2104(
+3.60,
+1.00), (
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00), (
+3.20,
-1.00), (
+3.60,
-1.00)
2105(
+4.00,
+1.00), (
+3.60,
+1.00), (
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00), (
+3.20,
-1.00)
2106(
+4.40,
+1.00), (
+4.00,
+1.00), (
+3.60,
+1.00), (
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00), (
+2.80,
-1.00)
2107(
+4.80,
+1.00), (
+4.40,
+1.00), (
+4.00,
+1.00), (
+3.60,
+1.00), (
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00), (
+2.40,
-1.00)
2108(
+5.20,
+1.00), (
+4.80,
+1.00), (
+4.40,
+1.00), (
+4.00,
+1.00), (
+3.60,
+1.00), (
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00)
2110if (info
/= 0)
error stop 'LUP factorization failed.'
2112(
+5.20,
+1.00), (
+4.80,
+1.00), (
+4.40,
+1.00), (
+4.00,
+1.00), (
+3.60,
+1.00), (
+3.20,
+1.00), (
+2.80,
+1.00), (
+2.40,
+1.00), (
+2.00,
+1.00)
2113(
+0.41,
+0.11), (
+0.56,
-1.95), (
+1.13,
-1.91), (
+1.69,
-1.86), (
+2.25,
-1.82), (
+2.81,
-1.77), (
+3.38,
-1.73), (
+3.94,
-1.68), (
+4.50,
-1.63)
2114(
+0.48,
+0.10), (
-0.05,
-0.09), (
+0.61,
-1.91), (
+1.23,
-1.81), (
+1.84,
-1.72), (
+2.46,
-1.62), (
+3.07,
-1.53), (
+3.69,
-1.43), (
+4.30,
-1.34)
2115(
+0.55,
+0.09), (
-0.04,
-0.08), (
-0.06,
-0.06), (
+0.67,
-1.87), (
+1.35,
-1.75), (
+2.02,
-1.62), (
+2.69,
-1.50), (
+3.37,
-1.37), (
+4.04,
-1.25)
2116(
+0.63,
+0.07), (
-0.03,
-0.07), (
-0.05,
-0.05), (
-0.07,
-0.03), (
+0.74,
-1.87), (
+1.47,
-1.73), (
+2.21,
-1.60), (
+2.94,
-1.46), (
+3.68,
-1.33)
2117(
+0.70,
+0.06), (
-0.03,
-0.05), (
-0.04,
-0.04), (
-0.05,
-0.03), (
-0.06,
-0.00), (
+0.79,
-1.88), (
+1.58,
-1.76), (
+2.38,
-1.64), (
+3.17,
-1.53)
2118(
+0.78,
+0.04), (
-0.02,
-0.04), (
-0.03,
-0.03), (
-0.04,
-0.02), (
-0.04,
-0.00), (
-0.04,
+0.01), (
+0.83,
-1.92), (
+1.66,
-1.83), (
+2.49,
-1.75)
2119(
+0.85,
+0.03), (
-0.01,
-0.03), (
-0.02,
-0.02), (
-0.03,
-0.01), (
-0.03,
-0.00), (
-0.03,
+0.01), (
-0.02,
+0.02), (
+0.84,
-1.96), (
+1.68,
-1.92)
2120(
+0.93,
+0.01), (
-0.01,
-0.01), (
-0.01,
-0.01), (
-0.01,
-0.01), (
-0.01,
-0.00), (
-0.01,
+0.00), (
-0.01,
+0.01), (
-0.01,
+0.01), (
+0.82,
-1.99)
2123(
-0.17,
-0.42), (
-0.12,
+0.13), (
-0.06,
+0.15), (
+0.00,
+0.15), (
+0.05,
+0.13), (
+0.09,
+0.09), (
+0.11,
+0.05), (
+0.11,
+0.00), (
+0.04,
-0.28)
2124(
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.06,
-0.03), (
-0.06,
-0.01), (
-0.05,
+0.01), (
-0.04,
+0.03), (
-0.03,
+0.04), (
-0.01,
+0.04), (
+0.11,
+0.00)
2125(
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.06,
-0.03), (
-0.06,
-0.01), (
-0.05,
+0.01), (
-0.04,
+0.03), (
-0.03,
+0.04), (
+0.11,
+0.05)
2126(
+0.01,
+0.00), (
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.06,
-0.03), (
-0.06,
-0.01), (
-0.05,
+0.01), (
-0.04,
+0.03), (
+0.09,
+0.09)
2127(
+0.00,
+0.01), (
+0.01,
+0.00), (
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.06,
-0.03), (
-0.06,
-0.01), (
-0.05,
+0.01), (
+0.05,
+0.13)
2128(
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.01,
+0.00), (
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.06,
-0.03), (
-0.06,
-0.01), (
+0.00,
+0.15)
2129(
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.01,
+0.00), (
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.06,
-0.03), (
-0.06,
+0.15)
2130(
-0.00,
+0.01), (
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.01,
+0.00), (
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.04,
-0.55), (
-0.12,
+0.13)
2131(
-0.01,
+0.01), (
-0.00,
+0.01), (
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.00,
+0.01), (
+0.01,
+0.00), (
+0.01,
+0.00), (
+0.18,
+0.43), (
-0.17,
-0.42)
2132mul
= matmul(mat, inv)
2134(
+1.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00)
2135(
-0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
+0.00)
2136(
-0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
2137(
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
-0.00), (
+1.00,
-0.00), (
-0.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
2138(
-0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
-0.00), (
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
2139(
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
2140(
+0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
+0.00)
2141(
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
2142(
-0.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
2155 +1.00,
-0.91,
+0.82,
+0.33
2156 -0.91,
+1.00,
-0.93,
-0.46
2157 +0.82,
-0.93,
+1.00,
+0.14
2158 +0.33,
-0.46,
+0.14,
+1.00
2161 +1.00,
+0.00,
+0.00,
+0.00
2162 -0.91,
+0.41,
+0.00,
+0.00
2163 +0.82,
-0.45,
+0.37,
+0.00
2164 +0.33,
-0.39,
-0.83,
+0.21
2166call setMatInv(inv, chol, auxil
= choLow)
2168 +25.49,
+79.28,
+49.81,
+21.28
2169 +79.28,
+285.91,
+189.55,
+79.55
2170 +49.81,
+189.55,
+128.85,
+53.21
2171 +21.28,
+79.55,
+53.21,
+23.32
2172mul
= matmul(mat, inv)
2174 +1.00,
+0.00,
+0.00,
-0.00
2175 +0.00,
+1.00,
-0.00,
+0.00
2176 -0.00,
-0.00,
+1.00,
-0.00
2177 -0.00,
+0.00,
+0.00,
+1.00
2185 +1.00,
+0.08,
+0.25,
-0.49
2186 +0.08,
+1.00,
-0.88,
-0.76
2187 +0.25,
-0.88,
+1.00,
+0.57
2188 -0.49,
-0.76,
+0.57,
+1.00
2191 +1.00,
+0.00,
+0.00,
+0.00
2192 +0.08,
+1.00,
+0.00,
+0.00
2193 +0.25,
-0.90,
+0.36,
+0.00
2194 -0.49,
-0.72,
+0.13,
+0.48
2196call setMatInv(inv, chol, auxil
= choLow)
2198 +3.11,
-1.37,
-3.36,
+2.39
2199 -1.37,
+7.94,
+6.32,
+1.72
2200 -3.36,
+6.32,
+8.30,
-1.60
2201 +2.39,
+1.72,
-1.60,
+4.37
2202mul
= matmul(mat, inv)
2204 +1.00,
+0.00,
+0.00,
+0.00
2205 -0.00,
+1.00,
+0.00,
+0.00
2206 +0.00,
+0.00,
+1.00,
+0.00
2207 +0.00,
+0.00,
-0.00,
+1.00
2220call setMatInv(inv, chol, auxil
= choLow)
2223mul
= matmul(mat, inv)
2233 +1.00,
+0.24,
+0.60,
+0.39
2234 +0.24,
+1.00,
+0.83,
+0.64
2235 +0.60,
+0.83,
+1.00,
+0.79
2236 +0.39,
+0.64,
+0.79,
+1.00
2239 +1.00,
+0.00,
+0.00,
+0.00
2240 +0.24,
+0.97,
+0.00,
+0.00
2241 +0.60,
+0.71,
+0.37,
+0.00
2242 +0.39,
+0.56,
+0.42,
+0.60
2244call setMatInv(inv, chol, auxil
= choLow)
2246 +2.50,
+2.16,
-3.79,
+0.64
2247 +2.16,
+5.14,
-6.13,
+0.72
2248 -3.79,
-6.13,
+10.85,
-3.16
2249 +0.64,
+0.72,
-3.16,
+2.78
2250mul
= matmul(mat, inv)
2252 +1.00,
+0.00,
-0.00,
+0.00
2253 +0.00,
+1.00,
-0.00,
+0.00
2254 +0.00,
+0.00,
+1.00,
+0.00
2255 +0.00,
+0.00,
-0.00,
+1.00
2272call setMatInv(inv, chol, auxil
= choLow)
2275 -2.88,
+11.63,
-8.77
2277mul
= matmul(mat, inv)
2289 +1.00,
-0.53,
-0.80,
-0.74
2290 -0.53,
+1.00,
+0.07,
+0.87
2291 -0.80,
+0.07,
+1.00,
+0.21
2292 -0.74,
+0.87,
+0.21,
+1.00
2295 +1.00,
+0.00,
+0.00,
+0.00
2296 -0.53,
+0.85,
+0.00,
+0.00
2297 -0.80,
-0.42,
+0.43,
+0.00
2298 -0.74,
+0.56,
-0.37,
+0.00
2300call setMatInv(inv, chol, auxil
= choLow)
2302+1198849.75,
-216577.14,
+786819.25,
+915451.44
2303 -216577.14,
+39129.80,
-142141.64,
-165383.86
2304 +786819.25,
-142141.64,
+516399.94,
+600821.00
2305 +915451.44,
-165383.86,
+600821.00,
+699050.62
2306mul
= matmul(mat, inv)
2308 +0.94,
-0.01,
-0.12,
-0.03
2309 +0.06,
+0.98,
+0.09,
+0.06
2310 -0.03,
+0.01,
+1.02,
+0.00
2311 +0.00,
+0.02,
+0.06,
+1.00
2328call setMatInv(inv, chol, auxil
= choLow)
2333mul
= matmul(mat, inv)
2350call setMatInv(inv, chol, auxil
= choLow)
2353mul
= matmul(mat, inv)
2370call setMatInv(inv, chol, auxil
= choLow)
2374mul
= matmul(mat, inv)
2394call setMatInv(inv, chol, auxil
= choLow)
2399mul
= matmul(mat, inv)
2412(
+9.00,
+0.00), (
+3.00,
+3.00), (
+3.00,
-3.00)
2413(
+3.00,
-3.00), (
+18.00,
+0.00), (
+8.00,
-6.00)
2414(
+3.00,
+3.00), (
+8.00,
+6.00), (
+43.00,
+0.00)
2417(
+3.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
2418(
+1.00,
-1.00), (
+4.00,
+0.00), (
+0.00,
+0.00)
2419(
+1.00,
+1.00), (
+2.00,
+1.00), (
+6.00,
+0.00)
2421call setMatInv(inv, chol, auxil
= choLow)
2423(
+0.13,
-0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
2424(
-0.02,
+0.03), (
+0.07,
-0.00), (
-0.01,
+0.01)
2425(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
-0.00)
2426mul
= matmul(mat, inv)
2428(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
2429(
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
2430(
+0.00,
+0.00), (
-0.00,
-0.00), (
+1.00,
-0.00)
2434(
+9.00,
+0.00), (
+3.00,
+3.00), (
+3.00,
-3.00)
2435(
+3.00,
-3.00), (
+18.00,
+0.00), (
+8.00,
-6.00)
2436(
+3.00,
+3.00), (
+8.00,
+6.00), (
+43.00,
+0.00)
2439(
+3.00,
+0.00), (
+1.00,
+1.00), (
+1.00,
-1.00)
2440(
+0.00,
+0.00), (
+4.00,
+0.00), (
+2.00,
-1.00)
2441(
+0.00,
+0.00), (
+0.00,
+0.00), (
+6.00,
+0.00)
2443call setMatInv(inv, chol, auxil
= choUpp)
2445(
+0.13,
-0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
2446(
-0.02,
+0.03), (
+0.07,
-0.00), (
-0.01,
+0.01)
2447(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
-0.00)
2448mul
= matmul(mat, inv)
2450(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
2451(
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
2452(
+0.00,
+0.00), (
-0.00,
-0.00), (
+1.00,
-0.00)
2472call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2475call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2478mul
= matmul(mat, inv)
2483call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2486call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2489mul
= matmul(mat, inv)
2499 +1.00,
-0.70,
+0.82,
+0.01,
+0.26
2500 -0.70,
+1.00,
-0.24,
-0.58,
-0.30
2501 +0.82,
-0.24,
+1.00,
-0.34,
+0.38
2502 +0.01,
-0.58,
-0.34,
+1.00,
+0.33
2503 +0.26,
-0.30,
+0.38,
+0.33,
+1.00
2506 +1.00,
+0.00,
+0.00,
+0.00,
+0.00
2507 -0.70,
+0.72,
+0.00,
+0.00,
+0.00
2508 +0.82,
+0.47,
+0.34,
+0.00,
+0.00
2509 +0.01,
-0.81,
+0.08,
+0.58,
+0.00
2510 +0.26,
-0.17,
+0.74,
+0.23,
+0.55
2514call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2516 +35.64,
+20.92,
-26.81,
+0.33,
+7.25
2517 +0.00,
+14.73,
-14.20,
+2.35,
+3.68
2518 +0.00,
+0.00,
+22.79,
+1.98,
-6.74
2519 +0.00,
+0.00,
+0.00,
+3.48,
-1.30
2520 +0.00,
+0.00,
+0.00,
+0.00,
+3.25
2521call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2523 +35.64,
+20.92,
-26.81,
+0.33,
+7.25
2524 +20.92,
+14.73,
-14.20,
+2.35,
+3.68
2525 -26.81,
-14.20,
+22.79,
+1.98,
-6.74
2526 +0.33,
+2.35,
+1.98,
+3.48,
-1.30
2527 +7.25,
+3.68,
-6.74,
-1.30,
+3.25
2528mul
= matmul(mat, inv)
2530 +1.00,
+0.00,
-0.00,
-0.00,
+0.00
2531 -0.00,
+1.00,
+0.00,
-0.00,
-0.00
2532 +0.00,
+0.00,
+1.00,
-0.00,
+0.00
2533 -0.00,
-0.00,
+0.00,
+1.00,
-0.00
2534 +0.00,
-0.00,
-0.00,
-0.00,
+1.00
2537call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2539 +35.64,
+0.00,
+0.00,
+0.00,
+0.00
2540 +20.92,
+14.73,
+0.00,
+0.00,
+0.00
2541 -26.81,
-14.20,
+22.79,
+0.00,
+0.00
2542 +0.33,
+2.35,
+1.98,
+3.48,
+0.00
2543 +7.25,
+3.68,
-6.74,
-1.30,
+3.25
2544call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2546 +35.64,
+20.92,
-26.81,
+0.33,
+7.25
2547 +20.92,
+14.73,
-14.20,
+2.35,
+3.68
2548 -26.81,
-14.20,
+22.79,
+1.98,
-6.74
2549 +0.33,
+2.35,
+1.98,
+3.48,
-1.30
2550 +7.25,
+3.68,
-6.74,
-1.30,
+3.25
2551mul
= matmul(mat, inv)
2553 +1.00,
+0.00,
-0.00,
-0.00,
+0.00
2554 -0.00,
+1.00,
+0.00,
-0.00,
-0.00
2555 +0.00,
+0.00,
+1.00,
-0.00,
+0.00
2556 -0.00,
-0.00,
+0.00,
+1.00,
-0.00
2557 +0.00,
-0.00,
-0.00,
-0.00,
+1.00
2576call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2578 +153.95,
-105.27,
-54.66
2579 +0.00,
+75.09,
+34.81
2580 +0.00,
+0.00,
+22.51
2581call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2583 +153.95,
-105.27,
-54.66
2584 -105.27,
+75.09,
+34.81
2585 -54.66,
+34.81,
+22.51
2586mul
= matmul(mat, inv)
2593call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2595 +153.95,
+0.00,
+0.00
2596 -105.27,
+75.09,
+0.00
2597 -54.66,
+34.81,
+22.51
2598call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2600 +153.95,
-105.27,
-54.66
2601 -105.27,
+75.09,
+34.81
2602 -54.66,
+34.81,
+22.51
2603mul
= matmul(mat, inv)
2615 +1.00,
-0.49,
+0.23,
-0.57,
+0.38
2616 -0.49,
+1.00,
-0.23,
+0.75,
-0.41
2617 +0.23,
-0.23,
+1.00,
+0.17,
+0.87
2618 -0.57,
+0.75,
+0.17,
+1.00,
-0.03
2619 +0.38,
-0.41,
+0.87,
-0.03,
+1.00
2622 +1.00,
+0.00,
+0.00,
+0.00,
+0.00
2623 -0.49,
+0.87,
+0.00,
+0.00,
+0.00
2624 +0.23,
-0.13,
+0.96,
+0.00,
+0.00
2625 -0.57,
+0.54,
+0.38,
+0.49,
+0.00
2626 +0.38,
-0.26,
+0.77,
+0.06,
+0.43
2630call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2632 +1.90,
-0.51,
-0.09,
+1.45,
-0.81
2633 +0.00,
+3.60,
+0.28,
-2.99,
+1.33
2634 +0.00,
+0.00,
+4.73,
-1.18,
-3.99
2635 +0.00,
+0.00,
+0.00,
+4.24,
-0.62
2636 +0.00,
+0.00,
+0.00,
+0.00,
+5.29
2637call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2639 +1.90,
-0.51,
-0.09,
+1.45,
-0.81
2640 -0.51,
+3.60,
+0.28,
-2.99,
+1.33
2641 -0.09,
+0.28,
+4.73,
-1.18,
-3.99
2642 +1.45,
-2.99,
-1.18,
+4.24,
-0.62
2643 -0.81,
+1.33,
-3.99,
-0.62,
+5.29
2644mul
= matmul(mat, inv)
2646 +1.00,
+0.00,
+0.00,
+0.00,
-0.00
2647 -0.00,
+1.00,
-0.00,
+0.00,
+0.00
2648 -0.00,
+0.00,
+1.00,
-0.00,
-0.00
2649 +0.00,
+0.00,
+0.00,
+1.00,
+0.00
2650 -0.00,
+0.00,
+0.00,
-0.00,
+1.00
2653call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2655 +1.90,
+0.00,
+0.00,
+0.00,
+0.00
2656 -0.51,
+3.60,
+0.00,
+0.00,
+0.00
2657 -0.09,
+0.28,
+4.73,
+0.00,
+0.00
2658 +1.45,
-2.99,
-1.18,
+4.24,
+0.00
2659 -0.81,
+1.33,
-3.99,
-0.62,
+5.29
2660call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2662 +1.90,
-0.51,
-0.09,
+1.45,
-0.81
2663 -0.51,
+3.60,
+0.28,
-2.99,
+1.33
2664 -0.09,
+0.28,
+4.73,
-1.18,
-3.99
2665 +1.45,
-2.99,
-1.18,
+4.24,
-0.62
2666 -0.81,
+1.33,
-3.99,
-0.62,
+5.29
2667mul
= matmul(mat, inv)
2669 +1.00,
+0.00,
+0.00,
+0.00,
-0.00
2670 -0.00,
+1.00,
-0.00,
+0.00,
+0.00
2671 -0.00,
+0.00,
+1.00,
-0.00,
-0.00
2672 +0.00,
+0.00,
+0.00,
+1.00,
+0.00
2673 -0.00,
+0.00,
+0.00,
-0.00,
+1.00
2681 +1.00,
-0.43,
-0.16,
+0.81,
-0.90
2682 -0.43,
+1.00,
-0.37,
-0.09,
+0.66
2683 -0.16,
-0.37,
+1.00,
-0.46,
+0.15
2684 +0.81,
-0.09,
-0.46,
+1.00,
-0.61
2685 -0.90,
+0.66,
+0.15,
-0.61,
+1.00
2688 +1.00,
+0.00,
+0.00,
+0.00,
+0.00
2689 -0.43,
+0.90,
+0.00,
+0.00,
+0.00
2690 -0.16,
-0.48,
+0.86,
+0.00,
+0.00
2691 +0.81,
+0.29,
-0.22,
+0.46,
+0.00
2692 -0.90,
+0.30,
+0.18,
+0.18,
+0.18
2696call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2698 +35.04,
-9.76,
-9.88,
-15.69,
+30.04
2699 +0.00,
+6.04,
+4.01,
+3.48,
-11.30
2700 +0.00,
+0.00,
+4.56,
+4.83,
-9.32
2701 +0.00,
+0.00,
+0.00,
+9.28,
-11.55
2702 +0.00,
+0.00,
+0.00,
+0.00,
+29.97
2703call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2705 +35.04,
-9.76,
-9.88,
-15.69,
+30.04
2706 -9.76,
+6.04,
+4.01,
+3.48,
-11.30
2707 -9.88,
+4.01,
+4.56,
+4.83,
-9.32
2708 -15.69,
+3.48,
+4.83,
+9.28,
-11.55
2709 +30.04,
-11.30,
-9.32,
-11.55,
+29.97
2710mul
= matmul(mat, inv)
2712 +1.00,
+0.00,
-0.00,
+0.00,
-0.00
2713 +0.00,
+1.00,
+0.00,
+0.00,
+0.00
2714 -0.00,
-0.00,
+1.00,
-0.00,
+0.00
2715 +0.00,
+0.00,
+0.00,
+1.00,
-0.00
2716 +0.00,
-0.00,
+0.00,
-0.00,
+1.00
2719call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2721 +35.04,
+0.00,
+0.00,
+0.00,
+0.00
2722 -9.76,
+6.04,
+0.00,
+0.00,
+0.00
2723 -9.88,
+4.01,
+4.56,
+0.00,
+0.00
2724 -15.69,
+3.48,
+4.83,
+9.28,
+0.00
2725 +30.04,
-11.30,
-9.32,
-11.55,
+29.97
2726call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2728 +35.04,
-9.76,
-9.88,
-15.69,
+30.04
2729 -9.76,
+6.04,
+4.01,
+3.48,
-11.30
2730 -9.88,
+4.01,
+4.56,
+4.83,
-9.32
2731 -15.69,
+3.48,
+4.83,
+9.28,
-11.55
2732 +30.04,
-11.30,
-9.32,
-11.55,
+29.97
2733mul
= matmul(mat, inv)
2735 +1.00,
+0.00,
-0.00,
+0.00,
-0.00
2736 +0.00,
+1.00,
+0.00,
+0.00,
+0.00
2737 -0.00,
-0.00,
+1.00,
-0.00,
+0.00
2738 +0.00,
+0.00,
+0.00,
+1.00,
-0.00
2739 +0.00,
-0.00,
+0.00,
-0.00,
+1.00
2747 +1.00,
-0.62,
-0.87,
-0.53,
-0.57,
-0.32
2748 -0.62,
+1.00,
+0.28,
+0.53,
-0.13,
+0.40
2749 -0.87,
+0.28,
+1.00,
+0.19,
+0.67,
+0.02
2750 -0.53,
+0.53,
+0.19,
+1.00,
-0.13,
+0.16
2751 -0.57,
-0.13,
+0.67,
-0.13,
+1.00,
+0.40
2752 -0.32,
+0.40,
+0.02,
+0.16,
+0.40,
+1.00
2755 +1.00,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
2756 -0.62,
+0.79,
+0.00,
+0.00,
+0.00,
+0.00
2757 -0.87,
-0.33,
+0.38,
+0.00,
+0.00,
+0.00
2758 -0.53,
+0.26,
-0.49,
+0.64,
+0.00,
+0.00
2759 -0.57,
-0.62,
-0.06,
-0.48,
+0.26,
+0.00
2760 -0.32,
+0.26,
-0.45,
-0.46,
+0.51,
+0.40
2764call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2766 +149.61,
+68.18,
+56.56,
+44.01,
+68.09,
-14.84
2767 +0.00,
+35.89,
+20.89,
+19.78,
+36.49,
-10.71
2768 +0.00,
+0.00,
+29.21,
+15.84,
+17.67,
-0.49
2769 +0.00,
+0.00,
+0.00,
+14.43,
+20.90,
-4.82
2770 +0.00,
+0.00,
+0.00,
+0.00,
+40.64,
-12.74
2771 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+6.40
2772call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2774 +149.61,
+68.18,
+56.56,
+44.01,
+68.09,
-14.84
2775 +68.18,
+35.89,
+20.89,
+19.78,
+36.49,
-10.71
2776 +56.56,
+20.89,
+29.21,
+15.84,
+17.67,
-0.49
2777 +44.01,
+19.78,
+15.84,
+14.43,
+20.90,
-4.82
2778 +68.09,
+36.49,
+17.67,
+20.90,
+40.64,
-12.74
2779 -14.84,
-10.71,
-0.49,
-4.82,
-12.74,
+6.40
2780mul
= matmul(mat, inv)
2782 +1.00,
+0.00,
+0.00,
+0.00,
+0.00,
-0.00
2783 -0.00,
+1.00,
-0.00,
-0.00,
-0.00,
+0.00
2784 -0.00,
-0.00,
+1.00,
-0.00,
-0.00,
-0.00
2785 +0.00,
-0.00,
-0.00,
+1.00,
-0.00,
+0.00
2786 -0.00,
-0.00,
-0.00,
-0.00,
+1.00,
-0.00
2787 +0.00,
-0.00,
-0.00,
-0.00,
-0.00,
+1.00
2790call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2792 +149.61,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
2793 +68.18,
+35.89,
+0.00,
+0.00,
+0.00,
+0.00
2794 +56.56,
+20.89,
+29.21,
+0.00,
+0.00,
+0.00
2795 +44.01,
+19.78,
+15.84,
+14.43,
+0.00,
+0.00
2796 +68.09,
+36.49,
+17.67,
+20.90,
+40.64,
+0.00
2797 -14.84,
-10.71,
-0.49,
-4.82,
-12.74,
+6.40
2798call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2800 +149.61,
+68.18,
+56.56,
+44.01,
+68.09,
-14.84
2801 +68.18,
+35.89,
+20.89,
+19.78,
+36.49,
-10.71
2802 +56.56,
+20.89,
+29.21,
+15.84,
+17.67,
-0.49
2803 +44.01,
+19.78,
+15.84,
+14.43,
+20.90,
-4.82
2804 +68.09,
+36.49,
+17.67,
+20.90,
+40.64,
-12.74
2805 -14.84,
-10.71,
-0.49,
-4.82,
-12.74,
+6.40
2806mul
= matmul(mat, inv)
2808 +1.00,
+0.00,
+0.00,
+0.00,
+0.00,
-0.00
2809 -0.00,
+1.00,
-0.00,
-0.00,
-0.00,
+0.00
2810 -0.00,
-0.00,
+1.00,
-0.00,
-0.00,
-0.00
2811 +0.00,
-0.00,
-0.00,
+1.00,
-0.00,
+0.00
2812 -0.00,
-0.00,
-0.00,
-0.00,
+1.00,
-0.00
2813 +0.00,
-0.00,
-0.00,
-0.00,
-0.00,
+1.00
2821 +1.00,
+0.94,
-0.38,
-0.35,
+0.12,
+0.29
2822 +0.94,
+1.00,
-0.16,
-0.53,
+0.05,
+0.38
2823 -0.38,
-0.16,
+1.00,
-0.53,
-0.63,
-0.10
2824 -0.35,
-0.53,
-0.53,
+1.00,
+0.41,
-0.46
2825 +0.12,
+0.05,
-0.63,
+0.41,
+1.00,
+0.51
2826 +0.29,
+0.38,
-0.10,
-0.46,
+0.51,
+1.00
2829 +1.00,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
2830 +0.94,
+0.33,
+0.00,
+0.00,
+0.00,
+0.00
2831 -0.38,
+0.62,
+0.69,
+0.00,
+0.00,
+0.00
2832 -0.35,
-0.61,
-0.42,
+0.58,
+0.00,
+0.00
2833 +0.12,
-0.20,
-0.67,
+0.08,
+0.70,
+0.00
2834 +0.29,
+0.33,
-0.29,
-0.47,
+0.55,
+0.44
2838call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2840 +24.17,
-21.82,
+7.49,
+1.23,
+1.35,
+2.02
2841 +0.00,
+22.09,
-5.35,
+1.68,
-2.06,
-0.89
2842 +0.00,
+0.00,
+5.09,
+2.61,
+0.94,
+1.12
2843 +0.00,
+0.00,
+0.00,
+7.66,
-4.22,
+4.93
2844 +0.00,
+0.00,
+0.00,
+0.00,
+5.39,
-4.19
2845 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+5.27
2846call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2848 +24.17,
-21.82,
+7.49,
+1.23,
+1.35,
+2.02
2849 -21.82,
+22.09,
-5.35,
+1.68,
-2.06,
-0.89
2850 +7.49,
-5.35,
+5.09,
+2.61,
+0.94,
+1.12
2851 +1.23,
+1.68,
+2.61,
+7.66,
-4.22,
+4.93
2852 +1.35,
-2.06,
+0.94,
-4.22,
+5.39,
-4.19
2853 +2.02,
-0.89,
+1.12,
+4.93,
-4.19,
+5.27
2854mul
= matmul(mat, inv)
2856 +1.00,
+0.00,
+0.00,
+0.00,
-0.00,
+0.00
2857 +0.00,
+1.00,
+0.00,
+0.00,
-0.00,
+0.00
2858 -0.00,
-0.00,
+1.00,
-0.00,
+0.00,
-0.00
2859 +0.00,
-0.00,
+0.00,
+1.00,
+0.00,
-0.00
2860 +0.00,
-0.00,
+0.00,
+0.00,
+1.00,
+0.00
2861 +0.00,
-0.00,
+0.00,
+0.00,
-0.00,
+1.00
2864call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2866 +24.17,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
2867 -21.82,
+22.09,
+0.00,
+0.00,
+0.00,
+0.00
2868 +7.49,
-5.35,
+5.09,
+0.00,
+0.00,
+0.00
2869 +1.23,
+1.68,
+2.61,
+7.66,
+0.00,
+0.00
2870 +1.35,
-2.06,
+0.94,
-4.22,
+5.39,
+0.00
2871 +2.02,
-0.89,
+1.12,
+4.93,
-4.19,
+5.27
2872call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2874 +24.17,
-21.82,
+7.49,
+1.23,
+1.35,
+2.02
2875 -21.82,
+22.09,
-5.35,
+1.68,
-2.06,
-0.89
2876 +7.49,
-5.35,
+5.09,
+2.61,
+0.94,
+1.12
2877 +1.23,
+1.68,
+2.61,
+7.66,
-4.22,
+4.93
2878 +1.35,
-2.06,
+0.94,
-4.22,
+5.39,
-4.19
2879 +2.02,
-0.89,
+1.12,
+4.93,
-4.19,
+5.27
2880mul
= matmul(mat, inv)
2882 +1.00,
+0.00,
+0.00,
+0.00,
-0.00,
+0.00
2883 +0.00,
+1.00,
+0.00,
+0.00,
-0.00,
+0.00
2884 -0.00,
-0.00,
+1.00,
-0.00,
+0.00,
-0.00
2885 +0.00,
-0.00,
+0.00,
+1.00,
+0.00,
-0.00
2886 +0.00,
-0.00,
+0.00,
+0.00,
+1.00,
+0.00
2887 +0.00,
-0.00,
+0.00,
+0.00,
-0.00,
+1.00
2904call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2908call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2912mul
= matmul(mat, inv)
2918call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2922call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2926mul
= matmul(mat, inv)
2944call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2947call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2950mul
= matmul(mat, inv)
2955call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
2958call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
2961mul
= matmul(mat, inv)
2971 +1.00,
-0.42,
+0.35,
-0.17,
+0.64
2972 -0.42,
+1.00,
-0.95,
+0.65,
-0.13
2973 +0.35,
-0.95,
+1.00,
-0.82,
-0.11
2974 -0.17,
+0.65,
-0.82,
+1.00,
+0.47
2975 +0.64,
-0.13,
-0.11,
+0.47,
+1.00
2978 +1.00,
+0.00,
+0.00,
+0.00,
+0.00
2979 -0.42,
+0.91,
+0.00,
+0.00,
+0.00
2980 +0.35,
-0.89,
+0.29,
+0.00,
+0.00
2981 -0.17,
+0.64,
-0.66,
+0.36,
+0.00
2982 +0.64,
+0.15,
-0.66,
+0.11,
+0.33
2986call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
2988 +4.22,
-6.49,
-6.68,
+1.85,
-5.17
2989 +0.00,
+52.16,
+59.34,
+7.07,
+14.26
2990 +0.00,
+0.00,
+72.06,
+12.96,
+13.97
2991 +0.00,
+0.00,
+0.00,
+8.71,
-2.92
2992 +0.00,
+0.00,
+0.00,
+0.00,
+9.11
2993call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
2995 +4.22,
-6.49,
-6.68,
+1.85,
-5.17
2996 -6.49,
+52.16,
+59.34,
+7.07,
+14.26
2997 -6.68,
+59.34,
+72.06,
+12.96,
+13.97
2998 +1.85,
+7.07,
+12.96,
+8.71,
-2.92
2999 -5.17,
+14.26,
+13.97,
-2.92,
+9.11
3000mul
= matmul(mat, inv)
3002 +1.00,
+0.00,
+0.00,
+0.00,
-0.00
3003 +0.00,
+1.00,
-0.00,
-0.00,
-0.00
3004 -0.00,
+0.00,
+1.00,
+0.00,
+0.00
3005 +0.00,
-0.00,
-0.00,
+1.00,
-0.00
3006 +0.00,
-0.00,
-0.00,
-0.00,
+1.00
3009call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
3011 +4.22,
+0.00,
+0.00,
+0.00,
+0.00
3012 -6.49,
+52.16,
+0.00,
+0.00,
+0.00
3013 -6.68,
+59.34,
+72.06,
+0.00,
+0.00
3014 +1.85,
+7.07,
+12.96,
+8.71,
+0.00
3015 -5.17,
+14.26,
+13.97,
-2.92,
+9.11
3016call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3018 +4.22,
-6.49,
-6.68,
+1.85,
-5.17
3019 -6.49,
+52.16,
+59.34,
+7.07,
+14.26
3020 -6.68,
+59.34,
+72.06,
+12.96,
+13.97
3021 +1.85,
+7.07,
+12.96,
+8.71,
-2.92
3022 -5.17,
+14.26,
+13.97,
-2.92,
+9.11
3023mul
= matmul(mat, inv)
3025 +1.00,
+0.00,
+0.00,
+0.00,
-0.00
3026 +0.00,
+1.00,
-0.00,
-0.00,
-0.00
3027 -0.00,
+0.00,
+1.00,
+0.00,
+0.00
3028 +0.00,
-0.00,
-0.00,
+1.00,
-0.00
3029 +0.00,
-0.00,
-0.00,
-0.00,
+1.00
3037 +1.00,
+0.93,
-0.85,
+0.83,
-0.56,
-0.41
3038 +0.93,
+1.00,
-0.73,
+0.82,
-0.52,
-0.26
3039 -0.85,
-0.73,
+1.00,
-0.87,
+0.31,
+0.28
3040 +0.83,
+0.82,
-0.87,
+1.00,
-0.60,
-0.42
3041 -0.56,
-0.52,
+0.31,
-0.60,
+1.00,
+0.81
3042 -0.41,
-0.26,
+0.28,
-0.42,
+0.81,
+1.00
3045 +1.00,
+0.93,
-0.85,
+0.83,
-0.56,
-0.41
3046 +0.00,
+0.36,
+0.18,
+0.12,
+0.02,
+0.37
3047 +0.00,
+0.00,
+0.50,
-0.37,
-0.35,
-0.27
3048 +0.00,
+0.00,
+0.00,
+0.40,
-0.68,
-0.55
3049 +0.00,
+0.00,
+0.00,
+0.00,
+0.31,
+0.31
3050 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+0.47
3054call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3056 +73.94,
-52.67,
+55.33,
+44.08,
+22.38,
+1.78
3057 +0.00,
+41.25,
-37.70,
-32.34,
-13.72,
-2.97
3058 +0.00,
+0.00,
+49.69,
+40.87,
+22.40,
-1.84
3059 +0.00,
+0.00,
+0.00,
+37.35,
+19.15,
-1.48
3060 +0.00,
+0.00,
+0.00,
+0.00,
+14.76,
-4.51
3061 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+4.52
3062call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3064 +73.94,
-52.67,
+55.33,
+44.08,
+22.38,
+1.78
3065 -52.67,
+41.25,
-37.70,
-32.34,
-13.72,
-2.97
3066 +55.33,
-37.70,
+49.69,
+40.87,
+22.40,
-1.84
3067 +44.08,
-32.34,
+40.87,
+37.35,
+19.15,
-1.48
3068 +22.38,
-13.72,
+22.40,
+19.15,
+14.76,
-4.51
3069 +1.78,
-2.97,
-1.84,
-1.48,
-4.51,
+4.52
3070mul
= matmul(mat, inv)
3072 +1.00,
+0.00,
-0.00,
+0.00,
-0.00,
-0.00
3073 -0.00,
+1.00,
-0.00,
+0.00,
-0.00,
-0.00
3074 +0.00,
-0.00,
+1.00,
-0.00,
+0.00,
+0.00
3075 -0.00,
-0.00,
+0.00,
+1.00,
-0.00,
-0.00
3076 +0.00,
-0.00,
+0.00,
-0.00,
+1.00,
+0.00
3077 +0.00,
-0.00,
+0.00,
-0.00,
+0.00,
+1.00
3080call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3082 +73.94,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
3083 -52.67,
+41.25,
+0.00,
+0.00,
+0.00,
+0.00
3084 +55.33,
-37.70,
+49.69,
+0.00,
+0.00,
+0.00
3085 +44.08,
-32.34,
+40.87,
+37.35,
+0.00,
+0.00
3086 +22.38,
-13.72,
+22.40,
+19.15,
+14.76,
+0.00
3087 +1.78,
-2.97,
-1.84,
-1.48,
-4.51,
+4.52
3088call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3090 +73.94,
-52.67,
+55.33,
+44.08,
+22.38,
+1.78
3091 -52.67,
+41.25,
-37.70,
-32.34,
-13.72,
-2.97
3092 +55.33,
-37.70,
+49.69,
+40.87,
+22.40,
-1.84
3093 +44.08,
-32.34,
+40.87,
+37.35,
+19.15,
-1.48
3094 +22.38,
-13.72,
+22.40,
+19.15,
+14.76,
-4.51
3095 +1.78,
-2.97,
-1.84,
-1.48,
-4.51,
+4.52
3096mul
= matmul(mat, inv)
3098 +1.00,
+0.00,
-0.00,
+0.00,
-0.00,
-0.00
3099 -0.00,
+1.00,
-0.00,
+0.00,
-0.00,
-0.00
3100 +0.00,
-0.00,
+1.00,
-0.00,
+0.00,
+0.00
3101 -0.00,
-0.00,
+0.00,
+1.00,
-0.00,
-0.00
3102 +0.00,
-0.00,
+0.00,
-0.00,
+1.00,
+0.00
3103 +0.00,
-0.00,
+0.00,
-0.00,
+0.00,
+1.00
3118call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3121call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3124mul
= matmul(mat, inv)
3129call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3132call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3135mul
= matmul(mat, inv)
3156call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3158 +63.31,
+61.69,
-1.21
3159 +0.00,
+66.98,
+5.17
3161call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3163 +63.31,
+61.69,
-1.21
3164 +61.69,
+66.98,
+5.17
3166mul
= matmul(mat, inv)
3173call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3175 +63.31,
+0.00,
+0.00
3176 +61.69,
+66.98,
+0.00
3178call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3180 +63.31,
+61.69,
-1.21
3181 +61.69,
+66.98,
+5.17
3183mul
= matmul(mat, inv)
3202call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3205call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3208mul
= matmul(mat, inv)
3213call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3216call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3219mul
= matmul(mat, inv)
3240call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3245call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3250mul
= matmul(mat, inv)
3257call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3262call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3267mul
= matmul(mat, inv)
3290call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3292 +15.01,
+16.90,
-12.42
3293 +0.00,
+20.46,
-14.76
3294 +0.00,
+0.00,
+11.69
3295call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3297 +15.01,
+16.90,
-12.42
3298 +16.90,
+20.46,
-14.76
3299 -12.42,
-14.76,
+11.69
3300mul
= matmul(mat, inv)
3307call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3309 +15.01,
+0.00,
+0.00
3310 +16.90,
+20.46,
+0.00
3311 -12.42,
-14.76,
+11.69
3312call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3314 +15.01,
+16.90,
-12.42
3315 +16.90,
+20.46,
-14.76
3316 -12.42,
-14.76,
+11.69
3317mul
= matmul(mat, inv)
3338call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3342call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3346mul
= matmul(mat, inv)
3352call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3356call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3360mul
= matmul(mat, inv)
3380call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3384call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3388mul
= matmul(mat, inv)
3394call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3398call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3402mul
= matmul(mat, inv)
3420call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3423call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3426mul
= matmul(mat, inv)
3431call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3434call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3437mul
= matmul(mat, inv)
3447 +1.00,
+0.63,
+0.89,
-0.38,
+0.09,
+0.34
3448 +0.63,
+1.00,
+0.24,
-0.67,
+0.37,
-0.20
3449 +0.89,
+0.24,
+1.00,
-0.01,
-0.06,
+0.49
3450 -0.38,
-0.67,
-0.01,
+1.00,
-0.04,
+0.10
3451 +0.09,
+0.37,
-0.06,
-0.04,
+1.00,
-0.71
3452 +0.34,
-0.20,
+0.49,
+0.10,
-0.71,
+1.00
3455 +1.00,
+0.63,
+0.89,
-0.38,
+0.09,
+0.34
3456 +0.00,
+0.78,
-0.41,
-0.55,
+0.41,
-0.54
3457 +0.00,
+0.00,
+0.18,
+0.61,
+0.16,
-0.18
3458 +0.00,
+0.00,
+0.00,
+0.43,
+0.28,
+0.09
3459 +0.00,
+0.00,
+0.00,
+0.00,
+0.85,
-0.60
3460 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+0.44
3464call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3466 +219.38,
-70.11,
-170.01,
+37.19,
-15.98,
-19.77
3467 +0.00,
+25.15,
+53.56,
-10.10,
+4.20,
+6.44
3468 +0.00,
+0.00,
+133.66,
-29.20,
+11.67,
+13.76
3469 +0.00,
+0.00,
+0.00,
+8.48,
-3.42,
-3.51
3470 +0.00,
+0.00,
+0.00,
+0.00,
+4.01,
+3.68
3471 +0.00,
+0.00,
+0.00,
+0.00,
+0.00,
+5.17
3472call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3474 +219.38,
-70.11,
-170.01,
+37.19,
-15.98,
-19.77
3475 -70.11,
+25.15,
+53.56,
-10.10,
+4.20,
+6.44
3476 -170.01,
+53.56,
+133.66,
-29.20,
+11.67,
+13.76
3477 +37.19,
-10.10,
-29.20,
+8.48,
-3.42,
-3.51
3478 -15.98,
+4.20,
+11.67,
-3.42,
+4.01,
+3.68
3479 -19.77,
+6.44,
+13.76,
-3.51,
+3.68,
+5.17
3480mul
= matmul(mat, inv)
3482 +1.00,
+0.00,
+0.00,
-0.00,
-0.00,
-0.00
3483 +0.00,
+1.00,
-0.00,
+0.00,
-0.00,
-0.00
3484 +0.00,
-0.00,
+1.00,
-0.00,
-0.00,
-0.00
3485 -0.00,
-0.00,
-0.00,
+1.00,
+0.00,
+0.00
3486 +0.00,
+0.00,
+0.00,
-0.00,
+1.00,
+0.00
3487 +0.00,
-0.00,
+0.00,
-0.00,
+0.00,
+1.00
3490call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3492 +219.38,
+0.00,
+0.00,
+0.00,
+0.00,
+0.00
3493 -70.11,
+25.15,
+0.00,
+0.00,
+0.00,
+0.00
3494 -170.01,
+53.56,
+133.66,
+0.00,
+0.00,
+0.00
3495 +37.19,
-10.10,
-29.20,
+8.48,
+0.00,
+0.00
3496 -15.98,
+4.20,
+11.67,
-3.42,
+4.01,
+0.00
3497 -19.77,
+6.44,
+13.76,
-3.51,
+3.68,
+5.17
3498call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3500 +219.38,
-70.11,
-170.01,
+37.19,
-15.98,
-19.77
3501 -70.11,
+25.15,
+53.56,
-10.10,
+4.20,
+6.44
3502 -170.01,
+53.56,
+133.66,
-29.20,
+11.67,
+13.76
3503 +37.19,
-10.10,
-29.20,
+8.48,
-3.42,
-3.51
3504 -15.98,
+4.20,
+11.67,
-3.42,
+4.01,
+3.68
3505 -19.77,
+6.44,
+13.76,
-3.51,
+3.68,
+5.17
3506mul
= matmul(mat, inv)
3508 +1.00,
+0.00,
+0.00,
-0.00,
-0.00,
-0.00
3509 +0.00,
+1.00,
-0.00,
+0.00,
-0.00,
-0.00
3510 +0.00,
-0.00,
+1.00,
-0.00,
-0.00,
-0.00
3511 -0.00,
-0.00,
-0.00,
+1.00,
+0.00,
+0.00
3512 +0.00,
+0.00,
+0.00,
-0.00,
+1.00,
+0.00
3513 +0.00,
-0.00,
+0.00,
-0.00,
+0.00,
+1.00
3522(
+9.00,
+0.00), (
+3.00,
+3.00), (
+3.00,
-3.00)
3523(
+3.00,
-3.00), (
+18.00,
+0.00), (
+8.00,
-6.00)
3524(
+3.00,
+3.00), (
+8.00,
+6.00), (
+43.00,
+0.00)
3527(
+3.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3528(
+1.00,
-1.00), (
+4.00,
+0.00), (
+0.00,
+0.00)
3529(
+1.00,
+1.00), (
+2.00,
+1.00), (
+6.00,
+0.00)
3532call setMatInv(inv, chol, auxil
= choLow, subset
= uppDia)
3534(
+0.13,
+0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
3535(
+0.00,
+0.00), (
+0.07,
+0.00), (
-0.01,
+0.01)
3536(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.03,
+0.00)
3537call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3539(
+0.13,
+0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
3540(
-0.02,
+0.03), (
+0.07,
+0.00), (
-0.01,
+0.01)
3541(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
+0.00)
3542mul
= matmul(mat, inv)
3544(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
3545(
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
3546(
+0.00,
+0.00), (
-0.00,
-0.00), (
+1.00,
-0.00)
3549call setMatInv(inv, chol, auxil
= choLow, subset
= lowDia)
3551(
+0.13,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3552(
-0.02,
+0.03), (
+0.07,
+0.00), (
+0.00,
+0.00)
3553(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
+0.00)
3554call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3556(
+0.13,
+0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
3557(
-0.02,
+0.03), (
+0.07,
+0.00), (
-0.01,
+0.01)
3558(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
+0.00)
3559mul
= matmul(mat, inv)
3561(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
3562(
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
3563(
+0.00,
+0.00), (
-0.00,
-0.00), (
+1.00,
-0.00)
3567(
+9.00,
+0.00), (
+3.00,
+3.00), (
+3.00,
-3.00)
3568(
+3.00,
-3.00), (
+18.00,
+0.00), (
+8.00,
-6.00)
3569(
+3.00,
+3.00), (
+8.00,
+6.00), (
+43.00,
+0.00)
3572(
+3.00,
+0.00), (
+1.00,
+1.00), (
+1.00,
-1.00)
3573(
+0.00,
+0.00), (
+4.00,
+0.00), (
+2.00,
-1.00)
3574(
+0.00,
+0.00), (
+0.00,
+0.00), (
+6.00,
+0.00)
3577call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3579(
+0.13,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3580(
-0.02,
+0.03), (
+0.07,
+0.00), (
+0.00,
+0.00)
3581(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
+0.00)
3582call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3584(
+0.13,
+0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
3585(
-0.02,
+0.03), (
+0.07,
+0.00), (
-0.01,
+0.01)
3586(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
+0.00)
3587mul
= matmul(mat, inv)
3589(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
3590(
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
3591(
+0.00,
+0.00), (
-0.00,
-0.00), (
+1.00,
-0.00)
3594call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3596(
+0.13,
+0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
3597(
+0.00,
+0.00), (
+0.07,
+0.00), (
-0.01,
+0.01)
3598(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.03,
+0.00)
3599call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3601(
+0.13,
+0.00), (
-0.02,
-0.03), (
-0.00,
+0.01)
3602(
-0.02,
+0.03), (
+0.07,
+0.00), (
-0.01,
+0.01)
3603(
-0.00,
-0.01), (
-0.01,
-0.01), (
+0.03,
+0.00)
3604mul
= matmul(mat, inv)
3606(
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
3607(
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
3608(
+0.00,
+0.00), (
-0.00,
-0.00), (
+1.00,
-0.00)
3621call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3624call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3627mul
= matmul(mat, inv)
3632call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3635call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3638mul
= matmul(mat, inv)
3647(
+1.76,
+0.00), (
+1.47,
-0.73), (
+1.77,
-0.83)
3648(
+1.47,
+0.73), (
+3.16,
+0.00), (
+3.14,
-1.05)
3649(
+1.77,
+0.83), (
+3.14,
+1.05), (
+5.03,
+0.00)
3652(
+1.33,
+0.00), (
+1.11,
-0.55), (
+1.34,
-0.63)
3653(
+0.00,
+0.00), (
+1.28,
+0.00), (
+1.03,
-0.85)
3654(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.04,
+0.00)
3657call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3659(
+1.60,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3660(
-0.30,
-0.94), (
+1.62,
+0.00), (
+0.00,
+0.00)
3661(
-0.57,
+0.38), (
-0.75,
-0.62), (
+0.93,
+0.00)
3662call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3664(
+1.60,
+0.00), (
-0.30,
+0.94), (
-0.57,
-0.38)
3665(
-0.30,
-0.94), (
+1.62,
+0.00), (
-0.75,
+0.62)
3666(
-0.57,
+0.38), (
-0.75,
-0.62), (
+0.93,
+0.00)
3667mul
= matmul(mat, inv)
3669(
+1.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
+0.00)
3670(
-0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
3671(
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
3674call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3676(
+1.60,
+0.00), (
-0.30,
+0.94), (
-0.57,
-0.38)
3677(
+0.00,
+0.00), (
+1.62,
+0.00), (
-0.75,
+0.62)
3678(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.93,
+0.00)
3679call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3681(
+1.60,
+0.00), (
-0.30,
+0.94), (
-0.57,
-0.38)
3682(
-0.30,
-0.94), (
+1.62,
+0.00), (
-0.75,
+0.62)
3683(
-0.57,
+0.38), (
-0.75,
-0.62), (
+0.93,
+0.00)
3684mul
= matmul(mat, inv)
3686(
+1.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
+0.00)
3687(
-0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
3688(
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
3695(
+2.12,
+0.00), (
+1.95,
+0.42), (
+2.01,
-1.20), (
+2.19,
+0.16)
3696(
+1.95,
-0.42), (
+5.43,
+0.00), (
+3.63,
-3.20), (
+4.45,
+0.83)
3697(
+2.01,
+1.20), (
+3.63,
+3.20), (
+5.55,
+0.00), (
+4.49,
+3.85)
3698(
+2.19,
-0.16), (
+4.45,
-0.83), (
+4.49,
-3.85), (
+8.94,
+0.00)
3701(
+1.46,
+0.00), (
+1.34,
+0.29), (
+1.38,
-0.82), (
+1.50,
+0.11)
3702(
+0.00,
+0.00), (
+1.89,
+0.00), (
+1.07,
-0.91), (
+1.28,
+0.59)
3703(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+1.67,
+0.68)
3704(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.21,
+0.00)
3707call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3709(
+0.94,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3710(
+0.05,
-0.16), (
+1.44,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3711(
-0.69,
-0.14), (
-0.91,
-1.59), (
+3.23,
+0.00), (
+0.00,
+0.00)
3712(
+0.17,
-0.12), (
+0.41,
+0.50), (
-1.15,
+0.47), (
+0.69,
+0.00)
3713call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3715(
+0.94,
+0.00), (
+0.05,
+0.16), (
-0.69,
+0.14), (
+0.17,
+0.12)
3716(
+0.05,
-0.16), (
+1.44,
+0.00), (
-0.91,
+1.59), (
+0.41,
-0.50)
3717(
-0.69,
-0.14), (
-0.91,
-1.59), (
+3.23,
+0.00), (
-1.15,
-0.47)
3718(
+0.17,
-0.12), (
+0.41,
+0.50), (
-1.15,
+0.47), (
+0.69,
+0.00)
3719mul
= matmul(mat, inv)
3721(
+1.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
3722(
-0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
3723(
-0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
-0.00), (
-0.00,
+0.00)
3724(
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00), (
+1.00,
+0.00)
3727call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3729(
+0.94,
+0.00), (
+0.05,
+0.16), (
-0.69,
+0.14), (
+0.17,
+0.12)
3730(
+0.00,
+0.00), (
+1.44,
+0.00), (
-0.91,
+1.59), (
+0.41,
-0.50)
3731(
+0.00,
+0.00), (
+0.00,
+0.00), (
+3.23,
+0.00), (
-1.15,
-0.47)
3732(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.69,
+0.00)
3733call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3735(
+0.94,
+0.00), (
+0.05,
+0.16), (
-0.69,
+0.14), (
+0.17,
+0.12)
3736(
+0.05,
-0.16), (
+1.44,
+0.00), (
-0.91,
+1.59), (
+0.41,
-0.50)
3737(
-0.69,
-0.14), (
-0.91,
-1.59), (
+3.23,
+0.00), (
-1.15,
-0.47)
3738(
+0.17,
-0.12), (
+0.41,
+0.50), (
-1.15,
+0.47), (
+0.69,
+0.00)
3739mul
= matmul(mat, inv)
3741(
+1.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
3742(
-0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
3743(
-0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
-0.00), (
-0.00,
+0.00)
3744(
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
-0.00), (
+1.00,
+0.00)
3757call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3760call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3763mul
= matmul(mat, inv)
3768call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3771call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3774mul
= matmul(mat, inv)
3783(
+1.12,
+0.00), (
+1.77,
+0.56), (
+1.45,
-0.17)
3784(
+1.77,
-0.56), (
+6.26,
+0.00), (
+5.66,
-0.72)
3785(
+1.45,
+0.17), (
+5.66,
+0.72), (
+6.74,
+0.00)
3788(
+1.06,
+0.00), (
+1.68,
+0.52), (
+1.37,
-0.16)
3789(
+0.00,
+0.00), (
+1.78,
+0.00), (
+1.93,
+0.15)
3790(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.03,
+0.00)
3793call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3795(
+2.53,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3796(
-0.94,
+0.97), (
+1.42,
+0.00), (
+0.00,
+0.00)
3797(
+0.35,
-0.78), (
-1.02,
+0.08), (
+0.94,
+0.00)
3798call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3800(
+2.53,
+0.00), (
-0.94,
-0.97), (
+0.35,
+0.78)
3801(
-0.94,
+0.97), (
+1.42,
+0.00), (
-1.02,
-0.08)
3802(
+0.35,
-0.78), (
-1.02,
+0.08), (
+0.94,
+0.00)
3803mul
= matmul(mat, inv)
3805(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3806(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
3807(
-0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
3810call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3812(
+2.53,
+0.00), (
-0.94,
-0.97), (
+0.35,
+0.78)
3813(
+0.00,
+0.00), (
+1.42,
+0.00), (
-1.02,
-0.08)
3814(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.94,
+0.00)
3815call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3817(
+2.53,
+0.00), (
-0.94,
-0.97), (
+0.35,
+0.78)
3818(
-0.94,
+0.97), (
+1.42,
+0.00), (
-1.02,
-0.08)
3819(
+0.35,
-0.78), (
-1.02,
+0.08), (
+0.94,
+0.00)
3820mul
= matmul(mat, inv)
3822(
+1.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3823(
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
+0.00)
3824(
-0.00,
+0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
3831(
+3.30,
+0.00), (
+3.47,
+1.56), (
+2.85,
+1.35), (
+3.57,
+1.12)
3832(
+3.47,
-1.56), (
+6.36,
+0.00), (
+5.48,
-0.75), (
+6.80,
-0.45)
3833(
+2.85,
-1.35), (
+5.48,
+0.75), (
+7.08,
+0.00), (
+7.44,
+1.21)
3834(
+3.57,
-1.12), (
+6.80,
+0.45), (
+7.44,
-1.21), (
+10.01,
+0.00)
3837(
+1.82,
+0.00), (
+1.91,
+0.86), (
+1.57,
+0.74), (
+1.97,
+0.61)
3838(
+0.00,
+0.00), (
+1.40,
+0.00), (
+1.31,
-0.59), (
+1.79,
+0.04)
3839(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.41,
+0.00), (
+1.11,
+0.42)
3840(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.07,
+0.00)
3843call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3845(
+1.36,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3846(
-0.76,
+0.31), (
+1.19,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3847(
+0.11,
+0.59), (
-0.16,
-0.16), (
+1.11,
+0.00), (
+0.00,
+0.00)
3848(
-0.11,
-0.45), (
-0.36,
+0.07), (
-0.68,
+0.26), (
+0.87,
+0.00)
3849call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3851(
+1.36,
+0.00), (
-0.76,
-0.31), (
+0.11,
-0.59), (
-0.11,
+0.45)
3852(
-0.76,
+0.31), (
+1.19,
+0.00), (
-0.16,
+0.16), (
-0.36,
-0.07)
3853(
+0.11,
+0.59), (
-0.16,
-0.16), (
+1.11,
+0.00), (
-0.68,
-0.26)
3854(
-0.11,
-0.45), (
-0.36,
+0.07), (
-0.68,
+0.26), (
+0.87,
+0.00)
3855mul
= matmul(mat, inv)
3857(
+1.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00)
3858(
-0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
-0.00)
3859(
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
-0.00)
3860(
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
+1.00,
-0.00)
3863call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3865(
+1.36,
+0.00), (
-0.76,
-0.31), (
+0.11,
-0.59), (
-0.11,
+0.45)
3866(
+0.00,
+0.00), (
+1.19,
+0.00), (
-0.16,
+0.16), (
-0.36,
-0.07)
3867(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.11,
+0.00), (
-0.68,
-0.26)
3868(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.87,
+0.00)
3869call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3871(
+1.36,
+0.00), (
-0.76,
-0.31), (
+0.11,
-0.59), (
-0.11,
+0.45)
3872(
-0.76,
+0.31), (
+1.19,
+0.00), (
-0.16,
+0.16), (
-0.36,
-0.07)
3873(
+0.11,
+0.59), (
-0.16,
-0.16), (
+1.11,
+0.00), (
-0.68,
-0.26)
3874(
-0.11,
-0.45), (
-0.36,
+0.07), (
-0.68,
+0.26), (
+0.87,
+0.00)
3875mul
= matmul(mat, inv)
3877(
+1.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00)
3878(
-0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
-0.00)
3879(
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00), (
+0.00,
-0.00)
3880(
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00), (
+1.00,
-0.00)
3893call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3896call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3899mul
= matmul(mat, inv)
3904call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3907call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3910mul
= matmul(mat, inv)
3919(
+1.10,
+0.00), (
+1.96,
-0.08), (
+1.92,
+0.96), (
+1.18,
+1.04), (
+1.81,
+0.67)
3920(
+1.96,
+0.08), (
+5.46,
+0.00), (
+4.79,
+0.67), (
+4.08,
+1.40), (
+5.86,
+0.17)
3921(
+1.92,
-0.96), (
+4.79,
-0.67), (
+9.60,
+0.00), (
+8.43,
+2.10), (
+9.43,
-0.81)
3922(
+1.18,
-1.04), (
+4.08,
-1.40), (
+8.43,
-2.10), (
+11.51,
+0.00), (
+11.15,
-4.08)
3923(
+1.81,
-0.67), (
+5.86,
-0.17), (
+9.43,
+0.81), (
+11.15,
+4.08), (
+16.59,
+0.00)
3926(
+1.05,
+0.00), (
+1.87,
-0.08), (
+1.83,
+0.91), (
+1.13,
+0.99), (
+1.73,
+0.64)
3927(
+0.00,
+0.00), (
+1.39,
+0.00), (
+1.03,
-0.85), (
+1.47,
-0.39), (
+1.92,
-0.83)
3928(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.91,
+0.00), (
+1.90,
+0.24), (
+1.57,
-0.62)
3929(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.81,
+0.00), (
+1.43,
-0.58)
3930(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.89,
+0.00)
3933call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
3935(
+4.39,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3936(
-1.44,
-0.64), (
+0.94,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3937(
-0.47,
+0.71), (
-0.12,
-0.25), (
+0.59,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
3938(
+0.08,
-0.17), (
+0.08,
-0.00), (
-0.30,
+0.06), (
+0.51,
+0.00), (
+0.00,
+0.00)
3939(
+0.24,
+0.11), (
-0.20,
+0.01), (
+0.00,
-0.03), (
-0.22,
-0.09), (
+0.28,
+0.00)
3940call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
3942(
+4.39,
+0.00), (
-1.44,
+0.64), (
-0.47,
-0.71), (
+0.08,
+0.17), (
+0.24,
-0.11)
3943(
-1.44,
-0.64), (
+0.94,
+0.00), (
-0.12,
+0.25), (
+0.08,
+0.00), (
-0.20,
-0.01)
3944(
-0.47,
+0.71), (
-0.12,
-0.25), (
+0.59,
+0.00), (
-0.30,
-0.06), (
+0.00,
+0.03)
3945(
+0.08,
-0.17), (
+0.08,
-0.00), (
-0.30,
+0.06), (
+0.51,
+0.00), (
-0.22,
+0.09)
3946(
+0.24,
+0.11), (
-0.20,
+0.01), (
+0.00,
-0.03), (
-0.22,
-0.09), (
+0.28,
+0.00)
3947mul
= matmul(mat, inv)
3949(
+1.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
+0.00)
3950(
+0.00,
-0.00), (
+1.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
+0.00)
3951(
+0.00,
-0.00), (
-0.00,
+0.00), (
+1.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
3952(
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
+0.00), (
+0.00,
-0.00)
3953(
+0.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
3956call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
3958(
+4.39,
+0.00), (
-1.44,
+0.64), (
-0.47,
-0.71), (
+0.08,
+0.17), (
+0.24,
-0.11)
3959(
+0.00,
+0.00), (
+0.94,
+0.00), (
-0.12,
+0.25), (
+0.08,
+0.00), (
-0.20,
-0.01)
3960(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.59,
+0.00), (
-0.30,
-0.06), (
+0.00,
+0.03)
3961(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.51,
+0.00), (
-0.22,
+0.09)
3962(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.28,
+0.00)
3963call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
3965(
+4.39,
+0.00), (
-1.44,
+0.64), (
-0.47,
-0.71), (
+0.08,
+0.17), (
+0.24,
-0.11)
3966(
-1.44,
-0.64), (
+0.94,
+0.00), (
-0.12,
+0.25), (
+0.08,
+0.00), (
-0.20,
-0.01)
3967(
-0.47,
+0.71), (
-0.12,
-0.25), (
+0.59,
+0.00), (
-0.30,
-0.06), (
+0.00,
+0.03)
3968(
+0.08,
-0.17), (
+0.08,
-0.00), (
-0.30,
+0.06), (
+0.51,
+0.00), (
-0.22,
+0.09)
3969(
+0.24,
+0.11), (
-0.20,
+0.01), (
+0.00,
-0.03), (
-0.22,
-0.09), (
+0.28,
+0.00)
3970mul
= matmul(mat, inv)
3972(
+1.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
+0.00)
3973(
+0.00,
-0.00), (
+1.00,
+0.00), (
+0.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
+0.00)
3974(
+0.00,
-0.00), (
-0.00,
+0.00), (
+1.00,
-0.00), (
-0.00,
-0.00), (
+0.00,
-0.00)
3975(
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
+0.00), (
+0.00,
-0.00)
3976(
+0.00,
-0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
3983(
+2.42,
+0.00), (
+2.19,
-0.41), (
+2.54,
+1.18), (
+1.65,
+0.59), (
+2.11,
-0.71), (
+2.12,
-1.47)
3984(
+2.19,
+0.41), (
+4.49,
+0.00), (
+4.96,
+0.14), (
+3.19,
+1.28), (
+4.50,
+0.39), (
+4.97,
-1.42)
3985(
+2.54,
-1.18), (
+4.96,
-0.14), (
+9.55,
+0.00), (
+5.89,
+1.71), (
+6.88,
-0.38), (
+7.98,
-2.94)
3986(
+1.65,
-0.59), (
+3.19,
-1.28), (
+5.89,
-1.71), (
+6.03,
+0.00), (
+6.86,
-2.84), (
+7.89,
-4.66)
3987(
+2.11,
+0.71), (
+4.50,
-0.39), (
+6.88,
+0.38), (
+6.86,
+2.84), (
+11.06,
+0.00), (
+13.37,
-1.69)
3988(
+2.12,
+1.47), (
+4.97,
+1.42), (
+7.98,
+2.94), (
+7.89,
+4.66), (
+13.37,
+1.69), (
+19.56,
+0.00)
3991(
+1.56,
+0.00), (
+1.40,
-0.26), (
+1.63,
+0.76), (
+1.06,
+0.38), (
+1.36,
-0.46), (
+1.36,
-0.95)
3992(
+0.00,
+0.00), (
+1.57,
+0.00), (
+1.83,
-0.86), (
+1.15,
+0.30), (
+1.58,
+0.43), (
+1.79,
-0.29)
3993(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.49,
+0.00), (
+1.36,
+0.24), (
+1.68,
-0.50), (
+1.98,
-0.93)
3994(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.20,
+0.00), (
+1.24,
-0.65), (
+1.96,
-0.45)
3995(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.15,
+0.00), (
+1.64,
+0.04)
3996(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.42,
+0.00)
3999call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
4001(
+2.77,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
4002(
-1.53,
-1.02), (
+1.87,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
4003(
-0.19,
+1.45), (
-0.68,
-0.85), (
+1.27,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
4004(
+0.18,
-0.17), (
+0.04,
+0.03), (
-0.69,
+0.04), (
+1.82,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00)
4005(
-0.24,
-1.06), (
+0.25,
+0.75), (
-0.34,
+0.38), (
-0.68,
-0.68), (
+1.78,
+0.00), (
+0.00,
+0.00)
4006(
+0.28,
+0.37), (
-0.20,
-0.24), (
+0.18,
-0.27), (
-0.07,
+0.19), (
-0.71,
+0.02), (
+0.50,
+0.00)
4007call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
4009(
+2.77,
+0.00), (
-1.53,
+1.02), (
-0.19,
-1.45), (
+0.18,
+0.17), (
-0.24,
+1.06), (
+0.28,
-0.37)
4010(
-1.53,
-1.02), (
+1.87,
+0.00), (
-0.68,
+0.85), (
+0.04,
-0.03), (
+0.25,
-0.75), (
-0.20,
+0.24)
4011(
-0.19,
+1.45), (
-0.68,
-0.85), (
+1.27,
+0.00), (
-0.69,
-0.04), (
-0.34,
-0.38), (
+0.18,
+0.27)
4012(
+0.18,
-0.17), (
+0.04,
+0.03), (
-0.69,
+0.04), (
+1.82,
+0.00), (
-0.68,
+0.68), (
-0.07,
-0.19)
4013(
-0.24,
-1.06), (
+0.25,
+0.75), (
-0.34,
+0.38), (
-0.68,
-0.68), (
+1.78,
+0.00), (
-0.71,
-0.02)
4014(
+0.28,
+0.37), (
-0.20,
-0.24), (
+0.18,
-0.27), (
-0.07,
+0.19), (
-0.71,
+0.02), (
+0.50,
+0.00)
4015mul
= matmul(mat, inv)
4017(
+1.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
-0.00)
4018(
-0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
4019(
-0.00,
-0.00), (
-0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
4020(
-0.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
4021(
+0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
4022(
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
4025call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
4027(
+2.77,
+0.00), (
-1.53,
+1.02), (
-0.19,
-1.45), (
+0.18,
+0.17), (
-0.24,
+1.06), (
+0.28,
-0.37)
4028(
+0.00,
+0.00), (
+1.87,
+0.00), (
-0.68,
+0.85), (
+0.04,
-0.03), (
+0.25,
-0.75), (
-0.20,
+0.24)
4029(
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.27,
+0.00), (
-0.69,
-0.04), (
-0.34,
-0.38), (
+0.18,
+0.27)
4030(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.82,
+0.00), (
-0.68,
+0.68), (
-0.07,
-0.19)
4031(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+1.78,
+0.00), (
-0.71,
-0.02)
4032(
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
+0.50,
+0.00)
4033call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
4035(
+2.77,
+0.00), (
-1.53,
+1.02), (
-0.19,
-1.45), (
+0.18,
+0.17), (
-0.24,
+1.06), (
+0.28,
-0.37)
4036(
-1.53,
-1.02), (
+1.87,
+0.00), (
-0.68,
+0.85), (
+0.04,
-0.03), (
+0.25,
-0.75), (
-0.20,
+0.24)
4037(
-0.19,
+1.45), (
-0.68,
-0.85), (
+1.27,
+0.00), (
-0.69,
-0.04), (
-0.34,
-0.38), (
+0.18,
+0.27)
4038(
+0.18,
-0.17), (
+0.04,
+0.03), (
-0.69,
+0.04), (
+1.82,
+0.00), (
-0.68,
+0.68), (
-0.07,
-0.19)
4039(
-0.24,
-1.06), (
+0.25,
+0.75), (
-0.34,
+0.38), (
-0.68,
-0.68), (
+1.78,
+0.00), (
-0.71,
-0.02)
4040(
+0.28,
+0.37), (
-0.20,
-0.24), (
+0.18,
-0.27), (
-0.07,
+0.19), (
-0.71,
+0.02), (
+0.50,
+0.00)
4041mul
= matmul(mat, inv)
4043(
+1.00,
-0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
-0.00)
4044(
-0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
4045(
-0.00,
-0.00), (
-0.00,
+0.00), (
+1.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
-0.00,
+0.00)
4046(
-0.00,
-0.00), (
+0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
-0.00), (
+0.00,
-0.00), (
-0.00,
+0.00)
4047(
+0.00,
+0.00), (
-0.00,
+0.00), (
-0.00,
+0.00), (
+0.00,
-0.00), (
+1.00,
+0.00), (
-0.00,
+0.00)
4048(
+0.00,
-0.00), (
+0.00,
+0.00), (
+0.00,
+0.00), (
-0.00,
-0.00), (
+0.00,
+0.00), (
+1.00,
+0.00)
4061call setMatInv(inv, chol, auxil
= choUpp, subset
= lowDia)
4064call setMatCopy(inv(
1:,
2:), rdpack, inv(
2:,
1:), rdpack, lowDia, transHerm)
4067mul
= matmul(mat, inv)
4072call setMatInv(inv, chol, auxil
= choUpp, subset
= uppDia)
4075call setMatCopy(inv(
2:,
1:), rdpack, inv(
1:,
2:), rdpack, uppDia, transHerm)
4078mul
= matmul(mat, inv)