| /phasta/phSolver/compressible/ |
| H A D | e3dc.f | 3 & rmi, stiff, A0DC) argument
38 & ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1)),
344 rmi( :,1) = rmi( :,1) + DC(:,intp) * gAgyi( :,1)
346 rmi( :,2) = rmi( :,2) + DC(:,intp) * gAgyi( :,2)
348 rmi( :,3) = rmi( :,3) + DC(:,intp) * gAgyi( :,3)
350 rmi( :,4) = rmi( :,4) + DC(:,intp) * gAgyi( :,4)
352 rmi( :,5) = rmi( :,5) + DC(:,intp) * gAgyi( :,5)
355 rmi( :,6) = rmi( :,6) + DC(:,intp) * gAgyi( :,6)
357 rmi( :,7) = rmi( :,7) + DC(:,intp) * gAgyi( :,7)
359 rmi( :,8) = rmi( :,8) + DC(:,intp) * gAgyi( :,8)
[all …]
|
| H A D | e3massr.f | 2 & rmi, A0) argument
29 & rmi(npro,nflow*(nsd+1)), A0(npro,nflow,nflow)
79 rmi(:,16) = rmi(:,16) + fct1*dui(:,1)
80 rmi(:,17) = rmi(:,17) + fct1*dui(:,2)
81 rmi(:,18) = rmi(:,18) + fct1*dui(:,3)
82 rmi(:,19) = rmi(:,19) + fct1*dui(:,4)
83 rmi(:,20) = rmi(:,20) + fct1*dui(:,5)
|
| H A D | e3wmlt.f | 2 & ri, rmi, rl, argument
56 & rmi(npro,nflow*(nsd+1)), stiff(npro,3*nflow,3*nflow),
100 & shg(:,i,1) * rmi(:, 1) + shg(:,i,2) * rmi(:, 6)
101 & + shg(:,i,3) * rmi(:,11)
102 & + shp(:,i) * rmi(:,16) )
104 & shg(:,i,1) * rmi(:, 2) + shg(:,i,2) * rmi(:, 7)
105 & + shg(:,i,3) * rmi(:,12)
106 & + shp(:,i) * rmi(:,17) )
108 & shg(:,i,1) * rmi(:, 3) + shg(:,i,2) * rmi(:, 8)
109 & + shg(:,i,3) * rmi(:,13)
[all …]
|
| H A D | e3ls.f | 6 & ri, rmi, rk, argument
70 & ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1)),
235 rmi(:,1) =
241 & + rmi(:,1)
242 rmi(:,2) =
248 & + rmi(:,2)
249 rmi(:,3) =
255 & + rmi(:,3)
256 rmi(:,4) =
262 & + rmi(:,4)
[all …]
|
| H A D | e3visc.f | 5 & ri, rmi, stiff,
47 & ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1))
278 rmi(:,2) = rlm2mu * g1yi(:,2)
282 rmi(:,3) = rmu * g1yi(:,3)
285 rmi(:,4) = rmu * g1yi(:,4)
288 rmi(:,5) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
295 ri (:,2:5) = ri (:,2:5) + rmi(:,2:5)
303 rmi(:, 7) = rmu * g1yi(:,3)
306 rmi(:, 8) = rlm * g1yi(:,2)
310 rmi(:, 9) = rmu * g2yi(:,4)
[all …]
|
| H A D | e3source.f | 1 subroutine e3source (ri, rmi, rlyi, argument
28 dimension ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1)),
213 rmi (:,16) = rmi (:,16) - src(:,1)
214 rmi (:,17) = rmi (:,17) - src(:,2)
215 rmi (:,18) = rmi (:,18) - src(:,3)
216 rmi (:,19) = rmi (:,19) - src(:,4)
217 rmi (:,20) = rmi (:,20) - src(:,5)
|
| H A D | e3.f | 85 dimension ri(npro,nflow*(nsd+1)), rmi(npro,nflow*(nsd+1)),
122 rmi = zero
168 & ri, rmi, EGmass,
181 & ri, rmi, stiff,
189 call e3source (ri, rmi, rlyi,
206 & ri, rmi, rk,
221 & ri, rmi, stiff, A0DC)
248 call e3massr (aci, dui, ri, rmi, A0)
293 & ri, rmi, rl,
|
| H A D | e3conv.f | 6 & ri, rmi, EGmass, argument
58 & rmi(npro,nflow*(nsd+1)), EGmass(npro,nedof,nedof),
184 & rmi(:,16:20) = rLyi ! modified residual uses non i.b.p form of conv.
|