def _opbend_transform_mean(rs, fn_low, deriv=<NUM_LIT:0>):
v = <NUM_LIT:0.0><EOL>d = np.zeros((<NUM_LIT:4>,<NUM_LIT:3>), float)<EOL>dd = np.zeros((<NUM_LIT:4>,<NUM_LIT:3>,<NUM_LIT:4>,<NUM_LIT:3>), float)<EOL>for p in np.array([[<NUM_LIT:0>,<NUM_LIT:1>,<NUM_LIT:2>], [<NUM_LIT:2>,<NUM_LIT:0>,<NUM_LIT:1>], [<NUM_LIT:1>,<NUM_LIT:2>,<NUM_LIT:0>]]):<EOL><INDENT>opbend = _opbend_transform([rs[p[<NUM_LIT:0>]], rs[p[<NUM_LIT:1>]], rs[p[<NUM_LIT:2>]], rs[<NUM_LIT:3>]], fn_low, deriv)<EOL>v += opbend[<NUM_LIT:0>]/<NUM_LIT:3><EOL>index0 = np.where(p==<NUM_LIT:0>)[<NUM_LIT:0>][<NUM_LIT:0>] <EOL>index1 = np.where(p==<NUM_LIT:1>)[<NUM_LIT:0>][<NUM_LIT:0>]<EOL>index2 = np.where(p==<NUM_LIT:2>)[<NUM_LIT:0>][<NUM_LIT:0>]<EOL>index3 = <NUM_LIT:3><EOL>if deriv><NUM_LIT:0>:<EOL><INDENT>d[<NUM_LIT:0>] += opbend[<NUM_LIT:1>][index0]/<NUM_LIT:3><EOL>d[<NUM_LIT:1>] += opbend[<NUM_LIT:1>][index1]/<NUM_LIT:3><EOL>d[<NUM_LIT:2>] += opbend[<NUM_LIT:1>][index2]/<NUM_LIT:3><EOL>d[<NUM_LIT:3>] += opbend[<NUM_LIT:1>][index3]/<NUM_LIT:3><EOL><DEDENT>if deriv><NUM_LIT:1>:<EOL><INDENT>dd[<NUM_LIT:0>, :, <NUM_LIT:0>, :] += opbend[<NUM_LIT:2>][index0, :, index0, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:0>, :, <NUM_LIT:1>, :] += opbend[<NUM_LIT:2>][index0, :, index1, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:0>, :, <NUM_LIT:2>, :] += opbend[<NUM_LIT:2>][index0, :, index2, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:0>, :, <NUM_LIT:3>, :] += opbend[<NUM_LIT:2>][index0, :, index3, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:1>, :, <NUM_LIT:0>, :] += opbend[<NUM_LIT:2>][index1, :, index0, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:1>, :, <NUM_LIT:1>, :] += opbend[<NUM_LIT:2>][index1, :, index1, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:1>, :, <NUM_LIT:2>, :] += opbend[<NUM_LIT:2>][index1, :, index2, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:1>, :, <NUM_LIT:3>, :] += opbend[<NUM_LIT:2>][index1, :, index3, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:2>, :, <NUM_LIT:0>, :] += opbend[<NUM_LIT:2>][index2, :, index0, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:2>, :, <NUM_LIT:1>, :] += opbend[<NUM_LIT:2>][index2, :, index1, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:2>, :, <NUM_LIT:2>, :] += opbend[<NUM_LIT:2>][index2, :, index2, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:2>, :, <NUM_LIT:3>, :] += opbend[<NUM_LIT:2>][index2, :, index3, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:3>, :, <NUM_LIT:0>, :] += opbend[<NUM_LIT:2>][index3, :, index0, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:3>, :, <NUM_LIT:1>, :] += opbend[<NUM_LIT:2>][index3, :, index1, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:3>, :, <NUM_LIT:2>, :] += opbend[<NUM_LIT:2>][index3, :, index2, :]/<NUM_LIT:3><EOL>dd[<NUM_LIT:3>, :, <NUM_LIT:3>, :] += opbend[<NUM_LIT:2>][index3, :, index3, :]/<NUM_LIT:3><EOL><DEDENT><DEDENT>if deriv==<NUM_LIT:0>:<EOL><INDENT>return v,<EOL><DEDENT>elif deriv==<NUM_LIT:1>:<EOL><INDENT>return v, d<EOL><DEDENT>elif deriv==<NUM_LIT:2>:<EOL><INDENT>return v, d, dd<EOL><DEDENT>else:<EOL><INDENT>raise ValueError("<STR_LIT>")<EOL><DEDENT>
Compute the mean of the 3 opbends