Home
last modified time | relevance | path

Searched refs:wdetJ (Results 1 – 16 of 16) sorted by relevance

/libCEED/examples/fluids/qfunctions/
H A Dturb_spanstats.h23 const CeedScalar wdetJ = q_data[0][i] * delta_t; in ChildStatsCollection() local
28 v[TURB_MEAN_DENSITY][i] = wdetJ * s.U.density; in ChildStatsCollection()
29 v[TURB_MEAN_PRESSURE][i] = wdetJ * s.Y.pressure; in ChildStatsCollection()
30 v[TURB_MEAN_PRESSURE_SQUARED][i] = wdetJ * Square(s.Y.pressure); in ChildStatsCollection()
31 v[TURB_MEAN_PRESSURE_VELOCITY_X][i] = wdetJ * s.Y.pressure * s.Y.velocity[0]; in ChildStatsCollection()
32 v[TURB_MEAN_PRESSURE_VELOCITY_Y][i] = wdetJ * s.Y.pressure * s.Y.velocity[1]; in ChildStatsCollection()
33 v[TURB_MEAN_PRESSURE_VELOCITY_Z][i] = wdetJ * s.Y.pressure * s.Y.velocity[2]; in ChildStatsCollection()
34 v[TURB_MEAN_DENSITY_TEMPERATURE][i] = wdetJ * s.U.density * s.Y.temperature; in ChildStatsCollection()
35 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
36 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
[all …]
H A Ddifferential_filter.h44 const CeedScalar wdetJ = q_data[0][i]; in DifferentialFilter_RHS() local
47 v0[DIFF_FILTER_PRESSURE][i] = wdetJ * s.Y.pressure; in DifferentialFilter_RHS()
48 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
50 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
51 v0[DIFF_FILTER_TEMPERATURE][i] = wdetJ * s.Y.temperature; in DifferentialFilter_RHS()
52 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
53 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
54 v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2]; in DifferentialFilter_RHS()
55 v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2]; in DifferentialFilter_RHS()
[all …]
H A Dsetupgeo2d.h53 const CeedScalar wdetJ = w[i] * detJ; in Setup2d() local
55 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data); in Setup2d()
94 const CeedScalar wdetJ = w[i] * detJb; in SetupBoundary2d() local
96 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur); in SetupBoundary2d()
148 const CeedScalar wdetJ = w[i] * detJb; in Setup2D_3Dcoords() local
150 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur); in Setup2D_3Dcoords()
H A Dadvection.h196 …ND(CeedInt N, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar *dXdx)… in QdataUnpack_ND() argument
200 StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataUnpack_ND()
204 StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataUnpack_ND()
210 …_ND(CeedInt N, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar *dXdx, in QdataBoundaryUnpack_ND() argument
215 if (wdetJ) StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataBoundaryUnpack_ND()
219 if (wdetJ) StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataBoundaryUnpack_ND()
300 CeedScalar wdetJ, dXdx[9]; in IFunction_AdvectionGeneric() local
301 QdataUnpack_ND(dim, Q, i, q_data, &wdetJ, dXdx); in IFunction_AdvectionGeneric()
309 v[f][i] = wdetJ * q_dot[f][i]; // K Mass/transient term in IFunction_AdvectionGeneric()
321 v[4][i] = wdetJ * q_dot[4][i]; // transient part (ALWAYS) in IFunction_AdvectionGeneric()
[all …]
H A Dsetupgeo.h57 const CeedScalar wdetJ = w[i] * detJ; in Setup() local
59 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data); in Setup()
117 const CeedScalar wdetJ = w[i] * detJb; in SetupBoundary() local
119 StoredValuesPack(Q, i, 0, 1, &wdetJ, q_data_sur); in SetupBoundary()
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()
80 v[f][i] = wdetJ * qi_dot[f]; 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()
186 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] = wdetJ * (dXdx[k][0] * Flux[j][0] + dXdx[k][1] * … in RHSFunction_Newtonian()
190 for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j]; in RHSFunction_Newtonian()
196 for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j]; in RHSFunction_Newtonian()
205 …for (CeedInt k = 0; k < 3; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] *… in RHSFunction_Newtonian()
[all …]
H A Dutils.h245 …dataUnpack_3D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar dXdx[3… in QdataUnpack_3D() argument
246 StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataUnpack_3D()
263 …daryUnpack_3D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar dXdx[2… in QdataBoundaryUnpack_3D() argument
265 if (wdetJ) StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataBoundaryUnpack_3D()
282 …dataUnpack_2D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar dXdx[2… in QdataUnpack_2D() argument
283 StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataUnpack_2D()
299 …daryUnpack_2D(CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedScalar normal… in QdataBoundaryUnpack_2D() argument
300 if (wdetJ) StoredValuesUnpack(Q, i, 0, 1, q_data, wdetJ); in QdataBoundaryUnpack_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()
36 v[j * 3 + k][i] = wdetJ * grad_velocity[j][k]; in VelocityGradientProjectionRHS()
H A Dshocktube.h248 CeedScalar wdetJ, dXdx[3][3]; in EulerShockTube() local
249 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in EulerShockTube()
280 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd… in EulerShockTube()
285 …dv[k][j + 1][i] += wdetJ * ((rho * u[j] * u[0] + (j == 0 ? P : 0)) * dXdx[k][0] + (rho * u[j] * u[… in EulerShockTube()
291 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[… in EulerShockTube()
316 for (CeedInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * nu_shock * drhodx[j]; in EulerShockTube()
319 for (CeedInt j = 0; j < 3; j++) dv[j][k][i] -= wdetJ * nu_shock * du[k][j]; in EulerShockTube()
322 for (CeedInt j = 0; j < 3; j++) dv[j][4][i] -= wdetJ * nu_shock * dEdx[j]; in EulerShockTube()
366 …for (CeedInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXd… in EulerShockTube()
H A Dgrid_anisotropy_tensor.h43 const CeedScalar wdetJ = q_data[0][i]; in AnisotropyTensorProjection() local
55 for (CeedInt j = 0; j < 6; j++) v[j][i] = wdetJ * km_A_ij[j]; in AnisotropyTensorProjection()
56 v[6][i] = wdetJ * delta; in AnisotropyTensorProjection()
H A Deulervortex.h293 CeedScalar wdetJ, dXdx[3][3]; in Euler() local
294 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in Euler()
323 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] += wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd… in Euler()
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()
334 …for (CeedInt j = 0; j < 3; j++) dv[j][4][i] += wdetJ * (E + P) * (u[0] * dXdx[j][0] + u[1] * dXdx[… 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()
416 CeedScalar wdetJ, dXdx[3][3]; in IFunction_Euler() local
417 QdataUnpack_3D(Q, i, q_data, &wdetJ, dXdx); in IFunction_Euler()
443 for (CeedInt j = 0; j < 5; j++) v[j][i] += wdetJ * q_dot[j][i]; in IFunction_Euler()
447 …for (CeedInt j = 0; j < 3; j++) dv[j][0][i] -= wdetJ * (rho * u[0] * dXdx[j][0] + rho * u[1] * dXd… in IFunction_Euler()
[all …]
/libCEED/examples/solids/qfunctions/
H A Dconstant-force.h41 CeedScalar wdetJ = q_data[i]; in SetupConstantForce() local
45 force[i + 0 * Q] = forcing_vector[0] * wdetJ; in SetupConstantForce()
48 force[i + 1 * Q] = forcing_vector[1] * wdetJ; in SetupConstantForce()
51 force[i + 2 * Q] = forcing_vector[2] * wdetJ; in SetupConstantForce()
H A Dmanufactured-force.h44 CeedScalar wdetJ = q_data[i]; in SetupMMSForce() local
55 wdetJ / 1e8; in SetupMMSForce()
65 wdetJ / 1e8; in SetupMMSForce()
75 wdetJ / 1e8; in SetupMMSForce()
H A Dlinear.h51 const CeedScalar wdetJ = q_data[0][i]; in ElasResidual_Linear() local
113 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * sigma[j][m] * wdetJ; in ElasResidual_Linear()
148 const CeedScalar wdetJ = q_data[0][i]; in ElasJacobian_Linear() local
208 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dsigma[j][m] * wdetJ; in ElasJacobian_Linear()
246 const CeedScalar wdetJ = q_data[0][i]; in ElasEnergy_Linear() local
276 … + strain_vol * mu + (e[0][1] * e[0][1] + e[0][2] * e[0][2] + e[1][2] * e[1][2]) * 2 * mu) * wdetJ; in ElasEnergy_Linear()
H A Dfinite-strain-neo-hookean.h188 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_NH() local
242 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * P[j][m] * wdetJ; in ElasFSResidual_NH()
283 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_NH() local
387 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * deltaP[j][m] * wdetJ; in ElasFSJacobian_NH()
423 const CeedScalar wdetJ = q_data[0][i]; in ElasFSEnergy_NH() local
458 …[i] = (lambda * logJ * logJ / 2. - mu * logJ + mu * (E2[0][0] + E2[1][1] + E2[2][2]) / 2.) * wdetJ; in ElasFSEnergy_NH()
H A Dfinite-strain-mooney-rivlin.h190 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_MR() local
245 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * P[j][m] * wdetJ; in ElasFSResidual_MR()
281 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_MR() local
408 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dP[j][m] * wdetJ; in ElasFSJacobian_MR()
441 const CeedScalar wdetJ = q_data[0][i]; in ElasFSEnergy_MR() local
500 …) * (logJ) - (mu_1 + 2 * mu_2) * logJ + (mu_1 / 2.) * (I_1 - 3) + (mu_2 / 2.) * (I_2 - 3)) * wdetJ; in ElasFSEnergy_MR()