Home
last modified time | relevance | path

Searched full:dim (Results 1 – 25 of 38) sorted by relevance

12

/honee/qfunctions/
H A Dadvection.h64 CEED_QFUNCTION_HELPER int Exact_AdvectionGeneric(CeedInt dim, CeedScalar time, const CeedScalar X[]… in Exact_AdvectionGeneric() argument
69 const CeedScalar lz = dim == 2 ? 0. : context->lz; in Exact_AdvectionGeneric()
73 const CeedScalar theta = dim == 2 ? M_PI / 3 : M_PI; in Exact_AdvectionGeneric()
76 const CeedScalar x = X[0], y = X[1], z = dim == 2 ? 0. : X[2]; in Exact_AdvectionGeneric()
89 q[3] = dim == 2 ? 0. : wind[2]; in Exact_AdvectionGeneric()
147 CeedScalar theta = context->wave_frequency * DotN(X, wind, dim) + context->wave_phase; in Exact_AdvectionGeneric()
240 …LPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXdx, CeedInt dim) { in Tau() argument
245 MatVecNM(dXdx, s.Y.velocity, dim, dim, CEED_NOTRANSPOSE, uX); in Tau()
246 return context->CtauS / sqrt(DotN(uX, uX, dim)); in Tau()
251 MatMatN(dXdx, dXdx, dim, CEED_TRANSPOSE, CEED_NOTRANSPOSE, gijd_mat); in Tau()
[all …]
H A Ddiff_flux_projection.h9 …neric(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out, const CeedInt dim, in ComputeDivDiffusiveFluxGeneric() argument
18 QdataUnpack_ND(dim, Q, i, q_data, NULL, dXdx); in ComputeDivDiffusiveFluxGeneric()
22 … // Get gradient into dim x dim matrix form, with orientation [flux_direction][gradient_direction] in ComputeDivDiffusiveFluxGeneric()
24 const CeedInt offset = Q * n * dim; // offset to reach nth component flux gradients in ComputeDivDiffusiveFluxGeneric()
25 for (CeedInt g = 0; g < dim; g++) { in ComputeDivDiffusiveFluxGeneric()
26 for (CeedInt f = 0; f < dim; f++) { in ComputeDivDiffusiveFluxGeneric()
27 grad_qn[f * dim + g] = grad_q[offset + (Q * num_comps * dim) * g + Q * f + i]; in ComputeDivDiffusiveFluxGeneric()
31 DivergenceND(grad_qn, dXdx, dim, &v[n][i]); in ComputeDivDiffusiveFluxGeneric()
H A Dutils.h338 @param[in] dim Dimension of the element
343 @param[out] dXdx Inverse of the mapping Jacobian (shape [dim][dim]), or `NULL`
347 CEED_QFUNCTION_HELPER int QdataUnpack_ND(CeedInt dim, CeedInt Q, CeedInt i, const CeedScalar *q_dat… in QdataUnpack_ND() argument
348 switch (dim) { in QdataUnpack_ND()
364 @param[in] dim Dimension of the element
369 @param[out] dXdx Inverse of the mapping Jacobian (shape [dim - 1][dim]), or `NULL`
370 @param[out] normal Components of the normal vector (shape [dim]), or `NULL`
374 CEED_QFUNCTION_HELPER int QdataBoundaryUnpack_ND(CeedInt dim, CeedInt Q, CeedInt i, const CeedScala… in QdataBoundaryUnpack_ND() argument
376 switch (dim) { in QdataBoundaryUnpack_ND()
393 @param[in] dim Dimension of the element
[all …]
H A Dmonitor_cfl.h9 CeedInt dim) { in MonitorCFL() argument
21 switch (dim) { in MonitorCFL()
28 ScaleN((CeedScalar *)gijd_mat, 0.25, Square(dim)); in MonitorCFL()
38 ScaleN((CeedScalar *)gijd_mat, 0.25, Square(dim)); in MonitorCFL()
H A Dadvection_types.h22 ADVDIF_IC_BUBBLE_SPHERE = 0, // dim=3
23 ADVDIF_IC_BUBBLE_CYLINDER = 1, // dim=2
24 ADVDIF_IC_COSINE_HILL = 2, // dim=2
/honee/src/
H A Dqdata.c95 PetscInt num_comp_x, dim; in QDataGetNumComponents() local
98 PetscCall(DMGetDimension(dm, &dim)); in QDataGetNumComponents()
101 switch (dim) { in QDataGetNumComponents()
117 …DM of dimension %" PetscInt_FMT " and coordinates with dimension %" PetscInt_FMT, dim, num_comp_x); in QDataGetNumComponents()
138 PetscInt dim, height = 0, num_comp_x; in QDataGet() local
145 PetscCall(DMGetDimension(dm, &dim)); in QDataGet()
146 switch (dim) { in QDataGet()
161 …(PetscObject)dm), PETSC_ERR_SUP, "%s not valid for DM of dimension %" PetscInt_FMT, __func__, dim); in QDataGet()
165 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_setup, "dx", num_comp_x * (dim - height), CEED_EVAL_G… in QDataGet()
198 PetscInt dim; in QDataBoundaryGetNumComponents() local
[all …]
H A Ddm-utils.c212 PetscInt num_dual_basis_vectors, dim, num_comp; in DMPlexCeedElemRestrictionStridedCreate() local
214 PetscCall(PetscFEGetSpatialDimension(fe, &dim)); in DMPlexCeedElemRestrictionStridedCreate()
285 @param[in] dim Spatial dimension
294 static PetscErrorCode PetscDTUniformTensorQuadrature(PetscInt dim, PetscInt num_comp, PetscInt num_… in PetscDTUniformTensorQuadrature() argument
297 …PetscInt num_total_points = dim > 1 ? (dim > 2 ? (num_points * num_points * num_points) : (n… in PetscDTUniformTensorQuadrature()
301 PetscCall(PetscMalloc1(num_total_points * dim, &coords)); in PetscDTUniformTensorQuadrature()
304 switch (dim) { in PetscDTUniformTensorQuadrature()
332 coords[(i * num_points + j) * dim + 0] = coords_1d[i]; in PetscDTUniformTensorQuadrature()
333 coords[(i * num_points + j) * dim + 1] = coords_1d[j]; in PetscDTUniformTensorQuadrature()
348 coords[((i * num_points + j) * num_points + k) * dim + 0] = coords_1d[i]; in PetscDTUniformTensorQuadrature()
[all …]
H A Dboundary_condition.c128 PetscInt dim, dim_sur, height = 1, num_comp_x, num_comp_q; in HoneeBCCreateIFunctionQF() local
138 PetscCall(DMGetDimension(dm, &dim)); in HoneeBCCreateIFunctionQF()
139 dim_sur = dim - height; in HoneeBCCreateIFunctionQF()
176 PetscInt dim, dim_sur, height = 1, num_comp_x, num_comp_q; in HoneeBCCreateIJacobianQF() local
186 PetscCall(DMGetDimension(dm, &dim)); in HoneeBCCreateIJacobianQF()
187 dim_sur = dim - height; in HoneeBCCreateIJacobianQF()
224 PetscInt dim, height = 1, num_comp_x, num_comp_q; in HoneeBCAddIFunctionOp() local
237 PetscCall(DMGetDimension(dm, &dim)); in HoneeBCAddIFunctionOp()
293 PetscInt dim, height = 1, num_comp_x, num_comp_q; in HoneeBCAddIJacobianOp() local
306 PetscCall(DMGetDimension(dm, &dim)); in HoneeBCAddIJacobianOp()
H A Dsetuplibceed.c94 PetscInt dim, num_comp_x; in SetupLibceed() local
102 PetscCall(DMGetDimension(dm, &dim)); in SetupLibceed()
147 PetscCallCeed(ceed, CeedQFunctionAddInput(qf_ics, "dx", num_comp_x * dim, CEED_EVAL_GRAD)); in SetupLibceed()
169 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_vol, "Grad_q", num_comp_q * dim, CEED_EVAL_GRAD)); in SetupLibceed()
175 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_vol, "Grad_v", num_comp_q * dim, CEED_EVAL_GRAD)… in SetupLibceed()
198 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_ifunction_vol, "Grad_q", num_comp_q * dim, CEED_EVAL_… in SetupLibceed()
205 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_ifunction_vol, "Grad_v", num_comp_q * dim, CEED_EVAL… in SetupLibceed()
231 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_ijacobian_vol, "Grad_dq", num_comp_q * dim, CEED_EVAL… in SetupLibceed()
235 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_ijacobian_vol, "Grad_v", num_comp_q * dim, CEED_EVAL… in SetupLibceed()
H A Dsetupts.c75 PetscInt dim = 3; in Surface_Forces_NS() local
102 for (PetscInt j = 0; j < dim; j++) { in Surface_Forces_NS()
103 reaction_force[w * dim + j] -= r[node].momentum[j]; in Surface_Forces_NS()
110 …PetscCallMPI(MPIU_Allreduce(MPI_IN_PLACE, reaction_force, dim * num_walls, MPIU_SCALAR, MPI_SUM, c… in Surface_Forces_NS()
257 PetscInt num_wall = honee->app_ctx->wall_forces.num_wall, dim = 3; in TSMonitor_WallForce() local
267 PetscCall(PetscCalloc1(num_wall * dim, &reaction_force)); in TSMonitor_WallForce()
282 … reaction_force[w * dim + 0], reaction_force[w * dim + 1], reaction_force[w * dim + 2])); in TSMonitor_WallForce()
286 … reaction_force[w * dim + 0], reaction_force[w * dim + 1], reaction_force[w * dim + 2])); in TSMonitor_WallForce()
H A Dhonee-meshtransform.c61 PetscInt narr, ncoords, dim; in HoneeMeshTransform_PlateMesh() local
67 PetscCall(DMGetDimension(dm, &dim)); in HoneeMeshTransform_PlateMesh()
78 PetscScalar(*coords)[dim] = (PetscScalar(*)[dim])arr_coords; in HoneeMeshTransform_PlateMesh()
79 ncoords = narr / dim; in HoneeMeshTransform_PlateMesh()
H A Dmonitor_cfl.c35 PetscInt num_comp_cfl = 1, dim, tab_level; in SetupMontiorCfl() local
44 PetscCall(DMGetDimension(honee->dm, &dim)); in SetupMontiorCfl()
63 switch (dim) { in SetupMontiorCfl()
92 …"Could not create CFL monitor QFunction for dim %" PetscInt_FMT " and state variable %s", dim, Sta… in SetupMontiorCfl()
H A Ddiff_flux_projection.c24 …PetscInt height = 0, dm_field = 0, dim, degree = honee->app_ctx->degree, q_extra … in DivDiffFluxProjectionCreate() local
43 PetscCall(DMGetDimension(projection->dm, &dim)); in DivDiffFluxProjectionCreate()
77 projection->num_comp = diff_flux_proj_->num_diff_flux_comps * dim; in DivDiffFluxProjectionCreate()
257 PetscInt dim; in DivDiffFluxProjectionSetup_Indirect() local
259 PetscCall(DMGetDimension(projection->dm, &dim)); in DivDiffFluxProjectionSetup_Indirect()
262 switch (dim) { in DivDiffFluxProjectionSetup_Indirect()
288 dim, diff_flux_proj->num_diff_flux_comps); in DivDiffFluxProjectionSetup_Indirect()
290 …dQFunctionAddInput(qf_calc_divergence, "Grad F_diff", projection->num_comp * dim, CEED_EVAL_GRAD)); in DivDiffFluxProjectionSetup_Indirect()
H A Dvelocity_gradient_projection.c45 PetscInt dim, height = 0, dm_field = 0, num_comp_input; in VelocityGradientProjectionSetup() local
53 PetscCall(DMGetDimension(grad_velo_proj->dm, &dim)); in VelocityGradientProjectionSetup()
87 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_assemble, "Grad_q", num_comp_input * dim, CEED_EV… in VelocityGradientProjectionSetup()
H A Dstrong_boundary_conditions.c21 PetscInt dim; in SetupStrongSTG_Ceed() local
26 PetscCall(DMGetDimension(dm, &dim)); in SetupStrongSTG_Ceed()
27 num_comp_x = dim; in SetupStrongSTG_Ceed()
28 dXdx_size = num_comp_x * (dim - height_cell); in SetupStrongSTG_Ceed()
/honee/problems/
H A Dadvection.c26 PetscInt dim; in PRINT_ADVECTION() local
29 PetscCall(DMGetDimension(honee->dm, &dim)); in PRINT_ADVECTION()
43 switch (dim) { in PRINT_ADVECTION()
84 PetscInt dim; in CreateKSPMassOperator_AdvectionStabilized() local
87 PetscCall(DMGetDimension(honee->dm, &dim)); in CreateKSPMassOperator_AdvectionStabilized()
105 switch (dim) { in CreateKSPMassOperator_AdvectionStabilized()
120 PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_mass, "Grad_v", 5 * dim, CEED_EVAL_GRAD)); in CreateKSPMassOperator_AdvectionStabilized()
149 PetscInt dim, num_comp_q; in DivDiffFluxProjectionCreateRHS_Direct_AdvDif() local
154 PetscCall(DMGetDimension(projection->dm, &dim)); in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
178 switch (dim) { in DivDiffFluxProjectionCreateRHS_Direct_AdvDif()
[all …]
H A Ddensitycurrent.c21 PetscInt dim; in NS_DENSITY_CURRENT() local
38 PetscCall(DMGetDimension(dm, &dim)); in NS_DENSITY_CURRENT()
39 for (PetscInt i = 0; i < dim; i++) domain_size[i] = domain_max[i] - domain_min[i]; in NS_DENSITY_CURRENT()
50 for (PetscInt i = 0; i < dim; i++) center[i] = .5 * domain_size[i]; in NS_DENSITY_CURRENT()
51 PetscInt n = dim; in NS_DENSITY_CURRENT()
53 n = dim; in NS_DENSITY_CURRENT()
61 for (PetscInt i = 0; i < dim; i++) dc_axis[i] /= norm; in NS_DENSITY_CURRENT()
74 for (PetscInt i = 0; i < dim; i++) center[i] *= units->meter; in NS_DENSITY_CURRENT()
H A Dchannel.c150 PetscInt narr, ncoords, dim; in DivDiffFluxVerifyMesh() local
156 PetscCall(DMGetDimension(dm, &dim)); in DivDiffFluxVerifyMesh()
166 PetscScalar(*coords)[dim] = (PetscScalar(*)[dim])arr_coords; in DivDiffFluxVerifyMesh()
167 ncoords = narr / dim; in DivDiffFluxVerifyMesh()
H A Deulervortex.c63 PetscInt dim; in NS_EULER_VORTEX() local
93 PetscCall(DMGetDimension(dm, &dim)); in NS_EULER_VORTEX()
94 for (PetscInt i = 0; i < dim; i++) domain_size[i] = domain_max[i] - domain_min[i]; in NS_EULER_VORTEX()
102 PetscInt n = dim; in NS_EULER_VORTEX()
105 …for (PetscInt i = 0; i < dim; i++) center[i] = .5 * domain_size[i] / units->meter; // Redimension… in NS_EULER_VORTEX()
106 n = dim; in NS_EULER_VORTEX()
H A Dnewtonian.c46 PetscInt dim = 3; in CreateKSPMassOperator_NewtonianStabilized() local
74 PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_mass, "Grad_v", 5 * dim, CEED_EVAL_GRAD)); in CreateKSPMassOperator_NewtonianStabilized()
103 PetscInt dim, num_comp_q; in DivDiffFluxProjectionCreateRHS_Direct_NS() local
110 PetscCall(DMGetDimension(projection->dm, &dim)); in DivDiffFluxProjectionCreateRHS_Direct_NS()
150 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_volume, "Grad_q", num_comp_q * dim, CEED_EVAL_GRA… in DivDiffFluxProjectionCreateRHS_Direct_NS()
152 …unctionAddOutput(qf_rhs_volume, "diffusive flux RHS", projection->num_comp * dim, CEED_EVAL_GRAD)); in DivDiffFluxProjectionCreateRHS_Direct_NS()
197 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs_boundary, "Grad_q", num_comp_q * dim, CEED_EVAL_G… in DivDiffFluxProjectionCreateRHS_Direct_NS()
271 PetscInt dim, num_comp_q; in DivDiffFluxProjectionCreateRHS_Indirect_NS() local
277 PetscCall(DMGetDimension(projection->dm, &dim)); in DivDiffFluxProjectionCreateRHS_Indirect_NS()
307 PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs, "Grad_q", num_comp_q * dim, CEED_EVAL_GRAD)); in DivDiffFluxProjectionCreateRHS_Indirect_NS()
[all …]
/honee/qfunctions/spanstats/
H A Dcflpe.h19 StateVariable state_var, CeedInt dim) { in ChildStatsCollection_CflPe() argument
33 switch (dim) { in ChildStatsCollection_CflPe()
42 ScaleN((CeedScalar *)gijd_mat, 0.25, Square(dim)); in ChildStatsCollection_CflPe()
55 ScaleN((CeedScalar *)gijd_mat, 0.25, Square(dim)); in ChildStatsCollection_CflPe()
/honee/src/spanstats/
H A Dcflpe.c21 PetscInt dim, num_comp_q, num_comp_x; in CreateStatisticCollectionOperator() local
31 PetscCall(DMGetDimension(honee->dm, &dim)); in CreateStatisticCollectionOperator()
36 switch (dim) { in CreateStatisticCollectionOperator()
71 …ot create CFL/Pe spanstats collection QFunction for dim %" PetscInt_FMT " and state variable %s", … in CreateStatisticCollectionOperator()
/honee/examples/
H A Dadvection.yaml6 dim: 3
H A Dadvection_bl.yaml26 dim: 2
/honee/tests/
H A Dchannel_divdiff_verify.yaml30 dim: 3

12