Home
last modified time | relevance | path

Searched refs:v (Results 1 – 25 of 31) sorted by relevance

12

/honee/qfunctions/spanstats/
H A Dturbulence.h44 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ChildStatsCollection() local
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()
63v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
64v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection()
[all …]
H A Dcflpe.h22 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ChildStatsCollection_CflPe() local
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/qfunctions/
H A Dutils.h43 CEED_QFUNCTION_HELPER CeedScalar DotN(const CeedScalar *u, const CeedScalar *v, const CeedInt N) { in DotN() argument
45 CeedPragmaSIMD for (CeedInt i = 0; i < N; i++) output += u[i] * v[i]; in DotN()
55 …dScalar Dot3(const CeedScalar *u, const CeedScalar *v) { return u[0] * v[0] + u[1] * v[1] + u[2] *… in Dot3() argument
58 …PER CeedScalar Dot2(const CeedScalar *u, const CeedScalar *v) { return u[0] * v[0] + u[1] * v[1]; } in Dot2() argument
67 CEED_QFUNCTION_HELPER void Cross3(const CeedScalar u[3], const CeedScalar v[3], CeedScalar w[3]) { in Cross3()
68 w[0] = (u[1] * v[2]) - (u[2] * v[1]); in Cross3()
69 w[1] = (u[2] * v[0]) - (u[0] * v[2]); in Cross3()
70 w[2] = (u[0] * v[1]) - (u[1] * v[0]); in Cross3()
74 CEED_QFUNCTION_HELPER void Curl3(const CeedScalar gradient[3][3], CeedScalar v[3]) { in Curl3()
75 v[0] = gradient[2][1] - gradient[1][2]; in Curl3()
[all …]
H A Dsgs_dd_training.h30 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ComputeSGS_DDAnisotropicTrainingDataNodal() local
50 for (int j = 0; j < 6; j++) v[j][i] = inv_multiplicity[i] * inputs[j]; in ComputeSGS_DDAnisotropicTrainingDataNodal()
52v[0 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XX][i] - Square(s.Y.velocity[0])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
53v[1 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YY][i] - Square(s.Y.velocity[1])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
54v[2 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] - Square(s.Y.velocity[2])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal()
55v[3 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] - s.Y.velocity[1] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
56v[4 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] - s.Y.velocity[0] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
57v[5 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XY][i] - s.Y.velocity[0] * s.Y.velocity[1]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
H A Dchannel.h82 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Channel_Inflow() local
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()
132 v[4][i] -= wdetJb * u_normal * (E + P); in Channel_Inflow()
143 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Channel_Outflow() local
159 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in Channel_Outflow()
166 v[0][i] -= wdetJb * rho * u_normal; in Channel_Outflow()
169 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Channel_Outflow()
172 v[4][i] -= wdetJb * u_normal * (E + P); in Channel_Outflow()
H A Dmonitor_totalkineticenergy.h13 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in MonitorTotalKineticEnergy() local
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()
H A Dgrid_anisotropy_tensor.h36 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in AnisotropyTensorProjection() local
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()
60 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in AnisotropyTensorCollocate() local
74 for (CeedInt j = 0; j < 6; j++) v[j][i] = km_A_ij[j]; in AnisotropyTensorCollocate()
75 v[6][i] = delta; in AnisotropyTensorCollocate()
H A Dutils_eigensolver_jacobi.h216 for (CeedInt v = 0; v < N; v++) { in ApplyRotLeft() local
217 CeedScalar Aiv = A[i * N + v]; in ApplyRotLeft()
218 A[i * N + v] = rotmat_cst[0] * A[i * N + v] - rotmat_cst[1] * A[j * N + v]; in ApplyRotLeft()
219 A[j * N + v] = rotmat_cst[1] * Aiv + rotmat_cst[0] * A[j * N + v]; in ApplyRotLeft()
H A Dadvection.h275 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in IFunction_AdvectionGeneric() local
293 v[f][i] = wdetJ * q_dot[f][i]; // K Mass/transient term in IFunction_AdvectionGeneric()
306 v[4][i] = wdetJ * q_dot[4][i]; // transient part (ALWAYS) in IFunction_AdvectionGeneric()
308 v[4][i] += wdetJ * strong_conv; in IFunction_AdvectionGeneric()
352 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in MassFunction_AdvectionGeneric() local
366 v[f][i] = wdetJ * q_dot[f][i]; // K Mass/transient term in MassFunction_AdvectionGeneric()
370 v[4][i] = wdetJ * q_dot[4][i]; in MassFunction_AdvectionGeneric()
410 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in RHSFunction_AdvectionGeneric() local
428 v[f][i] = 0.; in RHSFunction_AdvectionGeneric()
443 v[4][i] = -wdetJ * strong_conv; in RHSFunction_AdvectionGeneric()
[all …]
H A Deulervortex.h267 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Euler() local
311 v[j][i] = 0.; in Euler()
388 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in IFunction_Euler() local
431 v[j][i] = 0.; in IFunction_Euler()
435 for (CeedInt j = 0; j < 5; j++) v[j][i] += wdetJ * q_dot[j][i]; in IFunction_Euler()
521 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in TravelingVortex_Inflow() local
550 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in TravelingVortex_Inflow()
562 v[0][i] -= wdetJb * rho_inlet * face_normal; in TravelingVortex_Inflow()
565 …for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho_inlet * face_normal * mean_velocity[j… in TravelingVortex_Inflow()
568 v[4][i] -= wdetJb * face_normal * (E_inlet + P_inlet); in TravelingVortex_Inflow()
[all …]
H A Ddiff_flux_projection.h13 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ComputeDivDiffusiveFluxGeneric() local
30 v[n][i] = 0; in ComputeDivDiffusiveFluxGeneric()
31 DivergenceND(grad_qn, dXdx, dim, &v[n][i]); in ComputeDivDiffusiveFluxGeneric()
H A Dnewtonian.h47 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in MassFunction_Newtonian() local
63 v[f][i] = wdetJ * qi_dot[f]; in MassFunction_Newtonian()
196 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in RHSFunction_Newtonian() local
231 for (int j = 0; j < 5; j++) v[j][i] = wdetJ * body_force[j]; in RHSFunction_Newtonian()
238 for (int j = 0; j < 5; j++) v[j][i] -= wdetJ * idl_residual[j]; in RHSFunction_Newtonian()
318 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in IFunction_Newtonian() local
345 for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * v_i[j]; in IFunction_Newtonian()
427 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in IJacobian_Newtonian() local
453 for (int j = 0; j < 5; j++) v[j][i] += wdetJ * v_i[j]; in IJacobian_Newtonian()
456 for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * v_i[j]; in IJacobian_Newtonian()
[all …]
H A Dbc_outflow.h42 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in RiemannOutflow() local
79 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in RiemannOutflow()
111 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in RiemannOutflow_Jacobian() local
161 for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in RiemannOutflow_Jacobian()
190 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in PressureOutflow() local
220 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in PressureOutflow()
252 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in PressureOutflow_Jacobian() local
289 for (int j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in PressureOutflow_Jacobian()
H A Dmonitor_cfl.h13 CeedScalar(*v) = out[0]; in MonitorCFL()
30 v[i] = CalculateCFL_2D(s.Y.velocity, 1, gijd_mat); in MonitorCFL()
40 v[i] = CalculateCFL_3D(s.Y.velocity, 1, gijd_mat); in MonitorCFL()
H A Dblasius.h129 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Blasius_Inflow() local
171 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in Blasius_Inflow()
181 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Blasius_Inflow_Jacobian() local
221 v[0][i] = -wdetJb * drho * u_normal; in Blasius_Inflow_Jacobian()
223 v[j + 1][i] = -wdetJb * (drho * u_normal * s.Y.velocity[j] + normal[j] * dP); in Blasius_Inflow_Jacobian()
225 v[4][i] = -wdetJb * u_normal * (dE + dP); in Blasius_Inflow_Jacobian()
H A Dbc_slip.h15 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Slip() local
38 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in Slip()
62 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Slip_Jacobian() local
94 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in Slip_Jacobian()
H A Dbc_freestream.h19 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Freestream() local
44 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j]; in Freestream()
81 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Freestream_Jacobian() local
110 for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j]; in Freestream_Jacobian()
H A Dvelocity_gradient_projection.h14 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in VelocityGradientProjectionRHS() local
32 v[j * 3 + k][i] = wdetJ * grad_velocity[j][k]; in VelocityGradientProjectionRHS()
H A Dstg_shur14.h337 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in StgShur14Inflow() local
394 for (CeedInt j = 0; j < 5; j++) v[j][i] = 0.; in StgShur14Inflow()
398 v[0][i] -= wdetJb * rho * u_normal; in StgShur14Inflow()
401 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + normal[j] * P); in StgShur14Inflow()
404 v[4][i] -= wdetJb * u_normal * (E + P); in StgShur14Inflow()
417 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in StgShur14Inflow_Jacobian() local
460 v[0][i] = -wdetJb * drho * u_normal; in StgShur14Inflow_Jacobian()
461 …for (int j = 0; j < 3; j++) v[j + 1][i] = -wdetJb * (drho * u_normal * velocity[j] + normal[j] * d… in StgShur14Inflow_Jacobian()
462 v[4][i] = -wdetJb * u_normal * (dE + dP); in StgShur14Inflow_Jacobian()
H A Dmass.h22 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Mass_N() local
24 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) CeedPragmaSIMD for (CeedInt j = 0; j < N; j++) v[j]… in Mass_N()
H A Ddifferential_filter.h80 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in DifferentialFilter_LHS_N() local
126 v[j][i] = wdetJ * q[j][i]; in DifferentialFilter_LHS_N()
164 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in DifferentialFilter_MMS_RHS() local
168 v[0][i] = wdetJ * q[0][i]; in DifferentialFilter_MMS_RHS()
/honee/doc/
H A Dtheory.md85 …e first multiply the strong form {eq}`eq-vector-ns` by a test function $\bm v \in H^1(\Omega)$ and…
88 …} \bm v \cdot \left(\frac{\partial \bm{q}_N}{\partial t} + \nabla \cdot \bm{F}(\bm{q}_N) - \bm{S}(…
91 with $\mathcal{V}_p = \{ \bm v(\bm x) \in H^{1}(\Omega_e) \,|\, \bm v(\bm x_e(\bm X)) \in P_p(\bm{I…
97 \int_{\Omega} \bm v \cdot \left( \frac{\partial \bm{q}_N}{\partial t} - \bm{S}(\bm{q}_N) \right) \…
98 - \int_{\Omega} \nabla \bm v \!:\! \bm{F}(\bm{q}_N)\,dV & \\
99 + \int_{\partial \Omega} \bm v \cdot \bm{F}(\bm q_N) \cdot \widehat{\bm{n}} \,dS
100 &= 0 \, , \; \forall \bm v \in \mathcal{V}_p \,,
107v \!:\! \bm F$ represents contraction over both fields and spatial dimensions while a single dot r…
183 …\int_{\Omega} \bm v \cdot \left( \frac{\partial \bm{q}_N}{\partial t} - \bm{S}(\bm{q}_N) \right) …
184 - \int_{\Omega} \nabla \bm v \!:\! \bm{F}(\bm{q}_N)\,dV & \\
[all …]
/honee/src/
H A Dmisc.c229 Vec u, v; in HoneeCalculateDomainSize() local
245 PetscCall(DMCreateGlobalVector(dm_coord, &v)); in HoneeCalculateDomainSize()
247 PetscCall(ApplyCeedOperatorLocalToGlobal(u, v, op_mass_ctx)); in HoneeCalculateDomainSize()
248 PetscCall(VecSum(v, volume)); in HoneeCalculateDomainSize()
252 PetscCall(VecDestroy(&v)); in HoneeCalculateDomainSize()
/honee/tests/
H A Djunit_common.py43 default = [self.enum_type(v.lower()) for v in default]
55 values = [self.enum_type(v.lower()) for v in values]
370 for k, v in other.items():
545 for k, v in spec.key_values.items():
/honee/
H A DMakefile385 $(foreach v,$(CONFIG_VARS),$(if $(call needs_save,$(v)),$(file >> $(CONFIG),$(v) := $($(v)))))

12