Lines Matching refs:k
188 for (CeedInt k = 0; k < 3; k++) { // Columns of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
189 dF[i][0][k + 1] = ((i == k) ? 1. : 0.); 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()
191 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
301 for (CeedInt k = 0; k < 3; k++) { 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()
318 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in Euler() local
327 for (CeedInt k = 0; k < 3; k++) { 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()
346 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in Euler() local
352 for (CeedInt k = 0; k < 5; k++) { 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()
370 for (CeedInt k = 0; k < 5; k++) { 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()
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() local
424 for (CeedInt k = 0; k < 3; k++) { 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()
440 for (CeedInt k = 0; k < 3; k++) dv[k][j][i] = 0.; in IFunction_Euler() local
451 for (CeedInt k = 0; k < 3; k++) { 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()
470 for (CeedInt k = 0; k < 3; k++) dqdx[k + 1][j] = dUdx[k][j]; in IFunction_Euler() local
476 for (CeedInt k = 0; k < 5; k++) { 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()
498 for (CeedInt k = 0; k < 5; k++) { 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()
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() local
509 for (CeedInt k = 0; k < 5; k++) { 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()
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() local