12 type(display_type) :: disp
14 real(RKG) ,
parameter :: dummyr
= -huge(
0._RKG)
15 complex(CKG),
parameter :: dumm_cmplx_value
= cmplx(
-huge(
0._CKG),
-huge(
0._CKG), CKG)
21 real(RKG),
allocatable :: RefA(:,:), mat(:,:), vecA(:), vecB(:)
23 mat
= reshape( [
1._RKG,
2._RKG,
3._RKG &
24 ,
2._RKG,
2._RKG,
4._RKG &
25 ,
3._RKG,
2._RKG,
2._RKG &
26 ,
4._RKG,
2._RKG,
1._RKG &
27 ], shape
= [
4,
3], order
= [
2,
1])
28 RefA
= reshape( [
4._RKG,
8._RKG,
12._RKG &
29 ,
4._RKG,
6._RKG,
10._RKG &
30 ,
4._RKG,
4._RKG,
5._RKG &
31 ,
8._RKG,
10._RKG,
13._RKG &
32 ], shape
= [
4,
3], order
= [
2,
1])
34 call disp%show(
"dummyr ! some dummy value to illustrate functionality of `incA` and `incB` arguments.")
38 call disp%show(
"vecA = [3._RKG, 2._RKG, 1._RKG, 4._RKG]")
39 vecA
= [
3._RKG,
2._RKG,
1._RKG,
4._RKG]
40 call disp%show(
"vecB = [1._RKG, dummyr, 2._RKG, dummyr, 3._RKG]")
41 vecB
= [
1._RKG, dummyr,
2._RKG, dummyr,
3._RKG]
42 call disp%show(
"call setMatUpdateR1(mat, vecA, vecB, incA = 1_IK, incB = 2_IK, roff = 0_IK)")
43 call setMatUpdateR1(mat, vecA, vecB, incA
= 1_IK, incB
= 2_IK, roff
= 0_IK)
50 mat
= reshape( [
1._RKG,
2._RKG,
3._RKG &
51 ,
2._RKG,
2._RKG,
4._RKG &
52 ,
3._RKG,
2._RKG,
2._RKG &
53 ,
4._RKG,
2._RKG,
1._RKG &
54 , dummyr, dummyr, dummyr
&
55 , dummyr, dummyr, dummyr
&
56 , dummyr, dummyr, dummyr
&
57 , dummyr, dummyr, dummyr
&
58 , dummyr, dummyr, dummyr
&
59 , dummyr, dummyr, dummyr
&
60 ], shape
= [
10,
3], order
= [
2,
1])
62 call disp%show(
"dummyr ! some dummy value to illustrate functionality of `incA` and `incB` arguments.")
66 call disp%show(
"vecA = [3._RKG, 2._RKG, 1._RKG, 4._RKG]")
67 vecA
= [
3._RKG,
2._RKG,
1._RKG,
4._RKG]
68 call disp%show(
"vecB = [3._RKG, dummyr, 2._RKG, dummyr, 1._RKG]")
69 vecB
= [
3._RKG, dummyr,
2._RKG, dummyr,
1._RKG]
70 call disp%show(
"call setMatUpdateR1(mat, vecA, vecB, incA = 1_IK, incB = -2_IK, roff = 0_IK)")
71 call setMatUpdateR1(mat, vecA, vecB, incA
= 1_IK, incB
= -2_IK, roff
= 0_IK)
76 mat
= reshape( [ dummyr, dummyr, dummyr
&
77 , dummyr, dummyr, dummyr
&
78 , dummyr, dummyr, dummyr
&
79 ,
1._RKG,
2._RKG,
3._RKG &
80 ,
2._RKG,
2._RKG,
4._RKG &
81 ,
3._RKG,
2._RKG,
2._RKG &
82 ,
4._RKG,
2._RKG,
1._RKG &
83 , dummyr, dummyr, dummyr
&
84 , dummyr, dummyr, dummyr
&
85 , dummyr, dummyr, dummyr
&
86 ], shape
= [
10,
3], order
= [
2,
1])
88 call disp%show(
"dummyr ! some dummy value to illustrate functionality of `incA` and `incB` arguments.")
92 call disp%show(
"vecA = [3._RKG, dummyr, dummyr, 2._RKG, dummyr, dummyr, 1._RKG, dummyr, dummyr, 4._RKG]")
93 vecA
= [
3._RKG, dummyr, dummyr,
2._RKG, dummyr, dummyr,
1._RKG, dummyr, dummyr,
4._RKG]
94 call disp%show(
"vecB = [1._RKG, 2._RKG, 3._RKG]")
95 vecB
= [
1._RKG,
2._RKG,
3._RKG]
96 call disp%show(
"call setMatUpdateR1(mat, vecA, vecB, incA = 3_IK, incB = 1_IK, roff = 3_IK)")
97 call setMatUpdateR1(mat, vecA, vecB, incA
= 3_IK, incB
= 1_IK, roff
= 3_IK)
102 mat
= reshape( [
1._RKG,
2._RKG,
3._RKG &
103 ,
2._RKG,
2._RKG,
4._RKG &
104 ,
3._RKG,
2._RKG,
2._RKG &
105 ,
4._RKG,
2._RKG,
1._RKG &
106 ], shape
= [
4,
3], order
= [
2,
1])
107 RefA
= reshape( [
7._RKG,
14._RKG,
21._RKG &
108 ,
6._RKG,
10._RKG,
16._RKG &
109 ,
5._RKG,
6._RKG,
8._RKG &
110 ,
12._RKG,
18._RKG,
25._RKG &
111 ], shape
= [
4,
3], order
= [
2,
1])
113 call disp%show(
"dummyr ! some dummy value to illustrate functionality of `incA` and `incB` arguments.")
117 call disp%show(
"vecA = [3._RKG, 2._RKG, 1._RKG, 4._RKG]")
118 vecA
= [
3._RKG,
2._RKG,
1._RKG,
4._RKG]
119 call disp%show(
"vecB = [1._RKG, dummyr, 2._RKG, dummyr, 3._RKG]")
120 vecB
= [
1._RKG, dummyr,
2._RKG, dummyr,
3._RKG]
121 call disp%show(
"call setMatUpdateR1(mat, vecA, vecB, alpha = 2._RKG, incA = 1_IK, incB = 2_IK, roff = 0_IK)")
122 call setMatUpdateR1(mat, vecA, vecB, alpha
= 2._RKG, incA
= 1_IK, incB
= 2_IK, roff
= 0_IK)
132 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
133 call disp%show(
"!Complex matrix update of rank one.")
134 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
139 complex(CKG),
allocatable :: RefA(:,:), mat(:,:), vecA(:), vecB(:)
141 mat
= reshape( [ (
1._CKG,
2._CKG), (
3._CKG,
5._CKG), (
2._CKG,
0._CKG)
&
142 , (
2._CKG,
3._CKG), (
7._CKG,
9._CKG), (
4._CKG,
8._CKG)
&
143 , (
7._CKG,
4._CKG), (
1._CKG,
4._CKG), (
6._CKG,
0._CKG)
&
144 , (
8._CKG,
2._CKG), (
2._CKG,
5._CKG), (
8._CKG,
0._CKG)
&
145 , (
9._CKG,
1._CKG), (
3._CKG,
6._CKG), (
1._CKG,
0._CKG)
&
146 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
147 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
148 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
149 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
150 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
151 ], shape
= [
10,
3], order
= [
2,
1])
152 RefA
= reshape( [ (
-2._CKG,
6._CKG), (
7._CKG,
13._CKG), (
5._CKG,
1._CKG)
&
153 , (
6._CKG,
11._CKG), (
23._CKG,
9._CKG), (
8._CKG,
4._CKG)
&
154 , (
6._CKG,
7._CKG), (
5._CKG,
8._CKG), (
8._CKG,
0._CKG)
&
155 , (
3._CKG,
12._CKG), (
14._CKG,
21._CKG), (
15._CKG,
1._CKG)
&
156 , (
11._CKG,
5._CKG), (
11._CKG,
6._CKG), (
3._CKG,
-2._CKG)
&
157 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
158 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
159 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
160 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
161 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
162 ], shape
= [
10,
3], order
= [
2,
1])
165 call disp%show(
"dummyr ! some dummy value to illustrate functionality of `incA` and `incB` arguments.")
169 call disp%show(
"vecA = [(1._CKG, 2._CKG), (4._CKG, 0._CKG), (1._CKG, 1._CKG), (3._CKG, 4._CKG), (2._CKG, 0._CKG)]")
170 vecA
= [(
1._CKG,
2._CKG), (
4._CKG,
0._CKG), (
1._CKG,
1._CKG), (
3._CKG,
4._CKG), (
2._CKG,
0._CKG)]
171 call disp%show(
"vecB = [(1._CKG, 2._CKG), (4._CKG, 0._CKG), (1._CKG, -1._CKG)]")
172 vecB
= [(
1._CKG,
2._CKG), (
4._CKG,
0._CKG), (
1._CKG,
-1._CKG)]
173 call disp%show(
"call setMatUpdateR1(mat, vecA, vecB, incA = 1_IK, incB = 1_IK, roff = 0_IK)")
174 call setMatUpdateR1(mat, vecA, vecB, incA
= 1_IK, incB
= 1_IK, roff
= 0_IK)
180 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
181 call disp%show(
"!Complex matrix update of rank one with complex conjugate transpose of `Y`.")
182 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
186 mat
= reshape( [ (
1._CKG,
2._CKG), (
3._CKG,
5._CKG), (
2._CKG,
0._CKG)
&
187 , (
2._CKG,
3._CKG), (
7._CKG,
9._CKG), (
4._CKG,
8._CKG)
&
188 , (
7._CKG,
4._CKG), (
1._CKG,
4._CKG), (
6._CKG,
0._CKG)
&
189 , (
8._CKG,
2._CKG), (
2._CKG,
5._CKG), (
8._CKG,
0._CKG)
&
190 , (
9._CKG,
1._CKG), (
3._CKG,
6._CKG), (
1._CKG,
0._CKG)
&
191 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
192 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
193 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
194 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
195 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
196 ], shape
= [
10,
3], order
= [
2,
1])
197 RefA
= reshape( [ (
6._CKG,
2._CKG), (
7._CKG,
13._CKG), (
1._CKG,
3._CKG)
&
198 , (
6._CKG,
-5._CKG), (
23._CKG,
9._CKG), (
8._CKG,
12._CKG)
&
199 , (
10._CKG,
3._CKG), (
5._CKG,
8._CKG), (
6._CKG,
2._CKG)
&
200 , (
19._CKG,
0._CKG), (
14._CKG,
21._CKG), (
7._CKG,
7._CKG)
&
201 , (
11._CKG,
-3._CKG), (
11._CKG,
6._CKG), (
3._CKG,
2._CKG)
&
202 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
203 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
204 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
205 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
206 , dumm_cmplx_value, dumm_cmplx_value, dumm_cmplx_value
&
207 ], shape
= [
10,
3], order
= [
2,
1])
210 call disp%show(
"dummyr ! some dummy value to illustrate functionality of `incA` and `incB` arguments.")
214 call disp%show(
"vecA = [(1._CKG, 2._CKG), (4._CKG, 0._CKG), (1._CKG, 1._CKG), (3._CKG, 4._CKG), (2._CKG, 0._CKG)]")
215 vecA
= [(
1._CKG,
2._CKG), (
4._CKG,
0._CKG), (
1._CKG,
1._CKG), (
3._CKG,
4._CKG), (
2._CKG,
0._CKG)]
216 call disp%show(
"vecB = [(1._CKG, 2._CKG), (4._CKG, 0._CKG), (1._CKG, -1._CKG)]")
217 vecB
= [(
1._CKG,
2._CKG), (
4._CKG,
0._CKG), (
1._CKG,
-1._CKG)]
218 call disp%show(
"call setMatUpdateR1(mat, vecA, vecB, operationB = transHerm, incA = 1_IK, incB = 1_IK, roff = 0_IK)")
219 call setMatUpdateR1(mat, vecA, vecB, operationB
= transHerm, incA
= 1_IK, incB
= 1_IK, roff
= 0_IK)
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.