Home
last modified time | relevance | path

Searched full:s (Results 1 – 25 of 85) sorted by relevance

1234

/honee/qfunctions/
H A Dnewtonian_state.h24 CEED_QFUNCTION_HELPER void UnpackState_U(StateConservative s, CeedScalar U[5]) { in UnpackState_U() argument
25 U[0] = s.density; in UnpackState_U()
26 for (int i = 0; i < 3; i++) U[i + 1] = s.momentum[i]; in UnpackState_U()
27 U[4] = s.E_total; in UnpackState_U()
30 CEED_QFUNCTION_HELPER void UnpackState_Y(StatePrimitive s, CeedScalar Y[5]) { in UnpackState_Y() argument
31 Y[0] = s.pressure; in UnpackState_Y()
32 for (int i = 0; i < 3; i++) Y[i + 1] = s.velocity[i]; in UnpackState_Y()
33 Y[4] = s.temperature; in UnpackState_Y()
36 CEED_QFUNCTION_HELPER void UnpackState_V(StateEntropy s, CeedScalar V[5]) { in UnpackState_V() argument
37 V[0] = s.S_density; in UnpackState_V()
[all …]
H A Dsgs_dd_training.h6 …ess modeling in the strain-rate eigenframe for large eddy simulation' 2022 and 'S-frame discrepancy
44 const State s = StateFromQ(gas, qi, state_var); in ComputeSGS_DDAnisotropicTrainingDataNodal() local
48 …ComputeSgsDDInputs(grad_velo_aniso, km_A_ij, delta, gas.mu / s.U.density, eigenvectors, inputs, &g… 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 Dnewtonian.h24 State s = StateFromPrimitive(gas, context->reference); in ICsNewtonianIG() local
25 StateToQ(gas, s, q, state_var); in ICsNewtonianIG()
56 const State s = StateFromQ(gas, qi, state_var); in MassFunction_Newtonian() local
70 Tau_diagPrim(context->tau_coeffs, gas, s, dXdx, context->dt, Tau_d); in MassFunction_Newtonian()
71 Stabilization(context->stabilization, gas, s, Tau_d, grad_s, U_dot, body_force, divFdiff, stab); in MassFunction_Newtonian()
88 … InternalDampingLayer_Residual(const NewtonianIGProperties gas, const State s, const CeedScalar si… in InternalDampingLayer_Residual() argument
91 State damp_s = StateFromY_fwd(gas, s, damp_Y); in InternalDampingLayer_Residual()
104 @param[in] s Solution `State`
114 CEED_QFUNCTION_HELPER void InternalDampingLayer_IFunction_Integrand(const State s, const NewtonianI… in InternalDampingLayer_IFunction_Integrand() argument
118 CeedScalar damp_state[5] = {s.Y.pressure - pressure, 0, 0, 0, 0}, idl_residual[5] = {0.}; in InternalDampingLayer_IFunction_Integrand()
[all …]
H A Dbc_slip.h22 State s = StateFromQ(gas, qi, state_var); in Slip() local
29 const CeedScalar vel_normal = Dot3(s.Y.velocity, normal); in Slip()
30 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip()
31 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y… in Slip()
34 StateConservative flux = RiemannFlux_HLLC(gas, s, s_reflect, normal); in Slip()
75 State s = StateFromQ(gas, qi, state_var); in Slip_Jacobian() local
76 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in Slip_Jacobian()
79 const CeedScalar vel_normal = Dot3(s.Y.velocity, normal); in Slip_Jacobian()
80 for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * normal[j] * vel_normal; in Slip_Jacobian()
81 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y… in Slip_Jacobian()
[all …]
H A Ddifferential_filter.h41 const State s = StateFromQ(gas, qi, state_var); in DifferentialFilter_RHS() local
43 v0[DIFF_FILTER_PRESSURE][i] = wdetJ * s.Y.pressure; in DifferentialFilter_RHS()
44 v0[DIFF_FILTER_VELOCITY_X][i] = wdetJ * s.Y.velocity[0]; in DifferentialFilter_RHS()
45 v0[DIFF_FILTER_VELOCITY_Y][i] = wdetJ * s.Y.velocity[1]; in DifferentialFilter_RHS()
46 v0[DIFF_FILTER_VELOCITY_Z][i] = wdetJ * s.Y.velocity[2]; in DifferentialFilter_RHS()
47 v0[DIFF_FILTER_TEMPERATURE][i] = wdetJ * s.Y.temperature; in DifferentialFilter_RHS()
48 v1[DIFF_FILTER_VELOCITY_SQUARED_XX][i] = wdetJ * s.Y.velocity[0] * s.Y.velocity[0]; in DifferentialFilter_RHS()
49 v1[DIFF_FILTER_VELOCITY_SQUARED_YY][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[1]; in DifferentialFilter_RHS()
50 v1[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] = wdetJ * s.Y.velocity[2] * s.Y.velocity[2]; in DifferentialFilter_RHS()
51 v1[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] = wdetJ * s.Y.velocity[1] * s.Y.velocity[2]; in DifferentialFilter_RHS()
[all …]
H A Dstabilization.h24 CEED_QFUNCTION_HELPER void StabilizationMatrix(const NewtonianIGProperties gas, const State s, cons… in StabilizationMatrix() argument
31 State ds = StateFromY_fwd(gas, s, dY); in StabilizationMatrix()
32 FluxInviscid_fwd(gas, s, ds, dF); in StabilizationMatrix()
40 …n(StabilizationType stab_type, const NewtonianIGProperties gas, const State s, const CeedScalar Ta… in Stabilization() argument
49 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
52 FluxInviscidStrong(gas, s, ds, strong_residual); in Stabilization()
56 StabilizationMatrix(gas, s, Tau_d, strong_residual, stab); in Stabilization()
66 …u_diagPrim(TauDiagCoefficients tau_coeffs, NewtonianIGProperties gas, State s, const CeedScalar dX… in Tau_diagPrim() argument
76 const CeedScalar rho = s.U.density; in Tau_diagPrim()
89 MatVec3(gijd_mat, s.Y.velocity, CEED_NOTRANSPOSE, gij_uj); in Tau_diagPrim()
[all …]
H A Dmonitor_totalkineticenergy.h18 const State s = StateFromQ(gas, qi, state_var); in MonitorTotalKineticEnergy() local
23 StatePhysicalGradientFromReference(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in MonitorTotalKineticEnergy()
25 v[0][i] = wdetJ * 0.5 * s.U.density * Dot3(s.Y.velocity, s.Y.velocity); in MonitorTotalKineticEnergy()
32 v[3][i] = wdetJ * s.Y.pressure * div_u; in MonitorTotalKineticEnergy()
H A Dblasius.h114 State s = BlasiusSolution(context, x, x0, x_inflow, S_infty.U.density, &t12); in ICsBlasius() local
117 StateToQ(gas, s, q, context->newt_ctx.state_var); in ICsBlasius()
146 State s = BlasiusSolution(context, x, x0, context->x_inflow, rho_0, &t12); in Blasius_Inflow() local
153 s.U.density = s_int.U.density; in Blasius_Inflow()
154 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
156 s.U.E_total = s_int.U.E_total; in Blasius_Inflow()
157 s.Y = StatePrimitiveFromConservative(gas, s.U); in Blasius_Inflow()
161 FluxInviscid(gas, s, Flux_inviscid); in Blasius_Inflow()
201 State s = BlasiusSolution(context, x, x0, 0, rho_0, &t12); in Blasius_Inflow_Jacobian() local
209 CeedScalar dE_kinetic = .5 * drho * Dot3(s.Y.velocity, s.Y.velocity); in Blasius_Inflow_Jacobian()
[all …]
H A Dbc_freestream.h27 const State s = StateFromQ(gas, qi, state_var); in Freestream() local
36 flux = RiemannFlux_HLL(gas, s, context->S_infty, normal); in Freestream()
39 flux = RiemannFlux_HLLC(gas, s, context->S_infty, normal); in Freestream()
96 State s = StateFromQ(gas, qi, state_var); in Freestream_Jacobian() local
97 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in Freestream_Jacobian()
102 dflux = RiemannFlux_HLL_fwd(gas, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
105 dflux = RiemannFlux_HLLC_fwd(gas, s, ds, context->S_infty, dS_infty, normal); in Freestream_Jacobian()
H A Dadvection.h210 …ysicalGradientFromReference_ND(CeedInt N, CeedInt Q, CeedInt i, NewtonianIGProperties gas, State s, in StatePhysicalGradientFromReference_ND() argument
220 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
233 grad_s[k] = StateFromQ_fwd(gas, s, dqi, state_var); in StatePhysicalGradientFromReference_ND()
240 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd… in Tau() argument
245 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
254 MatVecNM(gijd_mat, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, gij_uj); in Tau()
255 …return 1 / sqrt(Square(2 * context->Ctau_t / context->dt) + DotN(s.Y.velocity, gij_uj, dim) * Squa… in Tau()
282 const State s = StateFromU(gas, qi); in IFunction_AdvectionGeneric() local
287 …StatePhysicalGradientFromReference_ND(dim, Q, i, gas, s, STATEVAR_CONSERVATIVE, grad_q, dXdx, grad… in IFunction_AdvectionGeneric()
303 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in IFunction_AdvectionGeneric()
[all …]
H A Dutils_eigensolver_jacobi.h61 /// A[i][j] = 0. The results will be stored in c, s, and t
90 /// angle θ. This function assumes that c=cos(θ). s=sin(θ), t=tan(θ)
104 /// | c ... s |
108 /// | -s ... c |
160 // A[i][i] = c*c*A[i][i] + s*s*A[j][j] - 2*s*c*A[i][j] in ApplyRot()
161 // A[j][j] = s*s*A[i][i] + c*c*A[j][j] + 2*s*c*A[i][j] in ApplyRot()
206 /// the arguments "s" and "c" refer to cos(θ) and sin(θ), respectively).
215 // Recall that c = cos(θ) and s = sin(θ) in ApplyRotLeft()
228 /// @param[in] s sort decreasing order?
H A Dbc_outflow.h198 State s = StateFromQ(gas, qi, state_var); in PressureOutflow() local
199 s.Y.pressure = outflow->pressure; in PressureOutflow()
206 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_q, dXdx, grad_s); in PressureOutflow()
212 ViscousEnergyFlux(gas, s.Y, grad_s, stress, Fe); in PressureOutflow()
215 FluxInviscid(gas, s, F_inviscid); in PressureOutflow()
268 State s = StateFromQ(gas, qi, state_var); in PressureOutflow_Jacobian() local
269 State ds = StateFromQ_fwd(gas, s, dqi, state_var); in PressureOutflow_Jacobian()
270 s.Y.pressure = outflow->pressure; in PressureOutflow_Jacobian()
274 StatePhysicalGradientFromReference_Boundary(Q, i, gas, s, state_var, Grad_dq, dXdx, grad_ds); in PressureOutflow_Jacobian()
281 ViscousEnergyFlux_fwd(gas, s.Y, ds.Y, grad_ds, stress, dstress, dFe); in PressureOutflow_Jacobian()
[all …]
H A Dmonitor_cfl.h19 const State s = StateFromQ(gas, qi, state_var); in MonitorCFL() local
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()
/honee/qfunctions/spanstats/
H A Dturbulence.h54 const State s = StateFromQ(gas, qi, state_var); 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 Dcflpe.h30 const State s = StateFromQ(gas, qi, state_var); in ChildStatsCollection_CflPe() local
44 cfl = CalculateCFL_2D(s.Y.velocity, context->timestep, gijd_mat); in ChildStatsCollection_CflPe()
45 Pe = CalculatePe_2D(s.Y.velocity, context->diffusion_coeff, gijd_mat); in ChildStatsCollection_CflPe()
57 cfl = CalculateCFL_3D(s.Y.velocity, context->timestep, gijd_mat); in ChildStatsCollection_CflPe()
58 Pe = CalculatePe_3D(s.Y.velocity, context->diffusion_coeff, gijd_mat); in ChildStatsCollection_CflPe()
/honee/
H A Dcommon.mk8 ansicolor = $(shell echo $(call last_path,$(1)) | cksum | cut -b1-2 | xargs -IS expr 2 \* S + 17)
9 emacs_out = @printf " %10s %s/%s\n" $(1) $(call rule_path,$(2)) $(call rule_file,$(2))
11 printf " %10s \033[38;5;%d;1m%s\033[m/%s\n" \
14 printf " %10s %s\n" $(1) $(2); fi
/honee/src/
H A Dmonitor_totalkineticenergy.c172 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# honee_git_commit: %s\n", buf_const)); in TSMonitor_TotalKineticEnergy()
174 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# petsc_version: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
176 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# petsc_directory: %s\n", buf_const)); in TSMonitor_TotalKineticEnergy()
178 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# petsc_arch: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
180 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# date: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
182 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# date_timezone: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
184 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# username: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
186 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# hostname: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
188 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# working_directory: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
191 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# original_file_path: %s\n", buf)); in TSMonitor_TotalKineticEnergy()
[all …]
H A Dmonitor_cfl.c92 … create CFL monitor QFunction for dim %" PetscInt_FMT " and state variable %s", dim, StateVariable… in SetupMontiorCfl()
171 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# petsc_version: %s\n", buf)); in TSMonitor_Cfl()
173 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# petsc_directory: %s\n", buf_const)); in TSMonitor_Cfl()
175 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# petsc_arch: %s\n", buf)); in TSMonitor_Cfl()
177 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# date: %s\n", buf)); in TSMonitor_Cfl()
179 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# date_timezone: %s\n", buf)); in TSMonitor_Cfl()
181 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# username: %s\n", buf)); in TSMonitor_Cfl()
183 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# hostname: %s\n", buf)); in TSMonitor_Cfl()
185 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# working_directory: %s\n", buf)); in TSMonitor_Cfl()
189 PetscCall(PetscViewerASCIIPrintf(ctx->viewer, "# original_file_path: %s\n", buf)); in TSMonitor_Cfl()
H A Dhonee-file.c12 Developer's note: Could instead use PetscStrendswith
25 PetscCheck(ext_len, comm, PETSC_ERR_ARG_WRONG, "Zero-size extension: %s", extension); in HoneeCheckFilenameExtension()
35 // @brief Read in binary int based on it's data type
104 …"CGNS file %s does not have FlowIterations array and solution name have 'FlowSolution' (found '%s'… in HoneeLoadInitialCondition()
112 …(PETSC_COMM_WORLD, PETSC_ERR_SUP, "File does not have a valid extension, recieved '%s'", filename); in HoneeLoadInitialCondition()
158 …comm, PETSC_ERR_FILE_UNEXPECTED, "Not a HONEE header token or a PETSc Vec in file '%s'", filename); in HoneeLoadBinaryVec()
162 …PetscCheck(length == N, comm, PETSC_ERR_ARG_INCOMP, "File Vec '%s' has length %" PetscInt_FMT " bu… in HoneeLoadBinaryVec()
195 …viewer, comm, PETSC_ERR_ARG_WRONGSTATE, "Viewer must be binary type; instead got %s", viewer_type); in HoneeWriteBinaryVec()
197 …e == FILE_MODE_WRITE, comm, PETSC_ERR_ARG_WRONGSTATE, "Viewer must be binary type; instead got %s", in HoneeWriteBinaryVec()
231 …ETSC_ERR_FILE_UNEXPECTED, "Found %d dimensions instead of 2 on the first line of %s", ndims, path); in PhastaDatFileOpen()
[all …]
H A Dmisc.c164 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Time integrator %s on time step %" PetscInt_FMT " with fi… in PostProcess()
278 " Host Name : %s\n" in PrintRunInfo()
292 " libCEED Backend : %s\n" in PrintRunInfo()
293 " libCEED Backend MemType : %s\n", in PrintRunInfo()
307 " Box Faces : %s\n" in PrintRunInfo()
308 " DM VecType : %s\n" in PrintRunInfo()
309 " Time Stepping Scheme : %s\n", in PrintRunInfo()
328 …PetscCall(PetscSNPrintf(pmat_coo_type_str, sizeof(pmat_coo_type_str), " (COO MatType: %s)", pmat_c… in PrintRunInfo()
333 " IJacobian A MatType : %s\n" in PrintRunInfo()
334 " IJacobian P MatType : %s\n", in PrintRunInfo()
[all …]
/honee/tests/
H A Dsmartsim_regression_framework.py29 s = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
30 s.bind(('', 0))
31 addr = s.getsockname()
32 s.close()
60 Check whether the S-frame-oriented vorticity vector's second component is just flipped.
207 '%Y-%m-%d %H:%M:%S %Z', time.localtime(start)),
/honee/problems/
H A Dbc_freestream.c181 // @brief Calculate relative error, (A - B) / S
182 // If S < threshold, then set S=1
183 static inline CeedScalar RelativeError(CeedScalar S, CeedScalar A, CeedScalar B, CeedScalar thresho… in RelativeError() argument
184 return (A - B) / (fabs(S) > threshold ? S : 1); in RelativeError()
203 printf("%s[0] error %g (expected %.10e, got %.10e)\n", name, relative_error[0], Q_s[0], Q_a[0]); in CheckQWithTolerance()
207 …printf("%s[%d] error %g (expected %.10e, got %.10e)\n", name, i, relative_error[i], Q_s[i], Q_a[i]… in CheckQWithTolerance()
211 printf("%s[4] error %g (expected %.10e, got %.10e)\n", name, relative_error[4], Q_s[4], Q_a[4]); in CheckQWithTolerance()
408 …if (fabs(ds_left_err) >= rtol) printf("%s ds_left error %g (expected %.10e, got %.10e)\n", buf, ds… in TestComputeHLLSpeeds_Roe_fwd()
409 …if (fabs(ds_right_err) >= rtol) printf("%s ds_right error %g (expected %.10e, got %.10e)\n", buf, … in TestComputeHLLSpeeds_Roe_fwd()
456 …if (fabs(dH_err) >= rtol) printf("%s dH error %g (expected %.10e, got %.10e)\n", buf, dH_err, dH_f… in TestTotalSpecificEnthalpy_fwd()
[all …]
H A Dsgs_model_torch_weak.c11 SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_SUP, "Must build with USE_TORCH set to run %s", __func__); in LoadModel_Torch()
17 SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_SUP, "Must build with USE_TORCH set to run %s", __func__); in ModelInference_Torch()
/honee/doc/
H A Dtheory.md31 \frac{\partial \bm{q}}{\partial t} + \nabla \cdot \bm{F}(\bm{q}) -S(\bm{q}) = 0 \, ,
65 S(\bm{q}) &=
88 …(\frac{\partial \bm{q}_N}{\partial t} + \nabla \cdot \bm{F}(\bm{q}_N) - \bm{S}(\bm{q}_N) \right) \…
97 \int_{\Omega} \bm v \cdot \left( \frac{\partial \bm{q}_N}{\partial t} - \bm{S}(\bm{q}_N) \right) \…
107 …ion in just one, which should be clear from context, e.g., $\bm v \cdot \bm S$ contracts over fiel…
111 <!-- TODO: This should be reframed in terms of PETSc TS's F(t, u, \dot u) = G(t, u) rather than spe…
120 \bm{q}_N^{n+1} = \bm{q}_N^n + \Delta t \sum_{i=1}^{s} b_i k_i \, ,
131 k_i &= f\left(t^n + c_i \Delta t, \bm{q}_N^n + \Delta t \sum_{j=1}^s a_{ij} k_j \right)\\
138 f(t^n, \bm{q}_N^n) = - [\nabla \cdot \bm{F}(\bm{q}_N)]^n + [S(\bm{q}_N)]^n \, .
168 More details of PETSc's time stepping solvers can be found in the [TS User Guide](https://petsc.org…
[all …]
H A Dexamples.md6 …014`. It is intended to test non-reflecting/Riemann boundary conditions. It's primarily intended f…
76 Given the force components $\bm F = (F_x, F_y, F_z)$ and surface area $S = \pi D L_z$ where $L_z$ i…
80 C_L &= \frac{2 F_y}{\rho_\infty u_\infty^2 S} \\
81 C_D &= \frac{2 F_x}{\rho_\infty u_\infty^2 S} \\
171 - `1/s`
207 - `m/s`
318 - `m/s`
410 - `m/s`
436 There is no perturbation in the entropy $S=P/\rho^\gamma$ ($\delta S=0)$.
477 This test problem is based on Sod's Shock Tube (from{cite}`sodshocktubewiki`), a canonical test cas…
[all …]

1234