honee/src/differential_filter.c

17   Ceed                ceed      = honee->ceed;  in DifferentialFilterCreateOperators()  local
29 …PetscCall(DMPlexCeedCoordinateCreateField(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE,…  in DifferentialFilterCreateOperators()
30 …PetscCall(QDataGet(ceed, dm_filter, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, &elem_restr_qd, &q_dat…  in DifferentialFilterCreateOperators()
39 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_RHS_Prim, Differential…  in DifferentialFilterCreateOperators()
42 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_RHS_Conserv, Different…  in DifferentialFilterCreateOperators()
45 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_RHS_Entropy, Different…  in DifferentialFilterCreateOperators()
49       PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_rhs));  in DifferentialFilterCreateOperators()
50 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_MMS_RHS, DifferentialF…  in DifferentialFilterCreateOperators()
52 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE,…  in DifferentialFilterCreateOperators()
53 …PetscCall(DMPlexCeedBasisCreate(ceed, honee->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, 0, 0, &ba…  in DifferentialFilterCreateOperators()
55     PetscCallCeed(ceed, CeedQFunctionSetContext(qf_rhs, diff_filter_qfctx));  in DifferentialFilterCreateOperators()
56     PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs, "q", num_comp_q, CEED_EVAL_INTERP));  in DifferentialFilterCreateOperators()
57     PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs, "qdata", q_data_size, CEED_EVAL_NONE));  in DifferentialFilterCreateOperators()
58     PetscCallCeed(ceed, CeedQFunctionAddInput(qf_rhs, "x", num_comp_x, CEED_EVAL_INTERP));  in DifferentialFilterCreateOperators()
63 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs, field_name, diff_filter->num_field_components[i…  in DifferentialFilterCreateOperators()
66     PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_rhs, NULL, NULL, &op_rhs));  in DifferentialFilterCreateOperators()
67 … PetscCallCeed(ceed, CeedOperatorSetField(op_rhs, "q", elem_restr_q, basis_q, CEED_VECTOR_ACTIVE));  in DifferentialFilterCreateOperators()
68 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs, "qdata", elem_restr_qd, CEED_BASIS_NONE, q_data));  in DifferentialFilterCreateOperators()
69     PetscCallCeed(ceed, CeedOperatorSetField(op_rhs, "x", elem_restr_x, basis_x, x_coord));  in DifferentialFilterCreateOperators()
75 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, dm_filter, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE,…  in DifferentialFilterCreateOperators()
76 …PetscCall(DMPlexCeedBasisCreate(ceed, dm_filter, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, height, d…  in DifferentialFilterCreateOperators()
79 …PetscCallCeed(ceed, CeedOperatorSetField(op_rhs, field_name, elem_restr_filter, basis_filter, CEED…  in DifferentialFilterCreateOperators()
81       PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_filter));  in DifferentialFilterCreateOperators()
82       PetscCallCeed(ceed, CeedBasisDestroy(&basis_filter));  in DifferentialFilterCreateOperators()
85 …PetscCall(OperatorApplyContextCreate(honee->dm, dm_filter, ceed, op_rhs, NULL, NULL, honee->Q_loc,…  in DifferentialFilterCreateOperators()
87     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q));  in DifferentialFilterCreateOperators()
88     PetscCallCeed(ceed, CeedBasisDestroy(&basis_q));  in DifferentialFilterCreateOperators()
89     PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_rhs));  in DifferentialFilterCreateOperators()
90     PetscCallCeed(ceed, CeedOperatorDestroy(&op_rhs));  in DifferentialFilterCreateOperators()
101 …PetscCall(GridAnisotropyTensorCalculateCollocatedVector(ceed, honee, &elem_restr_grid_aniso, &grid…  in DifferentialFilterCreateOperators()
103     PetscCallCeed(ceed, CeedOperatorCreateComposite(ceed, &op_lhs));  in DifferentialFilterCreateOperators()
113 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_LHS_1, DifferentialFil…  in DifferentialFilterCreateOperators()
116 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_LHS_5, DifferentialFil…  in DifferentialFilterCreateOperators()
119 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_LHS_6, DifferentialFil…  in DifferentialFilterCreateOperators()
122 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, DifferentialFilter_LHS_11, DifferentialFi…  in DifferentialFilterCreateOperators()
129       PetscCallCeed(ceed, CeedQFunctionSetContext(qf_lhs, diff_filter_qfctx));  in DifferentialFilterCreateOperators()
130       PetscCallCeed(ceed, CeedQFunctionSetUserFlopsEstimate(qf_lhs, 0));  in DifferentialFilterCreateOperators()
131       PetscCallCeed(ceed, CeedQFunctionAddInput(qf_lhs, "q", num_comp_filter, CEED_EVAL_INTERP));  in DifferentialFilterCreateOperators()
132 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_lhs, "Grad_q", num_comp_filter * dim, CEED_EVAL_GRAD)…  in DifferentialFilterCreateOperators()
133 …PetscCallCeed(ceed, CeedQFunctionAddInput(qf_lhs, "anisotropy tensor", num_comp_grid_aniso, CEED_E…  in DifferentialFilterCreateOperators()
134       PetscCallCeed(ceed, CeedQFunctionAddInput(qf_lhs, "x", num_comp_x, CEED_EVAL_INTERP));  in DifferentialFilterCreateOperators()
135       PetscCallCeed(ceed, CeedQFunctionAddInput(qf_lhs, "qdata", q_data_size, CEED_EVAL_NONE));  in DifferentialFilterCreateOperators()
136       PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_lhs, "v", num_comp_filter, CEED_EVAL_INTERP));  in DifferentialFilterCreateOperators()
137 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_lhs, "Grad_v", num_comp_filter * dim, CEED_EVAL_GRAD…  in DifferentialFilterCreateOperators()
144 …PetscCallCeed(ceed, CeedOperatorGetFieldByName(diff_filter->op_rhs_ctx->op, field_name, &op_field)…  in DifferentialFilterCreateOperators()
145 …PetscCallCeed(ceed, CeedOperatorFieldGetData(op_field, NULL, &elem_restr_filter, &basis_filter, NU…  in DifferentialFilterCreateOperators()
148       PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_lhs, NULL, NULL, &op_lhs_sub));  in DifferentialFilterCreateOperators()
149 …PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "q", elem_restr_filter, basis_filter, CEED_VE…  in DifferentialFilterCreateOperators()
150 …PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "Grad_q", elem_restr_filter, basis_filter, CE…  in DifferentialFilterCreateOperators()
151 …PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "anisotropy tensor", elem_restr_grid_aniso, C…  in DifferentialFilterCreateOperators()
152       PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "x", elem_restr_x, basis_x, x_coord));  in DifferentialFilterCreateOperators()
153 …PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "qdata", elem_restr_qd, CEED_BASIS_NONE, q_da…  in DifferentialFilterCreateOperators()
154 …PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "v", elem_restr_filter, basis_filter, CEED_VE…  in DifferentialFilterCreateOperators()
155 …PetscCallCeed(ceed, CeedOperatorSetField(op_lhs_sub, "Grad_v", elem_restr_filter, basis_filter, CE…  in DifferentialFilterCreateOperators()
157       PetscCallCeed(ceed, CeedOperatorCompositeAddSub(op_lhs, op_lhs_sub));  in DifferentialFilterCreateOperators()
158       PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_filter));  in DifferentialFilterCreateOperators()
159       PetscCallCeed(ceed, CeedBasisDestroy(&basis_filter));  in DifferentialFilterCreateOperators()
160       PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_lhs));  in DifferentialFilterCreateOperators()
161       PetscCallCeed(ceed, CeedOperatorDestroy(&op_lhs_sub));  in DifferentialFilterCreateOperators()
163     PetscCallCeed(ceed, CeedVectorDestroy(&grid_aniso_ceed));  in DifferentialFilterCreateOperators()
164     PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_grid_aniso));  in DifferentialFilterCreateOperators()
166 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_lhs, "filter width scaling", &diff_filter-…  in DifferentialFilterCreateOperators()
183     PetscCallCeed(ceed, CeedOperatorDestroy(&op_lhs));  in DifferentialFilterCreateOperators()
186   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_x));  in DifferentialFilterCreateOperators()
187   PetscCallCeed(ceed, CeedBasisDestroy(&basis_x));  in DifferentialFilterCreateOperators()
188   PetscCallCeed(ceed, CeedVectorDestroy(&x_coord));  in DifferentialFilterCreateOperators()
189   PetscCallCeed(ceed, CeedVectorDestroy(&q_data));  in DifferentialFilterCreateOperators()
190   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_qd));  in DifferentialFilterCreateOperators()
197   Ceed                      ceed    = honee->ceed;  in DifferentialFilterSetup()  local
285 …PetscCallCeed(ceed, CeedQFunctionContextGetDataRead(problem->apply_vol_ifunction.qfctx, CEED_MEM_H…  in DifferentialFilterSetup()
287 …PetscCallCeed(ceed, CeedQFunctionContextRestoreDataRead(problem->apply_vol_ifunction.qfctx, &newt_…  in DifferentialFilterSetup()
289   PetscCallCeed(ceed, CeedQFunctionContextCreate(ceed, &diff_filter_qfctx));  in DifferentialFilterSetup()
290 …PetscCallCeed(ceed, CeedQFunctionContextSetData(diff_filter_qfctx, CEED_MEM_HOST, CEED_USE_POINTER…  in DifferentialFilterSetup()
291 …PetscCallCeed(ceed, CeedQFunctionContextSetDataDestroy(diff_filter_qfctx, CEED_MEM_HOST, FreeConte…  in DifferentialFilterSetup()
292   PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(diff_filter_qfctx, "filter width scaling",  in DifferentialFilterSetup()
301   PetscCallCeed(ceed, CeedQFunctionContextDestroy(&diff_filter_qfctx));  in DifferentialFilterSetup()