Lines Matching refs:j
180 for (CeedInt j = 0; j < 3; j++) { // Rows of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
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()
187 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
234 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsEuler() local
294 for (CeedInt j = 0; j < 3; j++) { in Euler() local
296 drhodx[j] += drho[k] * dXdx[k][j]; in Euler()
297 dEdx[j] += dE[k] * dXdx[k][j]; in Euler()
299 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in Euler()
300 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in Euler()
310 for (CeedInt j = 0; j < 5; j++) { in Euler() local
311 v[j][i] = 0.; in Euler()
312 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in Euler()
317 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd… in Euler() local
320 for (CeedInt j = 0; j < 3; j++) { in Euler() local
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 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[… in Euler() local
337 for (CeedInt j = 0; j < 3; j++) { in Euler() local
338 dqdx[0][j] = drhodx[j]; in Euler()
339 dqdx[4][j] = dEdx[j]; in Euler()
340 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in Euler()
345 for (CeedInt j = 0; j < 3; j++) { in Euler() local
347 for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j]; in Euler()
363 for (CeedInt j = 0; j < 3; j++) { in Euler() local
365 … for (CeedInt l = 0; l < 5; l++) stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in Euler()
369 for (CeedInt j = 0; j < 5; j++) { in Euler() local
370 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd… in Euler()
415 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
417 drhodx[j] += drho[k] * dXdx[k][j]; in IFunction_Euler()
418 dEdx[j] += dE[k] * dXdx[k][j]; in IFunction_Euler()
420 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in IFunction_Euler()
421 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in IFunction_Euler()
430 for (CeedInt j = 0; j < 5; j++) { in IFunction_Euler() local
431 v[j][i] = 0.; in IFunction_Euler()
432 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in IFunction_Euler()
435 for (CeedInt j = 0; j < 5; j++) v[j][i] += wdetJ * q_dot[j][i]; in IFunction_Euler() local
439 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd… in IFunction_Euler() local
442 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
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 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] -= wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[… in IFunction_Euler() local
459 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
460 dqdx[0][j] = drhodx[j]; in IFunction_Euler()
461 dqdx[4][j] = dEdx[j]; in IFunction_Euler()
462 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in IFunction_Euler()
467 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
469 for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j]; in IFunction_Euler()
475 for (CeedInt j = 0; j < 5; j++) strong_res[j] = q_dot[j][i] + strong_conv[j]; in IFunction_Euler() local
489 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
491 … for (CeedInt l = 0; l < 5; l++) stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in IFunction_Euler()
495 for (CeedInt j = 0; j < 5; j++) { in IFunction_Euler() local
496 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd… in IFunction_Euler()
500 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
502 … for (CeedInt l = 0; l < 5; l++) stab[k][j] = jacob_F_conv[j][k][l] * Tau_x[j] * strong_res[l]; in IFunction_Euler()
506 for (CeedInt j = 0; j < 5; j++) { in IFunction_Euler() local
507 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd… in IFunction_Euler()
550 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in TravelingVortex_Inflow() local
565 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_inlet * face_normal * mean_velocity[j… in TravelingVortex_Inflow() local
606 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0; in Euler_Outflow() local
618 … for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P); in Euler_Outflow() local