174     PetscCallCeed(ceed, CeedOperatorCompositeGetSubList(user->op_rhs_ctx->op, &sub_ops));  in CreateKSPMassOperator_NewtonianStabilized()
175     PetscCallCeed(ceed, CeedOperatorGetFieldByName(sub_ops[sub_op_index], "q", &field));  in CreateKSPMassOperator_NewtonianStabilized()
176     PetscCallCeed(ceed, CeedOperatorFieldGetData(field, NULL, &elem_restr_q, &basis_q, NULL));  in CreateKSPMassOperator_NewtonianStabilized()
178     PetscCallCeed(ceed, CeedOperatorGetFieldByName(sub_ops[sub_op_index], "qdata", &field));  in CreateKSPMassOperator_NewtonianStabilized()
179     PetscCallCeed(ceed, CeedOperatorFieldGetData(field, NULL, &elem_restr_qd_i, NULL, &q_data));  in CreateKSPMassOperator_NewtonianStabilized()
181     PetscCallCeed(ceed, CeedOperatorGetContext(sub_ops[sub_op_index], &qf_ctx));  in CreateKSPMassOperator_NewtonianStabilized()
184   PetscCallCeed(ceed, CeedElemRestrictionGetNumComponents(elem_restr_q, &num_comp_q));  in CreateKSPMassOperator_NewtonianStabilized()
185   PetscCallCeed(ceed, CeedElemRestrictionGetNumComponents(elem_restr_qd_i, &q_data_size));  in CreateKSPMassOperator_NewtonianStabilized()
187 …PetscCallCeed(ceed, CeedQFunctionCreateInterior(ceed, 1, MassFunction_Newtonian_Conserv, MassFunct…  in CreateKSPMassOperator_NewtonianStabilized()
189   PetscCallCeed(ceed, CeedQFunctionSetContext(qf_mass, qf_ctx));  in CreateKSPMassOperator_NewtonianStabilized()
190   PetscCallCeed(ceed, CeedQFunctionSetUserFlopsEstimate(qf_mass, 0));  in CreateKSPMassOperator_NewtonianStabilized()
191   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_mass, "q_dot", 5, CEED_EVAL_INTERP));  in CreateKSPMassOperator_NewtonianStabilized()
192   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_mass, "q", 5, CEED_EVAL_INTERP));  in CreateKSPMassOperator_NewtonianStabilized()
193   PetscCallCeed(ceed, CeedQFunctionAddInput(qf_mass, "qdata", q_data_size, CEED_EVAL_NONE));  in CreateKSPMassOperator_NewtonianStabilized()
194   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_mass, "v", 5, CEED_EVAL_INTERP));  in CreateKSPMassOperator_NewtonianStabilized()
195   PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_mass, "Grad_v", 5 * dim, CEED_EVAL_GRAD));  in CreateKSPMassOperator_NewtonianStabilized()
197   PetscCallCeed(ceed, CeedOperatorCreate(ceed, qf_mass, NULL, NULL, op_mass));  in CreateKSPMassOperator_NewtonianStabilized()
198 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "q_dot", elem_restr_q, basis_q, CEED_VECTOR_ACT…  in CreateKSPMassOperator_NewtonianStabilized()
199   PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "q", elem_restr_q, basis_q, user->q_ceed));  in CreateKSPMassOperator_NewtonianStabilized()
200 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "qdata", elem_restr_qd_i, CEED_BASIS_NONE, q_da…  in CreateKSPMassOperator_NewtonianStabilized()
201 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "v", elem_restr_q, basis_q, CEED_VECTOR_ACTIVE)…  in CreateKSPMassOperator_NewtonianStabilized()
202 …PetscCallCeed(ceed, CeedOperatorSetField(*op_mass, "Grad_v", elem_restr_q, basis_q, CEED_VECTOR_AC…  in CreateKSPMassOperator_NewtonianStabilized()
204   PetscCallCeed(ceed, CeedVectorDestroy(&q_data));  in CreateKSPMassOperator_NewtonianStabilized()
205   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_q));  in CreateKSPMassOperator_NewtonianStabilized()
206   PetscCallCeed(ceed, CeedElemRestrictionDestroy(&elem_restr_qd_i));  in CreateKSPMassOperator_NewtonianStabilized()
207   PetscCallCeed(ceed, CeedBasisDestroy(&basis_q));  in CreateKSPMassOperator_NewtonianStabilized()
208   PetscCallCeed(ceed, CeedQFunctionContextDestroy(&qf_ctx));  in CreateKSPMassOperator_NewtonianStabilized()
209   PetscCallCeed(ceed, CeedQFunctionDestroy(&qf_mass));  in CreateKSPMassOperator_NewtonianStabilized()
443   PetscCallCeed(ceed, CeedQFunctionContextCreate(user->ceed, &problem->ics.qfunction_context));  in NS_NEWTONIAN_IG()
444   PetscCallCeed(ceed,  in NS_NEWTONIAN_IG()
446 …PetscCallCeed(ceed, CeedQFunctionContextSetDataDestroy(problem->ics.qfunction_context, CEED_MEM_HO…  in NS_NEWTONIAN_IG()
447 …PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(problem->ics.qfunction_context, "evaluation…  in NS_NEWTONIAN_IG()
450   PetscCallCeed(ceed, CeedQFunctionContextCreate(user->ceed, &newtonian_ig_context));  in NS_NEWTONIAN_IG()
451   PetscCallCeed(ceed,  in NS_NEWTONIAN_IG()
453 …PetscCallCeed(ceed, CeedQFunctionContextSetDataDestroy(newtonian_ig_context, CEED_MEM_HOST, FreeCo…  in NS_NEWTONIAN_IG()
454 …PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(newtonian_ig_context, "timestep size", offs…  in NS_NEWTONIAN_IG()
456 …PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(newtonian_ig_context, "ijacobian time shift…  in NS_NEWTONIAN_IG()
459 …PetscCallCeed(ceed, CeedQFunctionContextRegisterDouble(newtonian_ig_context, "solution time", offs…  in NS_NEWTONIAN_IG()
463 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_vol_if…  in NS_NEWTONIAN_IG()
464 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_vol_ij…  in NS_NEWTONIAN_IG()
465 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_inflow…  in NS_NEWTONIAN_IG()
466 …PetscCallCeed(ceed, CeedQFunctionContextReferenceCopy(newtonian_ig_context, &problem->apply_inflow…  in NS_NEWTONIAN_IG()
484 …PetscCallCeed(ceed, CeedQFunctionContextGetData(problem->apply_vol_rhs.qfunction_context, CEED_MEM…  in PRINT_NEWTONIAN()
490 …PetscCallCeed(ceed, CeedQFunctionContextRestoreData(problem->apply_vol_rhs.qfunction_context, &new…  in PRINT_NEWTONIAN()