251     CeedInt             b_in, b_out, num_elem, num_nodes, num_qpts, num_comp;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()  local
258       b_in = b;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
260         if (active_bases_in[b_in] == active_bases_out[b_out]) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
269       for (b_in = 0; b_in < num_active_bases_in; b_in++) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
270         if (active_bases_in[b_in] == active_bases_out[b_out]) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
274       if (b_in == num_active_bases_in) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
278 …CeedCheck(active_elem_rstrs_in[b_in] == active_elem_rstrs_out[b_out], CeedOperatorReturnCeed(op), …  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
283 …CeedCall(CeedOperatorCreateActivePointBlockRestriction(active_elem_rstrs_in[b_in], &diag_elem_rstr…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
285       CeedCall(CeedElemRestrictionCreateUnsignedCopy(active_elem_rstrs_in[b_in], &diag_elem_rstr));  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
295     CeedCall(CeedBasisGetNumNodes(active_bases_in[b_in], &num_nodes));  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
296     CeedCall(CeedBasisGetNumComponents(active_bases_in[b_in], &num_comp));  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
297     if (active_bases_in[b_in] == CEED_BASIS_NONE) num_qpts = num_nodes;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
298     else CeedCall(CeedBasisGetNumQuadraturePoints(active_bases_in[b_in], &num_qpts));  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
301     for (CeedInt i = 0; i < num_eval_modes_in[b_in]; i++) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
302       has_eval_none = has_eval_none || (eval_modes_in[b_in][i] == CEED_EVAL_NONE);  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
332         for (CeedInt e_in = 0; e_in < num_eval_modes_in[b_in]; e_in++) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
335 …CeedCall(CeedOperatorGetBasisPointer(active_bases_in[b_in], eval_modes_in[b_in][e_in], identity, &…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
336 …CeedCall(CeedBasisGetNumQuadratureComponents(active_bases_in[b_in], eval_modes_in[b_in][e_in], &q_…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
338             if (e_in == 0 || eval_modes_in[b_in][e_in] != eval_mode_in_prev) d_in = 0;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
341           eval_mode_in_prev = eval_modes_in[b_in][e_in];  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
351 …(eval_mode_offsets_in[b_in][e_in] + c_in) * num_output_components + eval_mode_offsets_out[b_out][e…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
362 …(eval_mode_offsets_in[b_in][e_in] + c_out) * num_output_components + eval_mode_offsets_out[b_out][…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()