Lines Matching refs:rmu
3 & rmu, rlm, rlm2mu,
50 dimension rmu(npro), rlm(npro),
70 stiff(:, 3, 3) = rmu
71 stiff(:, 4, 4) = rmu
73 stiff(:, 5, 3) = rmu * u2
74 stiff(:, 5, 4) = rmu * u3
80 stiff(:, 3, 7) = rmu
81 stiff(:, 5, 7) = rmu * u2
87 stiff(:, 4,12) = rmu
88 stiff(:, 5,12) = rmu * u3
94 stiff(:, 7, 3) = rmu
97 stiff(:,10, 3) = rmu * u1
102 stiff(:, 7, 7) = rmu
104 stiff(:, 9, 9) = rmu
105 stiff(:,10, 7) = rmu * u1
107 stiff(:,10, 9) = rmu * u3
113 stiff(:, 9,13) = rmu
114 stiff(:,10,13) = rmu * u3
119 stiff(:,12, 4) = rmu
122 stiff(:,15, 4) = rmu * u1
126 stiff(:,13, 9) = rmu
129 stiff(:,15, 9) = rmu * u2
133 stiff(:,12,12) = rmu
134 stiff(:,13,13) = rmu
136 stiff(:,15,12) = rmu * u1
137 stiff(:,15,13) = rmu * u2
230 compK(:, 1) = f1 * T * rlm2mu + f3 * T * rmu
231 & + f6 * T * rmu
233 compK(:, 2) = f2 * T * (rlm + rmu)
234 compK(:, 3) = f1 * T * rmu + f3 * T * rlm2mu
235 & + f6 * T * rmu
237 compK(:, 4) = f4 * T * (rlm + rmu)
238 compK(:, 5) = f5 * T * (rlm + rmu)
239 compK(:, 6) = f1 * T * rmu + f3 * T * rmu
242 compK(:, 7) = f1 * T * rlm2mu * u1 + f2 * T * (rlm + rmu) * u2
243 & + f3 * T * rmu * u1 + f4 * T * (rlm + rmu) * u3
244 & + f6 * T * rmu * u1
245 compK(:, 8) = f1 * T * rmu * u2 + f2 * T * (rlm + rmu) * u1
246 & + f3 * T * rlm2mu * u2 + f5 * T * (rlm + rmu) * u3
247 & + f6 * T * rmu * u2
248 compK(:, 9) = f1 * T * rmu * u3 + f3 * T * rmu * u3
249 & + f4 * T * (rlm + rmu) * u1 + f5 * T * (rlm + rmu) * u2
254 compK(:,10) = f1 * T * (con * T + two * rmu * rk + (rlm +
255 & rmu) * u1**2) + f2 * T * (rlm + rmu) * two * u1 * u2
256 & + f3 * T * (con * T + two * rmu * rk + (rlm + rmu) *
257 & u2**2) + f4 * T * (rlm + rmu) * two * u1 * u3
258 & + f5 * T * (rlm + rmu) * two * u2 * u3 + f6 * T * (con
259 & * T + two * rmu * rk + (rlm + rmu) * u3**2)
282 rmi(:,3) = rmu * g1yi(:,3)
283 & + rmu * g2yi(:,2)
285 rmi(:,4) = rmu * g1yi(:,4)
286 & + rmu * g3yi(:,2)
288 rmi(:,5) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
289 & + rmu * u3 * g1yi(:,4)
290 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3)
291 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4)
303 rmi(:, 7) = rmu * g1yi(:,3)
304 & + rmu * g2yi(:,2)
310 rmi(:, 9) = rmu * g2yi(:,4)
311 & + rmu * g3yi(:,3)
313 rmi(:,10) = rlm * u2 * g1yi(:,2) + rmu * u1 * g1yi(:,3)
314 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3)
315 & + rmu * u3 * g2yi(:,4)
316 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4)
327 rmi(:,12) = rmu * g1yi(:,4)
328 & + rmu * g3yi(:,2)
330 rmi(:,13) = rmu * g2yi(:,4)
331 & + rmu * g3yi(:,3)
337 rmi(:,15) = rlm * u3 * g1yi(:,2) + rmu * u1 * g1yi(:,4)
338 & + rlm * u3 * g2yi(:,3) + rmu * u2 * g2yi(:,4)
339 & + rmu * u1 * g3yi(:,2) + rmu * u2 * g3yi(:,3)
362 & rmu, Sclr, rho,
395 & g3yti(npro), rmu(npro),
414 stifft(:,1,1)=rmu
416 stifft(:,1,1)=saSigmaInv*((rmu/rho)+max(zero,Sclr))