Home
last modified time | relevance | path

Searched refs:U (Results 1 – 25 of 43) sorted by relevance

12

/libCEED/examples/fluids/qfunctions/
H A Driemann_solver.h57 UnpackState_U(left.U, U_left); in Flux_HLL()
58 UnpackState_U(right.U, U_right); in Flux_HLL()
76 UnpackState_U(left.U, U_l); in Flux_HLL_fwd()
77 UnpackState_U(right.U, U_r); in Flux_HLL_fwd()
82 UnpackState_U(dleft.U, dU_l); in Flux_HLL_fwd()
83 UnpackState_U(dright.U, dU_r); in Flux_HLL_fwd()
105 RoeWeights r = RoeSetup(left.U.density, right.U.density); in ComputeHLLSpeeds_Roe()
126 RoeWeights r = RoeSetup(left.U.density, right.U.density); in ComputeHLLSpeeds_Roe_fwd()
127 RoeWeights dr = RoeSetup_fwd(left.U.density, right.U.density, dleft.U.density, dright.U.density); in ComputeHLLSpeeds_Roe_fwd()
221 CeedScalar fact = side.U.density * (s_side - u_side) / (s_side - s_star); in RiemannFlux_HLLC_Star()
[all …]
H A Dnewtonian_state.h27 StateConservative U; member
31 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U()
32 U[0] = s.density; in UnpackState_U()
33 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U()
34 U[4] = s.E_total; in UnpackState_U()
62 return e_internal + e_kinetic + s.Y.pressure / s.U.density; in TotalSpecificEnthalpy()
68 …turn de_internal + de_kinetic + ds.Y.pressure / s.U.density - s.Y.pressure / Square(s.U.density) *… in TotalSpecificEnthalpy_fwd()
71 … StatePrimitive StatePrimitiveFromConservative(NewtonianIdealGasContext gas, StateConservative U) { in StatePrimitiveFromConservative() argument
73 for (CeedInt i = 0; i < 3; i++) Y.velocity[i] = U.momentum[i] / U.density; in StatePrimitiveFromConservative()
75 CeedScalar e_total = U.E_total / U.density; in StatePrimitiveFromConservative()
[all …]
H A Dgaussianwave.h44 CeedScalar U[5] = {0.}, qi[5] = {0.}; in IC_GaussianWave() local
48 …const CeedScalar e_kinetic = 0.5 * S_infty.U.density * Dot3(S_infty.Y.velocity, S_infty.Y.velocity… in IC_GaussianWave()
52 U[0] = S_infty.U.density * perturbation; in IC_GaussianWave()
53 U[1] = S_infty.Y.velocity[0] * U[0]; in IC_GaussianWave()
54 U[2] = S_infty.Y.velocity[1] * U[0]; in IC_GaussianWave()
55 U[3] = S_infty.Y.velocity[2] * U[0]; in IC_GaussianWave()
56 U[4] = S_infty.Y.pressure / (gamma - 1) * perturbation + e_kinetic; in IC_GaussianWave()
58 State initCond = StateFromU(newt_ctx, U); in IC_GaussianWave()
H A Dturb_spanstats.h28 v[TURB_MEAN_DENSITY][i] = wdetJ * s.U.density; 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()
37 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Z][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
38 v[TURB_MEAN_MOMENTUM_X][i] = wdetJ * s.U.momentum[0]; in ChildStatsCollection()
39 v[TURB_MEAN_MOMENTUM_Y][i] = wdetJ * s.U.momentum[1]; in ChildStatsCollection()
40 v[TURB_MEAN_MOMENTUM_Z][i] = wdetJ * s.U.momentum[2]; in ChildStatsCollection()
41 v[TURB_MEAN_MOMENTUMFLUX_XX][i] = wdetJ * s.U.momentum[0] * s.Y.velocity[0]; in ChildStatsCollection()
42 v[TURB_MEAN_MOMENTUMFLUX_YY][i] = wdetJ * s.U.momentum[1] * s.Y.velocity[1]; in ChildStatsCollection()
[all …]
H A Dnewtonian.h26 CeedScalar U[5]; in InternalDampingLayer() local
27 UnpackState_U(damp_s.U, U); in InternalDampingLayer()
28 for (int i = 0; i < 5; i++) damp_residual[i] += U[i]; in InternalDampingLayer()
86 UnpackState_U(s_dot.U, U_dot); in MassFunction_Newtonian()
163 CeedScalar U[5], wdetJ, dXdx[3][3]; in RHSFunction_Newtonian() local
165 for (int j = 0; j < 5; j++) U[j] = q[j][i]; in RHSFunction_Newtonian()
167 State s = StateFromU(context, U); in RHSFunction_Newtonian()
189 …const CeedScalar body_force[5] = {0, s.U.density * g[0], s.U.density * g[1], s.U.density * g[2], D… in RHSFunction_Newtonian()
262 …const CeedScalar body_force[5] = {0, s.U.density * g[0], s.U.density * g[1], s.U.density * g[2], D… in IFunction_Newtonian()
268 UnpackState_U(s_dot.U, U_dot); in IFunction_Newtonian()
[all …]
H A Dblasius.h117 const CeedScalar x0 = U_infty * S_infty.U.density / (mu * 25 / Square(delta0)); in ICsBlasius()
121 State s = BlasiusSolution(context, x, x0, x_inflow, S_infty.U.density, &t12); in ICsBlasius()
142 const CeedScalar rho_0 = S_infty.U.density; in Blasius_Inflow()
162 s.U.density = s_int.U.density; in Blasius_Inflow()
163 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
165 s.U.E_total = s_int.U.E_total; in Blasius_Inflow()
166 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
199 const CeedScalar rho_0 = S_infty.U.density; in Blasius_Inflow_Jacobian()
H A Dadvection.h239 CeedScalar U[5] = {0.}; in StatePhysicalGradientFromReference_ND() local
240 grad_s[2] = StateFromU(gas, U); in StatePhysicalGradientFromReference_ND()
305 const CeedScalar Grad_E[3] = {grad_s[0].U.E_total, grad_s[1].U.E_total, grad_s[2].U.E_total}; in IFunction_AdvectionGeneric()
318 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in IFunction_AdvectionGeneric()
329 for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = -wdetJ * s.U.E_total * uX[j]; in IFunction_AdvectionGeneric()
335 for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total; in IFunction_AdvectionGeneric()
441 const CeedScalar Grad_E[3] = {grad_s[0].U.E_total, grad_s[1].U.E_total, grad_s[2].U.E_total}; in RHSFunction_AdvectionGeneric()
454 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in RHSFunction_AdvectionGeneric()
463 for (CeedInt j = 0; j < dim; j++) grad_v[j][4][i] = wdetJ * s.U.E_total * uX[j]; in RHSFunction_AdvectionGeneric()
470 for (CeedInt i = 0; i < dim; i++) Fe[i] = -context->diffusion_coeff * grad_s[i].U.E_total; in RHSFunction_AdvectionGeneric()
/libCEED/include/ceed/jit-source/hip/
H A Dhip-ref-basis-nontensor-templates.h20 const CeedScalar *U; in Contract() local
26 U = &d_U[elem * strides_elem_U + comp * strides_comp_U]; in Contract()
29 const CeedScalar val = U[i]; in Contract()
47 const CeedScalar *U; in ContractTranspose() local
55 U = &d_U[elem * strides_elem_U + comp * strides_comp_U + d * strides_q_comp_U]; in ContractTranspose()
56 for (CeedInt i = 0; i < Q; i++) r_V += d_B[t_id + i * P + d * P * Q] * U[i]; in ContractTranspose()
H A Dhip-shared-basis-tensor-templates.h20 inline __device__ void ContractX1d(SharedData_Hip &data, const CeedScalar *U, const CeedScalar *B, … in ContractX1d() argument
22 data.slice[data.t_id_x] = *U; in ContractX1d()
36 inline __device__ void ContractTransposeX1d(SharedData_Hip &data, const CeedScalar *U, const CeedSc… in ContractTransposeX1d() argument
38 data.slice[data.t_id_x] = *U; in ContractTransposeX1d()
129 inline __device__ void ContractX2d(SharedData_Hip &data, const CeedScalar *U, const CeedScalar *B, … in ContractX2d() argument
131 data.slice[data.t_id_x + data.t_id_y * T_1D] = *U; in ContractX2d()
145 inline __device__ void ContractY2d(SharedData_Hip &data, const CeedScalar *U, const CeedScalar *B, … in ContractY2d() argument
147 data.slice[data.t_id_x + data.t_id_y * T_1D] = *U; in ContractY2d()
161 inline __device__ void ContractTransposeY2d(SharedData_Hip &data, const CeedScalar *U, const CeedSc… in ContractTransposeY2d() argument
163 data.slice[data.t_id_x + data.t_id_y * T_1D] = *U; in ContractTransposeY2d()
[all …]
H A Dhip-shared-basis-tensor-flattened-templates.h20 …SharedData_Hip &data, const int t_id_x, const int t_id_y, const CeedScalar *U, const CeedScalar *B, in ContractX2dFlattened() argument
23 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractX2dFlattened()
37 …SharedData_Hip &data, const int t_id_x, const int t_id_y, const CeedScalar *U, const CeedScalar *B, in ContractY2dFlattened() argument
40 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractY2dFlattened()
54 …ransposeY2dFlattened(SharedData_Hip &data, const int t_id_x, const int t_id_y, const CeedScalar *U, in ContractTransposeY2dFlattened() argument
57 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractTransposeY2dFlattened()
71 …ransposeX2dFlattened(SharedData_Hip &data, const int t_id_x, const int t_id_y, const CeedScalar *U, in ContractTransposeX2dFlattened() argument
74 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractTransposeX2dFlattened()
88 …sposeAddX2dFlattened(SharedData_Hip &data, const int t_id_x, const int t_id_y, const CeedScalar *U, in ContractTransposeAddX2dFlattened() argument
91 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractTransposeAddX2dFlattened()
[all …]
H A Dhip-shared-basis-nontensor-templates.h16 inline __device__ void Contract1d(SharedData_Hip &data, const CeedScalar *U, const CeedScalar *B, C… in Contract1d() argument
17 data.slice[data.t_id_x] = *U; in Contract1d()
32 inline __device__ void ContractTranspose1d(SharedData_Hip &data, const CeedScalar *U, const CeedSca… in ContractTranspose1d() argument
33 data.slice[data.t_id_x] = *U; in ContractTranspose1d()
/libCEED/include/ceed/jit-source/cuda/
H A Dcuda-ref-basis-nontensor-templates.h20 const CeedScalar *U; in Contract() local
26 U = d_U + elem * strides_elem_U + comp * strides_comp_U; in Contract()
29 const CeedScalar val = U[i]; in Contract()
47 const CeedScalar *U; in ContractTranspose() local
55 U = &d_U[elem * strides_elem_U + comp * strides_comp_U + d * strides_q_comp_U]; in ContractTranspose()
56 for (CeedInt i = 0; i < Q; i++) r_V += d_B[t_id + i * P + d * P * Q] * U[i]; in ContractTranspose()
H A Dcuda-shared-basis-tensor-templates.h20 inline __device__ void ContractX1d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B,… in ContractX1d() argument
22 data.slice[data.t_id_x] = *U; in ContractX1d()
36 inline __device__ void ContractTransposeX1d(SharedData_Cuda &data, const CeedScalar *U, const CeedS… in ContractTransposeX1d() argument
38 data.slice[data.t_id_x] = *U; in ContractTransposeX1d()
129 inline __device__ void ContractX2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B,… in ContractX2d() argument
131 data.slice[data.t_id_x + data.t_id_y * T_1D] = *U; in ContractX2d()
145 inline __device__ void ContractY2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B,… in ContractY2d() argument
147 data.slice[data.t_id_x + data.t_id_y * T_1D] = *U; in ContractY2d()
161 inline __device__ void ContractTransposeY2d(SharedData_Cuda &data, const CeedScalar *U, const CeedS… in ContractTransposeY2d() argument
163 data.slice[data.t_id_x + data.t_id_y * T_1D] = *U; in ContractTransposeY2d()
[all …]
H A Dcuda-shared-basis-tensor-flattened-templates.h20 …haredData_Cuda &data, const int t_id_x, const int t_id_y, const CeedScalar *U, const CeedScalar *B, in ContractX2dFlattened() argument
23 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractX2dFlattened()
37 …haredData_Cuda &data, const int t_id_x, const int t_id_y, const CeedScalar *U, const CeedScalar *B, in ContractY2dFlattened() argument
40 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractY2dFlattened()
54 …ansposeY2dFlattened(SharedData_Cuda &data, const int t_id_x, const int t_id_y, const CeedScalar *U, in ContractTransposeY2dFlattened() argument
57 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractTransposeY2dFlattened()
71 …ansposeX2dFlattened(SharedData_Cuda &data, const int t_id_x, const int t_id_y, const CeedScalar *U, in ContractTransposeX2dFlattened() argument
74 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractTransposeX2dFlattened()
88 …poseAddX2dFlattened(SharedData_Cuda &data, const int t_id_x, const int t_id_y, const CeedScalar *U, in ContractTransposeAddX2dFlattened() argument
91 data.slice[t_id_x + t_id_y * T_1D] = *U; in ContractTransposeAddX2dFlattened()
[all …]
H A Dcuda-shared-basis-nontensor-templates.h16 inline __device__ void Contract1d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, … in Contract1d() argument
17 data.slice[data.t_id_x] = *U; in Contract1d()
32 inline __device__ void ContractTranspose1d(SharedData_Cuda &data, const CeedScalar *U, const CeedSc… in ContractTranspose1d() argument
33 data.slice[data.t_id_x] = *U; in ContractTranspose1d()
/libCEED/include/ceed/jit-source/sycl/
H A Dsycl-shared-basis-tensor-templates.h19 …(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *restrict U, local const CeedSca… in ContractX1d() argument
23 scratch[item_id_x] = *U; in ContractX1d()
38 …(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *restrict U, local const CeedSca… in ContractTransposeX1d() argument
42 scratch[item_id_x] = *U; in ContractTransposeX1d()
109 …(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *restrict U, local const CeedSca… in ContractX2d() argument
114 scratch[item_id_x + item_id_y * T_1D] = *U; in ContractX2d()
129 …(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *restrict U, local const CeedSca… in ContractY2d() argument
134 scratch[item_id_x + item_id_y * T_1D] = *U; in ContractY2d()
149 …(const CeedInt P_1D, const CeedInt Q_1D, private const CeedScalar *restrict U, local const CeedSca… in ContractTransposeY2d() argument
154 scratch[item_id_x + item_id_y * T_1D] = *U; in ContractTransposeY2d()
[all …]
/libCEED/examples/solids/include/
H A Dmisc.h29 PetscErrorCode FormJacobian(SNES snes, Vec U, Mat J, Mat J_pre, void *ctx);
34 PetscErrorCode ViewSolution(MPI_Comm comm, AppCtx app_ctx, Vec U, PetscInt increment, PetscScalar l…
36 …nosticQuantities(MPI_Comm comm, DM dm_U, UserMult user, AppCtx app_ctx, Vec U, CeedElemRestriction…
/libCEED/examples/petsc/
H A Darea.c68 Vec U, U_loc, V, V_loc; in main() local
138 PetscCall(DMCreateGlobalVector(dm, &U)); in main()
139 PetscCall(VecGetLocalSize(U, &l_size)); in main()
140 PetscCall(VecGetSize(U, &g_size)); in main()
143 PetscCall(VecDuplicate(U, &V)); in main()
212 PetscCall(VecDestroy(&U)); in main()
/libCEED/backends/ref/
H A Dceed-ref-qfunction.c17 static int CeedQFunctionApply_Ref(CeedQFunction qf, CeedInt Q, CeedVector *U, CeedVector *V) { in CeedQFunctionApply_Ref() argument
29 CeedCallBackend(CeedVectorGetArrayRead(U[i], CEED_MEM_HOST, &impl->inputs[i])); in CeedQFunctionApply_Ref()
38 CeedCallBackend(CeedVectorRestoreArrayRead(U[i], &impl->inputs[i])); in CeedQFunctionApply_Ref()
/libCEED/backends/memcheck/
H A Dceed-memcheck-qfunction.c19 static int CeedQFunctionApply_Memcheck(CeedQFunction qf, CeedInt Q, CeedVector *U, CeedVector *V) { in CeedQFunctionApply_Memcheck() argument
37 CeedCallBackend(CeedVectorGetArrayRead(U[i], CEED_MEM_HOST, &impl->inputs[i])); in CeedQFunctionApply_Memcheck()
39 CeedCallBackend(CeedVectorGetLength(U[i], &len)); in CeedQFunctionApply_Memcheck()
64 CeedCallBackend(CeedVectorRestoreArrayRead(U[i], &impl->inputs[i])); in CeedQFunctionApply_Memcheck()
/libCEED/examples/solids/src/
H A Dmisc.c78 PetscErrorCode FormJacobian(SNES snes, Vec U, Mat J, Mat J_pre, void *ctx) { in FormJacobian() argument
112 PetscErrorCode ViewSolution(MPI_Comm comm, AppCtx app_ctx, Vec U, PetscInt increment, PetscScalar l… in ViewSolution() argument
130 PetscCall(VecGetDM(U, &dm)); in ViewSolution()
135 PetscCall(VecView(U, viewer)); in ViewSolution()
144 PetscErrorCode ViewDiagnosticQuantities(MPI_Comm comm, DM dmU, UserMult user, AppCtx app_ctx, Vec U… in ViewDiagnosticQuantities() argument
170 PetscCall(DMGlobalToLocal(dmU, U, INSERT_VALUES, user->X_loc)); in ViewDiagnosticQuantities()
/libCEED/backends/cuda-ref/
H A Dceed-cuda-ref-qfunction.c21 static int CeedQFunctionApply_Cuda(CeedQFunction qf, CeedInt Q, CeedVector *U, CeedVector *V) { in CeedQFunctionApply_Cuda() argument
38 CeedCallBackend(CeedVectorGetArrayRead(U[i], CEED_MEM_DEVICE, &data->fields.inputs[i])); in CeedQFunctionApply_Cuda()
53 CeedCallBackend(CeedVectorRestoreArrayRead(U[i], &data->fields.inputs[i])); in CeedQFunctionApply_Cuda()
/libCEED/backends/hip-ref/
H A Dceed-hip-ref-qfunction.c21 static int CeedQFunctionApply_Hip(CeedQFunction qf, CeedInt Q, CeedVector *U, CeedVector *V) { in CeedQFunctionApply_Hip() argument
39 CeedCallBackend(CeedVectorGetArrayRead(U[i], CEED_MEM_DEVICE, &data->fields.inputs[i])); in CeedQFunctionApply_Hip()
55 CeedCallBackend(CeedVectorRestoreArrayRead(U[i], &data->fields.inputs[i])); in CeedQFunctionApply_Hip()
/libCEED/examples/solids/
H A Delasticity.c54 Vec U, *U_g, *U_loc; // U: solution, R: residual, F: forcing in main() local
203 PetscCall(VecDuplicate(U_g[fine_level], &U)); in main()
234 PetscCall(VecDuplicate(U, &neumann_bcs)); in main()
329 PetscCall(VecGetType(U, &vecType)); in main()
579 PetscCall(PetscObjectSetName((PetscObject)U, "")); in main()
580 PetscCall(VecSet(U, 0.0)); in main()
584 PetscCall(ViewSolution(comm, app_ctx, U, 0, 0.0)); in main()
618 PetscCall(SNESSolve(snes, F, U)); in main()
622 PetscCall(ViewSolution(comm, app_ctx, U, increment, load_increment)); in main()
639 … PetscCall(ComputeStrainEnergy(dm_energy, res_ctx, ceed_data[fine_level]->op_energy, U, &energy)); in main()
[all …]
/libCEED/backends/sycl-ref/
H A Dceed-sycl-ref-qfunction.sycl.cpp26 static int CeedQFunctionApply_Sycl(CeedQFunction qf, CeedInt Q, CeedVector *U, CeedVector *V) { in CeedQFunctionApply_Sycl() argument
46 const CeedVector *U_i = U; in CeedQFunctionApply_Sycl()
93 U_i = U; in CeedQFunctionApply_Sycl()

12