Make optimisers optimal

[config-i1]

Example. When bounds="a" in es() function, the optimal smoothing parameters differ from the case of bounds="u", and usually the value of CF in the former case is larger than in the latter case.
This means that the surface of CF is complicated and we need to find fast and reliable methods of getting to optimal for different types of bounds.