14 type(display_type) :: disp
16 real(RKG) ,
parameter :: DUM
= -huge(
0._RKG)
17 complex(CKG) ,
parameter :: CMPLX_DUMM
= cmplx(
-huge(
0._CKG),
-huge(
0._CKG), CKG)
18 integer(IK) :: ndim, ndum, roff, coff, roffA, coffA
19 character(:, SK),
allocatable :: cform
20 cform
= getFormat([
cmplx(
0.,
0., CKG)], ed
= SK_
'f', signed
= .true.)
26 real(RKG) :: alpha, beta
27 real(RKG),
allocatable :: ref(:,:), mat(:,:), matA(:,:), VecX(:), VecY(:)
29 matA
= reshape( [
real(RKG) ::
0.0,
8.0 &
38 ], shape
= [
9,
2], order
= [
2,
1])
39 mat
= reshape( [
real(RKG) ::
0.0,
1.0,
3.0,
6.0,
10.0,
15.0,
21.0,
28.0 &
40 , DUM,
2.0,
4.0,
7.0,
11.0,
16.0,
22.0,
29.0 &
41 , DUM, DUM,
5.0,
8.0,
12.0,
17.0,
23.0,
30.0 &
42 , DUM, DUM, DUM,
9.0,
13.0,
18.0,
24.0,
31.0 &
43 , DUM, DUM, DUM, DUM,
14.0,
19.0,
25.0,
32.0 &
44 , DUM, DUM, DUM, DUM, DUM,
20.0,
26.0,
33.0 &
45 , DUM, DUM, DUM, DUM, DUM, DUM,
27.0,
34.0 &
46 , DUM, DUM, DUM, DUM, DUM, DUM, DUM,
35.0 &
47 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
48 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
49 ], shape
= [
10,
8], order
= [
2,
1])
50 ref
= reshape( [
real(RKG) ::
64.0,
73.0,
83.0,
94.0,
106.0,
119.0,
133.0,
148.0 &
51 , DUM,
84.0,
96.0,
109.0,
123.0,
138.0,
154.0,
171.0 &
52 , DUM, DUM,
109.0,
124.0,
140.0,
157.0,
175.0,
194.0 &
53 , DUM, DUM, DUM,
139.0,
157.0,
176.0,
196.0,
217.0 &
54 , DUM, DUM, DUM, DUM,
174.0,
195.0,
217.0,
240.0 &
55 , DUM, DUM, DUM, DUM, DUM,
214.0,
238.0,
263.0 &
56 , DUM, DUM, DUM, DUM, DUM, DUM,
259.0,
286.0 &
57 , DUM, DUM, DUM, DUM, DUM, DUM, DUM,
309.0 &
58 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
59 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
60 ], shape
= [
10,
8], order
= [
2,
1])
66 call disp%show(
"alpha = 1._RKG; beta = 1._RKG; ndim = 8; ndum = 2; roff = 0; coff = 0; roffA = 0; coffA = 0;")
67 alpha
= 1._RKG; beta
= 1._RKG; ndim
= 8; ndum
= 2; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
68 call disp%show(
"call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)")
69 call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
76 matA
= reshape( [
real(RKG) ::
0.0,
3.0,
6.0,
9.0,
12.0,
15.0,
18.0,
21.0 &
77 ,
1.0,
4.0,
7.0,
10.0,
13.0,
16.0,
19.0,
22.0 &
78 ,
2.0,
5.0,
8.0,
11.0,
14.0,
17.0,
20.0,
23.0 &
79 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
80 ], shape
= [
4,
8], order
= [
2,
1])
81 mat
= reshape( [
real(RKG) ::
0.0, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
82 ,
1.0,
8.0, DUM, DUM, DUM, DUM, DUM, DUM
&
83 ,
2.0,
9.0,
15.0, DUM, DUM, DUM, DUM, DUM
&
84 ,
3.0,
10.0,
16.0,
21.0, DUM, DUM, DUM, DUM
&
85 ,
4.0,
11.0,
17.0,
22.0,
26.0, DUM, DUM, DUM
&
86 ,
5.0,
12.0,
18.0,
23.0,
27.0,
30.0, DUM, DUM
&
87 ,
6.0,
13.0,
19.0,
24.0,
28.0,
31.0,
33.0, DUM
&
88 ,
7.0,
14.0,
20.0,
25.0,
29.0,
32.0,
34.0,
35.0 &
89 ], shape
= [
8,
8], order
= [
2,
1])
90 ref
= reshape( [
real(RKG) ::
5.0, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
91 ,
15.0,
58.0, DUM, DUM, DUM, DUM, DUM, DUM
&
92 ,
25.0,
95.0,
164.0, DUM, DUM, DUM, DUM, DUM
&
93 ,
35.0,
132.0,
228.0,
323.0, DUM, DUM, DUM, DUM
&
94 ,
45.0,
169.0,
292.0,
414.0,
535.0, DUM, DUM, DUM
&
95 ,
55.0,
206.0,
356.0,
505.0,
653.0,
800.0, DUM, DUM
&
96 ,
65.0,
243.0,
420.0,
596.0,
771.0,
945.0,
1118.0, DUM
&
97 ,
75.0,
280.0,
484.0,
687.0,
889.0,
1090.0,
1290.0,
1489.0 &
98 ], shape
= [
8,
8], order
= [
2,
1])
104 call disp%show(
"alpha = 1._RKG; beta = 1._RKG; ndim = 8; ndum = 3; roff = 0; coff = 0; roffA = 0; coffA = 0;")
105 alpha
= 1._RKG; beta
= 1._RKG; ndim
= 8; ndum
= 3; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
106 call disp%show(
"call setMatUpdateTriang(mat, lowDia, matA, transSymm, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)")
107 call setMatUpdateTriang(mat, lowDia, matA, transSymm, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
117 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
118 call disp%show(
"!Complex Symmetric matrix update.")
119 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
124 complex(CKG) :: alpha, beta
125 complex(CKG),
allocatable :: ref(:,:), mat(:,:), matA(:,:), VecX(:), VecY(:)
127 matA
= reshape( [
complex(CKG) :: (
2.0,
0.0), (
3.0,
2.0), (
4.0,
1.0), (
1.0,
7.0), (
0.0,
0.0)
&
128 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0), (
2.0,
4.0), (
1.0,
2.0)
&
129 , (
1.0,
3.0), (
2.0,
1.0), (
6.0,
0.0), (
3.0,
2.0), (
2.0,
2.0)
&
130 ], shape
= [
3,
5], order
= [
2,
1])
131 mat
= reshape( [
complex(CKG) :: (
2.0,
1.0), (
1.0,
9.0), (
4.0,
5.0)
&
132 , CMPLX_DUMM, (
3.0,
1.0), (
6.0,
7.0)
&
133 , CMPLX_DUMM, CMPLX_DUMM, (
8.0,
1.0)
&
134 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
135 ], shape
= [
4,
3], order
= [
2,
1])
136 ref
= reshape( [
complex(CKG) :: (
-57.0,
13.0), (
-63.0,
79.0), (
-24.0,
70.0)
&
137 , CMPLX_DUMM, (
-28.0,
90.0), (
-55.0,
103.0)
&
138 , CMPLX_DUMM, CMPLX_DUMM, (
13.0,
75.0)
&
139 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
140 ], shape
= [
4,
3], order
= [
2,
1])
143 call disp%show( matA ,
format = cform )
145 call disp%show( mat ,
format = cform )
146 call disp%show(
"alpha = (1._CKG, 1._CKG); beta = (1._CKG, 1._CKG); ndim = 3; ndum = 5; roff = 0; coff = 0; roffA = 0; coffA = 0;")
147 alpha
= (
1._CKG,
1._CKG); beta
= (
1._CKG,
1._CKG); ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
148 call disp%show(
"call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)")
149 call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
151 call disp%show( mat ,
format = cform )
153 call disp%show( mat
- ref ,
format = cform )
159 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
160 call disp%show(
"!Complex hermitian matrix update.")
161 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
166 real(CKG) :: alpha, beta
167 complex(CKG),
allocatable :: ref(:,:), mat(:,:), matA(:,:), VecX(:), VecY(:)
169 matA
= reshape( [
complex(CKG) :: (
2.0,
0.0), (
3.0,
2.0), (
4.0,
1.0)
&
170 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0)
&
171 , (
1.0,
3.0), (
2.0,
1.0), (
6.0,
0.0)
&
172 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0)
&
173 , (
1.0,
9.0), (
3.0,
0.0), (
6.0,
7.0)
&
174 ], shape
= [
5,
3], order
= [
2,
1])
175 mat
= reshape( [
complex(CKG) :: (
6.0, DUM), CMPLX_DUMM, CMPLX_DUMM
&
176 , (
3.0,
4.0), (
10.0, DUM), CMPLX_DUMM
&
177 , (
9.0,
1.0), (
12.0,
2.0), (
3.0, DUM)
&
178 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
179 ], shape
= [
4,
3], order
= [
2,
1])
180 ref
= reshape( [
complex(CKG) :: (
138.0,
0.0), CMPLX_DUMM, CMPLX_DUMM
&
181 , (
65.0,
80.0), (
165.0,
0.0), CMPLX_DUMM
&
182 , (
134.0,
46.0), (
88.0,
-88.0), (
199.0,
0.0)
&
183 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
184 ], shape
= [
4,
3], order
= [
2,
1])
187 call disp%show( matA ,
format = cform )
189 call disp%show( mat ,
format = cform )
190 call disp%show(
"alpha = (1._CKG, 1._CKG); beta = (1._CKG, 1._CKG); ndim = 3; ndum = 5; roff = 0; coff = 0; roffA = 0; coffA = 0;")
191 alpha
= (
1._CKG,
1._CKG); beta
= (
1._CKG,
1._CKG); ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
192 call disp%show(
"call setMatUpdateTriang(mat, lowDia, matA, transHerm, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)")
193 call setMatUpdateTriang(mat, lowDia, matA, transHerm, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
195 call disp%show( mat ,
format = cform )
196 call disp%show(
"mat - ref ! Note the imaginary parts of the diagonals are set to zero on output.")
197 call disp%show( mat
- ref ,
format = cform )
201 matA
= reshape( [
complex(CKG) :: (
2.0,
0.0), (
3.0,
2.0), (
4.0,
1.0)
&
202 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0)
&
203 , (
1.0,
3.0), (
2.0,
1.0), (
6.0,
0.0)
&
204 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0)
&
205 , (
1.0,
9.0), (
3.0,
0.0), (
6.0,
7.0)
&
206 ], shape
= [
3,
5], order
= [
1,
2])
208 mat
= reshape( [
complex(CKG) :: (
6.0, DUM), CMPLX_DUMM, CMPLX_DUMM
&
209 , (
3.0,
4.0), (
10.0, DUM), CMPLX_DUMM
&
210 , (
9.0,
1.0), (
12.0,
2.0), (
3.0, DUM)
&
211 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
212 ], shape
= [
4,
3], order
= [
2,
1])
213 ref
= reshape( [
complex(CKG) :: (
138.0,
0.0), CMPLX_DUMM, CMPLX_DUMM
&
214 , (
65.0,
80.0), (
165.0,
0.0), CMPLX_DUMM
&
215 , (
134.0,
46.0), (
88.0,
-88.0), (
199.0,
0.0)
&
216 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
217 ], shape
= [
4,
3], order
= [
2,
1])
220 call disp%show( matA ,
format = cform )
222 call disp%show( mat ,
format = cform )
223 call disp%show(
"alpha = 1._CKG; beta = 1._CKG; ndim = 3; ndum = 5; roff = 0; coff = 0; roffA = 0; coffA = 0;")
224 alpha
= 1._CKG; beta
= 1._CKG; ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
225 call disp%show(
"call setMatUpdateTriang(mat, lowDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)")
226 call setMatUpdateTriang(mat, lowDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
228 call disp%show( mat ,
format = cform )
229 call disp%show(
"mat - ref ! Note the imaginary parts of the diagonals are set to zero on output.")
230 call disp%show( mat
- ref ,
format = cform )
234 matA
= reshape( [
complex(CKG) :: (
2.0,
0.0), (
3.0,
2.0), (
4.0,
1.0)
&
235 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0)
&
236 , (
1.0,
3.0), (
2.0,
1.0), (
6.0,
0.0)
&
237 , (
3.0,
3.0), (
8.0,
0.0), (
2.0,
5.0)
&
238 , (
1.0,
9.0), (
3.0,
0.0), (
6.0,
7.0)
&
241 mat
= reshape( [
complex(CKG) :: (
6.0, DUM), CMPLX_DUMM, CMPLX_DUMM
&
242 , (
3.0,
4.0), (
10.0, DUM), CMPLX_DUMM
&
243 , (
9.0,
1.0), (
12.0,
2.0), (
3.0, DUM)
&
244 , CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
246 ref
= reshape( [
complex(CKG) :: (
138., .
0), (
65.,
-72.), (
134.,
-44.), CMPLX_DUMM
&
247 , CMPLX_DUMM, (
165.,
0.), (
88.,
92.), CMPLX_DUMM
&
248 , CMPLX_DUMM, CMPLX_DUMM, (
199.,
0.), CMPLX_DUMM
&
249 ], shape
= [
3,
4], order
= [
2,
1])
252 call disp%show( matA ,
format = cform )
254 call disp%show( mat ,
format = cform )
255 call disp%show(
"alpha = 1._CKG; beta = 1._CKG; ndim = 3; ndum = 5; roff = 0; coff = 0; roffA = 0; coffA = 0;")
256 alpha
= 1._CKG; beta
= 1._CKG; ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
257 call disp%show(
"call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)")
258 call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
260 call disp%show( mat ,
format = cform )
261 call disp%show(
"mat - ref ! Note the imaginary parts of the diagonals are set to zero on output.")
262 call disp%show( mat
- ref ,
format = cform )
This is a generic method of the derived type display_type with pass attribute.
This is a generic method of the derived type display_type with pass attribute.
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 CKS
The single-precision complex kind in Fortran mode. On most platforms, this is a 32-bit real 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.
Generate and return an object of type display_type.
3+0.00000000,
+8.00000000
4+1.00000000,
+9.00000000
5+2.00000000,
+10.0000000
6+3.00000000,
+11.0000000
7+4.00000000,
+12.0000000
8+5.00000000,
+13.0000000
9+6.00000000,
+14.0000000
10+7.00000000,
+15.0000000
11-0.340282347E+39,
-0.340282347E+39
13+0.00000000,
+1.00000000,
+3.00000000,
+6.00000000,
+10.0000000,
+15.0000000,
+21.0000000,
+28.0000000
14-0.340282347E+39,
+2.00000000,
+4.00000000,
+7.00000000,
+11.0000000,
+16.0000000,
+22.0000000,
+29.0000000
15-0.340282347E+39,
-0.340282347E+39,
+5.00000000,
+8.00000000,
+12.0000000,
+17.0000000,
+23.0000000,
+30.0000000
16-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+9.00000000,
+13.0000000,
+18.0000000,
+24.0000000,
+31.0000000
17-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+14.0000000,
+19.0000000,
+25.0000000,
+32.0000000
18-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+20.0000000,
+26.0000000,
+33.0000000
19-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+27.0000000,
+34.0000000
20-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+35.0000000
21-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
22-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
23alpha
= 1._RKG; beta
= 1._RKG; ndim
= 8; ndum
= 2; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
24call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
26+64.0000000,
+73.0000000,
+83.0000000,
+94.0000000,
+106.000000,
+119.000000,
+133.000000,
+148.000000
27-0.340282347E+39,
+84.0000000,
+96.0000000,
+109.000000,
+123.000000,
+138.000000,
+154.000000,
+171.000000
28-0.340282347E+39,
-0.340282347E+39,
+109.000000,
+124.000000,
+140.000000,
+157.000000,
+175.000000,
+194.000000
29-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+139.000000,
+157.000000,
+176.000000,
+196.000000,
+217.000000
30-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+174.000000,
+195.000000,
+217.000000,
+240.000000
31-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+214.000000,
+238.000000,
+263.000000
32-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+259.000000,
+286.000000
33-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
+309.000000
34-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
35-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
37+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
38+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
39+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
40+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
41+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
42+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
43+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
44+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
45+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
46+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
50+0.00000000,
+3.00000000,
+6.00000000,
+9.00000000,
+12.0000000,
+15.0000000,
+18.0000000,
+21.0000000
51+1.00000000,
+4.00000000,
+7.00000000,
+10.0000000,
+13.0000000,
+16.0000000,
+19.0000000,
+22.0000000
52+2.00000000,
+5.00000000,
+8.00000000,
+11.0000000,
+14.0000000,
+17.0000000,
+20.0000000,
+23.0000000
53-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
55+0.00000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
56+1.00000000,
+8.00000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
57+2.00000000,
+9.00000000,
+15.0000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
58+3.00000000,
+10.0000000,
+16.0000000,
+21.0000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
59+4.00000000,
+11.0000000,
+17.0000000,
+22.0000000,
+26.0000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
60+5.00000000,
+12.0000000,
+18.0000000,
+23.0000000,
+27.0000000,
+30.0000000,
-0.340282347E+39,
-0.340282347E+39
61+6.00000000,
+13.0000000,
+19.0000000,
+24.0000000,
+28.0000000,
+31.0000000,
+33.0000000,
-0.340282347E+39
62+7.00000000,
+14.0000000,
+20.0000000,
+25.0000000,
+29.0000000,
+32.0000000,
+34.0000000,
+35.0000000
63alpha
= 1._RKG; beta
= 1._RKG; ndim
= 8; ndum
= 3; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
64call setMatUpdateTriang(mat, lowDia, matA, transSymm, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
66+5.00000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
67+15.0000000,
+58.0000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
68+25.0000000,
+95.0000000,
+164.000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
69+35.0000000,
+132.000000,
+228.000000,
+323.000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
70+45.0000000,
+169.000000,
+292.000000,
+414.000000,
+535.000000,
-0.340282347E+39,
-0.340282347E+39,
-0.340282347E+39
71+55.0000000,
+206.000000,
+356.000000,
+505.000000,
+653.000000,
+800.000000,
-0.340282347E+39,
-0.340282347E+39
72+65.0000000,
+243.000000,
+420.000000,
+596.000000,
+771.000000,
+945.000000,
+1118.00000,
-0.340282347E+39
73+75.0000000,
+280.000000,
+484.000000,
+687.000000,
+889.000000,
+1090.00000,
+1290.00000,
+1489.00000
75+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
76+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
77+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
78+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
79+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
80+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
81+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
82+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000,
+0.00000000
91(
+2.000000,
+0.000000), (
+3.000000,
+2.000000), (
+4.000000,
+1.000000), (
+1.000000,
+7.000000), (
+0.000000,
+0.000000)
92(
+3.000000,
+3.000000), (
+8.000000,
+0.000000), (
+2.000000,
+5.000000), (
+2.000000,
+4.000000), (
+1.000000,
+2.000000)
93(
+1.000000,
+3.000000), (
+2.000000,
+1.000000), (
+6.000000,
+0.000000), (
+3.000000,
+2.000000), (
+2.000000,
+2.000000)
95(
+2.000000,
+1.000000), (
+1.000000,
+9.000000), (
+4.000000,
+5.000000)
96(
***************,
***************), (
+3.000000,
+1.000000), (
+6.000000,
+7.000000)
97(
***************,
***************), (
***************,
***************), (
+8.000000,
+1.000000)
98(
***************,
***************), (
***************,
***************), (
***************,
***************)
99alpha
= (
1._CKG,
1._CKG); beta
= (
1._CKG,
1._CKG); ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
100call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
102(
-57.000000,
+13.000000), (
-63.000000,
+79.000000), (
-24.000000,
+70.000000)
103(
***************,
***************), (
-28.000000,
+90.000000), (
-55.000000,
+103.000000)
104(
***************,
***************), (
***************,
***************), (
+13.000000,
+75.000000)
105(
***************,
***************), (
***************,
***************), (
***************,
***************)
107(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
108(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
109(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
110(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
119(
+2.000000,
+0.000000), (
+3.000000,
+2.000000), (
+4.000000,
+1.000000)
120(
+3.000000,
+3.000000), (
+8.000000,
+0.000000), (
+2.000000,
+5.000000)
121(
+1.000000,
+3.000000), (
+2.000000,
+1.000000), (
+6.000000,
+0.000000)
122(
+3.000000,
+3.000000), (
+8.000000,
+0.000000), (
+2.000000,
+5.000000)
123(
+1.000000,
+9.000000), (
+3.000000,
+0.000000), (
+6.000000,
+7.000000)
125(
+6.000000,
***************), (
***************,
***************), (
***************,
***************)
126(
+3.000000,
+4.000000), (
+10.000000,
***************), (
***************,
***************)
127(
+9.000000,
+1.000000), (
+12.000000,
+2.000000), (
+3.000000,
***************)
128(
***************,
***************), (
***************,
***************), (
***************,
***************)
129alpha
= (
1._CKG,
1._CKG); beta
= (
1._CKG,
1._CKG); ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
130call setMatUpdateTriang(mat, lowDia, matA, transHerm, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
132(
+138.000000,
+0.000000), (
***************,
***************), (
***************,
***************)
133(
+65.000000,
+80.000000), (
+165.000000,
+0.000000), (
***************,
***************)
134(
+134.000000,
+46.000000), (
+88.000000,
-88.000000), (
+199.000000,
+0.000000)
135(
***************,
***************), (
***************,
***************), (
***************,
***************)
137(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
138(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
139(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
140(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
144(
+2.000000,
-0.000000), (
+3.000000,
-3.000000), (
+1.000000,
-3.000000), (
+3.000000,
-3.000000), (
+1.000000,
-9.000000)
145(
+3.000000,
-2.000000), (
+8.000000,
-0.000000), (
+2.000000,
-1.000000), (
+8.000000,
-0.000000), (
+3.000000,
-0.000000)
146(
+4.000000,
-1.000000), (
+2.000000,
-5.000000), (
+6.000000,
-0.000000), (
+2.000000,
-5.000000), (
+6.000000,
-7.000000)
148(
+6.000000,
***************), (
***************,
***************), (
***************,
***************)
149(
+3.000000,
+4.000000), (
+10.000000,
***************), (
***************,
***************)
150(
+9.000000,
+1.000000), (
+12.000000,
+2.000000), (
+3.000000,
***************)
151(
***************,
***************), (
***************,
***************), (
***************,
***************)
152alpha
= 1._CKG; beta
= 1._CKG; ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
153call setMatUpdateTriang(mat, lowDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
155(
+138.000000,
+0.000000), (
***************,
***************), (
***************,
***************)
156(
+65.000000,
+80.000000), (
+165.000000,
+0.000000), (
***************,
***************)
157(
+134.000000,
+46.000000), (
+88.000000,
-88.000000), (
+199.000000,
+0.000000)
158(
***************,
***************), (
***************,
***************), (
***************,
***************)
160(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
161(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
162(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
163(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
167(
+2.000000,
-0.000000), (
+3.000000,
-3.000000), (
+1.000000,
-3.000000), (
+3.000000,
-3.000000), (
+1.000000,
-9.000000)
168(
+3.000000,
-2.000000), (
+8.000000,
-0.000000), (
+2.000000,
-1.000000), (
+8.000000,
-0.000000), (
+3.000000,
-0.000000)
169(
+4.000000,
-1.000000), (
+2.000000,
-5.000000), (
+6.000000,
-0.000000), (
+2.000000,
-5.000000), (
+6.000000,
-7.000000)
171(
+6.000000,
***************), (
+3.000000,
+4.000000), (
+9.000000,
+1.000000), (
***************,
***************)
172(
***************,
***************), (
+10.000000,
***************), (
+12.000000,
+2.000000), (
***************,
***************)
173(
***************,
***************), (
***************,
***************), (
+3.000000,
***************), (
***************,
***************)
174alpha
= 1._CKG; beta
= 1._CKG; ndim
= 3; ndum
= 5; roff
= 0; coff
= 0; roffA
= 0; coffA
= 0;
175call setMatUpdateTriang(mat, uppDia, matA, alpha, beta, ndim, ndum, roff, coff, roffA, coffA)
177(
+138.000000,
+0.000000), (
+65.000000,
-72.000000), (
+134.000000,
-44.000000), (
***************,
***************)
178(
***************,
***************), (
+165.000000,
+0.000000), (
+88.000000,
+92.000000), (
***************,
***************)
179(
***************,
***************), (
***************,
***************), (
+199.000000,
+0.000000), (
***************,
***************)
181(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
182(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)
183(
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000), (
+0.000000,
+0.000000)