Lines Matching refs:rho
81 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
90 rho = pow(T / S_vortex, 1 / (gamma - 1.)); in Exact_Euler()
91 P = rho * T; in Exact_Euler()
96 q[0] = rho; in Exact_Euler()
97 q[1] = rho * u[0]; in Exact_Euler()
98 q[2] = rho * u[1]; in Exact_Euler()
99 q[3] = rho * u[2]; in Exact_Euler()
100 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
103 rho = 1.; in Exact_Euler()
107 q[0] = rho; in Exact_Euler()
108 q[1] = rho * u[0]; in Exact_Euler()
109 q[2] = rho * u[1]; in Exact_Euler()
110 q[3] = rho * u[2]; in Exact_Euler()
114 rho = 1.; in Exact_Euler()
120 q[0] = rho; in Exact_Euler()
121 q[1] = rho * u[0]; in Exact_Euler()
122 q[2] = rho * u[1]; in Exact_Euler()
123 q[3] = rho * u[2]; in Exact_Euler()
131 rho = P / (R * T); in Exact_Euler()
134 q[0] = rho; in Exact_Euler()
135 q[1] = rho * u[0]; in Exact_Euler()
136 q[2] = rho * u[1]; in Exact_Euler()
137 q[3] = rho * u[2]; in Exact_Euler()
138 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler()
144 rho = P / (R * T); in Exact_Euler()
149 q[0] = rho; in Exact_Euler()
150 q[1] = rho * u[0]; in Exact_Euler()
151 q[2] = rho * u[1]; in Exact_Euler()
152 q[3] = rho * u[2]; in Exact_Euler()
153 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler()
158 rho = P / (R * T); in Exact_Euler()
163 q[0] = rho; in Exact_Euler()
164 q[1] = rho * u[0]; in Exact_Euler()
165 q[2] = rho * u[1]; in Exact_Euler()
166 q[3] = rho * u[2]; in Exact_Euler()
167 q[4] = rho * (cv * T + (u[0] * u[0] + u[1] * u[1]) / 2.); in Exact_Euler()
176 …void ConvectiveFluxJacobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u… in ConvectiveFluxJacobian_Euler() argument
185 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
189 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
277 const CeedScalar rho = q[0][i]; in Euler() local
278 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler()
305 …const CeedScalar E_kinetic = 0.5 * rho * (u[0] * u[0] + u[1] * u[1] + u[2] * u[2]), E_internal = E… in Euler()
317 …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()
322 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0.)) * dXdx[k][0] + (rho * u[j] * u… in Euler()
323 (rho * u[j] * u[2] + (j == 2 ? P : 0.)) * dXdx[k][2]); in Euler()
333 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in Euler()
353 const CeedScalar sound_speed = sqrt(gamma * P / rho); in Euler()
398 const CeedScalar rho = q[0][i]; in IFunction_Euler() local
399 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in IFunction_Euler()
425 …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()
439 …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()
444 …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()
445 (rho * u[j] * u[2] + (j == 2 ? P : 0.)) * dXdx[k][2]); in IFunction_Euler()
455 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in IFunction_Euler()
479 const CeedScalar sound_speed = sqrt(gamma * P / rho); in IFunction_Euler()
594 const CeedScalar rho = q[0][i]; in Euler_Outflow() local
595 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler_Outflow()
611 const CeedScalar P = (E - E_kinetic * rho) * (gamma - 1.); // pressure in Euler_Outflow()
615 v[0][i] -= wdetJb * rho * u_normal; in Euler_Outflow()
618 … for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P); in Euler_Outflow()