\frac{x \cdot \left(\frac{x}{y} + 1\right)}{x + 1}\left(x \cdot \left(\sqrt[3]{\frac{x}{y} + 1} \cdot \frac{\sqrt[3]{\frac{x}{y} + 1}}{\sqrt[3]{x + 1} \cdot \sqrt[3]{x + 1}}\right)\right) \cdot \frac{\sqrt[3]{\frac{x}{y} + 1}}{\sqrt[3]{x + 1}}double code(double x, double y) {
return (((double) (x * ((double) ((x / y) + 1.0)))) / ((double) (x + 1.0)));
}
double code(double x, double y) {
return ((double) (((double) (x * ((double) (((double) cbrt(((double) ((x / y) + 1.0)))) * (((double) cbrt(((double) ((x / y) + 1.0)))) / ((double) (((double) cbrt(((double) (x + 1.0)))) * ((double) cbrt(((double) (x + 1.0))))))))))) * (((double) cbrt(((double) ((x / y) + 1.0)))) / ((double) cbrt(((double) (x + 1.0)))))));
}




Bits error versus x




Bits error versus y
Results
| Original | 9.2 |
|---|---|
| Target | 0.1 |
| Herbie | 0.9 |
Initial program 9.2
Simplified0.1
rmApplied add-cube-cbrt0.6
Applied add-cube-cbrt0.9
Applied times-frac0.9
Applied associate-*r*0.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2020196
(FPCore (x y)
:name "Codec.Picture.Types:toneMapping from JuicyPixels-3.2.6.1"
:precision binary64
:herbie-target
(* (/ x 1.0) (/ (+ (/ x y) 1.0) (+ x 1.0)))
(/ (* x (+ (/ x y) 1.0)) (+ x 1.0)))