optionally (and by default) use bernstein expansion to compute bounds