| /honee/qfunctions/spanstats/ |
| H A D | turbulence.h | 44 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() 63 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_X][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection() 64 …v[TURB_MEAN_DENSITY_TEMPERATURE_FLUX_Y][i] = wdetJ * s.U.density * s.Y.temperature * s.Y.velocity[… in ChildStatsCollection() [all …]
|
| H A D | cflpe.h | 22 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 D | utils.h | 43 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 D | sgs_dd_training.h | 30 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() 52 …v[0 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XX][i] - Square(s.Y.velocity[0])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal() 53 …v[1 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YY][i] - Square(s.Y.velocity[1])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal() 54 …v[2 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] - Square(s.Y.velocity[2])) * inv_mult… in ComputeSGS_DDAnisotropicTrainingDataNodal() 55 …v[3 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] - s.Y.velocity[1] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal() 56 …v[4 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] - s.Y.velocity[0] * s.Y.velocity[2]) … in ComputeSGS_DDAnisotropicTrainingDataNodal() 57 …v[5 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XY][i] - s.Y.velocity[0] * s.Y.velocity[1]) … in ComputeSGS_DDAnisotropicTrainingDataNodal()
|
| H A D | channel.h | 82 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 D | monitor_totalkineticenergy.h | 13 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 D | grid_anisotropy_tensor.h | 36 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 D | utils_eigensolver_jacobi.h | 216 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 D | advection.h | 275 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 D | eulervortex.h | 267 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 D | diff_flux_projection.h | 13 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 D | newtonian.h | 47 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 D | bc_outflow.h | 42 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 D | monitor_cfl.h | 13 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 D | blasius.h | 129 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 D | bc_slip.h | 15 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 D | bc_freestream.h | 19 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 D | velocity_gradient_projection.h | 14 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 D | stg_shur14.h | 337 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 D | mass.h | 22 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 D | differential_filter.h | 80 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 D | theory.md | 85 …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 \,, 107 …v \!:\! \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 D | misc.c | 229 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 D | junit_common.py | 43 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 D | Makefile | 385 $(foreach v,$(CONFIG_VARS),$(if $(call needs_save,$(v)),$(file >> $(CONFIG),$(v) := $($(v)))))
|