Lines Matching +full:- +full:f
1 // SPDX-FileCopyrightText: Copyright (c) 2017-2024, HONEE contributors.
2 // SPDX-License-Identifier: Apache-2.0 OR BSD-2-Clause
15 MPI_Comm comm = honee->comm; in PRINT_EULER_VORTEX()
16 Ceed ceed = honee->ceed; in PRINT_EULER_VORTEX()
20 PetscCallCeed(ceed, CeedQFunctionContextGetData(problem->ics.qfctx, CEED_MEM_HOST, &euler_ctx)); in PRINT_EULER_VORTEX()
25 " Background Velocity : %f,%f,%f\n" in PRINT_EULER_VORTEX()
27 …app_ctx->problem_name, EulerTestTypes[euler_ctx->euler_test], euler_ctx->mean_velocity[0], euler_c… in PRINT_EULER_VORTEX()
28 euler_ctx->mean_velocity[2], StabilizationTypes[euler_ctx->stabilization])); in PRINT_EULER_VORTEX()
30 PetscCallCeed(ceed, CeedQFunctionContextRestoreData(problem->ics.qfctx, &euler_ctx)); in PRINT_EULER_VORTEX()
39 …PetscCall(HoneeBCCreateIFunctionQF(bc_def, Euler_Outflow, Euler_Outflow_loc, honee_bc->qfctx, qf)); in EulerVortexOutflowBCSetup_CreateIFunctionQF()
48 …reateIFunctionQF(bc_def, TravelingVortex_Inflow, TravelingVortex_Inflow_loc, honee_bc->qfctx, qf)); in EulerVortexInflowBCSetup_CreateIFunctionQF()
58 MPI_Comm comm = honee->comm; in NS_EULER_VORTEX()
59 Ceed ceed = honee->ceed; in NS_EULER_VORTEX()
68 // ------------------------------------------------------ in NS_EULER_VORTEX()
70 // ------------------------------------------------------ in NS_EULER_VORTEX()
71 …problem->ics = (HoneeQFSpec){.qf_func_ptr = ICsEuler, .qf_loc = ICsEuler_… in NS_EULER_VORTEX()
72 problem->apply_vol_rhs = (HoneeQFSpec){.qf_func_ptr = Euler, .qf_loc = Euler_loc}; in NS_EULER_VORTEX()
73 …problem->apply_vol_ifunction = (HoneeQFSpec){.qf_func_ptr = IFunction_Euler, .qf_loc = IF… in NS_EULER_VORTEX()
74 problem->num_comps_jac_data = 0; in NS_EULER_VORTEX()
75 problem->compute_exact_solution_error = PETSC_TRUE; in NS_EULER_VORTEX()
76 problem->print_info = PRINT_EULER_VORTEX; in NS_EULER_VORTEX()
78 problem->num_components = 5; in NS_EULER_VORTEX()
79 PetscCall(PetscMalloc1(problem->num_components, &problem->component_names)); in NS_EULER_VORTEX()
80 …for (PetscInt i = 0; i < 5; i++) PetscCall(PetscStrallocpy(component_names[i], &problem->component… in NS_EULER_VORTEX()
82 // ------------------------------------------------------ in NS_EULER_VORTEX()
84 // ------------------------------------------------------ in NS_EULER_VORTEX()
85 Units units = honee->units; in NS_EULER_VORTEX()
86 CeedScalar vortex_strength = 5.; // - in NS_EULER_VORTEX()
87 CeedScalar c_tau = 0.5; // - 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()
96 // ------------------------------------------------------ in NS_EULER_VORTEX()
98 // ------------------------------------------------------ in NS_EULER_VORTEX()
100 // -- Physics in NS_EULER_VORTEX()
101 …PetscCall(PetscOptionsScalar("-vortex_strength", "Strength of Vortex", NULL, vortex_strength, &vor… in NS_EULER_VORTEX()
104 …PetscCall(PetscOptionsRealArray("-mean_velocity", "Background velocity vector", NULL, mean_velocit… 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()
107 PetscCall(PetscOptionsRealArray("-center", "Location of vortex center", NULL, center, &n, NULL)); in NS_EULER_VORTEX()
108 …PetscCall(PetscOptionsBool("-implicit", "Use implicit (IFunction) formulation", NULL, implicit = P… in NS_EULER_VORTEX()
109 …PetscCall(PetscOptionsEnum("-euler_test", "Euler test option", NULL, EulerTestTypes, (PetscEnum)(e… in NS_EULER_VORTEX()
111 …PetscCall(PetscOptionsEnum("-stab", "Stabilization method", NULL, StabilizationTypes, (PetscEnum)(… in NS_EULER_VORTEX()
112 PetscCall(PetscOptionsScalar("-c_tau", "Stabilization constant", NULL, c_tau, &c_tau, NULL)); in NS_EULER_VORTEX()
114 // -- Warnings in NS_EULER_VORTEX()
116 PetscCall(PetscPrintf(comm, "Warning! Use -stab supg only with -implicit\n")); in NS_EULER_VORTEX()
119 …PetscCall(PetscPrintf(comm, "Warning! Background velocity vector for -euler_test t1 and -euler_tes… in NS_EULER_VORTEX()
124 // ------------------------------------------------------ in NS_EULER_VORTEX()
126 // ------------------------------------------------------ in NS_EULER_VORTEX()
127 // -- Scale variables to desired units in NS_EULER_VORTEX()
129 center[i] *= units->meter; in NS_EULER_VORTEX()
130 mean_velocity[i] *= (units->meter / units->second); in NS_EULER_VORTEX()
133 // -- QFunction Context in NS_EULER_VORTEX()
134 honee->phys->implicit = implicit; in NS_EULER_VORTEX()
135 euler_ctx->curr_time = 0.; in NS_EULER_VORTEX()
136 euler_ctx->implicit = implicit; in NS_EULER_VORTEX()
137 euler_ctx->euler_test = euler_test; in NS_EULER_VORTEX()
138 euler_ctx->center[0] = center[0]; in NS_EULER_VORTEX()
139 euler_ctx->center[1] = center[1]; in NS_EULER_VORTEX()
140 euler_ctx->center[2] = center[2]; in NS_EULER_VORTEX()
141 euler_ctx->vortex_strength = vortex_strength; in NS_EULER_VORTEX()
142 euler_ctx->c_tau = c_tau; in NS_EULER_VORTEX()
143 euler_ctx->mean_velocity[0] = mean_velocity[0]; in NS_EULER_VORTEX()
144 euler_ctx->mean_velocity[1] = mean_velocity[1]; in NS_EULER_VORTEX()
145 euler_ctx->mean_velocity[2] = mean_velocity[2]; in NS_EULER_VORTEX()
146 euler_ctx->stabilization = stab; in NS_EULER_VORTEX()
148 PetscCallCeed(ceed, CeedQFunctionContextCreate(honee->ceed, &euler_qfctx)); in NS_EULER_VORTEX()
153 PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(euler_qfctx, &problem->ics.qfctx)); in NS_EULER_VORTEX()
154 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(euler_qfctx, &problem->apply_vol_rhs.qfctx)); in NS_EULER_VORTEX()
155 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(euler_qfctx, &problem->apply_vol_ifunction.q… in NS_EULER_VORTEX()
157 for (PetscCount b = 0; b < problem->num_bc_defs; b++) { in NS_EULER_VORTEX()
158 BCDefinition bc_def = problem->bc_defs[b]; in NS_EULER_VORTEX()
166 PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(euler_qfctx, &honee_bc->qfctx)); in NS_EULER_VORTEX()
167 honee_bc->honee = honee; in NS_EULER_VORTEX()
168 honee_bc->num_comps_jac_data = 0; in NS_EULER_VORTEX()
177 PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(euler_qfctx, &honee_bc->qfctx)); in NS_EULER_VORTEX()
178 honee_bc->honee = honee; in NS_EULER_VORTEX()
179 honee_bc->num_comps_jac_data = 0; in NS_EULER_VORTEX()