Lines Matching refs:u
87 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
98 u[0] = mean_velocity[0] - C * y0; in Exact_Euler()
99 u[1] = mean_velocity[1] + C * x0; 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()
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()
122 u[0] = mean_velocity[0]; in Exact_Euler()
123 u[1] = mean_velocity[1]; 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()
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()
151 u[0] = mean_velocity[0]; in Exact_Euler()
152 u[1] = mean_velocity[1]; 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()
165 u[0] = mean_velocity[0]; in Exact_Euler()
166 u[1] = mean_velocity[1]; 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 …cobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u[3], const CeedScalar… in ConvectiveFluxJacobian_Euler()
184 CeedScalar u_sq = u[0] * u[0] + u[1] * u[1] + u[2] * u[2]; // Velocity square in ConvectiveFluxJacobian_Euler()
187 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; in ConvectiveFluxJacobian_Euler()
190 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
191 … = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.) * u[i] * u[k]; in ConvectiveFluxJacobian_Euler()
195 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
196 dF[i][4][4] = u[i] * gamma; in ConvectiveFluxJacobian_Euler()
215 …_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const CeedScalar u[3], const CeedScalar… in Tau_spatial()
221 CeedScalar fastest_wave = fabs(u[i]) + sound_speed; in Tau_spatial()
284 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler() local
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 … 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()
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()
334 …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()
339 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in Euler()
361 Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau); in Euler()
407 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in IFunction_Euler() local
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 … 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()
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()
458 …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()
463 ConvectiveFluxJacobian_Euler(jacob_F_conv, rho, u, E, gamma); in IFunction_Euler()
489 Tau_spatial(Tau_x, dXdx, u, sound_speed, c_tau); in IFunction_Euler()
604 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Euler_Outflow() local
619 const CeedScalar E_kinetic = (u[0] * u[0] + u[1] * u[1]) / 2.; in Euler_Outflow()
621 … const CeedScalar u_normal = norm[0] * u[0] + norm[1] * u[1] + norm[2] * u[2]; // Normal velocity 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()