Lines Matching refs:j
186 for (CeedInt j = 0; j < 3; j++) { // Rows of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
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()
193 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
240 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsEuler() local
300 for (CeedInt j = 0; j < 3; j++) { in Euler() local
302 drhodx[j] += drho[k] * dXdx[k][j]; in Euler()
303 dEdx[j] += dE[k] * dXdx[k][j]; in Euler()
305 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in Euler()
306 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in Euler()
316 for (CeedInt j = 0; j < 5; j++) { in Euler() local
317 v[j][i] = 0.; in Euler()
318 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in Euler()
323 …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
326 for (CeedInt j = 0; j < 3; j++) { in Euler() local
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 …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
343 for (CeedInt j = 0; j < 3; j++) { in Euler() local
344 dqdx[0][j] = drhodx[j]; in Euler()
345 dqdx[4][j] = dEdx[j]; in Euler()
346 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in Euler()
351 for (CeedInt j = 0; j < 3; j++) { in Euler() local
353 for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j]; in Euler()
369 for (CeedInt j = 0; j < 3; j++) { in Euler() local
371 … for (CeedInt l = 0; l < 5; l++) stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in Euler()
375 for (CeedInt j = 0; j < 5; j++) { in Euler() local
376 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd… in Euler()
423 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
425 drhodx[j] += drho[k] * dXdx[k][j]; in IFunction_Euler()
426 dEdx[j] += dE[k] * dXdx[k][j]; in IFunction_Euler()
428 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in IFunction_Euler()
429 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in IFunction_Euler()
438 for (CeedInt j = 0; j < 5; j++) { in IFunction_Euler() local
439 v[j][i] = 0.; in IFunction_Euler()
440 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in IFunction_Euler()
443 for (CeedInt j = 0; j < 5; j++) v[j][i] += wdetJ * q_dot[j][i]; in IFunction_Euler() local
447 …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
450 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
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 …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
467 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
468 dqdx[0][j] = drhodx[j]; in IFunction_Euler()
469 dqdx[4][j] = dEdx[j]; in IFunction_Euler()
470 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in IFunction_Euler()
475 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
477 for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j]; in IFunction_Euler()
483 for (CeedInt j = 0; j < 5; j++) strong_res[j] = q_dot[j][i] + strong_conv[j]; in IFunction_Euler() local
497 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
499 … 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()
503 for (CeedInt j = 0; j < 5; j++) { in IFunction_Euler() local
504 …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()
508 for (CeedInt j = 0; j < 3; j++) { in IFunction_Euler() local
510 … 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()
514 for (CeedInt j = 0; j < 5; j++) { in IFunction_Euler() local
515 …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()
559 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in TravelingVortex_Inflow() local
574 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_inlet * face_normal * mean_velocity[j… in TravelingVortex_Inflow() local
615 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0; in Euler_Outflow() local
627 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Euler_Outflow() local