Lines Matching refs:r
1 subroutine i3LU (Diag, r, code) argument
34 dimension Diag(nshg,nflow,nflow), r(nshg,nflow)
105 r(:,2) = r(:,2) - Diag(:,2,1) * r(:,1)
106 r(:,3) = r(:,3) - Diag(:,3,1) * r(:,1)
107 & - Diag(:,3,2) * r(:,2)
108 r(:,4) = r(:,4) - Diag(:,4,1) * r(:,1)
109 & - Diag(:,4,2) * r(:,2)
110 & - Diag(:,4,3) * r(:,3)
111 r(:,5) = r(:,5) - Diag(:,5,1) * r(:,1)
112 & - Diag(:,5,2) * r(:,2)
113 & - Diag(:,5,3) * r(:,3)
114 & - Diag(:,5,4) * r(:,4)
127 r(:,5) = Diag(:,5,5) * r(:,5)
128 r(:,4) = Diag(:,4,4) * ( r(:,4)
129 & - r(:,5) * Diag(:,4,5) )
130 r(:,3) = Diag(:,3,3) * ( r(:,3)
131 & - r(:,5) * Diag(:,3,5)
132 & - r(:,4) * Diag(:,3,4) )
133 r(:,2) = Diag(:,2,2) * ( r(:,2)
134 & - r(:,5) * Diag(:,2,5)
135 & - r(:,4) * Diag(:,2,4)
136 & - r(:,3) * Diag(:,2,3) )
137 r(:,1) = Diag(:,1,1) * ( r(:,1)
138 & - r(:,5) * Diag(:,1,5)
139 & - r(:,4) * Diag(:,1,4)
140 & - r(:,3) * Diag(:,1,3)
141 & - r(:,2) * Diag(:,1,2) )
154 r(:,1) = r(:,1) / Diag(:,1,1) + r(:,2) * Diag(:,1,2) +
155 & r(:,3) * Diag(:,1,3) + r(:,4) * Diag(:,1,4) +
156 & r(:,5) * Diag(:,1,5)
157 r(:,2) = r(:,2) / Diag(:,2,2) +
158 & r(:,3) * Diag(:,2,3) + r(:,4) * Diag(:,2,4) +
159 & r(:,5) * Diag(:,2,5)
160 r(:,3) = r(:,3) / Diag(:,3,3) +
161 & r(:,4) * Diag(:,3,4) +
162 & r(:,5) * Diag(:,3,5)
163 r(:,4) = r(:,4) / Diag(:,4,4) +
164 & r(:,5) * Diag(:,4,5)
165 r(:,5) = r(:,5) / Diag(:,5,5)