VandenBroeck and Keller, Equation (6)

Percentage Accurate: 76.9% → 99.3%
Time: 6.7s
Alternatives: 8
Speedup: 4.4×

Specification

?
\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot \ell\\ t\_0 - \frac{1}{F \cdot F} \cdot \tan t\_0 \end{array} \end{array} \]
(FPCore (F l)
 :precision binary64
 (let* ((t_0 (* (PI) l))) (- t_0 (* (/ 1.0 (* F F)) (tan t_0)))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \mathsf{PI}\left(\right) \cdot \ell\\
t\_0 - \frac{1}{F \cdot F} \cdot \tan t\_0
\end{array}
\end{array}

Sampling outcomes in binary64 precision:

Local Percentage Accuracy vs ?

The average percentage accuracy by input value. Horizontal axis shows value of an input variable; the variable is choosen in the title. Vertical axis is accuracy; higher is better. Red represent the original program, while blue represents Herbie's suggestion. These can be toggled with buttons below the plot. The line is an average while dots represent individual samples.

Accuracy vs Speed?

Herbie found 8 alternatives:

AlternativeAccuracySpeedup
The accuracy (vertical axis) and speed (horizontal axis) of each alternatives. Up and to the right is better. The red square shows the initial program, and each blue circle shows an alternative.The line shows the best available speed-accuracy tradeoffs.

Initial Program: 76.9% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot \ell\\ t\_0 - \frac{1}{F \cdot F} \cdot \tan t\_0 \end{array} \end{array} \]
(FPCore (F l)
 :precision binary64
 (let* ((t_0 (* (PI) l))) (- t_0 (* (/ 1.0 (* F F)) (tan t_0)))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \mathsf{PI}\left(\right) \cdot \ell\\
t\_0 - \frac{1}{F \cdot F} \cdot \tan t\_0
\end{array}
\end{array}

Alternative 1: 99.3% accurate, 0.6× speedup?

\[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ l\_s \cdot \begin{array}{l} \mathbf{if}\;l\_m \leq 10^{+15}:\\ \;\;\;\;\mathsf{PI}\left(\right) \cdot l\_m - \frac{\frac{\tan \left(l\_m \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}\\ \mathbf{else}:\\ \;\;\;\;{\mathsf{PI}\left(\right)}^{0.25} \cdot \left({\mathsf{PI}\left(\right)}^{0.75} \cdot l\_m\right)\\ \end{array} \end{array} \]
l\_m = (fabs.f64 l)
l\_s = (copysign.f64 #s(literal 1 binary64) l)
(FPCore (l_s F l_m)
 :precision binary64
 (*
  l_s
  (if (<= l_m 1e+15)
    (- (* (PI) l_m) (/ (/ (tan (* l_m (PI))) F) F))
    (* (pow (PI) 0.25) (* (pow (PI) 0.75) l_m)))))
\begin{array}{l}
l\_m = \left|\ell\right|
\\
l\_s = \mathsf{copysign}\left(1, \ell\right)

\\
l\_s \cdot \begin{array}{l}
\mathbf{if}\;l\_m \leq 10^{+15}:\\
\;\;\;\;\mathsf{PI}\left(\right) \cdot l\_m - \frac{\frac{\tan \left(l\_m \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}\\

\mathbf{else}:\\
\;\;\;\;{\mathsf{PI}\left(\right)}^{0.25} \cdot \left({\mathsf{PI}\left(\right)}^{0.75} \cdot l\_m\right)\\


\end{array}
\end{array}
Derivation
  1. Split input into 2 regimes
  2. if l < 1e15

    1. Initial program 83.7%

      \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
    2. Add Preprocessing
    3. Step-by-step derivation
      1. lift-*.f64N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
      2. *-commutativeN/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \frac{1}{F \cdot F}} \]
      3. lift-/.f64N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \color{blue}{\frac{1}{F \cdot F}} \]
      4. associate-*r/N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot 1}{F \cdot F}} \]
      5. *-rgt-identityN/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F \cdot F} \]
      6. lift-*.f64N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{\color{blue}{F \cdot F}} \]
      7. associate-/r*N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
      8. lower-/.f64N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
      9. lower-/.f6490.0

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}}{F} \]
      10. lift-*.f64N/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F}}{F} \]
      11. *-commutativeN/A

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
      12. lower-*.f6490.0

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
    4. Applied rewrites90.0%

      \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\ell \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}} \]

    if 1e15 < l

    1. Initial program 63.8%

      \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
    2. Add Preprocessing
    3. Taylor expanded in F around inf

      \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
    4. Step-by-step derivation
      1. *-commutativeN/A

        \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      2. lower-*.f64N/A

        \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      3. lower-PI.f6499.5

        \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
    5. Applied rewrites99.5%

      \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
    6. Step-by-step derivation
      1. Applied rewrites99.6%

        \[\leadsto {\mathsf{PI}\left(\right)}^{0.25} \cdot \color{blue}{\left({\mathsf{PI}\left(\right)}^{0.75} \cdot \ell\right)} \]
    7. Recombined 2 regimes into one program.
    8. Add Preprocessing

    Alternative 2: 99.4% accurate, 1.0× speedup?

    \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot l\_m\\ l\_s \cdot \begin{array}{l} \mathbf{if}\;l\_m \leq 10^{+15}:\\ \;\;\;\;t\_0 - \frac{\frac{\tan \left(l\_m \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \end{array} \]
    l\_m = (fabs.f64 l)
    l\_s = (copysign.f64 #s(literal 1 binary64) l)
    (FPCore (l_s F l_m)
     :precision binary64
     (let* ((t_0 (* (PI) l_m)))
       (* l_s (if (<= l_m 1e+15) (- t_0 (/ (/ (tan (* l_m (PI))) F) F)) t_0))))
    \begin{array}{l}
    l\_m = \left|\ell\right|
    \\
    l\_s = \mathsf{copysign}\left(1, \ell\right)
    
    \\
    \begin{array}{l}
    t_0 := \mathsf{PI}\left(\right) \cdot l\_m\\
    l\_s \cdot \begin{array}{l}
    \mathbf{if}\;l\_m \leq 10^{+15}:\\
    \;\;\;\;t\_0 - \frac{\frac{\tan \left(l\_m \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}\\
    
    \mathbf{else}:\\
    \;\;\;\;t\_0\\
    
    
    \end{array}
    \end{array}
    \end{array}
    
    Derivation
    1. Split input into 2 regimes
    2. if l < 1e15

      1. Initial program 83.7%

        \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
        2. *-commutativeN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \frac{1}{F \cdot F}} \]
        3. lift-/.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \color{blue}{\frac{1}{F \cdot F}} \]
        4. associate-*r/N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot 1}{F \cdot F}} \]
        5. *-rgt-identityN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F \cdot F} \]
        6. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{\color{blue}{F \cdot F}} \]
        7. associate-/r*N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
        8. lower-/.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
        9. lower-/.f6490.0

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}}{F} \]
        10. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F}}{F} \]
        11. *-commutativeN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
        12. lower-*.f6490.0

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
      4. Applied rewrites90.0%

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\ell \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}} \]

      if 1e15 < l

      1. Initial program 63.8%

        \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
      2. Add Preprocessing
      3. Taylor expanded in F around inf

        \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
      4. Step-by-step derivation
        1. *-commutativeN/A

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
        2. lower-*.f64N/A

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
        3. lower-PI.f6499.5

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
      5. Applied rewrites99.5%

        \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
    3. Recombined 2 regimes into one program.
    4. Add Preprocessing

    Alternative 3: 98.3% accurate, 1.0× speedup?

    \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ l\_s \cdot \begin{array}{l} \mathbf{if}\;l\_m \leq 0.5:\\ \;\;\;\;\mathsf{fma}\left({F}^{-1}, \mathsf{PI}\left(\right) \cdot \frac{l\_m}{-F}, l\_m \cdot \mathsf{PI}\left(\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{PI}\left(\right) \cdot l\_m\\ \end{array} \end{array} \]
    l\_m = (fabs.f64 l)
    l\_s = (copysign.f64 #s(literal 1 binary64) l)
    (FPCore (l_s F l_m)
     :precision binary64
     (*
      l_s
      (if (<= l_m 0.5)
        (fma (pow F -1.0) (* (PI) (/ l_m (- F))) (* l_m (PI)))
        (* (PI) l_m))))
    \begin{array}{l}
    l\_m = \left|\ell\right|
    \\
    l\_s = \mathsf{copysign}\left(1, \ell\right)
    
    \\
    l\_s \cdot \begin{array}{l}
    \mathbf{if}\;l\_m \leq 0.5:\\
    \;\;\;\;\mathsf{fma}\left({F}^{-1}, \mathsf{PI}\left(\right) \cdot \frac{l\_m}{-F}, l\_m \cdot \mathsf{PI}\left(\right)\right)\\
    
    \mathbf{else}:\\
    \;\;\;\;\mathsf{PI}\left(\right) \cdot l\_m\\
    
    
    \end{array}
    \end{array}
    
    Derivation
    1. Split input into 2 regimes
    2. if l < 0.5

      1. Initial program 83.6%

        \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
        2. *-commutativeN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \frac{1}{F \cdot F}} \]
        3. lift-/.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \color{blue}{\frac{1}{F \cdot F}} \]
        4. associate-*r/N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot 1}{F \cdot F}} \]
        5. *-rgt-identityN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F \cdot F} \]
        6. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{\color{blue}{F \cdot F}} \]
        7. associate-/r*N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
        8. lower-/.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
        9. lower-/.f6490.0

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}}{F} \]
        10. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F}}{F} \]
        11. *-commutativeN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
        12. lower-*.f6490.0

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
      4. Applied rewrites90.0%

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\ell \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}} \]
      5. Applied rewrites90.0%

        \[\leadsto \color{blue}{\mathsf{fma}\left({F}^{-1}, \frac{\tan \left(\ell \cdot \mathsf{PI}\left(\right)\right)}{-F}, \ell \cdot \mathsf{PI}\left(\right)\right)} \]
      6. Taylor expanded in l around 0

        \[\leadsto \mathsf{fma}\left({F}^{-1}, \color{blue}{-1 \cdot \frac{\ell \cdot \mathsf{PI}\left(\right)}{F}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
      7. Step-by-step derivation
        1. mul-1-negN/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \color{blue}{\mathsf{neg}\left(\frac{\ell \cdot \mathsf{PI}\left(\right)}{F}\right)}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        2. associate-/l*N/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \mathsf{neg}\left(\color{blue}{\ell \cdot \frac{\mathsf{PI}\left(\right)}{F}}\right), \ell \cdot \mathsf{PI}\left(\right)\right) \]
        3. distribute-rgt-neg-inN/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \color{blue}{\ell \cdot \left(\mathsf{neg}\left(\frac{\mathsf{PI}\left(\right)}{F}\right)\right)}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        4. mul-1-negN/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \color{blue}{\left(-1 \cdot \frac{\mathsf{PI}\left(\right)}{F}\right)}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        5. lower-*.f64N/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \color{blue}{\ell \cdot \left(-1 \cdot \frac{\mathsf{PI}\left(\right)}{F}\right)}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        6. mul-1-negN/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \color{blue}{\left(\mathsf{neg}\left(\frac{\mathsf{PI}\left(\right)}{F}\right)\right)}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        7. distribute-neg-frac2N/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \color{blue}{\frac{\mathsf{PI}\left(\right)}{\mathsf{neg}\left(F\right)}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        8. mul-1-negN/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \frac{\mathsf{PI}\left(\right)}{\color{blue}{-1 \cdot F}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        9. lower-/.f64N/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \color{blue}{\frac{\mathsf{PI}\left(\right)}{-1 \cdot F}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        10. lower-PI.f64N/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \frac{\color{blue}{\mathsf{PI}\left(\right)}}{-1 \cdot F}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        11. mul-1-negN/A

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \frac{\mathsf{PI}\left(\right)}{\color{blue}{\mathsf{neg}\left(F\right)}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
        12. lower-neg.f6483.3

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \ell \cdot \frac{\mathsf{PI}\left(\right)}{\color{blue}{-F}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
      8. Applied rewrites83.3%

        \[\leadsto \mathsf{fma}\left({F}^{-1}, \color{blue}{\ell \cdot \frac{\mathsf{PI}\left(\right)}{-F}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]
      9. Step-by-step derivation
        1. Applied rewrites83.3%

          \[\leadsto \mathsf{fma}\left({F}^{-1}, \mathsf{PI}\left(\right) \cdot \color{blue}{\frac{\ell}{-F}}, \ell \cdot \mathsf{PI}\left(\right)\right) \]

        if 0.5 < l

        1. Initial program 65.4%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Taylor expanded in F around inf

          \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
        4. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          2. lower-*.f64N/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          3. lower-PI.f6498.1

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
        5. Applied rewrites98.1%

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      10. Recombined 2 regimes into one program.
      11. Add Preprocessing

      Alternative 4: 98.3% accurate, 3.2× speedup?

      \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot l\_m\\ l\_s \cdot \begin{array}{l} \mathbf{if}\;l\_m \leq 0.5:\\ \;\;\;\;t\_0 - \frac{\frac{\mathsf{PI}\left(\right)}{F} \cdot l\_m}{F}\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \end{array} \]
      l\_m = (fabs.f64 l)
      l\_s = (copysign.f64 #s(literal 1 binary64) l)
      (FPCore (l_s F l_m)
       :precision binary64
       (let* ((t_0 (* (PI) l_m)))
         (* l_s (if (<= l_m 0.5) (- t_0 (/ (* (/ (PI) F) l_m) F)) t_0))))
      \begin{array}{l}
      l\_m = \left|\ell\right|
      \\
      l\_s = \mathsf{copysign}\left(1, \ell\right)
      
      \\
      \begin{array}{l}
      t_0 := \mathsf{PI}\left(\right) \cdot l\_m\\
      l\_s \cdot \begin{array}{l}
      \mathbf{if}\;l\_m \leq 0.5:\\
      \;\;\;\;t\_0 - \frac{\frac{\mathsf{PI}\left(\right)}{F} \cdot l\_m}{F}\\
      
      \mathbf{else}:\\
      \;\;\;\;t\_0\\
      
      
      \end{array}
      \end{array}
      \end{array}
      
      Derivation
      1. Split input into 2 regimes
      2. if l < 0.5

        1. Initial program 83.6%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Step-by-step derivation
          1. lift-*.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
          2. *-commutativeN/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \frac{1}{F \cdot F}} \]
          3. lift-/.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \color{blue}{\frac{1}{F \cdot F}} \]
          4. associate-*r/N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot 1}{F \cdot F}} \]
          5. *-rgt-identityN/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F \cdot F} \]
          6. lift-*.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{\color{blue}{F \cdot F}} \]
          7. associate-/r*N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
          8. lower-/.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}{F}} \]
          9. lower-/.f6490.0

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{F}}}{F} \]
          10. lift-*.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\mathsf{PI}\left(\right) \cdot \ell\right)}}{F}}{F} \]
          11. *-commutativeN/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
          12. lower-*.f6490.0

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\tan \color{blue}{\left(\ell \cdot \mathsf{PI}\left(\right)\right)}}{F}}{F} \]
        4. Applied rewrites90.0%

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\frac{\tan \left(\ell \cdot \mathsf{PI}\left(\right)\right)}{F}}{F}} \]
        5. Taylor expanded in l around 0

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\ell \cdot \mathsf{PI}\left(\right)}{F}}}{F} \]
        6. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot \ell}}{F}}{F} \]
          2. associate-*l/N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\mathsf{PI}\left(\right)}{F} \cdot \ell}}{F} \]
          3. lower-*.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\mathsf{PI}\left(\right)}{F} \cdot \ell}}{F} \]
          4. lower-/.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\mathsf{PI}\left(\right)}{F}} \cdot \ell}{F} \]
          5. lower-PI.f6483.3

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\frac{\color{blue}{\mathsf{PI}\left(\right)}}{F} \cdot \ell}{F} \]
        7. Applied rewrites83.3%

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\frac{\mathsf{PI}\left(\right)}{F} \cdot \ell}}{F} \]

        if 0.5 < l

        1. Initial program 65.4%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Taylor expanded in F around inf

          \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
        4. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          2. lower-*.f64N/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          3. lower-PI.f6498.1

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
        5. Applied rewrites98.1%

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      3. Recombined 2 regimes into one program.
      4. Add Preprocessing

      Alternative 5: 92.6% accurate, 3.7× speedup?

      \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ \begin{array}{l} t_0 := \mathsf{PI}\left(\right) \cdot l\_m\\ l\_s \cdot \begin{array}{l} \mathbf{if}\;l\_m \leq 0.5:\\ \;\;\;\;t\_0 - \frac{\mathsf{PI}\left(\right)}{F \cdot F} \cdot l\_m\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \end{array} \]
      l\_m = (fabs.f64 l)
      l\_s = (copysign.f64 #s(literal 1 binary64) l)
      (FPCore (l_s F l_m)
       :precision binary64
       (let* ((t_0 (* (PI) l_m)))
         (* l_s (if (<= l_m 0.5) (- t_0 (* (/ (PI) (* F F)) l_m)) t_0))))
      \begin{array}{l}
      l\_m = \left|\ell\right|
      \\
      l\_s = \mathsf{copysign}\left(1, \ell\right)
      
      \\
      \begin{array}{l}
      t_0 := \mathsf{PI}\left(\right) \cdot l\_m\\
      l\_s \cdot \begin{array}{l}
      \mathbf{if}\;l\_m \leq 0.5:\\
      \;\;\;\;t\_0 - \frac{\mathsf{PI}\left(\right)}{F \cdot F} \cdot l\_m\\
      
      \mathbf{else}:\\
      \;\;\;\;t\_0\\
      
      
      \end{array}
      \end{array}
      \end{array}
      
      Derivation
      1. Split input into 2 regimes
      2. if l < 0.5

        1. Initial program 83.6%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Taylor expanded in l around 0

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\ell \cdot \mathsf{PI}\left(\right)}{{F}^{2}}} \]
        4. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\mathsf{PI}\left(\right) \cdot \ell}}{{F}^{2}} \]
          2. associate-*l/N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{PI}\left(\right)}{{F}^{2}} \cdot \ell} \]
          3. lower-*.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{PI}\left(\right)}{{F}^{2}} \cdot \ell} \]
          4. lower-/.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{PI}\left(\right)}{{F}^{2}}} \cdot \ell \]
          5. lower-PI.f64N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{\mathsf{PI}\left(\right)}}{{F}^{2}} \cdot \ell \]
          6. unpow2N/A

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\mathsf{PI}\left(\right)}{\color{blue}{F \cdot F}} \cdot \ell \]
          7. lower-*.f6477.0

            \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\mathsf{PI}\left(\right)}{\color{blue}{F \cdot F}} \cdot \ell \]
        5. Applied rewrites77.0%

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{PI}\left(\right)}{F \cdot F} \cdot \ell} \]

        if 0.5 < l

        1. Initial program 65.4%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Taylor expanded in F around inf

          \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
        4. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          2. lower-*.f64N/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          3. lower-PI.f6498.1

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
        5. Applied rewrites98.1%

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      3. Recombined 2 regimes into one program.
      4. Add Preprocessing

      Alternative 6: 92.6% accurate, 4.4× speedup?

      \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ l\_s \cdot \begin{array}{l} \mathbf{if}\;l\_m \leq 0.5:\\ \;\;\;\;\left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{F \cdot F}\right) \cdot l\_m\\ \mathbf{else}:\\ \;\;\;\;\mathsf{PI}\left(\right) \cdot l\_m\\ \end{array} \end{array} \]
      l\_m = (fabs.f64 l)
      l\_s = (copysign.f64 #s(literal 1 binary64) l)
      (FPCore (l_s F l_m)
       :precision binary64
       (* l_s (if (<= l_m 0.5) (* (- (PI) (/ (PI) (* F F))) l_m) (* (PI) l_m))))
      \begin{array}{l}
      l\_m = \left|\ell\right|
      \\
      l\_s = \mathsf{copysign}\left(1, \ell\right)
      
      \\
      l\_s \cdot \begin{array}{l}
      \mathbf{if}\;l\_m \leq 0.5:\\
      \;\;\;\;\left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{F \cdot F}\right) \cdot l\_m\\
      
      \mathbf{else}:\\
      \;\;\;\;\mathsf{PI}\left(\right) \cdot l\_m\\
      
      
      \end{array}
      \end{array}
      
      Derivation
      1. Split input into 2 regimes
      2. if l < 0.5

        1. Initial program 83.6%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Taylor expanded in l around 0

          \[\leadsto \color{blue}{\ell \cdot \left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{{F}^{2}}\right)} \]
        4. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \color{blue}{\left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{{F}^{2}}\right) \cdot \ell} \]
          2. lower-*.f64N/A

            \[\leadsto \color{blue}{\left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{{F}^{2}}\right) \cdot \ell} \]
          3. lower--.f64N/A

            \[\leadsto \color{blue}{\left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{{F}^{2}}\right)} \cdot \ell \]
          4. lower-PI.f64N/A

            \[\leadsto \left(\color{blue}{\mathsf{PI}\left(\right)} - \frac{\mathsf{PI}\left(\right)}{{F}^{2}}\right) \cdot \ell \]
          5. lower-/.f64N/A

            \[\leadsto \left(\mathsf{PI}\left(\right) - \color{blue}{\frac{\mathsf{PI}\left(\right)}{{F}^{2}}}\right) \cdot \ell \]
          6. lower-PI.f64N/A

            \[\leadsto \left(\mathsf{PI}\left(\right) - \frac{\color{blue}{\mathsf{PI}\left(\right)}}{{F}^{2}}\right) \cdot \ell \]
          7. unpow2N/A

            \[\leadsto \left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{\color{blue}{F \cdot F}}\right) \cdot \ell \]
          8. lower-*.f6477.0

            \[\leadsto \left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{\color{blue}{F \cdot F}}\right) \cdot \ell \]
        5. Applied rewrites77.0%

          \[\leadsto \color{blue}{\left(\mathsf{PI}\left(\right) - \frac{\mathsf{PI}\left(\right)}{F \cdot F}\right) \cdot \ell} \]

        if 0.5 < l

        1. Initial program 65.4%

          \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        2. Add Preprocessing
        3. Taylor expanded in F around inf

          \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
        4. Step-by-step derivation
          1. *-commutativeN/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          2. lower-*.f64N/A

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
          3. lower-PI.f6498.1

            \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
        5. Applied rewrites98.1%

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      3. Recombined 2 regimes into one program.
      4. Add Preprocessing

      Alternative 7: 73.8% accurate, 22.5× speedup?

      \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ l\_s \cdot \left(\mathsf{PI}\left(\right) \cdot l\_m\right) \end{array} \]
      l\_m = (fabs.f64 l)
      l\_s = (copysign.f64 #s(literal 1 binary64) l)
      (FPCore (l_s F l_m) :precision binary64 (* l_s (* (PI) l_m)))
      \begin{array}{l}
      l\_m = \left|\ell\right|
      \\
      l\_s = \mathsf{copysign}\left(1, \ell\right)
      
      \\
      l\_s \cdot \left(\mathsf{PI}\left(\right) \cdot l\_m\right)
      \end{array}
      
      Derivation
      1. Initial program 78.9%

        \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
      2. Add Preprocessing
      3. Taylor expanded in F around inf

        \[\leadsto \color{blue}{\ell \cdot \mathsf{PI}\left(\right)} \]
      4. Step-by-step derivation
        1. *-commutativeN/A

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
        2. lower-*.f64N/A

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
        3. lower-PI.f6472.0

          \[\leadsto \color{blue}{\mathsf{PI}\left(\right)} \cdot \ell \]
      5. Applied rewrites72.0%

        \[\leadsto \color{blue}{\mathsf{PI}\left(\right) \cdot \ell} \]
      6. Add Preprocessing

      Alternative 8: 3.1% accurate, 135.0× speedup?

      \[\begin{array}{l} l\_m = \left|\ell\right| \\ l\_s = \mathsf{copysign}\left(1, \ell\right) \\ l\_s \cdot 0 \end{array} \]
      l\_m = (fabs.f64 l)
      l\_s = (copysign.f64 #s(literal 1 binary64) l)
      (FPCore (l_s F l_m) :precision binary64 (* l_s 0.0))
      l\_m = fabs(l);
      l\_s = copysign(1.0, l);
      double code(double l_s, double F, double l_m) {
      	return l_s * 0.0;
      }
      
      l\_m =     private
      l\_s =     private
      module fmin_fmax_functions
          implicit none
          private
          public fmax
          public fmin
      
          interface fmax
              module procedure fmax88
              module procedure fmax44
              module procedure fmax84
              module procedure fmax48
          end interface
          interface fmin
              module procedure fmin88
              module procedure fmin44
              module procedure fmin84
              module procedure fmin48
          end interface
      contains
          real(8) function fmax88(x, y) result (res)
              real(8), intent (in) :: x
              real(8), intent (in) :: y
              res = merge(y, merge(x, max(x, y), y /= y), x /= x)
          end function
          real(4) function fmax44(x, y) result (res)
              real(4), intent (in) :: x
              real(4), intent (in) :: y
              res = merge(y, merge(x, max(x, y), y /= y), x /= x)
          end function
          real(8) function fmax84(x, y) result(res)
              real(8), intent (in) :: x
              real(4), intent (in) :: y
              res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
          end function
          real(8) function fmax48(x, y) result(res)
              real(4), intent (in) :: x
              real(8), intent (in) :: y
              res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
          end function
          real(8) function fmin88(x, y) result (res)
              real(8), intent (in) :: x
              real(8), intent (in) :: y
              res = merge(y, merge(x, min(x, y), y /= y), x /= x)
          end function
          real(4) function fmin44(x, y) result (res)
              real(4), intent (in) :: x
              real(4), intent (in) :: y
              res = merge(y, merge(x, min(x, y), y /= y), x /= x)
          end function
          real(8) function fmin84(x, y) result(res)
              real(8), intent (in) :: x
              real(4), intent (in) :: y
              res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
          end function
          real(8) function fmin48(x, y) result(res)
              real(4), intent (in) :: x
              real(8), intent (in) :: y
              res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
          end function
      end module
      
      real(8) function code(l_s, f, l_m)
      use fmin_fmax_functions
          real(8), intent (in) :: l_s
          real(8), intent (in) :: f
          real(8), intent (in) :: l_m
          code = l_s * 0.0d0
      end function
      
      l\_m = Math.abs(l);
      l\_s = Math.copySign(1.0, l);
      public static double code(double l_s, double F, double l_m) {
      	return l_s * 0.0;
      }
      
      l\_m = math.fabs(l)
      l\_s = math.copysign(1.0, l)
      def code(l_s, F, l_m):
      	return l_s * 0.0
      
      l\_m = abs(l)
      l\_s = copysign(1.0, l)
      function code(l_s, F, l_m)
      	return Float64(l_s * 0.0)
      end
      
      l\_m = abs(l);
      l\_s = sign(l) * abs(1.0);
      function tmp = code(l_s, F, l_m)
      	tmp = l_s * 0.0;
      end
      
      l\_m = N[Abs[l], $MachinePrecision]
      l\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[l]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
      code[l$95$s_, F_, l$95$m_] := N[(l$95$s * 0.0), $MachinePrecision]
      
      \begin{array}{l}
      l\_m = \left|\ell\right|
      \\
      l\_s = \mathsf{copysign}\left(1, \ell\right)
      
      \\
      l\_s \cdot 0
      \end{array}
      
      Derivation
      1. Initial program 78.9%

        \[\mathsf{PI}\left(\right) \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-*.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{1}{F \cdot F} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
        2. lift-/.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{1}{F \cdot F}} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        3. frac-2negN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{neg}\left(1\right)}{\mathsf{neg}\left(F \cdot F\right)}} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        4. metadata-evalN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\color{blue}{-1}}{\mathsf{neg}\left(F \cdot F\right)} \cdot \tan \left(\mathsf{PI}\left(\right) \cdot \ell\right) \]
        5. lift-tan.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{-1}{\mathsf{neg}\left(F \cdot F\right)} \cdot \color{blue}{\tan \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
        6. tan-quotN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{-1}{\mathsf{neg}\left(F \cdot F\right)} \cdot \color{blue}{\frac{\sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{\cos \left(\mathsf{PI}\left(\right) \cdot \ell\right)}} \]
        7. frac-timesN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{-1 \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)}{\left(\mathsf{neg}\left(F \cdot F\right)\right) \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \ell\right)}} \]
        8. frac-2negN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{neg}\left(-1 \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)\right)}{\mathsf{neg}\left(\left(\mathsf{neg}\left(F \cdot F\right)\right) \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \ell\right)\right)}} \]
        9. distribute-lft-neg-outN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\mathsf{neg}\left(-1 \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)\right)}{\color{blue}{\left(\mathsf{neg}\left(\left(\mathsf{neg}\left(F \cdot F\right)\right)\right)\right) \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \ell\right)}} \]
        10. remove-double-negN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\mathsf{neg}\left(-1 \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)\right)}{\color{blue}{\left(F \cdot F\right)} \cdot \cos \left(\mathsf{PI}\left(\right) \cdot \ell\right)} \]
        11. *-commutativeN/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \frac{\mathsf{neg}\left(-1 \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)\right)}{\color{blue}{\cos \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \left(F \cdot F\right)}} \]
        12. lower-/.f64N/A

          \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\mathsf{neg}\left(-1 \cdot \sin \left(\mathsf{PI}\left(\right) \cdot \ell\right)\right)}{\cos \left(\mathsf{PI}\left(\right) \cdot \ell\right) \cdot \left(F \cdot F\right)}} \]
      4. Applied rewrites59.3%

        \[\leadsto \mathsf{PI}\left(\right) \cdot \ell - \color{blue}{\frac{\sin \left(\mathsf{fma}\left(-\mathsf{PI}\left(\right), \ell, \mathsf{PI}\left(\right)\right)\right)}{\cos \left(\ell \cdot \mathsf{PI}\left(\right)\right) \cdot \left(F \cdot F\right)}} \]
      5. Taylor expanded in l around 0

        \[\leadsto \color{blue}{-1 \cdot \frac{\sin \mathsf{PI}\left(\right)}{{F}^{2}}} \]
      6. Step-by-step derivation
        1. associate-*r/N/A

          \[\leadsto \color{blue}{\frac{-1 \cdot \sin \mathsf{PI}\left(\right)}{{F}^{2}}} \]
        2. sin-PIN/A

          \[\leadsto \frac{-1 \cdot \color{blue}{0}}{{F}^{2}} \]
        3. metadata-evalN/A

          \[\leadsto \frac{\color{blue}{0}}{{F}^{2}} \]
        4. sin-PIN/A

          \[\leadsto \frac{\color{blue}{\sin \mathsf{PI}\left(\right)}}{{F}^{2}} \]
        5. sin-PIN/A

          \[\leadsto \frac{\color{blue}{0}}{{F}^{2}} \]
        6. div03.1

          \[\leadsto \color{blue}{0} \]
      7. Applied rewrites3.1%

        \[\leadsto \color{blue}{0} \]
      8. Add Preprocessing

      Reproduce

      ?
      herbie shell --seed 2024351 
      (FPCore (F l)
        :name "VandenBroeck and Keller, Equation (6)"
        :precision binary64
        (- (* (PI) l) (* (/ 1.0 (* F F)) (tan (* (PI) l)))))