Lines Matching refs:c
1 subroutine e3eig1 (rho, T, cp, gamb, c, argument
42 & c(npro), u1(npro),
173 fact = (one / sqrt( two * rho * T )) / c
175 Q(:,1,5) = fact * ((c + u) * c - eb1 * gamb)
176 Q(:,2,5) = -fact * (c + u * gamb) * Rcs(:,1,1)
177 Q(:,3,5) = -fact * (c + u * gamb) * Rcs(:,1,2)
178 Q(:,4,5) = -fact * (c + u * gamb) * Rcs(:,1,3)
181 Q(:,1,1) = fact * ((c - u) * c - eb1 * gamb)
182 Q(:,2,1) = fact * (c - u * gamb) * Rcs(:,1,1)
183 Q(:,3,1) = fact * (c - u * gamb) * Rcs(:,1,2)
184 Q(:,4,1) = fact * (c - u * gamb) * Rcs(:,1,3)
188 Q(:,2,3) = -fact * c * sqt2 * Rcs(:,2,1) ! + in original Jo/Sh code
189 Q(:,3,3) = -fact * c * sqt2 * Rcs(:,2,2) ! could be error here,
190 Q(:,4,3) = -fact * c * sqt2 * Rcs(:,2,3) ! but unlikely
194 Q(:,2,2) = fact * c * sqt2 * Rcs(:,3,1)
195 Q(:,3,2) = fact * c * sqt2 * Rcs(:,3,2)
196 Q(:,4,2) = fact * c * sqt2 * Rcs(:,3,3)
224 subroutine e3eig2 (u, c, AR1, AR2, AR3, argument
263 dimension u(npro), c(npro),
273 offd = (AR2 + AR3) * pt5 * c**2
275 rlam(:,1) = AR1 * (u + c)**2 + offd
276 rlam(:,2) = AR1 * u**2 + AR3 * c**2
277 rlam(:,3) = AR1 * u**2 + AR2 * c**2
279 rlam(:,5) = AR1 * (u - c)**2 + offd