1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
6 // This file is part of CEED:  http://github.com/ceed
10 #include <ceed.h>
16 PetscErrorCode SetupStrongSTG_Ceed(Ceed ceed, CeedData ceed_data, DM dm, ProblemData problem, Simpl…  in SetupStrongSTG_Ceed()  argument
36     PetscCall(CreateBasisFromPlex(ceed, dm, label, label_value, height, dm_field, &basis_q_sur));  in SetupStrongSTG_Ceed()
37 …PetscCall(CreateBasisFromPlex(ceed, dm_coord, label, label_value, height, dm_field, &basis_x_sur));  in SetupStrongSTG_Ceed()
39     PetscCallCeed(ceed, CeedBasisCreateProjection(basis_x_sur, basis_q_sur, &basis_x_to_q_sur));  in SetupStrongSTG_Ceed()
41     PetscCallCeed(ceed, CeedBasisDestroy(&basis_q_sur));  in SetupStrongSTG_Ceed()
42     PetscCallCeed(ceed, CeedBasisDestroy(&basis_x_sur));  in SetupStrongSTG_Ceed()
46 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, SetupStrongBC, SetupStrongBC_loc, &qf_set…  in SetupStrongSTG_Ceed()
47   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_setup, "x", num_comp_x, CEED_EVAL_INTERP));  in SetupStrongSTG_Ceed()
48 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_setup, "dxdX", num_comp_x * dim_boundary, CEED_EVAL_G…  in SetupStrongSTG_Ceed()
49   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_setup, "multiplicity", num_comp_q, CEED_EVAL_NONE));  in SetupStrongSTG_Ceed()
50   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_setup, "x stored", num_comp_x, CEED_EVAL_NONE));  in SetupStrongSTG_Ceed()
51   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_setup, "scale", 1, CEED_EVAL_NONE));  in SetupStrongSTG_Ceed()
52   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_setup, "dXdx", dXdx_size, CEED_EVAL_NONE));  in SetupStrongSTG_Ceed()
55 …PetscCall(SetupStrongStg_PreProcessing(ceed, problem, num_comp_x, stg_data_size, dXdx_size, &qf_st…  in SetupStrongSTG_Ceed()
56 …PetscCall(SetupStrongStg_QF(ceed, problem, num_comp_x, num_comp_q, stg_data_size, dXdx_size, &qf_s…  in SetupStrongSTG_Ceed()
61 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, dm, domain_label, bc->inflows[i], height, dm_field…  in SetupStrongSTG_Ceed()
62 …PetscCall(DMPlexCeedElemRestrictionCoordinateCreate(ceed, dm, domain_label, bc->inflows[i], height…  in SetupStrongSTG_Ceed()
63     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_q_sur, &multiplicity, NULL));  in SetupStrongSTG_Ceed()
64     PetscCallCeed(ceed, CeedElemRestrictionGetMultiplicity(elem_restr_q_sur, multiplicity));  in SetupStrongSTG_Ceed()
66 …PetscCall(DMPlexCeedElemRestrictionCollocatedCreate(ceed, dm, domain_label, bc->inflows[i], height…  in SetupStrongSTG_Ceed()
67 …PetscCall(DMPlexCeedElemRestrictionCollocatedCreate(ceed, dm, domain_label, bc->inflows[i], height…  in SetupStrongSTG_Ceed()
68 …PetscCall(DMPlexCeedElemRestrictionCollocatedCreate(ceed, dm, domain_label, bc->inflows[i], height…  in SetupStrongSTG_Ceed()
69 …PetscCall(DMPlexCeedElemRestrictionCollocatedCreate(ceed, dm, domain_label, bc->inflows[i], height…  in SetupStrongSTG_Ceed()
70     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_x_stored, &x_stored, NULL));  in SetupStrongSTG_Ceed()
71     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_scale, &scale_stored, NULL));  in SetupStrongSTG_Ceed()
72     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_stgdata, &stg_data, NULL));  in SetupStrongSTG_Ceed()
73     PetscCallCeed(ceed, CeedElemRestrictionCreateVector(elem_restr_dXdx, &dXdx, NULL));  in SetupStrongSTG_Ceed()
76     PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_setup, NULL, NULL, &op_setup));  in SetupStrongSTG_Ceed()
77 …PetscCallCeed(ceed, CeedOperatorSetName(op_setup, "Precomputed data for strong boundary conditions…  in SetupStrongSTG_Ceed()
78 …PetscCallCeed(ceed, CeedOperatorSetField(op_setup, "x", elem_restr_x_sur, basis_x_to_q_sur, CEED_V…  in SetupStrongSTG_Ceed()
79 …PetscCallCeed(ceed, CeedOperatorSetField(op_setup, "dxdX", elem_restr_x_sur, basis_x_to_q_sur, CEE…  in SetupStrongSTG_Ceed()
80 …PetscCallCeed(ceed, CeedOperatorSetField(op_setup, "multiplicity", elem_restr_q_sur, CEED_BASIS_NO…  in SetupStrongSTG_Ceed()
81 …PetscCallCeed(ceed, CeedOperatorSetField(op_setup, "x stored", elem_restr_x_stored, CEED_BASIS_NON…  in SetupStrongSTG_Ceed()
82 …PetscCallCeed(ceed, CeedOperatorSetField(op_setup, "scale", elem_restr_scale, CEED_BASIS_NONE, sca…  in SetupStrongSTG_Ceed()
83 …PetscCallCeed(ceed, CeedOperatorSetField(op_setup, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXdx)…  in SetupStrongSTG_Ceed()
86 …PetscCallCeed(ceed, CeedOperatorApply(op_setup, ceed_data->x_coord, CEED_VECTOR_NONE, CEED_REQUEST…  in SetupStrongSTG_Ceed()
89     PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_stgdata, NULL, NULL, &op_stgdata));  in SetupStrongSTG_Ceed()
90 …PetscCallCeed(ceed, CeedOperatorSetField(op_stgdata, "dXdx", elem_restr_dXdx, CEED_BASIS_NONE, dXd…  in SetupStrongSTG_Ceed()
91 …PetscCallCeed(ceed, CeedOperatorSetField(op_stgdata, "x", elem_restr_x_stored, CEED_BASIS_NONE, x_…  in SetupStrongSTG_Ceed()
92 …PetscCallCeed(ceed, CeedOperatorSetField(op_stgdata, "stg data", elem_restr_stgdata, CEED_BASIS_NO…  in SetupStrongSTG_Ceed()
94 …PetscCallCeed(ceed, CeedOperatorApply(op_stgdata, CEED_VECTOR_NONE, stg_data, CEED_REQUEST_IMMEDIA…  in SetupStrongSTG_Ceed()
97     PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_strongbc, NULL, NULL, &op_strong_bc_sub));  in SetupStrongSTG_Ceed()
98     PetscCallCeed(ceed, CeedOperatorSetName(op_strong_bc_sub, "Strong STG"));  in SetupStrongSTG_Ceed()
100 …PetscCallCeed(ceed, CeedOperatorSetField(op_strong_bc_sub, "dXdx", elem_restr_dXdx, CEED_BASIS_NON…  in SetupStrongSTG_Ceed()
101 …PetscCallCeed(ceed, CeedOperatorSetField(op_strong_bc_sub, "x", elem_restr_x_stored, CEED_BASIS_NO…  in SetupStrongSTG_Ceed()
102 …PetscCallCeed(ceed, CeedOperatorSetField(op_strong_bc_sub, "scale", elem_restr_scale, CEED_BASIS_N…  in SetupStrongSTG_Ceed()
103 …PetscCallCeed(ceed, CeedOperatorSetField(op_strong_bc_sub, "stg data", elem_restr_stgdata, CEED_BA…  in SetupStrongSTG_Ceed()
104 …PetscCallCeed(ceed, CeedOperatorSetField(op_strong_bc_sub, "q", elem_restr_q_sur, CEED_BASIS_NONE,…  in SetupStrongSTG_Ceed()
107     PetscCallCeed(ceed, CeedOperatorCompositeAddSub(op_strong_bc, op_strong_bc_sub));  in SetupStrongSTG_Ceed()
109     PetscCallCeed(ceed, CeedVectorDestroy(&multiplicity));  in SetupStrongSTG_Ceed()
110     PetscCallCeed(ceed, CeedVectorDestroy(&x_stored));  in SetupStrongSTG_Ceed()
111     PetscCallCeed(ceed, CeedVectorDestroy(&scale_stored));  in SetupStrongSTG_Ceed()
112     PetscCallCeed(ceed, CeedVectorDestroy(&stg_data));  in SetupStrongSTG_Ceed()
113     PetscCallCeed(ceed, CeedVectorDestroy(&dXdx));  in SetupStrongSTG_Ceed()
114     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_x_sur));  in SetupStrongSTG_Ceed()
115     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q_sur));  in SetupStrongSTG_Ceed()
116     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_x_stored));  in SetupStrongSTG_Ceed()
117     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_scale));  in SetupStrongSTG_Ceed()
118     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_stgdata));  in SetupStrongSTG_Ceed()
119     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_dXdx));  in SetupStrongSTG_Ceed()
120     PetscCallCeed(ceed, CeedOperatorDestroy(&op_strong_bc_sub));  in SetupStrongSTG_Ceed()
121     PetscCallCeed(ceed, CeedOperatorDestroy(&op_setup));  in SetupStrongSTG_Ceed()
122     PetscCallCeed(ceed, CeedOperatorDestroy(&op_stgdata));  in SetupStrongSTG_Ceed()
125 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_strong_bc, "solution time", &phys->stg_sol…  in SetupStrongSTG_Ceed()
127   PetscCallCeed(ceed, CeedBasisDestroy(&basis_x_to_q_sur));  in SetupStrongSTG_Ceed()
128   PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_strongbc));  in SetupStrongSTG_Ceed()
129   PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_stgdata));  in SetupStrongSTG_Ceed()
130   PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_setup));  in SetupStrongSTG_Ceed()
143 …PetscCallCeed(user->ceed, CeedOperatorSetContextDouble(user->op_strong_bc_ctx->op, user->phys->stg…  in DMPlexInsertBoundaryValues_StrongBCCeed()
155 PetscErrorCode SetupStrongBC_Ceed(Ceed ceed, CeedData ceed_data, DM dm, User user, ProblemData prob…  in SetupStrongBC_Ceed()  argument
171   PetscCallCeed(ceed, CeedOperatorCreateComposite(ceed, &op_strong_bc));  in SetupStrongBC_Ceed()
177       PetscCall(SetupStrongSTG_Ceed(ceed, ceed_data, dm, problem, bc, user->phys, op_strong_bc));  in SetupStrongBC_Ceed()
181 …PetscCall(OperatorApplyContextCreate(NULL, NULL, ceed, op_strong_bc, CEED_VECTOR_NONE, NULL, NULL,…  in SetupStrongBC_Ceed()
184   PetscCallCeed(ceed, CeedOperatorDestroy(&op_strong_bc));  in SetupStrongBC_Ceed()