def onp_fun(*args):
    args = [(x if (jnp.result_type(x) != jnp.bfloat16) else onp.asarray(x, onp.float32)) for x in args]
    return getattr(onp, op)(*args, axis=axis, keepdims=keepdims, interpolation=interpolation)
def onp_fun(*args):
    args = [(x if (jnp.result_type(x) != jnp.bfloat16) else onp.asarray(x, onp.float32)) for x in args]
    return getattr(onp, op)(*args, axis=axis, keepdims=keepdims, interpolation=interpolation)