Lines Matching refs:u
81 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
92 u[0] = mean_velocity[0] - C * y0; in Exact_Euler()
93 u[1] = mean_velocity[1] + C * x0; 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()
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()
116 u[0] = mean_velocity[0]; in Exact_Euler()
117 u[1] = mean_velocity[1]; 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()
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()
145 u[0] = mean_velocity[0]; in Exact_Euler()
146 u[1] = mean_velocity[1]; 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()
159 u[0] = mean_velocity[0]; in Exact_Euler()
160 u[1] = mean_velocity[1]; 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 …cobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar… in ConvectiveFluxJacobian_Euler()
178 CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2]; // Velocity square in ConvectiveFluxJacobian_Euler()
181 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; in ConvectiveFluxJacobian_Euler()
184 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
185 … = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.) * u[i] * u[k]; in ConvectiveFluxJacobian_Euler()
189 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
190 dF[i][4][4] = u[i] * gamma; in ConvectiveFluxJacobian_Euler()
209 …_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const CeedScalar u[3], const CeedScalar… in Tau_spatial()
215 CeedScalar fastest_wave = fabs(u[i]) + sound_speed; in Tau_spatial()
278 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler() local
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 … j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rh… 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()
328 …nt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[j][1] + u[… in Euler()
333 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in Euler()
355 Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau); in Euler()
399 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in IFunction_Euler() local
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 … j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rh… 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()
450 …nt j = 0; j < 3; j++) dv[j][4][i] -= wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[j][1] + u[… in IFunction_Euler()
455 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in IFunction_Euler()
481 Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau); in IFunction_Euler()
595 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler_Outflow() local
610 const CeedScalar E_kinetic = (u[0] * u[0] + u[1] * u[1]) / 2.; in Euler_Outflow()
612 const CeedScalar u_normal = Dot3(normal, u); // Normal velocity 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()