I’ve written a large numeric simulation that benefits substantially from starting Julia with `--math-mode=fast`

(79% speed-up). I know this flag can be dangerous, but up until now, everything has worked well.

I’d now like to add Bayesian Optimization to tune a few parameters. Unfortunately, the Gaussian Process MLE breaks when I use explicit parameter bounds with fast-math enabled. Performance is not important for the Gaussian Process – I just need it to work.

Is there a way to disable fast-math for *only* the call to the Bayesian Optimization function?

I’ve tried switching to explicit `@fastmath`

in the obvious parts of my code, but I still get a 43% speed-up by using the global `--math-mode=fast`

flag.