Home
last modified time | relevance | path

Searched refs:j (Results 1 – 25 of 41) sorted by relevance

12

/honee/qfunctions/
H A Deulervortex.h180 for (CeedInt j = 0; j < 3; j++) { // Rows of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
181 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; 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()
187 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
234 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsEuler() local
294 for (CeedInt j = 0; j < 3; j++) { 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()
[all …]
H A Dshocktube.h113 for (CeedInt j = 0; j < 3; j++) { // Rows of each Jacobian matrix in ConvectiveFluxJacobian_Euler() local
114 dF[i][j + 1][0] = ((i == j) ? ((gamma - 1.) * (u_sq / 2.)) : 0.) - u[i] * u[j]; 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()
120 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
180 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsShockTube() local
244 for (CeedInt j = 0; j < 3; j++) { 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()
[all …]
H A Dutils_eigensolver_jacobi.h32 for (CeedInt j = i + 2; j < N; j++) in MaxEntryRow() local
33 if (fabs(A[i * N + j]) > fabs(A[i * N + j_max])) j_max = j; in MaxEntryRow()
49 CeedInt j = max_idx_row[i]; in MaxEntry() local
50 if (fabs(A[i * N + j]) > max_entry) { in MaxEntry()
51 max_entry = fabs(A[i * N + j]); in MaxEntry()
53 *j_max = j; in MaxEntry()
67 CEED_QFUNCTION_HELPER void CalcRot(const CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedScalar… in CalcRot() argument
69 CeedScalar A_jj_ii = (A[j * N + j] - A[i * N + i]); in CalcRot()
74 CeedScalar A_ij = A[i * N + j]; in CalcRot()
155 CEED_QFUNCTION_HELPER void ApplyRot(CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedInt *max_id… in ApplyRot() argument
[all …]
H A Dbc_slip.h30 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip() local
38 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in Slip() local
74 for (int j = 0; j < 5; j++) dqi[j] = dq[j][i]; in Slip_Jacobian() local
80 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip_Jacobian() local
86 … for (CeedInt j = 0; j < 3; j++) dvel_reflect[j] = ds.Y.velocity[j] - 2. * normal[j] * dvel_normal; in Slip_Jacobian() local
94 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in Slip_Jacobian() local
H A Dnewtonian.h26 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsNewtonianIG() local
74 for (CeedInt j = 0; j < 5; j++) { in MassFunction_Newtonian() local
76 …Grad_v[k][j][i] = wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] * dXdx[k… in MassFunction_Newtonian()
206 for (int j = 0; j < 5; j++) U[j] = q[j][i]; in RHSFunction_Newtonian() local
226 for (CeedInt j = 0; j < 5; j++) { in RHSFunction_Newtonian() local
227 …nt k = 0; k < 3; k++) Grad_v[k][j][i] = wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * Flux[j][1]… in RHSFunction_Newtonian()
231 for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j]; in RHSFunction_Newtonian() local
238 for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j]; in RHSFunction_Newtonian() local
243 for (int j = 1; j < 5; j++) divFdiff_i[j] = divFdiff[j - 1][i]; in RHSFunction_Newtonian() local
250 for (CeedInt j = 0; j < 5; j++) { in RHSFunction_Newtonian() local
[all …]
H A Dchannel.h70 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsChannel() local
108 for (CeedInt j = 0; j < 5; j++) q_inside[j] = q[j][i]; in Channel_Inflow() local
120 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in Channel_Inflow() local
129 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_in * u_normal * s_exact.Y.velocity[j]… in Channel_Inflow() local
159 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in Channel_Outflow() local
169 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Channel_Outflow() local
H A Dadvection.h187 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection() local
205 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection2d() local
217 for (CeedInt j = 0; j < 5; j++) { in StatePhysicalGradientFromReference_ND() local
218 …dqi[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXd… in StatePhysicalGradientFromReference_ND()
229 for (CeedInt j = 0; j < 5; j++) { in StatePhysicalGradientFromReference_ND() local
230 …dqi[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0 * N + k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXd… in StatePhysicalGradientFromReference_ND()
231 grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k]; in StatePhysicalGradientFromReference_ND()
292 for (CeedInt j = 0; j < dim; j++) grad_v[j][f][i] = 0; // No Change in density or momentum in IFunction_AdvectionGeneric() local
297 for (CeedInt j = 0; j < dim; j++) { in IFunction_AdvectionGeneric() local
299 div_u += grad_s[k].Y.velocity[j]; in IFunction_AdvectionGeneric()
[all …]
H A Dutils.h88 …dInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[i] += A[j * M + i] * x[… in MatVecNM() local
109 … i = 0; i < N; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < M; j++) B[i * M + j] += D[i] * A[i * … in MatDiagNM() local
112 …dInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) B[i * N + j] += D[i] * A[… in MatDiagNM() local
130 CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) { in MatMatN() local
131 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[k * N + j]; in MatMatN()
137 CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) { in MatMatN() local
138 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[i * N + k] * B[j * N + k]; in MatMatN()
148 CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) { in MatMatN() local
149 … CeedPragmaSIMD for (CeedInt k = 0; k < N; k++) C[i * N + j] += A[k * N + i] * B[k * N + j]; in MatMatN()
155 CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) { in MatMatN() local
[all …]
H A Dblasius.h50 for (int j = 1; j < 4; j++) { in ChebyshevEval() local
52 table[j][i % 3] = i * (2 * table[j - 1][(i - 1) % 3] + table[j][(i - 2) % 3] / (i - 2)); in ChebyshevEval()
54 for (int j = 0; j < 4; j++) { in ChebyshevEval() local
55 f[j] += table[j][i % 3] * Tf[i]; in ChebyshevEval()
60 for (int j = 0; j < i; j++) f[i] *= dX_deta; in ChebyshevEval() local
118 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsBlasius() local
148 for (CeedInt j = 0; j < 5; j++) qi[j] = q[j][i]; in Blasius_Inflow() local
171 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in Blasius_Inflow() local
222 for (int j = 0; j < 3; j++) { in Blasius_Inflow_Jacobian() local
223 v[j + 1][i] = -wdetJb * (drho * u_normal * s.Y.velocity[j] + normal[j] * dP); in Blasius_Inflow_Jacobian()
H A Dbc_outflow.h60 for (CeedInt j = 0; j < 3; j++) { in RiemannOutflow() local
61 …y_ext.velocity[j] += normal[j] * proj; // (I - n n^T) projects into the plane tangent to the norm… in RiemannOutflow()
79 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in RiemannOutflow() local
125 for (int j = 0; j < 5; j++) dqi[j] = dq[j][i]; in RiemannOutflow_Jacobian() local
138 for (CeedInt j = 0; j < 3; j++) { in RiemannOutflow_Jacobian() local
139 y_ext.velocity[j] += normal[j] * proj; in RiemannOutflow_Jacobian()
140 dy_ext.velocity[j] += normal[j] * dproj; in RiemannOutflow_Jacobian()
161 for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in RiemannOutflow_Jacobian() local
220 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in PressureOutflow() local
266 for (int j = 0; j < 5; j++) dqi[j] = dq[j][i]; in PressureOutflow_Jacobian() local
[all …]
H A Dstg_shur14.h235 …for (CeedInt j = 0; j < 3; j++) lengths[j] = 2 / sqrt(Square(dXdx[0][j]) + Square(dXdx[1][j]) + Sq… in CalculateElementLengths() local
313 for (CeedInt j = 0; j < 3; j++) u[j] = ubar[j]; in ICsStg() local
319 for (CeedInt j = 0; j < 5; j++) { in ICsStg() local
320 q0[j][i] = q[j]; in ICsStg()
364 … for (CeedInt j = 1; j < 3; j++) h_node_sep[j] = 2 / sqrt(Square(dXdx[0][j]) + Square(dXdx[1][j])); in StgShur14Inflow() local
372 for (CeedInt j = 0; j < 3; j++) u[j] = ubar[j]; in StgShur14Inflow() local
394 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in StgShur14Inflow() local
401 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P); in StgShur14Inflow() local
439 for (CeedInt j = 0; j < 3; j++) velocity[j] = jac_data_sur[5 + j][i]; in StgShur14Inflow_Jacobian() local
461 …for (int j = 0; j < 3; j++) v[j + 1][i] = -wdetJb * (drho * u_normal * velocity[j] + normal[j] * d… in StgShur14Inflow_Jacobian() local
[all …]
H A Dnewtonian_state.h267 for (int j = 0; j < 3; j++) R.momentum[j] += a[i] * X[i].momentum[j]; in StateConservativeMult() local
428 …for (CeedInt j = 0; j < 3; j++) Flux[i].momentum[j] = s.U.momentum[i] * s.Y.velocity[j] + s.Y.pres… in FluxInviscid() local
436 for (CeedInt j = 0; j < 3; j++) { in FluxInviscid_fwd() local
437 …dFlux[i].momentum[j] = ds.U.momentum[i] * s.Y.velocity[j] + s.U.momentum[i] * ds.Y.velocity[j] + d… in FluxInviscid_fwd()
448 for (CeedInt j = 0; j < 3; j++) Flux_dot_n.momentum[j] += Flux[i].momentum[j] * normal[i]; in FluxInviscidDotNormal() local
460 for (CeedInt j = 0; j < 3; j++) Flux_dot_n.momentum[j] += dFlux[i].momentum[j] * normal[i]; in FluxInviscidDotNormal_fwd() local
473 for (CeedInt j = 0; j < 5; j++) strong_conv[j] += dF_i[j]; in FluxInviscidStrong() local
479 for (CeedInt j = 0; j < 3; j++) { in FluxTotal() local
480 Flux[0][j] = F_inviscid[j].density; in FluxTotal()
481 for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = F_inviscid[j].momentum[k] - stress[k][j]; in FluxTotal()
[all …]
H A Dstabilization.h36 for (CeedInt j = 0; j < 5; j++) stab[j][i] += dF_i[j]; in StabilizationMatrix() local
53 for (CeedInt j = 0; j < 5; j++) strong_residual[j] += U_dot[j] - body_force[j] + divFdiff[j]; in Stabilization() local
H A Dsetupgeo_helpers.h148 for (CeedInt j = 0; j < 2; j++) { in InvertBoundaryMappingJacobian_3D() local
150 for (CeedInt l = 0; l < 3; l++) dxdXTdxdX[j][k] += dxdX[l][j] * dxdX[l][k]; in InvertBoundaryMappingJacobian_3D()
164 for (CeedInt j = 0; j < 2; j++) { in InvertBoundaryMappingJacobian_3D() local
166 dXdx[j][k] = 0; in InvertBoundaryMappingJacobian_3D()
167 for (CeedInt l = 0; l < 2; l++) dXdx[j][k] += dxdXTdxdX_inv[l][j] * dxdX[k][l]; in InvertBoundaryMappingJacobian_3D()
H A Dsgs_dd_model.h92 for (int j = 0; j < 6; j++) v[j][i] = inv_multiplicity[i] * km_sgs[j]; in ComputeSgsDDNodal_Fused() local
197 for (int j = 0; j < 6; j++) kmsgs_stress[j][i] = inv_multiplicity[i] * kmsgs_stress_i[j]; in ComputeSgsDDNodal_Sequential_Outputs() local
207 for (CeedInt j = 0; j < 3; j++) { in FluxSubgridStress() local
208 Flux[0][j] = 0.; in FluxSubgridStress()
209 for (CeedInt k = 0; k < 3; k++) Flux[k + 1][j] = sgs[k][j]; in FluxSubgridStress()
210 Flux[4][j] = Y.velocity[0] * sgs[0][j] + Y.velocity[1] * sgs[1][j] + Y.velocity[2] * sgs[2][j]; in FluxSubgridStress()
235 for (CeedInt j = 0; j < 5; j++) { in IFunction_NodalSgs() local
236 …Grad_v[k][j][i] = -wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * Flux[j][1] + dXdx[k][2] * Flux[… in IFunction_NodalSgs()
H A Dgrid_anisotropy_tensor.h51 for (CeedInt j = 0; j < 6; j++) v[j][i] = wdetJ * km_A_ij[j]; in AnisotropyTensorProjection() local
74 for (CeedInt j = 0; j < 6; j++) v[j][i] = km_A_ij[j]; in AnisotropyTensorCollocate() local
H A Dbc_freestream.h44 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in Freestream() local
95 for (int j = 0; j < 5; j++) dqi[j] = dq[j][i]; in Freestream_Jacobian() local
110 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in Freestream_Jacobian() local
H A Dstrong_boundary_conditions.h16 for (int j = 0; j < 3; j++) coords_stored[j][i] = coords[j][i]; in SetupStrongBC() local
H A Ddifferential_filter.h86 CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) { in DifferentialFilter_LHS_N() local
126 v[j][i] = wdetJ * q[j][i]; in DifferentialFilter_LHS_N()
129 …dq[k] = Grad_q[0 * N + j][i] * dXdx[0][k] + Grad_q[1 * N + j][i] * dXdx[1][k] + Grad_q[2 * N + j][… in DifferentialFilter_LHS_N()
134 Grad_v[k * N + j][i] = wdetJ * dq_dXdx_a[k]; in DifferentialFilter_LHS_N()
214 for (CeedInt j = 0; j < 5; j++) q0[j][i] = phi; in DifferentialFilter_MMS_IC() local
H A Dvelocity_gradient_projection.h30 for (CeedInt j = 0; j < 3; j++) { in VelocityGradientProjectionRHS() local
32 v[j * 3 + k][i] = wdetJ * grad_velocity[j][k]; in VelocityGradientProjectionRHS()
H A Dtaylorgreen.h48 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsTaylorGreen() local
H A Dmass.h24 …eedInt i = 0; i < Q; i++) CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) v[j][i] = q_data[i] * u[j… in Mass_N() local
H A Dgaussianwave.h55 for (CeedInt j = 0; j < 5; j++) q0[j][i] = qi[j]; in IC_GaussianWave() local
H A Ddensitycurrent.h138 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsDC() local
H A Dsgs_dd_training.h50 for (int j = 0; j < 6; j++) v[j][i] = inv_multiplicity[i] * inputs[j]; in ComputeSGS_DDAnisotropicTrainingDataNodal() local

12