Lines Matching full:k
115 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
116 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); in ConvectiveFluxJacobian_Euler()
117 …dF[i][j + 1][k + 1] = ((j == k) ? u[i] : 0.) + ((i == k) ? u[j] : 0.) - ((i == j) ? u[k] : 0.) * (… in ConvectiveFluxJacobian_Euler()
118 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
245 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
246 du[j][k] = (dU[j][k] - drho[k] * u[j]) / rho; in EulerShockTube()
247 drhodx[j] += drho[k] * dXdx[k][j]; in EulerShockTube()
248 dEdx[j] += dE[k] * dXdx[k][j]; in EulerShockTube()
250 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in EulerShockTube()
251 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in EulerShockTube()
263 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0; in EulerShockTube() local
272 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
273 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[… in EulerShockTube()
274 (rho * u[j] * u[2] + (j == 2 ? P : 0)) * dXdx[k][2]); in EulerShockTube()
306 for (CeedInt k = 0; k < 3; k++) { in EulerShockTube() local
307 for (CeedInt j = 0; j < 3; j++) dv[j][k][i] -= wdetJ * nu_shock * du[k][j]; in EulerShockTube()
324 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in EulerShockTube() local
330 for (CeedInt k = 0; k < 5; k++) { in EulerShockTube() local
331 for (CeedInt l = 0; l < 5; l++) strong_conv[k] += jacob_F_conv[j][k][l] * dqdx[l][j]; in EulerShockTube()
347 for (CeedInt k = 0; k < 5; k++) { in EulerShockTube() local
349 stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in EulerShockTube()
354 …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