Home
last modified time | relevance | path

Searched refs:i (Results 1 – 25 of 74) sorted by relevance

123

/honee/qfunctions/
H A Dutils_eigensolver_jacobi.h30 CEED_QFUNCTION_HELPER CeedInt MaxEntryRow(const CeedScalar *A, CeedInt N, CeedInt i) { in MaxEntryRow() argument
31 CeedInt j_max = i + 1; in MaxEntryRow()
32 for (CeedInt j = i + 2; j < N; j++) in MaxEntryRow()
33 if (fabs(A[i * N + j]) > fabs(A[i * N + j_max])) j_max = j; in MaxEntryRow()
48 for (CeedInt i = 1; i < N - 1; i++) { in MaxEntry() local
49 CeedInt j = max_idx_row[i]; in MaxEntry()
50 if (fabs(A[i * N + j]) > max_entry) { in MaxEntry()
51 max_entry = fabs(A[i * N + j]); in MaxEntry()
52 *i_max = i; in MaxEntry()
67 CEED_QFUNCTION_HELPER void CalcRot(const CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedScalar… in CalcRot() argument
[all …]
H A Dsgs_dd_training.h35 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ComputeSGS_DDAnisotropicTrainingDataNodal() local
36 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ComputeSGS_DDAnisotropicTrainingDataNodal()
38 {grad_velo[0][0][i], grad_velo[0][1][i], grad_velo[0][2][i]}, in ComputeSGS_DDAnisotropicTrainingDataNodal()
39 {grad_velo[1][0][i], grad_velo[1][1][i], grad_velo[1][2][i]}, in ComputeSGS_DDAnisotropicTrainingDataNodal()
40 {grad_velo[2][0][i], grad_velo[2][1][i], grad_velo[2][2][i]} in ComputeSGS_DDAnisotropicTrainingDataNodal()
42 …ar km_A_ij[6] = {A_ij_delta[0][i], A_ij_delta[1][i], A_ij_delta[2][i], A_ij_delta[3][i], A_ij_delt… in ComputeSGS_DDAnisotropicTrainingDataNodal()
43 const CeedScalar delta = A_ij_delta[6][i]; in ComputeSGS_DDAnisotropicTrainingDataNodal()
50 for (int j = 0; j < 6; j++) v[j][i] = inv_multiplicity[i] * inputs[j]; in ComputeSGS_DDAnisotropicTrainingDataNodal()
52 …v[0 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XX][i] - Square(s.Y.velocity[0])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
53 …v[1 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YY][i] - Square(s.Y.velocity[1])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
[all …]
H A Dsgs_dd_model.h34 for (CeedInt i = 0; i < N; i++) x[i] *= (x[i] < 0 ? alpha : 1.); in LeakyReLU() local
78 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ComputeSgsDDNodal_Fused() local
79 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ComputeSgsDDNodal_Fused()
81 {grad_velo[0][0][i], grad_velo[0][1][i], grad_velo[0][2][i]}, in ComputeSgsDDNodal_Fused()
82 {grad_velo[1][0][i], grad_velo[1][1][i], grad_velo[1][2][i]}, in ComputeSgsDDNodal_Fused()
83 {grad_velo[2][0][i], grad_velo[2][1][i], grad_velo[2][2][i]} in ComputeSgsDDNodal_Fused()
85 …ar km_A_ij[6] = {A_ij_delta[0][i], A_ij_delta[1][i], A_ij_delta[2][i], A_ij_delta[3][i], A_ij_delt… in ComputeSgsDDNodal_Fused()
86 const CeedScalar delta = A_ij_delta[6][i]; in ComputeSgsDDNodal_Fused()
92 for (int j = 0; j < 6; j++) v[j][i] = inv_multiplicity[i] * km_sgs[j]; in ComputeSgsDDNodal_Fused()
122 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ComputeSgsDDNodal_Sequential_Inputs() local
[all …]
H A Dutils.h28 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] *= alpha; in ScaleN() local
33 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) u[i] = alpha; in SetValueN() local
37 …*x, CeedScalar *y, const CeedInt N) { CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) y[i] = x[i]; } in CopyN() local
45 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) output += u[i] * v[i]; in DotN() local
51 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) y[i] += alpha * x[i]; in AXPY() local
85 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) b[i] += DotN(&A[i * M], x, M); in MatVecNM() local
88 …CeedPragmaSIMD for (CeedInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) b[… in MatVecNM() local
109 …CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < M; j++) B[… in MatDiagNM() local
112 …dPragmaSIMD for (CeedInt i = 0; i < M; i++) { CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) B[i *… in MatDiagNM() local
129 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) { in MatMatN() local
[all …]
H A Dgrid_anisotropy_tensor.h21 for (int i = 0; i < 3; i++) evals[i] = 1 / sqrt(evals[i]); in AnisotropyTensor() local
38 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in AnisotropyTensorProjection() local
39 const CeedScalar wdetJ = q_data[0][i]; in AnisotropyTensorProjection()
41 {q_data[1][i], q_data[2][i], q_data[3][i]}, in AnisotropyTensorProjection()
42 {q_data[4][i], q_data[5][i], q_data[6][i]}, in AnisotropyTensorProjection()
43 {q_data[7][i], q_data[8][i], q_data[9][i]} in AnisotropyTensorProjection()
51 for (CeedInt j = 0; j < 6; j++) v[j][i] = wdetJ * km_A_ij[j]; in AnisotropyTensorProjection()
52 v[6][i] = wdetJ * delta; in AnisotropyTensorProjection()
62 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in AnisotropyTensorCollocate() local
64 {q_data[1][i], q_data[2][i], q_data[3][i]}, in AnisotropyTensorCollocate()
[all …]
H A Dnewtonian_state.h26 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U() local
32 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y() local
38 for (int i = 0; i < 3; i++) V[i + 1] = s.S_momentum[i]; in UnpackState_V() local
66 for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density; in StatePrimitiveFromConservative() local
77 for (CeedInt i = 0; i < 3; i++) { in StatePrimitiveFromConservative_fwd() local
78 dY.velocity[i] = (dU.momentum[i] - s.Y.velocity[i] * dU.density) / s.U.density; in StatePrimitiveFromConservative_fwd()
100 for (int i = 0; i < 3; i++) V.S_momentum[i] = rho_div_p * Y.velocity[i]; in StateEntropyFromPrimitive() local
118 …for (CeedInt i = 0; i < 3; i++) dV.S_momentum[i] = rho_div_p * dY.velocity[i] + drho_div_p * s.Y.v… in StateEntropyFromPrimitive_fwd() local
125 for (int i = 0; i < 3; i++) Y.velocity[i] = -V.S_momentum[i] / V.S_energy; in StatePrimitiveFromEntropy() local
138 …for (int i = 0; i < 3; i++) dY.velocity[i] = -(dV.S_momentum[i] - V.S_momentum[i] * dV.S_energy / … in StatePrimitiveFromEntropy_fwd() local
[all …]
H A Deulervortex.h179 for (CeedInt i = 0; i < 3; i++) { // Jacobian matrices for 3 directions 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()
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()
187 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
189 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
190 dF[i][4][4] = u[i] * gamma; in ConvectiveFluxJacobian_Euler()
211 for (CeedInt i = 0; i < 3; i++) { in Tau_spatial() local
213 CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i])); in Tau_spatial()
[all …]
H A Dshocktube.h112 for (CeedInt i = 0; i < 3; i++) { // Jacobian matrices for 3 directions 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()
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()
120 dF[i][j + 1][4] = ((i == j) ? (gamma - 1.) : 0.); in ConvectiveFluxJacobian_Euler()
122 dF[i][4][0] = u[i] * ((gamma - 1.) * u_sq - E * gamma / rho); in ConvectiveFluxJacobian_Euler()
123 dF[i][4][4] = u[i] * gamma; in ConvectiveFluxJacobian_Euler()
158 for (CeedInt i = 0; i < 3; i++) { in Tau_spatial() local
160 CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i])); in Tau_spatial()
[all …]
H A Dsetupgeo2d.h45 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in Setup2d() local
47 InvertMappingJacobian_2D(Q, i, J, dXdx, &detJ); in Setup2d()
48 const CeedScalar wdetJ = w[i] * detJ; in Setup2d()
50 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data); in Setup2d()
51 StoredValuesPack(Q, i, 1, 4, (const CeedScalar *)dXdx, q_data); in Setup2d()
86 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupBoundary2d() local
88 NormalVectorFromdxdX_2D(Q, i, J, normal, &detJb); in SetupBoundary2d()
89 const CeedScalar wdetJ = w[i] * detJb; in SetupBoundary2d()
91 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur); in SetupBoundary2d()
92 StoredValuesPack(Q, i, 1, 2, normal, q_data_sur); in SetupBoundary2d()
[all …]
H A Dadvection.h182 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsAdvection() local
183 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsAdvection()
187 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection()
200 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsAdvection2d() local
201 const CeedScalar x[] = {X[0][i], X[1][i]}; in ICsAdvection2d()
205 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsAdvection2d()
210 CEED_QFUNCTION_HELPER void StatePhysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, N… in StatePhysicalGradientFromReference_ND() argument
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()
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()
[all …]
H A Dsetupgeo.h47 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in Setup() local
49 InvertMappingJacobian_3D(Q, i, J, dXdx, &detJ); in Setup()
50 const CeedScalar wdetJ = w[i] * detJ; in Setup()
52 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data); in Setup()
53 StoredValuesPack(Q, i, 1, 9, (const CeedScalar *)dXdx, q_data); in Setup()
105 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupBoundary() local
108 NormalVectorFromdxdX_3D(Q, i, J, normal, &detJb); in SetupBoundary()
109 InvertBoundaryMappingJacobian_3D(Q, i, J, dXdx); in SetupBoundary()
110 const CeedScalar wdetJ = w[i] * detJb; in SetupBoundary()
112 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur); in SetupBoundary()
[all …]
H A Ddifferential_filter.h38 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in DifferentialFilter_RHS() local
39 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in DifferentialFilter_RHS()
40 const CeedScalar wdetJ = q_data[0][i]; in DifferentialFilter_RHS()
43 v0[DIFF_FILTER_PRESSURE][i] = wdetJ * s.Y.pressure; in DifferentialFilter_RHS()
44 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
45 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
46 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
47 v0[DIFF_FILTER_TEMPERATURE][i] = wdetJ * s.Y.temperature; in DifferentialFilter_RHS()
48 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
[all …]
H A Dblasius.h42 for (int i = 0; i < 4; i++) { in ChebyshevEval() local
44 f[i] = table[i][0] * Tf[0] + table[i][1] * Tf[1] + table[i][2] * Tf[2]; in ChebyshevEval()
46 for (int i = 3; i < N; i++) { in ChebyshevEval() local
48 table[0][i % 3] = 2 * x * table[0][(i - 1) % 3] - table[0][(i - 2) % 3]; in ChebyshevEval()
52 table[j][i % 3] = i * (2 * table[j - 1][(i - 1) % 3] + table[j][(i - 2) % 3] / (i - 2)); in ChebyshevEval()
55 f[j] += table[j][i % 3] * Tf[i]; in ChebyshevEval()
58 for (int i = 1; i < 4; i++) { in ChebyshevEval() local
60 for (int j = 0; j < i; j++) f[i] *= dX_deta; in ChebyshevEval()
112 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsBlasius() local
113 const CeedScalar x[3] = {X[0][i], X[1][i], X[2][i]}; in ICsBlasius()
[all …]
H A Dnewtonian.h22 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsNewtonianIG() local
26 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsNewtonianIG()
53 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in MassFunction_Newtonian() local
54 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MassFunction_Newtonian()
55 const CeedScalar qi_dot[5] = {q_dot[0][i], q_dot[1][i], q_dot[2][i], q_dot[3][i], q_dot[4][i]}; in MassFunction_Newtonian()
59 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in MassFunction_Newtonian()
63 v[f][i] = wdetJ * qi_dot[f]; in MassFunction_Newtonian()
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()
95 for (int i = 0; i < 5; i++) damp_residual[i] += U[i]; in InternalDampingLayer_Residual() local
203 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in RHSFunction_Newtonian() local
[all …]
H A Dchannel.h65 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ICsChannel() local
66 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsChannel()
70 for (CeedInt j = 0; j < 5; j++) q0[j][i] = q[j]; in ICsChannel()
89 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in Channel_Inflow() local
91 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, norm); in Channel_Inflow()
99 const CeedScalar x[3] = {0, X[1][i], X[2][i]}; in Channel_Inflow()
108 for (CeedInt j = 0; j < 5; j++) q_inside[j] = q[j][i]; in Channel_Inflow()
120 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in Channel_Inflow()
126 v[0][i] -= wdetJb * rho_in * u_normal; in Channel_Inflow()
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()
[all …]
H A Dstg_shur14.h46 for (CeedInt i = 0; i < nprofs; i++) { in InterpolateProfile() local
47 if (wall_dist < prof_wd[i]) { in InterpolateProfile()
48 idx = i; in InterpolateProfile()
165 for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i]; in StgShur14Calc() local
171 for (CeedInt i = 0; i < 3; i++) vp[i] *= 2 * sqrt(1.5); in StgShur14Calc() local
210 for (CeedInt i = 0; i < 3; i++) xdotd += d[i * nmodes + n] * xhat[i]; in StgShur14Calc_PrecompEktot() local
218 for (CeedInt i = 0; i < 3; i++) vp[i] *= 2 * sqrt(1.5); in StgShur14Calc_PrecompEktot() local
261 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in StgShur14Preprocess() local
262 const CeedScalar wall_dist = x[1][i]; in StgShur14Preprocess()
264 StoredValuesUnpack(Q, i, 0, 9, dXdx_q, (CeedScalar *)dXdx); in StgShur14Preprocess()
[all …]
H A Dbc_outflow.h48 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in RiemannOutflow() local
50 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in RiemannOutflow()
52 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in RiemannOutflow()
66 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_q, dXdx, grad_s); in RiemannOutflow()
79 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in RiemannOutflow()
83 StoredValuesPack(Q, i, 0, 5, qi, jac_data_sur); in RiemannOutflow()
84 StoredValuesPack(Q, i, 5, 6, kmstress, jac_data_sur); in RiemannOutflow()
117 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in RiemannOutflow_Jacobian() local
119 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in RiemannOutflow_Jacobian()
123 StoredValuesUnpack(Q, i, 0, 5, jac_data_sur, qi); in RiemannOutflow_Jacobian()
[all …]
H A Dmonitor_totalkineticenergy.h16 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in MonitorTotalKineticEnergy() local
17 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MonitorTotalKineticEnergy()
22 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in MonitorTotalKineticEnergy()
23 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in MonitorTotalKineticEnergy()
25 v[0][i] = wdetJ * 0.5 * s.U.density * Dot3(s.Y.velocity, s.Y.velocity); in MonitorTotalKineticEnergy()
30 v[1][i] = wdetJ * -2 * gas.mu * DotN((CeedScalar *)strain_rate, (CeedScalar *)strain_rate, 9); in MonitorTotalKineticEnergy()
31 v[2][i] = wdetJ * -gas.lambda * gas.mu * Square(div_u); in MonitorTotalKineticEnergy()
32 v[3][i] = wdetJ * s.Y.pressure * div_u; in MonitorTotalKineticEnergy()
35 v[4][i] = wdetJ * gas.mu * Dot3(vorticity, vorticity); in MonitorTotalKineticEnergy()
/honee/qfunctions/spanstats/
H A Dturbulence.h50 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ChildStatsCollection() local
51 const CeedScalar wdetJ = q_data[0][i] * delta_t; in ChildStatsCollection()
53 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ChildStatsCollection()
56 v[TURB_MEAN_DENSITY][i] = wdetJ * s.U.density; in ChildStatsCollection()
57 v[TURB_MEAN_PRESSURE][i] = wdetJ * s.Y.pressure; in ChildStatsCollection()
58 v[TURB_MEAN_PRESSURE_SQUARED][i] = wdetJ * Square(s.Y.pressure); in ChildStatsCollection()
59 v[TURB_MEAN_PRESSURE_VELOCITY_X][i] = wdetJ * s.Y.pressure * s.Y.velocity[0]; in ChildStatsCollection()
60 v[TURB_MEAN_PRESSURE_VELOCITY_Y][i] = wdetJ * s.Y.pressure * s.Y.velocity[1]; in ChildStatsCollection()
61 v[TURB_MEAN_PRESSURE_VELOCITY_Z][i] = wdetJ * s.Y.pressure * s.Y.velocity[2]; in ChildStatsCollection()
62 v[TURB_MEAN_DENSITY_TEMPERATURE][i] = wdetJ * s.U.density * s.Y.temperature; in ChildStatsCollection()
[all …]
H A Dcflpe.h28 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ChildStatsCollection_CflPe() local
29 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in ChildStatsCollection_CflPe()
37 QdataUnpack_2D(Q, i, q_data, &wdetJ, dXdx); in ChildStatsCollection_CflPe()
50 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in ChildStatsCollection_CflPe()
62 v[0][i] = wdetJ * cfl; in ChildStatsCollection_CflPe()
63 v[1][i] = wdetJ * Square(cfl); in ChildStatsCollection_CflPe()
64 v[2][i] = wdetJ * Cube(cfl); in ChildStatsCollection_CflPe()
65 v[3][i] = wdetJ * Pe; in ChildStatsCollection_CflPe()
66 v[4][i] = wdetJ * Square(Pe); in ChildStatsCollection_CflPe()
67 v[5][i] = wdetJ * Cube(Pe); in ChildStatsCollection_CflPe()
/honee/problems/
H A Dstg_shur14.c33 for (PetscInt i = 0; i < nprofs; i++) { in CalcCholeskyDecomp() local
34 Cij[0][i] = sqrt(Rij[0][i]); in CalcCholeskyDecomp()
35 Cij[3][i] = Rij[3][i] / Cij[0][i]; in CalcCholeskyDecomp()
36 Cij[1][i] = sqrt(Rij[1][i] - Square(Cij[3][i])); in CalcCholeskyDecomp()
37 Cij[4][i] = Rij[4][i] / Cij[0][i]; in CalcCholeskyDecomp()
38 Cij[5][i] = (Rij[5][i] - Cij[3][i] * Cij[4][i]) / Cij[1][i]; in CalcCholeskyDecomp()
39 Cij[2][i] = sqrt(Rij[2][i] - Square(Cij[4][i]) - Square(Cij[5][i])); in CalcCholeskyDecomp()
41 PetscCheck(!isnan(Cij[0][i]) && !isnan(Cij[1][i]) && !isnan(Cij[2][i]), comm, PETSC_ERR_FP, in CalcCholeskyDecomp()
42 … at profile point %" PetscInt_FMT ". Either STGInflow has non-SPD matrix or contains nan.", i + 1); in CalcCholeskyDecomp()
76 for (PetscInt i = 0; i < stg_ctx->nprofs; i++) { in ReadStgInflow() local
[all …]
H A Dbc_freestream.c126 …for (int i = 0; i < 3; i++) Y_inf.velocity[i] = reference->velocity[i] * units->second / units->me… in FreestreamBCSetup() local
138 for (int i = 0; i < 3; i++) Y_inf.velocity[i] *= units->meter / units->second; in FreestreamBCSetup() local
198 for (int i = 1; i < 4; i++) { in CheckQWithTolerance() local
199 relative_error[i] = RelativeError(u_magnitude, Q_a[i], Q_b[i], divisor_threshold); in CheckQWithTolerance()
205 for (int i = 1; i < 4; i++) { in CheckQWithTolerance() local
206 if (fabs(relative_error[i]) >= rtol_u) { in CheckQWithTolerance()
207 …rintf("%s[%d] error %g (expected %.10e, got %.10e)\n", name, i, relative_error[i], Q_s[i], Q_a[i]); in CheckQWithTolerance()
234 for (int i = 0; i < 10; i++) { in TestRiemannHLL_fwd() local
240 if (i < 5) { in TestRiemannHLL_fwd()
241 dY_left[i] = Y0_left[i]; in TestRiemannHLL_fwd()
[all …]
H A Ddensitycurrent.c39 for (PetscInt i = 0; i < dim; i++) domain_size[i] = domain_max[i] - domain_min[i]; in NS_DENSITY_CURRENT() local
50 for (PetscInt i = 0; i < dim; i++) center[i] = .5 * domain_size[i]; in NS_DENSITY_CURRENT() local
61 for (PetscInt i = 0; i < dim; i++) dc_axis[i] /= norm; in NS_DENSITY_CURRENT() local
74 for (PetscInt i = 0; i < dim; i++) center[i] *= units->meter; in NS_DENSITY_CURRENT() local
/honee/src/
H A Dhonee-meshtransform.c25 for (PetscInt i = 0; i < ndims; i++) dims[i] = atoi(array[i]); in GetYNodeLocs() local
32 for (PetscInt i = 0; i < dims[0]; i++) { in GetYNodeLocs() local
38 …FMT " of %s does not contain correct number of columns (%d instead of %d)", i, path, ndims, dims[1… in GetYNodeLocs()
40 node_locs[i] = (PetscReal)atof(array[0]); in GetYNodeLocs()
71 for (PetscInt i = 0; i < 3; i++) domain_size[i] = domain_max[i] - domain_min[i]; in HoneeMeshTransform_PlateMesh() local
98 for (PetscInt i = 0; i < ncoords; i++) { in HoneeMeshTransform_PlateMesh() local
99 PetscInt y_box_index = round(coords[i][1] / dybox); in HoneeMeshTransform_PlateMesh()
101 coords[i][1] = in HoneeMeshTransform_PlateMesh()
102 …(1 - (coords[i][0] - domain_min[0]) * angle_coeff / domain_max[1]) * dy1 * (pow(growth, coords[i][… in HoneeMeshTransform_PlateMesh()
105 …coords[i][1] = (1 - (coords[i][0] - domain_min[0]) * angle_coeff / domain_max[1]) * exp(log(refine… in HoneeMeshTransform_PlateMesh()
[all …]
H A Dsetupdm.c94 for (PetscInt i = 0; i < problem->num_bc_defs; i++) { in SetUpDM() local
95 BCDefinition bc_def = problem->bc_defs[i]; in SetUpDM()
100 PetscCall(BCDefinitionSetDM(problem->bc_defs[i], dm)); in SetUpDM()
116 for (PetscInt i = 0; i < problem->num_components; i++) { in SetUpDM() local
117 PetscCall(PetscSectionSetComponentName(section, 0, i, problem->component_names[i])); in SetUpDM()
131 for (PetscInt i = 0, d = honee->app_ctx->degree; i < honee->app_ctx->viz_refine; i++) { in VizRefineDM() local
133 PetscCall(DMRefine(dm_hierarchy[i], MPI_COMM_NULL, &dm_hierarchy[i + 1])); in VizRefineDM()
134 PetscCall(DMClearDS(dm_hierarchy[i + 1])); in VizRefineDM()
135 PetscCall(DMClearFields(dm_hierarchy[i + 1])); in VizRefineDM()
136 PetscCall(DMSetCoarseDM(dm_hierarchy[i + 1], dm_hierarchy[i])); in VizRefineDM()
[all …]

123