code: def _relu_complementarity(b, *output_index):
      return mpec.complements(
          b.z[output_index] - b.zhat[output_index] >= 0,
          b.z[output_index] >= 0,
      )
  
imports: import pyomo.environ
  import pyomo.mpec
label_desc: Use image node_id 2 for calling a global function with example usage: _relu_complementarity(b) and returns: mpec