Diagrams.Backend.Cairo.Internal:setTexture from diagrams-cairo-1.3.0.3

Average Error: 11.9 → 1.5

Time: 10.6s

Precision: 64

Math

\[\frac{x \cdot \left(y - z\right)}{y}\]

↓

\[\begin{array}{l} \mathbf{if}\;y \le -1.998165745453842247107774064380461371887 \cdot 10^{-24}:\\ \;\;\;\;x \cdot \frac{y - z}{y}\\ \mathbf{elif}\;y \le 9.436120878493814052713516163817953193206 \cdot 10^{-81}:\\ \;\;\;\;x - \frac{x \cdot z}{y}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{\frac{y}{y - z}}\\ \end{array}\]

C

double f(double x, double y, double z) {
        double r32980258 = x;
        double r32980259 = y;
        double r32980260 = z;
        double r32980261 = r32980259 - r32980260;
        double r32980262 = r32980258 * r32980261;
        double r32980263 = r32980262 / r32980259;
        return r32980263;
}

↓

double f(double x, double y, double z) {
        double r32980264 = y;
        double r32980265 = -1.9981657454538422e-24;
        bool r32980266 = r32980264 <= r32980265;
        double r32980267 = x;
        double r32980268 = z;
        double r32980269 = r32980264 - r32980268;
        double r32980270 = r32980269 / r32980264;
        double r32980271 = r32980267 * r32980270;
        double r32980272 = 9.436120878493814e-81;
        bool r32980273 = r32980264 <= r32980272;
        double r32980274 = r32980267 * r32980268;
        double r32980275 = r32980274 / r32980264;
        double r32980276 = r32980267 - r32980275;
        double r32980277 = r32980264 / r32980269;
        double r32980278 = r32980267 / r32980277;
        double r32980279 = r32980273 ? r32980276 : r32980278;
        double r32980280 = r32980266 ? r32980271 : r32980279;
        return r32980280;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Target

Original	11.9
Target	3.1
Herbie	1.5

\[\begin{array}{l} \mathbf{if}\;z \lt -2.060202331921739024383612783691266533098 \cdot 10^{104}:\\ \;\;\;\;x - \frac{z \cdot x}{y}\\ \mathbf{elif}\;z \lt 1.693976601382852594702773997610248441465 \cdot 10^{213}:\\ \;\;\;\;\frac{x}{\frac{y}{y - z}}\\ \mathbf{else}:\\ \;\;\;\;\left(y - z\right) \cdot \frac{x}{y}\\ \end{array}\]

Derivation

1. Split input into 3 regimes

2. if y < -1.9981657454538422e-24

   1. Initial program 14.8

      \[\frac{x \cdot \left(y - z\right)}{y}\]
   2. Using strategy rm

   3. Applied *-un-lft-identity 14.8

      \[\leadsto \frac{x \cdot \left(y - z\right)}{\color{blue}{1 \cdot y}}\]

   4. Applied times-frac 0.1

      \[\leadsto \color{blue}{\frac{x}{1} \cdot \frac{y - z}{y}}\]

   5. Simplified 0.1

      \[\leadsto \color{blue}{x} \cdot \frac{y - z}{y}\]

   if -1.9981657454538422e-24 < y < 9.436120878493814e-81

   1. Initial program 7.3

      \[\frac{x \cdot \left(y - z\right)}{y}\]
   2. Using strategy rm

   3. Applied associate-/l* 7.9

      \[\leadsto \color{blue}{\frac{x}{\frac{y}{y - z}}}\]

   4. Taylor expanded around 0 4.2

      \[\leadsto \color{blue}{x - \frac{x \cdot z}{y}}\]

   if 9.436120878493814e-81 < y

   1. Initial program 13.6

      \[\frac{x \cdot \left(y - z\right)}{y}\]
   2. Using strategy rm

   3. Applied associate-/l* 0.3

      \[\leadsto \color{blue}{\frac{x}{\frac{y}{y - z}}}\]
3. Recombined 3 regimes into one program.

4. Final simplification 1.5

   \[\leadsto \begin{array}{l} \mathbf{if}\;y \le -1.998165745453842247107774064380461371887 \cdot 10^{-24}:\\ \;\;\;\;x \cdot \frac{y - z}{y}\\ \mathbf{elif}\;y \le 9.436120878493814052713516163817953193206 \cdot 10^{-81}:\\ \;\;\;\;x - \frac{x \cdot z}{y}\\ \mathbf{else}:\\ \;\;\;\;\frac{x}{\frac{y}{y - z}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019179 +o rules:numerics
(FPCore (x y z)
  :name "Diagrams.Backend.Cairo.Internal:setTexture from diagrams-cairo-1.3.0.3"

  :herbie-target
  (if (< z -2.060202331921739e+104) (- x (/ (* z x) y)) (if (< z 1.6939766013828526e+213) (/ x (/ y (- y z))) (* (- y z) (/ x y))))

  (/ (* x (- y z)) y))