Lines Matching refs:Diag
1 subroutine i3LDU (Diag, r, code) argument
44 dimension Diag(nshg,nsymdf), r(nshg,nflow)
52 Diag(:, 1) = one / sqrt (Diag(:, 1))
54 Diag(:, 2) = Diag(:, 1) * Diag(:, 2)
55 Diag(:, 3) = Diag(:, 3) - Diag(:, 2) * Diag(:, 2)
56 Diag(:, 3) = one / sqrt (Diag(:, 3))
58 Diag(:, 4) = Diag(:, 1) * Diag(:, 4)
59 Diag(:, 5) = Diag(:, 3) * (Diag(:, 5)
60 & - Diag(:, 4) * Diag(:, 2))
61 Diag(:, 6) = Diag(:, 6) - Diag(:, 4) * Diag(:, 4)
62 & - Diag(:, 5) * Diag(:, 5)
63 Diag(:, 6) = one / sqrt (Diag(:, 6))
65 Diag(:, 7) = Diag(:, 1) * Diag(:, 7)
66 Diag(:, 8) = Diag(:, 3) * (Diag(:, 8)
67 & - Diag(:, 7) * Diag(:, 2))
68 Diag(:, 9) = Diag(:, 6) * (Diag(:, 9)
69 & - Diag(:, 7) * Diag(:, 4)
70 & - Diag(:, 8) * Diag(:, 5))
72 Diag(:,10) = Diag(:,10) - Diag(:, 7) * Diag(:, 7)
73 & - Diag(:, 8) * Diag(:, 8)
74 & - Diag(:, 9) * Diag(:, 9)
75 Diag(:,10) = one / sqrt (Diag(:,10))
77 Diag(:,11) = Diag(:, 1) * Diag(:,11)
78 Diag(:,12) = Diag(:, 3) * (Diag(:,12)
79 & - Diag(:,11) * Diag(:, 2))
80 Diag(:,13) = Diag(:, 6) * (Diag(:,13)
81 & - Diag(:,11) * Diag(:, 4)
82 & - Diag(:,12) * Diag(:, 5))
83 Diag(:,14) = Diag(:,10) * (Diag(:,14)
84 & - Diag(:,11) * Diag(:, 7)
85 & - Diag(:,12) * Diag(:, 8)
86 & - Diag(:,13) * Diag(:, 9))
88 Diag(:,15) = Diag(:,15) - Diag(:,11) * Diag(:,11)
89 & - Diag(:,12) * Diag(:,12)
90 & - Diag(:,13) * Diag(:,13)
91 & - Diag(:,14) * Diag(:,14)
92 Diag(:,15) = one / sqrt (Diag(:,15))
105 r(:,1) = Diag(:, 1) * r(:,1)
106 r(:,2) = Diag(:, 3) * ( r(:,2)
107 & - r(:,1) * Diag(:, 2) )
108 r(:,3) = Diag(:, 6) * ( r(:,3)
109 & - r(:,1) * Diag(:, 4)
110 & - r(:,2) * Diag(:, 5) )
111 r(:,4) = Diag(:,10) * ( r(:,4)
112 & - r(:,1) * Diag(:, 7)
113 & - r(:,2) * Diag(:, 8)
114 & - r(:,3) * Diag(:, 9) )
115 r(:,5) = Diag(:,15) * ( r(:,5)
116 & - r(:,1) * Diag(:,11)
117 & - r(:,2) * Diag(:,12)
118 & - r(:,3) * Diag(:,13)
119 & - r(:,4) * Diag(:,14) )
132 r(:,5) = Diag(:,15) * r(:,5)
133 r(:,4) = Diag(:,10) * ( r(:,4)
134 & - r(:,5) * Diag(:,14) )
135 r(:,3) = Diag(:, 6) * ( r(:,3)
136 & - r(:,5) * Diag(:,13)
137 & - r(:,4) * Diag(:, 9) )
138 r(:,2) = Diag(:, 3) * ( r(:,2)
139 & - r(:,5) * Diag(:,12)
140 & - r(:,4) * Diag(:, 8)
141 & - r(:,3) * Diag(:, 5) )
142 r(:,1) = Diag(:, 1) * ( r(:,1)
143 & - r(:,5) * Diag(:,11)
144 & - r(:,4) * Diag(:, 7)
145 & - r(:,3) * Diag(:, 4)
146 & - r(:,2) * Diag(:, 2) )
159 r(:,1) = r(:,1) / Diag(:, 1) + r(:,2) * Diag(:, 2) +
160 & r(:,3) * Diag(:, 4) + r(:,4) * Diag(:, 7) +
161 & r(:,5) * Diag(:,11)
162 r(:,2) = r(:,2) / Diag(:, 3) + r(:,3) * Diag(:, 5) +
163 & r(:,4) * Diag(:, 8) + r(:,5) * Diag(:,12)
164 r(:,3) = r(:,3) / Diag(:, 6) + r(:,4) * Diag(:, 9) +
165 & r(:,5) * Diag(:,13)
166 r(:,4) = r(:,4) / Diag(:,10) + r(:,5) * Diag(:,14)
167 r(:,5) = r(:,5) / Diag(:,15)