Home
last modified time | relevance | path

Searched refs:dXdx (Results 1 – 22 of 22) sorted by relevance

/honee/qfunctions/
H A Dshocktube.h134 …N_HELPER CeedScalar Covariant_length_along_vector(CeedScalar vec[3], const CeedScalar dXdx[3][3]) { in Covariant_length_along_vector()
137 MatVec3(dXdx, vec, CEED_TRANSPOSE, vec_dot_jacobian); in Covariant_length_along_vector()
156 CEED_QFUNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const Ceed… in Tau_spatial()
160 CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i])); in Tau_spatial()
236 CeedScalar wdetJ, dXdx[3][3]; in EulerShockTube() local
237 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); 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()
[all …]
H A Deulervortex.h209 CEED_QFUNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const Ceed… in Tau_spatial()
213 CeedScalar h = 2 / sqrt(Square(dXdx[0][i]) + Square(dXdx[1][i]) + Square(dXdx[2][i])); in Tau_spatial()
287 CeedScalar wdetJ, dXdx[3][3]; in Euler() local
288 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in Euler()
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()
317 … < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rho * u[2] … in Euler()
322 …J * ((rho * u[j] * u[0] + (j == 0 ? P : 0.)) * dXdx[k][0] + (rho * u[j] * u[1] + (j == 1 ? P : 0.)… in Euler()
[all …]
H A Dsetupgeo.h48 CeedScalar detJ, dXdx[3][3]; in Setup() local
49 InvertMappingJacobian_3D(Q, i, J, dXdx, &detJ); in Setup()
53 StoredValuesPack(Q, i, 1, 9, (const CeedScalar *)dXdx, q_data); in Setup()
106 CeedScalar detJb, normal[3], dXdx[2][3]; in SetupBoundary() local
109 InvertBoundaryMappingJacobian_3D(Q, i, J, dXdx); in SetupBoundary()
114 StoredValuesPack(Q, i, 4, 6, (const CeedScalar *)dXdx, q_data_sur); in SetupBoundary()
147 CeedScalar detJ_face, normal[3], dXdx[3][3]; in SetupBoundaryGradient() local
151 InvertMappingJacobian_3D(Q, i, J_cell, dXdx, NULL); in SetupBoundaryGradient()
154 StoredValuesPack(Q, i, 1, 9, (CeedScalar *)dXdx, q_data_sur); in SetupBoundaryGradient()
H A Dutils.h261 CEED_QFUNCTION_HELPER void KMMetricTensor(const CeedScalar dXdx[3][3], CeedScalar km_g_ij[6]) { in KMMetricTensor()
263 MatMat3(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, g_ij); in KMMetricTensor()
347 …CeedInt dim, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar *dXdx) { in QdataUnpack_ND() argument
351 if (dXdx) StoredValuesUnpack(Q, i, 1, 4, q_data, dXdx); in QdataUnpack_ND()
355 if (dXdx) StoredValuesUnpack(Q, i, 1, 9, q_data, dXdx); in QdataUnpack_ND()
374 …D(CeedInt dim, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar *dXdx, in QdataBoundaryUnpack_ND() argument
384 if (dXdx) StoredValuesUnpack(Q, i, 4, 6, q_data, (CeedScalar *)dXdx); in QdataBoundaryUnpack_ND()
404 CeedScalar *dXdx, CeedScalar *normal) { in QdataBoundaryGradientUnpack_ND() argument
408 if (dXdx) StoredValuesUnpack(Q, i, 1, 4, q_data, dXdx); in QdataBoundaryGradientUnpack_ND()
413 if (dXdx) StoredValuesUnpack(Q, i, 1, 9, q_data, dXdx); in QdataBoundaryGradientUnpack_ND()
[all …]
H A Dsetupgeo2d.h46 CeedScalar dXdx[2][2], detJ; in Setup2d() local
47 InvertMappingJacobian_2D(Q, i, J, dXdx, &detJ); in Setup2d()
51 StoredValuesPack(Q, i, 1, 4, (const CeedScalar *)dXdx, q_data); in Setup2d()
141 CeedScalar detJb, normal[3], dXdx[2][3]; in Setup2D_3Dcoords() local
144 InvertBoundaryMappingJacobian_3D(Q, i, J, dXdx); in Setup2D_3Dcoords()
148 StoredValuesPack(Q, i, 1, 6, (const CeedScalar *)dXdx, q_data_sur); in Setup2D_3Dcoords()
177 CeedScalar detJ_face, normal[2], dXdx[2][2]; in Setup2DBoundaryGradient() local
181 InvertMappingJacobian_2D(Q, i, J_cell, dXdx, NULL); in Setup2DBoundaryGradient()
184 StoredValuesPack(Q, i, 1, 4, (CeedScalar *)dXdx, q_data_sur); in Setup2DBoundaryGradient()
H A Dadvection.h211 … StateVariable state_var, const CeedScalar *grad_q, const CeedScalar *dXdx, 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()
240 …LPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXdx, CeedInt dim) { in Tau() argument
245 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
251 MatMatN(dXdx, dXdx, dim, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in Tau()
284 CeedScalar wdetJ, dXdx[9]; in IFunction_AdvectionGeneric() local
285 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in IFunction_AdvectionGeneric()
287 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in IFunction_AdvectionGeneric()
[all …]
H A Dnewtonian.h58 CeedScalar wdetJ, dXdx[3][3]; in MassFunction_Newtonian() local
59 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in MassFunction_Newtonian()
70 Tau_diagPrim(context->tau_coeffs, gas, s, dXdx, context->dt, Tau_d); 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()
204 CeedScalar U[5], wdetJ, dXdx[3][3]; in RHSFunction_Newtonian() local
207 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in RHSFunction_Newtonian()
211 StatePhysicalGradientFromReference(Q, i, gas, s, STATEVAR_CONSERVATIVE, Grad_q, dXdx, grad_s); in RHSFunction_Newtonian()
227 …t 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()
247 Tau_diagPrim(context->tau_coeffs, gas, s, dXdx, dt, Tau_d); in RHSFunction_Newtonian()
251 …; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] … in RHSFunction_Newtonian()
[all …]
H A Dmonitor_cfl.h23 CeedScalar wdetJ, dXdx[2][2], gijd_mat[2][2] = {{0.}}; in MonitorCFL() local
25 QdataUnpack_2D(Q, i, q_data, &wdetJ, dXdx); in MonitorCFL()
26 MatMat2(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in MonitorCFL()
33 CeedScalar wdetJ, dXdx[3][3], gijd_mat[3][3] = {{0.}}; in MonitorCFL() local
35 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in MonitorCFL()
36 MatMat3(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in MonitorCFL()
H A Dsetupgeo_helpers.h25 …Jacobian_3D(CeedInt Q, CeedInt i, const CeedScalar (*dxdX_q)[3][CEED_Q_VLA], CeedScalar dXdx[3][3], in InvertMappingJacobian_3D()
30 MatInv3(dxdX, dXdx, detJ_ptr); in InvertMappingJacobian_3D()
47 …Jacobian_2D(CeedInt Q, CeedInt i, const CeedScalar (*dxdX_q)[2][CEED_Q_VLA], CeedScalar dXdx[2][2], in InvertMappingJacobian_2D()
52 MatInv2(dxdX, dXdx, detJ_ptr); in InvertMappingJacobian_2D()
142 …cobian_3D(CeedInt Q, CeedInt i, const CeedScalar (*dxdX_q)[3][CEED_Q_VLA], CeedScalar dXdx[2][3]) { in InvertBoundaryMappingJacobian_3D()
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 Dbc_outflow.h49 CeedScalar wdetJb, dXdx[2][3], normal[3]; in RiemannOutflow() local
50 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in RiemannOutflow()
66 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_q, dXdx, grad_s); in RiemannOutflow()
118 CeedScalar wdetJb, dXdx[2][3], normal[3]; in RiemannOutflow_Jacobian() local
119 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in RiemannOutflow_Jacobian()
147 … StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_dq, dXdx, grad_ds); in RiemannOutflow_Jacobian()
201 CeedScalar wdetJb, dXdx[2][3], normal[3]; in PressureOutflow() local
202 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in PressureOutflow()
206 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in PressureOutflow()
259 CeedScalar wdetJb, dXdx[2][3], normal[3]; in PressureOutflow_Jacobian() local
[all …]
H A Dstrong_boundary_conditions.h18 CeedScalar dXdx[3][3]; in SetupStrongBC() local
19 InvertMappingJacobian_3D(Q, i, dxdX_q, dXdx, NULL); in SetupStrongBC()
21 StoredValuesPack(Q, i, 0, 9, (CeedScalar *)dXdx, dXdx_q); in SetupStrongBC()
H A Dstg_shur14.h234 CEED_QFUNCTION_HELPER void CalculateElementLengths(CeedScalar dXdx[3][3], CeedScalar scale, CeedSca… in CalculateElementLengths()
235 …dInt j = 0; j < 3; j++) lengths[j] = 2 / sqrt(Square(dXdx[0][j]) + Square(dXdx[1][j]) + Square(dXd… in CalculateElementLengths()
263 CeedScalar dXdx[3][3], h_node_sep[3]; in StgShur14Preprocess() local
264 StoredValuesUnpack(Q, i, 0, 9, dXdx_q, (CeedScalar *)dXdx); in StgShur14Preprocess()
266 CalculateElementLengths(dXdx, stg_ctx->h_scale_factor, h_node_sep); in StgShur14Preprocess()
301 CeedScalar dXdx[3][3]; in ICsStg() local
304 InvertMappingJacobian_3D(Q, i, J, dXdx, NULL); in ICsStg()
305 CalculateElementLengths(dXdx, stg_ctx->h_scale_factor, h_node_sep); in ICsStg()
358 CeedScalar wdetJb, dXdx[2][3], normal[3]; in StgShur14Inflow() local
359 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in StgShur14Inflow()
[all …]
H A Dgrid_anisotropy_tensor.h40 const CeedScalar dXdx[3][3] = { in AnisotropyTensorProjection() local
47 KMMetricTensor(dXdx, km_g_ij); in AnisotropyTensorProjection()
63 const CeedScalar dXdx[3][3] = { in AnisotropyTensorCollocate() local
70 KMMetricTensor(dXdx, km_g_ij); in AnisotropyTensorCollocate()
H A Ddiff_flux_projection.h16 CeedScalar dXdx[9]; in ComputeDivDiffusiveFluxGeneric() local
18 QdataUnpack_ND(dim, Q, i, q_data, NULL, dXdx); in ComputeDivDiffusiveFluxGeneric()
31 DivergenceND(grad_qn, dXdx, dim, &v[n][i]); in ComputeDivDiffusiveFluxGeneric()
H A Dvelocity_gradient_projection.h20 CeedScalar wdetJ, dXdx[3][3]; in VelocityGradientProjectionRHS() local
21 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in VelocityGradientProjectionRHS()
25 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, (CeedScalar *)Grad_q, dXdx, grad_s); in VelocityGradientProjectionRHS()
H A Dmonitor_totalkineticenergy.h19 CeedScalar wdetJ, dXdx[3][3], vorticity[3], kmstrain_rate[6], strain_rate[3][3]; in MonitorTotalKineticEnergy() local
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()
H A Ddifferential_filter.h88 CeedScalar wdetJ, dXdx[3][3]; in DifferentialFilter_LHS_N() local
89 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); 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()
131 MatVec3(dXdx, dq, CEED_NOTRANSPOSE, dq_dXdx); in DifferentialFilter_LHS_N()
H A Dstabilization.h66 …im(TauDiagCoefficients tau_coeffs, NewtonianIGProperties gas, State s, const CeedScalar dXdx[3][3], in Tau_diagPrim()
83 MatMat3(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in Tau_diagPrim()
H A Dsgs_dd_model.h227 CeedScalar wdetJ, dXdx[3][3]; in IFunction_NodalSgs() local
228 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in IFunction_NodalSgs()
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 Dnewtonian_state.h582 … const CeedScalar *grad_q, const CeedScalar dXdx[3][3], State grad_s[3]) { in StatePhysicalGradientFromReference()
586 MatMatNM((CeedScalar *)grad_qi, (CeedScalar *)dXdx, (CeedScalar *)dq, 5, 3, 3); in StatePhysicalGradientFromReference()
596 … StateVariable state_var, const CeedScalar *grad_q, const CeedScalar dXdx[2][3], in StatePhysicalGradientFromReference_Boundary()
601 …dqi[j] = grad_q[(Q * 5) * 0 + Q * j + i] * dXdx[0][k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1][… in StatePhysicalGradientFromReference_Boundary()
/honee/qfunctions/spanstats/
H A Dcflpe.h35 CeedScalar dXdx[2][2], gijd_mat[2][2] = {{0.}}; in ChildStatsCollection_CflPe() local
37 QdataUnpack_2D(Q, i, q_data, &wdetJ, dXdx); in ChildStatsCollection_CflPe()
40 MatMat2(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in ChildStatsCollection_CflPe()
48 CeedScalar dXdx[3][3], gijd_mat[3][3] = {{0.}}; in ChildStatsCollection_CflPe() local
50 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in ChildStatsCollection_CflPe()
53 MatMat3(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in ChildStatsCollection_CflPe()
/honee/src/
H A Dstrong_boundary_conditions.c14 CeedVector multiplicity, x_stored, scale_stored, stg_data, dXdx, x_coord; in SetupStrongSTG_Ceed() local
119 PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_dXdx, &dXdx, NULL)); in SetupStrongSTG_Ceed()
129 …etscCallCeed(ceed, CeedOperatorSetField(op_setup, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)); in SetupStrongSTG_Ceed()
136 …scCallCeed(ceed, CeedOperatorSetField(op_stgdata, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)); in SetupStrongSTG_Ceed()
146 …Ceed(ceed, CeedOperatorSetField(op_strong_bc_sub, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)); in SetupStrongSTG_Ceed()
159 PetscCallCeed(ceed, CeedVectorDestroy(&dXdx)); in SetupStrongSTG_Ceed()