42   // Check if NULL qf passed in  in CeedQFunctionCreateFallback()
43   if (!qf) return CEED_ERROR_SUCCESS;  in CeedQFunctionCreateFallback()
47   if (qf->source_path) {  in CeedQFunctionCreateFallback()
54   } else if (qf->user_source) {  in CeedQFunctionCreateFallback()
111   if (op->op_fallback) return CEED_ERROR_SUCCESS;  in CeedOperatorCreateFallback()
117   if (!ceed_fallback) return CEED_ERROR_SUCCESS;  in CeedOperatorCreateFallback()
123   if (is_composite) {  in CeedOperatorCreateFallback()
146     if (is_at_points) {  in CeedOperatorCreateFallback()
257     if (num_active_bases_in <= num_active_bases_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()
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()
281     // Assemble point block diagonal restriction, if needed  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
282     if (is_point_block) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
297     if (active_bases_in[b_in] == CEED_BASIS_NONE) num_qpts = num_nodes;  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
300     // Construct identity matrix for basis if required  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
307     if (has_eval_none) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
326         if (q_comp_out > 1) {  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()
337           if (q_comp_in > 1) {  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()
347               if (is_point_block) {  in CeedOperatorLinearAssembleAddDiagonalSingle_Mesh()
431     if (is_point_block) {  in CeedOperatorLinearAssembleAddDiagonalComposite()
489   if (elem_rstr_in != elem_rstr_out) {  in CeedOperatorAssembleSymbolicSingle()
544   if (elem_rstr_in != elem_rstr_out) {  in CeedOperatorAssembleSymbolicSingle()
561 …Note: If the value of `assembled` or `rstr` passed to this function are non-`NULL` , then it is as…
562 …    These objects will be destroyed if `*assembled` or `*rstr` is the only reference to the object.
582   // Determine if fallback parent or operator has implementation  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
584   if (op_fallback_parent && use_parent && op_fallback_parent->LinearAssembleQFunctionUpdate) {  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
585     // -- Backend version for op fallback parent is faster, if it exists  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
589   } else if (op->LinearAssembleQFunctionUpdate) {  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
596   if (LinearAssembleQFunctionUpdate) {  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
605     if (data_is_setup) {  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
610 …if (update_needed) CeedCall(LinearAssembleQFunctionUpdate(op_assemble, assembled_vec, assembled_rs…  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
628 …if (op_fallback) CeedCall(CeedOperatorLinearAssembleQFunctionBuildOrUpdate(op_fallback, assembled,…  in CeedOperatorLinearAssembleQFunctionBuildOrUpdate_Core()
641 …Note: If the value of `assembled` or `rstr` passed to this function are non-`NULL` , then it is as…
642 …    These objects will be destroyed if `*assembled` or `*rstr` is the only reference to the object.
682     if (num_elem == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorAssembleSingle()
685   if (op->LinearAssembleSingle) {  in CeedOperatorAssembleSingle()
695     if (op_fallback) {  in CeedOperatorAssembleSingle()
749   if (basis_in == CEED_BASIS_NONE) num_qpts_in = elem_size_in;  in CeedOperatorAssembleSingle()
753   if (elem_rstr_type_in == CEED_RESTRICTION_ORIENTED) {  in CeedOperatorAssembleSingle()
755   } else if (elem_rstr_type_in == CEED_RESTRICTION_CURL_ORIENTED) {  in CeedOperatorAssembleSingle()
759   if (elem_rstr_in != elem_rstr_out) {  in CeedOperatorAssembleSingle()
767     if (basis_out == CEED_BASIS_NONE) num_qpts_out = elem_size_out;  in CeedOperatorAssembleSingle()
775     if (elem_rstr_type_out == CEED_RESTRICTION_ORIENTED) {  in CeedOperatorAssembleSingle()
777     } else if (elem_rstr_type_out == CEED_RESTRICTION_CURL_ORIENTED) {  in CeedOperatorAssembleSingle()
799 …if (elem_rstr_curl_orients_in || elem_rstr_curl_orients_out) CeedCall(CeedCalloc(elem_size_out * e…  in CeedOperatorAssembleSingle()
825         if (contract) {  in CeedOperatorAssembleSingle()
836         // Transform the element matrix if required  in CeedOperatorAssembleSingle()
837         if (elem_rstr_orients_out) {  in CeedOperatorAssembleSingle()
847         } else if (elem_rstr_curl_orients_out) {  in CeedOperatorAssembleSingle()
860         if (elem_rstr_orients_in) {  in CeedOperatorAssembleSingle()
868         } else if (elem_rstr_curl_orients_in) {  in CeedOperatorAssembleSingle()
899   if (elem_rstr_type_in == CEED_RESTRICTION_ORIENTED) {  in CeedOperatorAssembleSingle()
901   } else if (elem_rstr_type_in == CEED_RESTRICTION_CURL_ORIENTED) {  in CeedOperatorAssembleSingle()
904   if (elem_rstr_in != elem_rstr_out) {  in CeedOperatorAssembleSingle()
905     if (elem_rstr_type_out == CEED_RESTRICTION_ORIENTED) {  in CeedOperatorAssembleSingle()
907     } else if (elem_rstr_type_out == CEED_RESTRICTION_CURL_ORIENTED) {  in CeedOperatorAssembleSingle()
940   if (rstr_in != rstr_out) {  in CeedOperatorAssemblyCountEntriesSingle()
975   if (is_composite) {  in CeedOperatorLinearAssembleGetNumEntries()
999 …@param[in]  p_mult_fine  L-vector multiplicity in parallel gather/scatter, or `NULL` if not creati…
1002 …@param[in]  basis_c_to_f `CeedBasis` for coarse to fine interpolation, or `NULL` if not creating p…
1032     if (is_at_points) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1055     if (vec == CEED_VECTOR_ACTIVE) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1058       if (!rstr_fine) CeedCall(CeedOperatorFieldGetElemRestriction(input_fields[i], &rstr_fine));  in CeedOperatorMultigridLevelCreateSingle_Core()
1063     if (dim == 0) CeedCall(CeedBasisGetDimension(basis, &dim));  in CeedOperatorMultigridLevelCreateSingle_Core()
1078     if (vec == CEED_VECTOR_ACTIVE) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1081       if (!rstr_fine) CeedCall(CeedOperatorFieldGetElemRestriction(output_fields[i], &rstr_fine));  in CeedOperatorMultigridLevelCreateSingle_Core()
1086     if (dim == 0) CeedCall(CeedBasisGetDimension(basis, &dim));  in CeedOperatorMultigridLevelCreateSingle_Core()
1104   if (op_restrict || op_prolong) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1131     if (use_scalar_mult) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1173   // Check that coarse to fine basis is provided if prolong/restrict operators are requested  in CeedOperatorMultigridLevelCreateSingle_Core()
1181   if (op_restrict) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1219   if (op_prolong) {  in CeedOperatorMultigridLevelCreateSingle_Core()
1377   if (comp_stride != 1) shift *= num_comp;  in CeedOperatorCreateActivePointBlockRestriction()
1404   if (!op->qf_assembled) {  in CeedOperatorGetQFunctionAssemblyData()
1465   @param[in]     needs_data_update Boolean flag indicating if update is needed or completed
1477   @brief Determine if `CeedQFunctionAssemblyData` needs update
1480   @param[out] is_update_needed Boolean flag indicating if re-assembly is required
1496 …Note: If the value of ` *data_copy` passed to this function is non-`NULL` , then it is assumed tha…
1497 …This `CeedQFunctionAssemblyData` will be destroyed if ` *data_copy` is the only reference to this …
1576   if (!*data || --(*data)->ref_count > 0) {  in CeedQFunctionAssemblyDataDestroy()
1599   if (!op->op_assembled) {  in CeedOperatorGetOperatorAssemblyData()
1654     if (vec == CEED_VECTOR_ACTIVE) {  in CeedOperatorAssemblyDataCreate()
1664         if ((*data)->active_bases_in[i] == basis_in) index = i;  in CeedOperatorAssemblyDataCreate()
1666       if (index == -1) {  in CeedOperatorAssemblyDataCreate()
1688       if (eval_mode != CEED_EVAL_WEIGHT) {  in CeedOperatorAssemblyDataCreate()
1689         // q_comp = 1 if CEED_EVAL_NONE, CEED_EVAL_WEIGHT caught by QF Assembly  in CeedOperatorAssemblyDataCreate()
1712     if (vec == CEED_VECTOR_ACTIVE) {  in CeedOperatorAssemblyDataCreate()
1722         if ((*data)->active_bases_out[i] == basis_out) index = i;  in CeedOperatorAssemblyDataCreate()
1724       if (index == -1) {  in CeedOperatorAssemblyDataCreate()
1746       if (eval_mode != CEED_EVAL_WEIGHT) {  in CeedOperatorAssemblyDataCreate()
1747         // q_comp = 1 if CEED_EVAL_NONE, CEED_EVAL_WEIGHT caught by QF Assembly  in CeedOperatorAssemblyDataCreate()
1799   if (num_active_bases_in) *num_active_bases_in = data->num_active_bases_in;  in CeedOperatorAssemblyDataGetEvalModes()
1800   if (num_eval_modes_in) *num_eval_modes_in = data->num_eval_modes_in;  in CeedOperatorAssemblyDataGetEvalModes()
1801   if (eval_modes_in) *eval_modes_in = (const CeedEvalMode **)data->eval_modes_in;  in CeedOperatorAssemblyDataGetEvalModes()
1802   if (eval_mode_offsets_in) *eval_mode_offsets_in = data->eval_mode_offsets_in;  in CeedOperatorAssemblyDataGetEvalModes()
1803   if (num_active_bases_out) *num_active_bases_out = data->num_active_bases_out;  in CeedOperatorAssemblyDataGetEvalModes()
1804   if (num_eval_modes_out) *num_eval_modes_out = data->num_eval_modes_out;  in CeedOperatorAssemblyDataGetEvalModes()
1805   if (eval_modes_out) *eval_modes_out = (const CeedEvalMode **)data->eval_modes_out;  in CeedOperatorAssemblyDataGetEvalModes()
1806   if (eval_mode_offsets_out) *eval_mode_offsets_out = data->eval_mode_offsets_out;  in CeedOperatorAssemblyDataGetEvalModes()
1807   if (num_output_components) *num_output_components = data->num_output_components;  in CeedOperatorAssemblyDataGetEvalModes()
1831   // Assemble B_in, B_out if needed  in CeedOperatorAssemblyDataGetBases()
1832   if (assembled_bases_in && !data->assembled_bases_in[0]) {  in CeedOperatorAssemblyDataGetBases()
1835 …if (data->active_bases_in[0] == CEED_BASIS_NONE) CeedCall(CeedElemRestrictionGetElementSize(data->…  in CeedOperatorAssemblyDataGetBases()
1848       if (has_eval_none) {  in CeedOperatorAssemblyDataGetBases()
1866             if (q_comp_in > 1) {  in CeedOperatorAssemblyDataGetBases()
1867               if (e_in == 0 || data->eval_modes_in[b][e_in] != eval_mode_in_prev) d_in = 0;  in CeedOperatorAssemblyDataGetBases()
1875       if (identity) CeedCall(CeedFree(&identity));  in CeedOperatorAssemblyDataGetBases()
1880   if (assembled_bases_out && !data->assembled_bases_out[0]) {  in CeedOperatorAssemblyDataGetBases()
1883 …if (data->active_bases_out[0] == CEED_BASIS_NONE) CeedCall(CeedElemRestrictionGetElementSize(data-…  in CeedOperatorAssemblyDataGetBases()
1896       if (has_eval_none) {  in CeedOperatorAssemblyDataGetBases()
1914             if (q_comp_out > 1) {  in CeedOperatorAssemblyDataGetBases()
1915               if (e_out == 0 || data->eval_modes_out[b][e_out] != eval_mode_out_prev) d_out = 0;  in CeedOperatorAssemblyDataGetBases()
1923       if (identity) CeedCall(CeedFree(&identity));  in CeedOperatorAssemblyDataGetBases()
1929   if (num_active_bases_in) *num_active_bases_in = data->num_active_bases_in;  in CeedOperatorAssemblyDataGetBases()
1930   if (active_bases_in) *active_bases_in = data->active_bases_in;  in CeedOperatorAssemblyDataGetBases()
1931   if (assembled_bases_in) *assembled_bases_in = (const CeedScalar **)data->assembled_bases_in;  in CeedOperatorAssemblyDataGetBases()
1932   if (num_active_bases_out) *num_active_bases_out = data->num_active_bases_out;  in CeedOperatorAssemblyDataGetBases()
1933   if (active_bases_out) *active_bases_out = data->active_bases_out;  in CeedOperatorAssemblyDataGetBases()
1934   if (assembled_bases_out) *assembled_bases_out = (const CeedScalar **)data->assembled_bases_out;  in CeedOperatorAssemblyDataGetBases()
1956   if (num_active_elem_rstrs_in) *num_active_elem_rstrs_in = data->num_active_bases_in;  in CeedOperatorAssemblyDataGetElemRestrictions()
1957   if (active_elem_rstrs_in) *active_elem_rstrs_in = data->active_elem_rstrs_in;  in CeedOperatorAssemblyDataGetElemRestrictions()
1958   if (num_active_elem_rstrs_out) *num_active_elem_rstrs_out = data->num_active_bases_out;  in CeedOperatorAssemblyDataGetElemRestrictions()
1959   if (active_elem_rstrs_out) *active_elem_rstrs_out = data->active_elem_rstrs_out;  in CeedOperatorAssemblyDataGetElemRestrictions()
1973   if (!*data) {  in CeedOperatorAssemblyDataDestroy()
2020   // Create if needed  in CeedOperatorGetFallback()
2021   if (!op->op_fallback) CeedCall(CeedOperatorCreateFallback(op));  in CeedOperatorGetFallback()
2022   if (op->op_fallback) {  in CeedOperatorGetFallback()
2028     if (is_debug) {  in CeedOperatorGetFallback()
2075 …if (op->op_fallback_parent) CeedCall(CeedReferenceCopy(CeedOperatorReturnCeed(op->op_fallback_pare…  in CeedOperatorGetFallbackParentCeed()
2111   if (op->LinearAssembleQFunction) {  in CeedOperatorLinearAssembleQFunction()
2120 …if (op_fallback) CeedCall(CeedOperatorLinearAssembleQFunction(op_fallback, assembled, rstr, reques…  in CeedOperatorLinearAssembleQFunction()
2133 …Note: If the value of `assembled` or `rstr` passed to this function are non-`NULL` , then it is as…
2134 …    These objects will be destroyed if `*assembled` or `*rstr` is the only reference to the object.
2177   if (!is_composite) {  in CeedOperatorLinearAssembleDiagonal()
2181     if (num_elem == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorLinearAssembleDiagonal()
2184   if (op->LinearAssembleDiagonal) {  in CeedOperatorLinearAssembleDiagonal()
2188   } else if (op->LinearAssembleAddDiagonal) {  in CeedOperatorLinearAssembleDiagonal()
2193   } else if (is_composite) {  in CeedOperatorLinearAssembleDiagonal()
2204     if (op_fallback) {  in CeedOperatorLinearAssembleDiagonal()
2243   if (!is_composite) {  in CeedOperatorLinearAssembleAddDiagonal()
2247     if (num_elem == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorLinearAssembleAddDiagonal()
2250   if (op->LinearAssembleAddDiagonal) {  in CeedOperatorLinearAssembleAddDiagonal()
2254   } else if (is_composite) {  in CeedOperatorLinearAssembleAddDiagonal()
2264     if (op_fallback) {  in CeedOperatorLinearAssembleAddDiagonal()
2305   if (is_composite) {  in CeedOperatorLinearAssemblePointBlockDiagonalSymbolic()
2408   if (!is_composite) {  in CeedOperatorLinearAssemblePointBlockDiagonal()
2412     if (num_elem == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorLinearAssemblePointBlockDiagonal()
2415   if (op->LinearAssemblePointBlockDiagonal) {  in CeedOperatorLinearAssemblePointBlockDiagonal()
2419   } else if (op->LinearAssembleAddPointBlockDiagonal) {  in CeedOperatorLinearAssemblePointBlockDiagonal()
2430     if (op_fallback) {  in CeedOperatorLinearAssemblePointBlockDiagonal()
2471   if (!is_composite) {  in CeedOperatorLinearAssembleAddPointBlockDiagonal()
2475     if (num_elem == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorLinearAssembleAddPointBlockDiagonal()
2478   if (op->LinearAssembleAddPointBlockDiagonal) {  in CeedOperatorLinearAssembleAddPointBlockDiagonal()
2488     if (op_fallback) {  in CeedOperatorLinearAssembleAddPointBlockDiagonal()
2494   if (is_composite) {  in CeedOperatorLinearAssembleAddPointBlockDiagonal()
2531   if (op->LinearAssembleSymbolic) {  in CeedOperatorLinearAssembleSymbolic()
2541     if (op_fallback) {  in CeedOperatorLinearAssembleSymbolic()
2555   if (is_composite) {  in CeedOperatorLinearAssembleSymbolic()
2597   if (!is_composite) {  in CeedOperatorLinearAssemble()
2601     if (num_elem == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorLinearAssemble()
2604   if (op->LinearAssemble) {  in CeedOperatorLinearAssemble()
2608   } else if (is_composite) {  in CeedOperatorLinearAssemble()
2619   } else if (op->LinearAssembleSingle) {  in CeedOperatorLinearAssemble()
2629     if (op_fallback) {  in CeedOperatorLinearAssemble()
2635   // Default to interface version if non-composite and no fallback  in CeedOperatorLinearAssemble()
2671   if (num_suboperators == 0) return CEED_ERROR_SUCCESS;  in CeedOperatorCompositeGetMultiplicity()
2689       if (skip_indices[j] == i) continue;  in CeedOperatorCompositeGetMultiplicity()
2703       if (sub_mult_array[j] > 0.0) mult_array[j] += 1.0;  in CeedOperatorCompositeGetMultiplicity()
2721 …@param[in]  p_mult_fine  L-vector multiplicity in parallel gather/scatter, or `NULL` if not creati…
2738   // Build prolongation matrix, if required  in CeedOperatorMultigridLevelCreate()
2739   if (op_prolong || op_restrict) {  in CeedOperatorMultigridLevelCreate()
2759 …@param[in]  p_mult_fine   L-vector multiplicity in parallel gather/scatter, or `NULL` if not creat…
2762 …@param[in]  interp_c_to_f Matrix for coarse to fine interpolation, or `NULL` if not creating prolo…
2790   // Create coarse to fine basis, if required  in CeedOperatorMultigridLevelCreateTensorH1()
2791   if (op_prolong || op_restrict) {  in CeedOperatorMultigridLevelCreateTensorH1()
2795     // Check if interpolation matrix is provided  in CeedOperatorMultigridLevelCreateTensorH1()
2826 …@param[in]  p_mult_fine   L-vector multiplicity in parallel gather/scatter, or `NULL` if not creat…
2829 …@param[in]  interp_c_to_f Matrix for coarse to fine interpolation, or `NULL` if not creating prolo…
2855   if (op_prolong || op_restrict) {  in CeedOperatorMultigridLevelCreateH1()
2860     // Check if interpolation matrix is provided  in CeedOperatorMultigridLevelCreateH1()
2920   if (op->CreateFDMElementInverse) {  in CeedOperatorCreateFDMElementInverse()
2930     if (op_fallback) {  in CeedOperatorCreateFDMElementInverse()
2948     if (vec == CEED_VECTOR_ACTIVE) {  in CeedOperatorCreateFDMElementInverse()
2954       if (!basis) CeedCall(CeedOperatorFieldGetBasis(op_fields[i], &basis));  in CeedOperatorCreateFDMElementInverse()
2955       if (!rstr) CeedCall(CeedOperatorFieldGetElemRestriction(op_fields[i], &rstr));  in CeedOperatorCreateFDMElementInverse()
3017 …       if (fabs(assembled_array[q * layout[0] + i * layout[1] + e * layout[2]]) > qf_value_bound) {  in CeedOperatorCreateFDMElementInverse()
3023       if (count) {  in CeedOperatorCreateFDMElementInverse()
3043         if (interp) fdm_diagonal[c * num_nodes + n] = 1.0;  in CeedOperatorCreateFDMElementInverse()
3044         if (grad) {  in CeedOperatorCreateFDMElementInverse()
3050 …if (fabs(fdm_diagonal[c * num_nodes + n]) < fdm_diagonal_bound) fdm_diagonal[c * num_nodes + n] = …  in CeedOperatorCreateFDMElementInverse()