251     CeedInt             b_in, b_out, num_elem, num_nodes, num_qpts, num_comp;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()  local
259       for (b_out = 0; b_out < num_active_bases_out; b_out++) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
260         if (active_bases_in[b_in] == active_bases_out[b_out]) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
264       if (b_out == num_active_bases_out) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
268       b_out = b;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
270         if (active_bases_in[b_in] == active_bases_out[b_out]) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
278 …CeedCheck(active_elem_rstrs_in[b_in] == active_elem_rstrs_out[b_out], CeedOperatorReturnCeed(op), …  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
304     for (CeedInt i = 0; i < num_eval_modes_out[b_out]; i++) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
305       has_eval_none = has_eval_none || (eval_modes_out[b_out][i] == CEED_EVAL_NONE);  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
319       for (CeedInt e_out = 0; e_out < num_eval_modes_out[b_out]; e_out++) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
324 …CeedCall(CeedOperatorGetBasisPointer(active_bases_out[b_out], eval_modes_out[b_out][e_out], identi…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
325 …CeedCall(CeedBasisGetNumQuadratureComponents(active_bases_out[b_out], eval_modes_out[b_out][e_out]…  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
327           if (e_out == 0 || eval_modes_out[b_out][e_out] != eval_mode_out_prev) d_out = 0;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
330         eval_mode_out_prev = eval_modes_out[b_out][e_out];  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
351 …fsets_in[b_in][e_in] + c_in) * num_output_components + eval_mode_offsets_out[b_out][e_out] + c_out;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
362 …sets_in[b_in][e_in] + c_out) * num_output_components + eval_mode_offsets_out[b_out][e_out] + c_out;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()