Lines Matching refs:rho
87 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
96 rho = pow(T / S_vortex, 1 / (gamma - 1.)); in Exact_Euler()
97 P = rho * T; in Exact_Euler()
102 q[0] = rho; in Exact_Euler()
103 q[1] = rho * u[0]; in Exact_Euler()
104 q[2] = rho * u[1]; in Exact_Euler()
105 q[3] = rho * u[2]; in Exact_Euler()
106 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
109 rho = 1.; in Exact_Euler()
113 q[0] = rho; in Exact_Euler()
114 q[1] = rho * u[0]; in Exact_Euler()
115 q[2] = rho * u[1]; in Exact_Euler()
116 q[3] = rho * u[2]; in Exact_Euler()
120 rho = 1.; in Exact_Euler()
126 q[0] = rho; in Exact_Euler()
127 q[1] = rho * u[0]; in Exact_Euler()
128 q[2] = rho * u[1]; in Exact_Euler()
129 q[3] = rho * u[2]; in Exact_Euler()
137 rho = P / (R * T); in Exact_Euler()
140 q[0] = rho; in Exact_Euler()
141 q[1] = rho * u[0]; in Exact_Euler()
142 q[2] = rho * u[1]; in Exact_Euler()
143 q[3] = rho * u[2]; in Exact_Euler()
144 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler()
150 rho = P / (R * T); in Exact_Euler()
155 q[0] = rho; in Exact_Euler()
156 q[1] = rho * u[0]; in Exact_Euler()
157 q[2] = rho * u[1]; in Exact_Euler()
158 q[3] = rho * u[2]; in Exact_Euler()
159 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler()
164 rho = P / (R * T); in Exact_Euler()
169 q[0] = rho; in Exact_Euler()
170 q[1] = rho * u[0]; in Exact_Euler()
171 q[2] = rho * u[1]; in Exact_Euler()
172 q[3] = rho * u[2]; in Exact_Euler()
173 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler()
182 …void ConvectiveFluxJacobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u… in ConvectiveFluxJacobian_Euler() argument
191 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
195 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
283 const CeedScalar rho = q[0][i]; in Euler() local
284 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler()
311 …const CeedScalar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E… in Euler()
323 …edInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1]… in Euler()
328 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0.)) * dXdx[k][0] + (rho * u[j] * u… in Euler()
329 (rho * u[j] * u[2] + (j == 2 ? P : 0.)) * dXdx[k][2]); in Euler()
339 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in Euler()
359 const CeedScalar sound_speed = sqrt(gamma * P / rho); in Euler()
406 const CeedScalar rho = q[0][i]; in IFunction_Euler() local
407 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in IFunction_Euler()
433 …const CeedScalar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E… in IFunction_Euler()
447 …edInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1]… in IFunction_Euler()
452 …dv[k][j + 1][i] -= wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0.)) * dXdx[k][0] + (rho * u[j] * u… in IFunction_Euler()
453 (rho * u[j] * u[2] + (j == 2 ? P : 0.)) * dXdx[k][2]); in IFunction_Euler()
463 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in IFunction_Euler()
487 const CeedScalar sound_speed = sqrt(gamma * P / rho); in IFunction_Euler()
603 const CeedScalar rho = q[0][i]; in Euler_Outflow() local
604 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler_Outflow()
620 const CeedScalar P = (E - E_kinetic * rho) * (gamma - 1.); // pressure in Euler_Outflow()
624 v[0][i] -= wdetJb * rho * u_normal; in Euler_Outflow()
627 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Euler_Outflow()