We can imrove convergence and rule out some of looping behaviors during Halley iterations by bracketing root like it's done in boost,