def METHOD_NAME(x, y, z):
   yfs = []
   if opts.twoD_example:
      yfs += [- x, - y, y - 1.0, - y - 0.5 + 0.5 * x]
   if opts.twoD_example_alternative:
      yfs += [y - 1.0, - y - 0.5 + 0.5 * x, - x, - y]
   if opts.dp:
      yfs += [y + x * dp_sinphi - dp_c * dp_cosphi, x - dp_t, -x - dp_tc]
   if opts.tensile:
      yfs += [x - opts.tensile_strength, y - opts.tensile_strength, z - opts.tensile_strength]
   if opts.mc:
      yfs += [0.5 * (x - z) + 0.5 * (x + z) * sinphi - cohcos,
              0.5 * (y - z) + 0.5 * (y + z) * sinphi - cohcos,
              0.5 * (x - y) + 0.5 * (x + y) * sinphi - cohcos,
              0.5 * (y - x) + 0.5 * (x + y) * sinphi - cohcos,
              0.5 * (z - y) + 0.5 * (y + z) * sinphi - cohcos,
              0.5 * (z - x) + 0.5 * (x + z) * sinphi - cohcos]
   yf = yfs[0]
   for i in range(1, len(yfs)):
      yf = yield_function_2(yf, yfs[i])
   return yf