Lines Matching refs:k
182 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
183 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); 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()
185 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
295 for (CeedInt k = 0; k < 3; k++) { 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()
312 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in Euler() local
321 for (CeedInt k = 0; k < 3; k++) { 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()
340 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in Euler() local
346 for (CeedInt k = 0; k < 5; k++) { 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()
364 for (CeedInt k = 0; k < 5; k++) { 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()
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() local
416 for (CeedInt k = 0; k < 3; k++) { 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()
432 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in IFunction_Euler() local
443 for (CeedInt k = 0; k < 3; k++) { 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()
462 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in IFunction_Euler() local
468 for (CeedInt k = 0; k < 5; k++) { 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()
490 for (CeedInt k = 0; k < 5; k++) { 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()
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() local
501 for (CeedInt k = 0; k < 5; k++) { 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()
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() local