| /honee/problems/ |
| H A D | shocktube.c | 34 PetscInt stab; in NS_SHOCKTUBE() local 68 …all(PetscOptionsEnum("-stab", "Stabilization method", NULL, StabilizationTypes, (PetscEnum)(stab =… in NS_SHOCKTUBE() 73 if (stab == STAB_SUPG) { in NS_SHOCKTUBE() 74 … PetscCall(PetscPrintf(comm, "Warning! -stab supg not implemented for the shocktube problem. \n")); in NS_SHOCKTUBE() 99 shocktube_ctx->stabilization = stab; in NS_SHOCKTUBE()
|
| H A D | eulervortex.c | 57 StabilizationType stab; in NS_EULER_VORTEX() local 111 …all(PetscOptionsEnum("-stab", "Stabilization method", NULL, StabilizationTypes, (PetscEnum)(stab =… in NS_EULER_VORTEX() 115 if (stab == STAB_SUPG && !implicit) { in NS_EULER_VORTEX() 116 PetscCall(PetscPrintf(comm, "Warning! Use -stab supg only with -implicit\n")); in NS_EULER_VORTEX() 146 euler_ctx->stabilization = stab; in NS_EULER_VORTEX()
|
| H A D | advection.c | 385 StabilizationType stab; in NS_ADVECTION() local 437 …all(PetscOptionsEnum("-stab", "Stabilization method", NULL, StabilizationTypes, (PetscEnum)(stab =… in NS_ADVECTION() 484 if (stab == STAB_NONE && CtauS != 0) { in NS_ADVECTION() 485 PetscCall(PetscPrintf(comm, "Warning! Use -CtauS only with -stab su or -stab supg\n")); in NS_ADVECTION() 489 if (stab == STAB_SUPG) problem->create_mass_operator = CreateKSPMassOperator_AdvectionStabilized; in NS_ADVECTION() 525 advection_ctx->stabilization = stab; in NS_ADVECTION()
|
| H A D | newtonian.c | 374 StabilizationType stab; in NS_NEWTONIAN_IG() local 431 …all(PetscOptionsEnum("-stab", "Stabilization method", NULL, StabilizationTypes, (PetscEnum)(stab =… in NS_NEWTONIAN_IG() 495 newtonian_ig_ctx->stabilization = stab; in NS_NEWTONIAN_IG() 570 …if (stab == STAB_SUPG && !implicit) problem->create_mass_operator = CreateKSPMassOperator_Newtonia… in NS_NEWTONIAN_IG()
|
| /honee/qfunctions/ |
| H A D | stabilization.h | 25 … const CeedScalar strong_residual[5], CeedScalar stab[5][3]) { in StabilizationMatrix() 28 // Zero stab so all future terms can safely sum into it in StabilizationMatrix() 29 SetValueN((CeedScalar *)stab, 0, 15); in StabilizationMatrix() 36 for (CeedInt j = 0; j < 5; j++) stab[j][i] += dF_i[j]; in StabilizationMatrix() 42 CeedScalar stab[5][3]) { in Stabilization() 56 StabilizationMatrix(gas, s, Tau_d, strong_residual, stab); in Stabilization()
|
| H A D | eulervortex.h | 358 CeedScalar stab[5][3] = {{0.}}; in Euler() local 365 … for (CeedInt l = 0; l < 5; l++) stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in Euler() 370 …dInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] … in Euler() 484 CeedScalar stab[5][3] = {{0.}}; in IFunction_Euler() local 491 … for (CeedInt l = 0; l < 5; l++) stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in IFunction_Euler() 496 …dInt k = 0; k < 3; k++) dv[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] … in IFunction_Euler() 502 … for (CeedInt l = 0; l < 5; l++) stab[k][j] = jacob_F_conv[j][k][l] * Tau_x[j] * strong_res[l]; in IFunction_Euler() 507 …dInt k = 0; k < 3; k++) dv[k][j][i] += wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] … in IFunction_Euler()
|
| H A D | newtonian.h | 68 CeedScalar Tau_d[3], stab[5][3], body_force[5] = {0.}, divFdiff[5] = {0.}, U_dot[5]; in MassFunction_Newtonian() local 71 Stabilization(context->stabilization, gas, s, Tau_d, grad_s, U_dot, body_force, divFdiff, stab); in MassFunction_Newtonian() 76 …Grad_v[k][j][i] = wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] + stab[j][2] * dXdx[k… in MassFunction_Newtonian() 246 CeedScalar Tau_d[3], stab[5][3], U_dot[5] = {0}; in RHSFunction_Newtonian() local 248 … Stabilization(context->stabilization, gas, s, Tau_d, grad_s, U_dot, body_force, divFdiff_i, stab); in RHSFunction_Newtonian() 251 … k = 0; k < 3; k++) Grad_v[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] … in RHSFunction_Newtonian() 301 CeedScalar stab[5][3]; in IFunction_Newtonian_Integrand() local 303 Stabilization(context->stabilization, gas, s, Tau_d, grad_s, U_dot, body_force, divFdiff_i, stab); in IFunction_Newtonian_Integrand() 304 AXPY(1, (CeedScalar *)stab, (CeedScalar *)grad_v_i, 15); in IFunction_Newtonian_Integrand()
|
| H A D | shocktube.h | 341 CeedScalar stab[5][3] = {0}; in EulerShockTube() local 349 stab[k][j] += jacob_F_conv[j][k][l] * Tau_x[j] * strong_conv[l]; in EulerShockTube() 354 …dInt k = 0; k < 3; k++) dv[k][j][i] -= wdetJ * (stab[j][0] * dXdx[k][0] + stab[j][1] * dXdx[k][1] … in EulerShockTube()
|
| /honee/examples/ |
| H A D | navierstokes.c | 200 …ex_box_faces 2,2,1 -ts_max_steps 5 -degree 3 -implicit false -ts_type rk -stab supg -state_var con… 201 …stab supg -Ctaus 0.5 -mass_ksp_type gmres -mass_pc_type vpbjacobi -compare_final_state_atol 1e-10 … 202 …0.83205029,0 -advection_ic_type skew -dm_plex_box_faces 2,1,1 -degree 2 -stab supg -stab_tau advd… 206 …type cylinder -dm_plex_box_faces 2,2,1 -diffusion_coeff 5e-3 -implicit 0 -stab supg -advection_ic_… 222 …-compare_final_state_atol 5E-4 -compare_final_state_filename tests/output/dc-implicit-stab-none.bin 223 …stab supg -degree 3 -dm_plex_box_faces 1,1,2 -dm_plex_box_lower 0,0,0 -dm_plex_box_upper 125,125,2… 224 …stab su -degree 3 -dm_plex_box_faces 1,1,2 -dm_plex_box_lower 0,0,0 -dm_plex_box_upper 125,125,250… 226 …stab supg -degree 3 -dm_plex_box_faces 1,1,2 -dm_plex_box_lower 0,0 -dm_plex_box_upper 125,125 -bc… 229 …plex_dim 3 -bc_symmetry_x 5,6 -bc_symmetry_y 3,4 -bc_symmetry_z 1,2 -yzb -stab su -ts_max_steps 10…
|
| H A D | advection_wave.yaml | 2 stab: supg
|
| H A D | advection.yaml | 3 stab: su
|
| H A D | advection_bl.yaml | 2 stab: supg
|
| H A D | taylor_green_vortex.yaml | 27 stab: 'supg'
|
| H A D | gaussianwave.yaml | 31 stab: supg
|
| H A D | flatplate_STG.yaml | 22 stab: 'supg'
|
| H A D | blasius.yaml | 31 stab: 'supg'
|
| H A D | vortexshedding.yaml | 5 stab: supg
|
| /honee/tests/ |
| H A D | blasius_test.yaml | 18 stab: 'none'
|
| H A D | blasius_stgtest.yaml | 17 stab: 'supg'
|
| H A D | channel_divdiff_verify.yaml | 22 stab: supg
|
| H A D | gaussianwave_cgns_load.yaml | 27 stab: supg
|
| H A D | stats_test.yaml | 17 stab: 'supg'
|
| /honee/doc/ |
| H A D | examples.md | 525 * - `-stab` 539 ./navierstokes -problem shocktube -yzb -stab su -bc_symmetry_z 3,4 -bc_symmetry_y 1,2 -bc_wall 5,6 … 632 * - `-stab` 721 …_box_upper 8000,8000,8000 -bc_wall 1,2,3,4,5,6 -wall_comps 4 -wind_type rotation -implicit -stab su 733 …dm_plex_box_upper 1000,1000 -bc_wall 1,2,3,4 -wall_comps 4 -wind_type rotation -implicit -stab supg
|
| H A D | runtime_options.md | 408 * - `-stab`
|
| H A D | theory.md | 192 This stabilization technique can be selected using the option `-stab supg`. 208 This stabilization technique can be selected using the option `-stab su`.
|