11 type(display_type) :: disp
13 character(:, SK),
allocatable :: cform, gform
14 real(TKG) ,
parameter :: DUM
= -huge(
0._TKG)
15 complex(TKG) ,
parameter :: CMPLX_DUMM
= cmplx(
-huge(
0._TKG),
-huge(
0._TKG), TKG)
16 cform
= getFormat([CMPLX_DUMM], ed
= SK_
'f', signed
= .true.)
17 gform
= getFormat([DUM], ed
= SK_
'f', signed
= .true.)
22 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
23 call disp%show(
"!Cholesky factorization of real positive-definite square matrix update.")
24 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
28 real(TKG),
allocatable :: cholref(:,:), cholmat(:,:)
30 cholmat
= reshape([
real(TKG) ::
1.0, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
31 ,
1.0,
2.0, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
32 ,
1.0,
2.0,
3.0, DUM, DUM, DUM, DUM, DUM, DUM
&
33 ,
1.0,
2.0,
3.0,
4.0, DUM, DUM, DUM, DUM, DUM
&
34 ,
1.0,
2.0,
3.0,
4.0,
5.0, DUM, DUM, DUM, DUM
&
35 ,
1.0,
2.0,
3.0,
4.0,
5.0,
6.0, DUM, DUM, DUM
&
36 ,
1.0,
2.0,
3.0,
4.0,
5.0,
6.0,
7.0, DUM, DUM
&
37 ,
1.0,
2.0,
3.0,
4.0,
5.0,
6.0,
7.0,
8.0, DUM
&
38 ,
1.0,
2.0,
3.0,
4.0,
5.0,
6.0,
7.0,
8.0,
9.0 &
39 ], shape
= [
9,
9], order
= [
2,
1])
40 cholref
= reshape([
real(TKG) ::
1.0, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
41 ,
1.0,
1.0, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
42 ,
1.0,
1.0,
1.0, DUM, DUM, DUM, DUM, DUM, DUM
&
43 ,
1.0,
1.0,
1.0,
1.0, DUM, DUM, DUM, DUM, DUM
&
44 ,
1.0,
1.0,
1.0,
1.0,
1.0, DUM, DUM, DUM, DUM
&
45 ,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0, DUM, DUM, DUM
&
46 ,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0, DUM, DUM
&
47 ,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0, DUM
&
48 ,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0 &
49 ], shape
= [
9,
9], order
= [
2,
1])
52 call disp%show( cholmat ,
format = gform )
53 call disp%show(
"cholmat = getMatChol(cholmat, lowDia)")
56 call disp%show( cholmat ,
format = gform )
58 call disp%show( cholref ,
format = gform )
62 cholmat
= reshape([
real(TKG) :: DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
63 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
64 , DUM,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0 &
65 , DUM, DUM,
2.0,
2.0,
2.0,
2.0,
2.0,
2.0,
2.0,
2.0 &
66 , DUM, DUM, DUM,
3.0,
3.0,
3.0,
3.0,
3.0,
3.0,
3.0 &
67 , DUM, DUM, DUM, DUM,
4.0,
4.0,
4.0,
4.0,
4.0,
4.0 &
68 , DUM, DUM, DUM, DUM, DUM,
5.0,
5.0,
5.0,
5.0,
5.0 &
69 , DUM, DUM, DUM, DUM, DUM, DUM,
6.0,
6.0,
6.0,
6.0 &
70 , DUM, DUM, DUM, DUM, DUM, DUM, DUM,
7.0,
7.0,
7.0 &
71 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM,
8.0,
8.0 &
72 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM,
9.0 &
73 ], shape
= [
11,
10], order
= [
2,
1])
74 cholref
= reshape([
real(TKG) :: DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
75 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM
&
76 , DUM,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0 &
77 , DUM, DUM,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0 &
78 , DUM, DUM, DUM,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0 &
79 , DUM, DUM, DUM, DUM,
1.0,
1.0,
1.0,
1.0,
1.0,
1.0 &
80 , DUM, DUM, DUM, DUM, DUM,
1.0,
1.0,
1.0,
1.0,
1.0 &
81 , DUM, DUM, DUM, DUM, DUM, DUM,
1.0,
1.0,
1.0,
1.0 &
82 , DUM, DUM, DUM, DUM, DUM, DUM, DUM,
1.0,
1.0,
1.0 &
83 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM,
1.0,
1.0 &
84 , DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM, DUM,
1.0 &
85 ], shape
= [
11,
10], order
= [
2,
1])
88 call disp%show( cholmat ,
format = gform )
89 call disp%show(
"cholmat = getMatChol(cholmat(3:,2:), uppDia)")
92 call disp%show( cholmat ,
format = gform )
94 call disp%show( cholref ,
format = gform )
100 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
101 call disp%show(
"!Cholesky factorization of complex positive-definite square matrix update.")
102 call disp%show(
"!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
107 complex(TKG),
allocatable :: cholref(:,:), cholmat(:,:)
109 cholmat
= reshape( [
complex(TKG) :: CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
110 , CMPLX_DUMM, CMPLX_DUMM, (
25.0, DUM), CMPLX_DUMM, CMPLX_DUMM
&
111 , CMPLX_DUMM, CMPLX_DUMM, (
-5.0,
5.0), (
51.0, DUM), CMPLX_DUMM
&
112 , CMPLX_DUMM, CMPLX_DUMM, (
10.0,
-5.0), (
4.0,
6.0), (
71.0, DUM)
&
113 ], shape
= [
4,
5], order
= [
2,
1])
114 cholref
= reshape( [
complex(TKG) :: CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
115 , CMPLX_DUMM, CMPLX_DUMM, (
5.0,
0.0), CMPLX_DUMM, CMPLX_DUMM
&
116 , CMPLX_DUMM, CMPLX_DUMM, (
-1.0,
1.0), (
7.0,
0.0), CMPLX_DUMM
&
117 , CMPLX_DUMM, CMPLX_DUMM, (
2.0,
-1.0), (
1.0,
1.0), (
8.0,
0.0)
&
118 ], shape
= [
4,
5], order
= [
2,
1])
121 call disp%show( cholmat ,
format = cform )
122 call disp%show(
"cholmat = getMatChol(cholmat(2:,3:), lowDia)")
125 call disp%show( cholmat ,
format = cform )
127 call disp%show( cholref ,
format = cform )
130 cholmat
= reshape( [
complex(TKG) :: CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
131 , (
9.0, DUM), (
3.0,
3.0), (
3.0,
-3.0), CMPLX_DUMM
&
132 , CMPLX_DUMM, (
18.0,DUM), (
8.0,
-6.0), CMPLX_DUMM
&
133 , CMPLX_DUMM, CMPLX_DUMM, (
43.0, DUM), CMPLX_DUMM
&
134 ], shape
= [
4,
4], order
= [
2,
1])
135 cholref
= reshape( [
complex(TKG) :: CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM, CMPLX_DUMM
&
136 , (
3.0,
0.0), (
1.0,
1.0), (
1.0,
-1.0), CMPLX_DUMM
&
137 , CMPLX_DUMM, (
4.0,
0.0), (
2.0,
-1.0), CMPLX_DUMM
&
138 , CMPLX_DUMM, CMPLX_DUMM, (
6.0,
0.0), CMPLX_DUMM
&
139 ], shape
= [
4,
4], order
= [
2,
1])
142 call disp%show( cholmat ,
format = cform )
143 call disp%show(
"cholmat = getMatChol(cholmat(2:,:3), uppDia)")
146 call disp%show( cholmat ,
format = cform )
148 call disp%show( cholref ,
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 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.