Home
last modified time | relevance | path

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

/libCEED/examples/fluids/qfunctions/
H A Dsetupgeo_helpers.h41 …Jacobian_3D(CeedInt Q, CeedInt i, const CeedScalar (*dxdX_q)[3][CEED_Q_VLA], CeedScalar dXdx[3][3], in InvertMappingJacobian_3D()
63 dXdx[0][0] = A11 / detJ; in InvertMappingJacobian_3D()
64 dXdx[0][1] = A12 / detJ; in InvertMappingJacobian_3D()
65 dXdx[0][2] = A13 / detJ; in InvertMappingJacobian_3D()
66 dXdx[1][0] = A21 / detJ; in InvertMappingJacobian_3D()
67 dXdx[1][1] = A22 / detJ; in InvertMappingJacobian_3D()
68 dXdx[1][2] = A23 / detJ; in InvertMappingJacobian_3D()
69 dXdx[2][0] = A31 / detJ; in InvertMappingJacobian_3D()
70 dXdx[2][1] = A32 / detJ; in InvertMappingJacobian_3D()
71 dXdx[2][2] = A33 / detJ; in InvertMappingJacobian_3D()
[all …]
H A Dstrong_boundary_conditions.h22 CeedScalar dXdx[2][3]; in SetupStrongBC() local
23 InvertBoundaryMappingJacobian_3D(Q, i, dxdX_q, dXdx); in SetupStrongBC()
24 dXdx_q[0][i] = dXdx[0][0]; in SetupStrongBC()
25 dXdx_q[1][i] = dXdx[0][1]; in SetupStrongBC()
26 dXdx_q[2][i] = dXdx[0][2]; in SetupStrongBC()
27 dXdx_q[3][i] = dXdx[1][0]; in SetupStrongBC()
28 dXdx_q[4][i] = dXdx[1][1]; in SetupStrongBC()
29 dXdx_q[5][i] = dXdx[1][2]; in SetupStrongBC()
H A Dshocktube.h139 …N_HELPER CeedScalar Covariant_length_along_vector(CeedScalar vec[3], const CeedScalar dXdx[3][3]) { in Covariant_length_along_vector()
144 vec_dot_jacobian[i] += dXdx[j][i] * vec[i]; in Covariant_length_along_vector()
168 CEED_QFUNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const Ceed… in Tau_spatial()
172 …CeedScalar h = 2 / sqrt(dXdx[0][i] * dXdx[0][i] + dXdx[1][i] * dXdx[1][i] + dXdx[2][i] * dXdx[2][i… in Tau_spatial()
248 CeedScalar wdetJ, dXdx[3][3]; in EulerShockTube() local
249 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in EulerShockTube()
259 drhodx[j] += drho[k] * dXdx[k][j]; in EulerShockTube()
260 dEdx[j] += dE[k] * dXdx[k][j]; in EulerShockTube()
262 dUdx[j][k] += dU[j][l] * dXdx[l][k]; in EulerShockTube()
263 dXdxdXdxT[j][k] += dXdx[j][l] * dXdx[k][l]; // dXdx_j,k * dXdx_k,j in EulerShockTube()
[all …]
H A Deulervortex.h215 CEED_QFUNCTION_HELPER void Tau_spatial(CeedScalar Tau_x[3], const CeedScalar dXdx[3][3], const Ceed… in Tau_spatial()
219 …CeedScalar h = 2 / sqrt(dXdx[0][i] * dXdx[0][i] + dXdx[1][i] * dXdx[1][i] + dXdx[2][i] * dXdx[2][i… in Tau_spatial()
293 CeedScalar wdetJ, dXdx[3][3]; in Euler() local
294 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in Euler()
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()
323 … < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXdx[j][1] + rho * u[2] … in Euler()
328 …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 Dnewtonian.h75 CeedScalar wdetJ, dXdx[3][3]; in MassFunction_Newtonian() local
76 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in MassFunction_Newtonian()
87 Tau_diagPrim(context, s, dXdx, context->dt, Tau_d); in MassFunction_Newtonian()
93 …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()
163 CeedScalar U[5], wdetJ, dXdx[3][3]; in RHSFunction_Newtonian() local
166 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in RHSFunction_Newtonian()
170 … StatePhysicalGradientFromReference(Q, i, context, s, STATEVAR_CONSERVATIVE, Grad_q, dXdx, grad_s); in RHSFunction_Newtonian()
186 …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()
201 Tau_diagPrim(context, s, dXdx, dt, Tau_d); in RHSFunction_Newtonian()
205 …; 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 Dsetupgeo.h55 CeedScalar detJ, dXdx[3][3]; in Setup() local
56 InvertMappingJacobian_3D(Q, i, J, dXdx, &detJ); in Setup()
60 StoredValuesPack(Q, i, 1, 9, (const CeedScalar *)dXdx, q_data); in Setup()
113 CeedScalar detJb, normal[3], dXdx[2][3]; in SetupBoundary() local
116 InvertBoundaryMappingJacobian_3D(Q, i, J, dXdx); in SetupBoundary()
121 StoredValuesPack(Q, i, 4, 6, (const CeedScalar *)dXdx, q_data_sur); in SetupBoundary()
H A Dadvection.h196 …D(CeedInt N, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar *dXdx) { in QdataUnpack_ND() argument
201 StoredValuesUnpack(Q, i, 1, 4, q_data, dXdx); in QdataUnpack_ND()
205 StoredValuesUnpack(Q, i, 1, 9, q_data, dXdx); in QdataUnpack_ND()
210 …_ND(CeedInt N, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar *dXdx, in QdataBoundaryUnpack_ND() argument
221 if (dXdx) StoredValuesUnpack(Q, i, 4, 6, q_data, (CeedScalar *)dXdx); in QdataBoundaryUnpack_ND()
228 … StateVariable state_var, const CeedScalar *grad_q, const CeedScalar *dXdx, in StatePhysicalGradientFromReference_ND() argument
235 …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()
247 …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()
248 grad_q[(Q * 5) * 2 + Q * j + i] * dXdx[2 * N + k]; in StatePhysicalGradientFromReference_ND()
257 …LPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXdx, CeedInt dim) { in Tau() argument
[all …]
H A Dsetupgeo2d.h51 CeedScalar dXdx[2][2], detJ; in Setup2d() local
52 InvertMappingJacobian_2D(Q, i, J, dXdx, &detJ); in Setup2d()
56 StoredValuesPack(Q, i, 1, 4, (const CeedScalar *)dXdx, q_data); in Setup2d()
144 CeedScalar detJb, normal[3], dXdx[2][3]; in Setup2D_3Dcoords() local
147 InvertBoundaryMappingJacobian_3D(Q, i, J, dXdx); in Setup2D_3Dcoords()
151 StoredValuesPack(Q, i, 1, 6, (const CeedScalar *)dXdx, q_data_sur); in Setup2D_3Dcoords()
H A Dgrid_anisotropy_tensor.h44 const CeedScalar dXdx[3][3] = { in AnisotropyTensorProjection() local
51 KMMetricTensor(dXdx, km_g_ij); in AnisotropyTensorProjection()
67 const CeedScalar dXdx[3][3] = { in AnisotropyTensorCollocate() local
74 KMMetricTensor(dXdx, km_g_ij); in AnisotropyTensorCollocate()
H A Dutils.h179 CEED_QFUNCTION_HELPER void KMMetricTensor(const CeedScalar dXdx[3][3], CeedScalar km_g_ij[6]) { in KMMetricTensor()
181 MatMat3(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, g_ij); in KMMetricTensor()
245 …pack_3D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar dXdx[3][3]) { in QdataUnpack_3D()
247 StoredValuesUnpack(Q, i, 1, 9, q_data, (CeedScalar *)dXdx); in QdataUnpack_3D()
263 …Unpack_3D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar dXdx[2][3], in QdataBoundaryUnpack_3D()
267 if (dXdx) StoredValuesUnpack(Q, i, 4, 6, q_data, (CeedScalar *)dXdx); in QdataBoundaryUnpack_3D()
282 …pack_2D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar dXdx[2][2]) { in QdataUnpack_2D()
284 StoredValuesUnpack(Q, i, 1, 4, q_data, (CeedScalar *)dXdx); in QdataUnpack_2D()
H A Dvelocity_gradient_projection.h24 CeedScalar wdetJ, dXdx[3][3]; in VelocityGradientProjectionRHS() local
25 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in VelocityGradientProjectionRHS()
29 …StatePhysicalGradientFromReference(Q, i, context, s, state_var, (CeedScalar *)Grad_q, dXdx, grad_s… in VelocityGradientProjectionRHS()
H A Dstg_shur14.h255 const CeedScalar dXdx[2][3] = { in StgShur14Preprocess() local
262 …r (CeedInt j = 1; j < 3; j++) h_node_sep[j] = 2 / sqrt(dXdx[0][j] * dXdx[0][j] + dXdx[1][j] * dXdx… in StgShur14Preprocess()
298 CeedScalar dXdx[3][3]; in ICsStg() local
299 InvertMappingJacobian_3D(Q, i, J, dXdx, NULL); in ICsStg()
302 …t j = 1; j < 3; j++) h_node_sep[j] = 2 / sqrt(Square(dXdx[0][j]) + Square(dXdx[1][j]) + Square(dXd… in ICsStg()
354 CeedScalar wdetJb, dXdx[2][3], norm[3]; in StgShur14Inflow() local
355 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, norm); in StgShur14Inflow()
360 … for (CeedInt j = 1; j < 3; j++) h_node_sep[j] = 2 / sqrt(Square(dXdx[0][j]) + Square(dXdx[1][j])); in StgShur14Inflow()
488 const CeedScalar dXdx[2][3] = { in StgShur14InflowStrongQF() local
495 … for (CeedInt j = 1; j < 3; j++) h_node_sep[j] = 2 / sqrt(Square(dXdx[0][j]) + Square(dXdx[1][j])); in StgShur14InflowStrongQF()
H A Dbc_freestream.h189 CeedScalar wdetJb, dXdx[2][3], normal[3]; in RiemannOutflow() local
190 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in RiemannOutflow()
206 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_q, dXdx, grad_s); in RiemannOutflow()
258 CeedScalar wdetJb, dXdx[2][3], normal[3]; in RiemannOutflow_Jacobian() local
259 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in RiemannOutflow_Jacobian()
287 … StatePhysicalGradientFromReference_Boundary(Q, i, gas, s_int, state_var, Grad_dq, dXdx, grad_ds); in RiemannOutflow_Jacobian()
341 CeedScalar wdetJb, dXdx[2][3], normal[3]; in PressureOutflow() local
342 QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, dXdx, normal); in PressureOutflow()
346 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in PressureOutflow()
399 CeedScalar wdetJb, dXdx[2][3], normal[3]; in PressureOutflow_Jacobian() local
[all …]
H A Dstabilization.h71 …void Tau_diagPrim(NewtonianIdealGasContext gas, State s, const CeedScalar dXdx[3][3], const CeedSc… in Tau_diagPrim()
88 MatMat3(dXdx, dXdx, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in Tau_diagPrim()
H A Ddifferential_filter.h92 CeedScalar wdetJ, dXdx[3][3]; in DifferentialFilter_LHS_N() local
93 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in DifferentialFilter_LHS_N()
133 …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()
135 MatVec3(dXdx, dq, CEED_NOTRANSPOSE, dq_dXdx); in DifferentialFilter_LHS_N()
H A Dnewtonian_state.h585 … const CeedScalar *grad_q, const CeedScalar dXdx[3][3], State grad_s[3]) { in StatePhysicalGradientFromReference()
590 …(Q * 5) * 0 + Q * j + i] * dXdx[0][k] + grad_q[(Q * 5) * 1 + Q * j + i] * dXdx[1][k] + grad_q[(Q *… in StatePhysicalGradientFromReference()
597 … StateVariable state_var, const CeedScalar *grad_q, const CeedScalar dXdx[2][3], in StatePhysicalGradientFromReference_Boundary()
602 …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()
/libCEED/examples/solids/qfunctions/
H A Dlinear.h52 const CeedScalar dXdx[3][3] = { in ElasResidual_Linear() local
65 for (CeedInt m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasResidual_Linear()
113 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * sigma[j][m] * wdetJ; in ElasResidual_Linear()
149 const CeedScalar dXdx[3][3] = { in ElasJacobian_Linear() local
162 for (CeedInt m = 0; m < 3; m++) graddeltau[j][k] += dXdx[m][k] * deltadu[j][m]; in ElasJacobian_Linear()
208 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dsigma[j][m] * wdetJ; in ElasJacobian_Linear()
247 const CeedScalar dXdx[3][3] = { in ElasEnergy_Linear() local
260 for (CeedInt m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasEnergy_Linear()
314 const CeedScalar dXdx[3][3] = { in ElasDiagnostic_Linear() local
327 for (CeedInt m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasDiagnostic_Linear()
H A Dfinite-strain-neo-hookean.h189 const CeedScalar dXdx[3][3] = { in ElasFSResidual_NH() local
201 for (CeedInt m = 0; m < 3; m++) grad_u[j][k][i] += dXdx[m][k] * du[j][m]; in ElasFSResidual_NH()
242 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * P[j][m] * wdetJ; in ElasFSResidual_NH()
284 const CeedScalar dXdx[3][3] = { in ElasFSJacobian_NH() local
297 for (CeedInt m = 0; m < 3; m++) graddeltau[j][k] += dXdx[m][k] * deltadu[j][m]; in ElasFSJacobian_NH()
387 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * deltaP[j][m] * wdetJ; in ElasFSJacobian_NH()
424 const CeedScalar dXdx[3][3] = { in ElasFSEnergy_NH() local
437 for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasFSEnergy_NH()
494 const CeedScalar dXdx[3][3] = { in ElasFSDiagnostic_NH() local
507 for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasFSDiagnostic_NH()
H A Dfinite-strain-mooney-rivlin.h191 const CeedScalar dXdx[3][3] = { in ElasFSResidual_MR() local
203 for (CeedInt m = 0; m < 3; m++) grad_u[j][k][i] += dXdx[m][k] * du[j][m]; in ElasFSResidual_MR()
245 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * P[j][m] * wdetJ; in ElasFSResidual_MR()
282 const CeedScalar dXdx[3][3] = { in ElasFSJacobian_MR() local
296 for (CeedInt m = 0; m < 3; m++) graddeltau[j][k] += dXdx[m][k] * deltadu[j][m]; in ElasFSJacobian_MR()
408 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dP[j][m] * wdetJ; in ElasFSJacobian_MR()
442 const CeedScalar dXdx[3][3] = { in ElasFSEnergy_MR() local
454 for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasFSEnergy_MR()
533 const CeedScalar dXdx[3][3] = { in ElasFSDiagnostic_MR() local
546 for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasFSDiagnostic_MR()
/libCEED/examples/fluids/src/
H A Dstrong_boundary_conditions.c18 CeedVector multiplicity, x_stored, scale_stored, stg_data, dXdx; in SetupStrongSTG_Ceed() local
73 PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_dXdx, &dXdx, NULL)); in SetupStrongSTG_Ceed()
83 …etscCallCeed(ceed, CeedOperatorSetField(op_setup, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)); in SetupStrongSTG_Ceed()
90 …scCallCeed(ceed, CeedOperatorSetField(op_stgdata, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)); in SetupStrongSTG_Ceed()
100 …Ceed(ceed, CeedOperatorSetField(op_strong_bc_sub, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)); in SetupStrongSTG_Ceed()
113 PetscCallCeed(ceed, CeedVectorDestroy(&dXdx)); in SetupStrongSTG_Ceed()