Lines Matching refs:k
120 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
121 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); in ConvectiveFluxJacobian_Euler()
122 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
123 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
257 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
258 du[j][k] = (dU[j][k] - drho[k] * u[j]) / rho; in EulerShockTube()
259 drhodx[j] += drho[k] * dXdx[k][j]; in EulerShockTube()
260 dEdx[j] += dE[k] * dXdx[k][j]; in EulerShockTube()
262 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in EulerShockTube()
263 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in EulerShockTube()
275 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0; in EulerShockTube() local
284 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
285 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[… in EulerShockTube()
286 (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]); in EulerShockTube()
318 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
319 for (CeedInt j = 0; j < 3; j++) dv[j][k][i] -= wdetJ * nu_shock * du[k][j]; in EulerShockTube()
336 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in EulerShockTube() local
342 for (CeedInt k = 0; k < 5; k++) { in EulerShockTube() local
343 for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j]; in EulerShockTube()
359 for (CeedInt k = 0; k < 5; k++) { in EulerShockTube() local
361 stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in EulerShockTube()
366 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd… in EulerShockTube() local