Migdal et al, Equation (51)

Percentage Accurate: 99.4% → 99.6%
Time: 9.8s
Alternatives: 9
Speedup: 1.1×

Specification

?
\[\begin{array}{l} \\ \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \end{array} \]
(FPCore (k n)
 :precision binary64
 (* (/ 1.0 (sqrt k)) (pow (* (* 2.0 (PI)) n) (/ (- 1.0 k) 2.0))))
\begin{array}{l}

\\
\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}
\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 9 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: 99.4% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \end{array} \]
(FPCore (k n)
 :precision binary64
 (* (/ 1.0 (sqrt k)) (pow (* (* 2.0 (PI)) n) (/ (- 1.0 k) 2.0))))
\begin{array}{l}

\\
\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}
\end{array}

Alternative 1: 99.6% accurate, 0.6× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\\ \frac{\sqrt{t\_0}}{{\left({t\_0}^{k}\right)}^{0.5} \cdot \sqrt{k}} \end{array} \end{array} \]
(FPCore (k n)
 :precision binary64
 (let* ((t_0 (* (* 2.0 (PI)) n)))
   (/ (sqrt t_0) (* (pow (pow t_0 k) 0.5) (sqrt k)))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\\
\frac{\sqrt{t\_0}}{{\left({t\_0}^{k}\right)}^{0.5} \cdot \sqrt{k}}
\end{array}
\end{array}
Derivation
  1. Initial program 99.5%

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

      \[\leadsto \color{blue}{\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}} \]
    2. *-commutativeN/A

      \[\leadsto \color{blue}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \cdot \frac{1}{\sqrt{k}}} \]
    3. lift-/.f64N/A

      \[\leadsto {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \cdot \color{blue}{\frac{1}{\sqrt{k}}} \]
    4. un-div-invN/A

      \[\leadsto \color{blue}{\frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}}} \]
    5. lift-pow.f64N/A

      \[\leadsto \frac{\color{blue}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}}}{\sqrt{k}} \]
    6. lift-/.f64N/A

      \[\leadsto \frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1 - k}{2}\right)}}}{\sqrt{k}} \]
    7. lift--.f64N/A

      \[\leadsto \frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{\color{blue}{1 - k}}{2}\right)}}{\sqrt{k}} \]
    8. div-subN/A

      \[\leadsto \frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1}{2} - \frac{k}{2}\right)}}}{\sqrt{k}} \]
    9. metadata-evalN/A

      \[\leadsto \frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\color{blue}{\frac{1}{2}} - \frac{k}{2}\right)}}{\sqrt{k}} \]
    10. pow-subN/A

      \[\leadsto \frac{\color{blue}{\frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{k}{2}\right)}}}}{\sqrt{k}} \]
    11. associate-/l/N/A

      \[\leadsto \color{blue}{\frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}}{\sqrt{k} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{k}{2}\right)}}} \]
    12. lower-/.f64N/A

      \[\leadsto \color{blue}{\frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}}{\sqrt{k} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{k}{2}\right)}}} \]
  4. Applied rewrites99.7%

    \[\leadsto \color{blue}{\frac{\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}}{\sqrt{k} \cdot {\left({\left(n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)\right)}^{k}\right)}^{0.5}}} \]
  5. Final simplification99.7%

    \[\leadsto \frac{\sqrt{\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n}}{{\left({\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{k}\right)}^{0.5} \cdot \sqrt{k}} \]
  6. Add Preprocessing

Alternative 2: 99.5% accurate, 0.9× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\\ \left({t\_0}^{\left(-0.5 \cdot k\right)} \cdot \sqrt{t\_0}\right) \cdot \frac{1}{\sqrt{k}} \end{array} \end{array} \]
(FPCore (k n)
 :precision binary64
 (let* ((t_0 (* (* 2.0 (PI)) n)))
   (* (* (pow t_0 (* -0.5 k)) (sqrt t_0)) (/ 1.0 (sqrt k)))))
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\\
\left({t\_0}^{\left(-0.5 \cdot k\right)} \cdot \sqrt{t\_0}\right) \cdot \frac{1}{\sqrt{k}}
\end{array}
\end{array}
Derivation
  1. Initial program 99.5%

    \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-pow.f64N/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}} \]
    2. lift-/.f64N/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1 - k}{2}\right)}} \]
    3. lift--.f64N/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{\color{blue}{1 - k}}{2}\right)} \]
    4. div-subN/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1}{2} - \frac{k}{2}\right)}} \]
    5. metadata-evalN/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\color{blue}{\frac{1}{2}} - \frac{k}{2}\right)} \]
    6. sub-negN/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1}{2} + \left(\mathsf{neg}\left(\frac{k}{2}\right)\right)\right)}} \]
    7. +-commutativeN/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\left(\mathsf{neg}\left(\frac{k}{2}\right)\right) + \frac{1}{2}\right)}} \]
    8. unpow-prod-upN/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{\left({\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\mathsf{neg}\left(\frac{k}{2}\right)\right)} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}\right)} \]
    9. lower-*.f64N/A

      \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{\left({\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\mathsf{neg}\left(\frac{k}{2}\right)\right)} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}\right)} \]
  4. Applied rewrites99.6%

    \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{\left({\left(n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)\right)}^{\left(k \cdot -0.5\right)} \cdot \sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)} \]
  5. Final simplification99.6%

    \[\leadsto \left({\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(-0.5 \cdot k\right)} \cdot \sqrt{\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n}\right) \cdot \frac{1}{\sqrt{k}} \]
  6. Add Preprocessing

Alternative 3: 98.3% accurate, 1.1× speedup?

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;k \leq 1:\\ \;\;\;\;\frac{\sqrt{\mathsf{PI}\left(\right) \cdot n}}{\sqrt{0.5 \cdot k}}\\ \mathbf{else}:\\ \;\;\;\;\frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(-0.5 \cdot k\right)}}{\sqrt{k}}\\ \end{array} \end{array} \]
(FPCore (k n)
 :precision binary64
 (if (<= k 1.0)
   (/ (sqrt (* (PI) n)) (sqrt (* 0.5 k)))
   (/ (pow (* (* 2.0 n) (PI)) (* -0.5 k)) (sqrt k))))
\begin{array}{l}

\\
\begin{array}{l}
\mathbf{if}\;k \leq 1:\\
\;\;\;\;\frac{\sqrt{\mathsf{PI}\left(\right) \cdot n}}{\sqrt{0.5 \cdot k}}\\

\mathbf{else}:\\
\;\;\;\;\frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(-0.5 \cdot k\right)}}{\sqrt{k}}\\


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

    1. Initial program 99.0%

      \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
    2. Add Preprocessing
    3. Taylor expanded in k around 0

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

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

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

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

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

        \[\leadsto \sqrt{2} \cdot \sqrt{\color{blue}{\frac{n \cdot \mathsf{PI}\left(\right)}{k}}} \]
      6. *-commutativeN/A

        \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
      7. lower-*.f64N/A

        \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
      8. lower-PI.f6474.7

        \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right)} \cdot n}{k}} \]
    5. Applied rewrites74.7%

      \[\leadsto \color{blue}{\sqrt{2} \cdot \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k}}} \]
    6. Step-by-step derivation
      1. Applied rewrites75.0%

        \[\leadsto \color{blue}{\sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2}} \]
      2. Step-by-step derivation
        1. Applied rewrites75.0%

          \[\leadsto \sqrt{\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot \frac{n}{k}} \]
        2. Step-by-step derivation
          1. Applied rewrites96.1%

            \[\leadsto \frac{\sqrt{n \cdot \mathsf{PI}\left(\right)}}{\color{blue}{\sqrt{0.5 \cdot k}}} \]

          if 1 < k

          1. Initial program 100.0%

            \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
          2. Add Preprocessing
          3. Taylor expanded in k around inf

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

              \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(k \cdot \frac{-1}{2}\right)}} \]
            2. lower-*.f64100.0

              \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(k \cdot -0.5\right)}} \]
          5. Applied rewrites100.0%

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(k \cdot -0.5\right)}} \]
          6. Step-by-step derivation
            1. lift-*.f64N/A

              \[\leadsto \color{blue}{\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(k \cdot \frac{-1}{2}\right)}} \]
            2. *-commutativeN/A

              \[\leadsto \color{blue}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(k \cdot \frac{-1}{2}\right)} \cdot \frac{1}{\sqrt{k}}} \]
            3. lift-/.f64N/A

              \[\leadsto {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(k \cdot \frac{-1}{2}\right)} \cdot \color{blue}{\frac{1}{\sqrt{k}}} \]
            4. un-div-invN/A

              \[\leadsto \color{blue}{\frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}}} \]
            5. lower-/.f64100.0

              \[\leadsto \color{blue}{\frac{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(k \cdot -0.5\right)}}{\sqrt{k}}} \]
            6. lift-*.f64N/A

              \[\leadsto \frac{{\color{blue}{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            7. *-commutativeN/A

              \[\leadsto \frac{{\color{blue}{\left(n \cdot \left(2 \cdot \mathsf{PI}\left(\right)\right)\right)}}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            8. lift-*.f64N/A

              \[\leadsto \frac{{\left(n \cdot \color{blue}{\left(2 \cdot \mathsf{PI}\left(\right)\right)}\right)}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            9. associate-*r*N/A

              \[\leadsto \frac{{\color{blue}{\left(\left(n \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right)}}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            10. lift-*.f64N/A

              \[\leadsto \frac{{\left(\color{blue}{\left(n \cdot 2\right)} \cdot \mathsf{PI}\left(\right)\right)}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            11. lower-*.f64100.0

              \[\leadsto \frac{{\color{blue}{\left(\left(n \cdot 2\right) \cdot \mathsf{PI}\left(\right)\right)}}^{\left(k \cdot -0.5\right)}}{\sqrt{k}} \]
            12. lift-*.f64N/A

              \[\leadsto \frac{{\left(\color{blue}{\left(n \cdot 2\right)} \cdot \mathsf{PI}\left(\right)\right)}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            13. *-commutativeN/A

              \[\leadsto \frac{{\left(\color{blue}{\left(2 \cdot n\right)} \cdot \mathsf{PI}\left(\right)\right)}^{\left(k \cdot \frac{-1}{2}\right)}}{\sqrt{k}} \]
            14. lower-*.f64100.0

              \[\leadsto \frac{{\left(\color{blue}{\left(2 \cdot n\right)} \cdot \mathsf{PI}\left(\right)\right)}^{\left(k \cdot -0.5\right)}}{\sqrt{k}} \]
          7. Applied rewrites100.0%

            \[\leadsto \color{blue}{\frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(-0.5 \cdot k\right)}}{\sqrt{k}}} \]
        3. Recombined 2 regimes into one program.
        4. Final simplification98.1%

          \[\leadsto \begin{array}{l} \mathbf{if}\;k \leq 1:\\ \;\;\;\;\frac{\sqrt{\mathsf{PI}\left(\right) \cdot n}}{\sqrt{0.5 \cdot k}}\\ \mathbf{else}:\\ \;\;\;\;\frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(-0.5 \cdot k\right)}}{\sqrt{k}}\\ \end{array} \]
        5. Add Preprocessing

        Alternative 4: 99.5% accurate, 1.1× speedup?

        \[\begin{array}{l} \\ \frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(\mathsf{fma}\left(k, -0.5, 0.5\right)\right)}}{\sqrt{k}} \end{array} \]
        (FPCore (k n)
         :precision binary64
         (/ (pow (* (* 2.0 n) (PI)) (fma k -0.5 0.5)) (sqrt k)))
        \begin{array}{l}
        
        \\
        \frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(\mathsf{fma}\left(k, -0.5, 0.5\right)\right)}}{\sqrt{k}}
        \end{array}
        
        Derivation
        1. Initial program 99.5%

          \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
        2. Add Preprocessing
        3. Step-by-step derivation
          1. lift-pow.f64N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}} \]
          2. lift-/.f64N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1 - k}{2}\right)}} \]
          3. clear-numN/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1}{\frac{2}{1 - k}}\right)}} \]
          4. associate-/r/N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\color{blue}{\left(\frac{1}{2} \cdot \left(1 - k\right)\right)}} \]
          5. metadata-evalN/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\color{blue}{\frac{1}{2}} \cdot \left(1 - k\right)\right)} \]
          6. pow-unpowN/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{{\left({\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}\right)}^{\left(1 - k\right)}} \]
          7. lower-pow.f64N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{{\left({\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\frac{1}{2}}\right)}^{\left(1 - k\right)}} \]
          8. unpow1/2N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\color{blue}{\left(\sqrt{\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n}\right)}}^{\left(1 - k\right)} \]
          9. lower-sqrt.f6499.5

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\color{blue}{\left(\sqrt{\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n}\right)}}^{\left(1 - k\right)} \]
          10. lift-*.f64N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\sqrt{\color{blue}{\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n}}\right)}^{\left(1 - k\right)} \]
          11. *-commutativeN/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\sqrt{\color{blue}{n \cdot \left(2 \cdot \mathsf{PI}\left(\right)\right)}}\right)}^{\left(1 - k\right)} \]
          12. lower-*.f6499.5

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\sqrt{\color{blue}{n \cdot \left(2 \cdot \mathsf{PI}\left(\right)\right)}}\right)}^{\left(1 - k\right)} \]
          13. lift-*.f64N/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\sqrt{n \cdot \color{blue}{\left(2 \cdot \mathsf{PI}\left(\right)\right)}}\right)}^{\left(1 - k\right)} \]
          14. *-commutativeN/A

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\sqrt{n \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot 2\right)}}\right)}^{\left(1 - k\right)} \]
          15. lower-*.f6499.5

            \[\leadsto \frac{1}{\sqrt{k}} \cdot {\left(\sqrt{n \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot 2\right)}}\right)}^{\left(1 - k\right)} \]
        4. Applied rewrites99.5%

          \[\leadsto \frac{1}{\sqrt{k}} \cdot \color{blue}{{\left(\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)}^{\left(1 - k\right)}} \]
        5. Step-by-step derivation
          1. lift-*.f64N/A

            \[\leadsto \color{blue}{\frac{1}{\sqrt{k}} \cdot {\left(\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)}^{\left(1 - k\right)}} \]
          2. *-commutativeN/A

            \[\leadsto \color{blue}{{\left(\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)}^{\left(1 - k\right)} \cdot \frac{1}{\sqrt{k}}} \]
          3. lift-/.f64N/A

            \[\leadsto {\left(\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)}^{\left(1 - k\right)} \cdot \color{blue}{\frac{1}{\sqrt{k}}} \]
          4. un-div-invN/A

            \[\leadsto \color{blue}{\frac{{\left(\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}}} \]
          5. lower-/.f6499.6

            \[\leadsto \color{blue}{\frac{{\left(\sqrt{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}}} \]
          6. lift-*.f64N/A

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{n \cdot \left(\mathsf{PI}\left(\right) \cdot 2\right)}}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          7. lift-*.f64N/A

            \[\leadsto \frac{{\left(\sqrt{n \cdot \color{blue}{\left(\mathsf{PI}\left(\right) \cdot 2\right)}}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          8. *-commutativeN/A

            \[\leadsto \frac{{\left(\sqrt{n \cdot \color{blue}{\left(2 \cdot \mathsf{PI}\left(\right)\right)}}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          9. associate-*r*N/A

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{\left(n \cdot 2\right) \cdot \mathsf{PI}\left(\right)}}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          10. lift-*.f64N/A

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{\left(n \cdot 2\right)} \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          11. lower-*.f6499.6

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{\left(n \cdot 2\right) \cdot \mathsf{PI}\left(\right)}}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          12. lift-*.f64N/A

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{\left(n \cdot 2\right)} \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          13. *-commutativeN/A

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{\left(2 \cdot n\right)} \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
          14. lower-*.f6499.6

            \[\leadsto \frac{{\left(\sqrt{\color{blue}{\left(2 \cdot n\right)} \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}} \]
        6. Applied rewrites99.6%

          \[\leadsto \color{blue}{\frac{{\left(\sqrt{\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}}} \]
        7. Step-by-step derivation
          1. lift-pow.f64N/A

            \[\leadsto \frac{\color{blue}{{\left(\sqrt{\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}}{\sqrt{k}} \]
          2. lift-sqrt.f64N/A

            \[\leadsto \frac{{\color{blue}{\left(\sqrt{\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)}\right)}}^{\left(1 - k\right)}}{\sqrt{k}} \]
          3. sqrt-pow2N/A

            \[\leadsto \frac{\color{blue}{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(\frac{1 - k}{2}\right)}}}{\sqrt{k}} \]
          4. lift--.f64N/A

            \[\leadsto \frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(\frac{\color{blue}{1 - k}}{2}\right)}}{\sqrt{k}} \]
          5. lower-pow.f64N/A

            \[\leadsto \frac{\color{blue}{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(\frac{1 - k}{2}\right)}}}{\sqrt{k}} \]
          6. lift-*.f64N/A

            \[\leadsto \frac{{\color{blue}{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \]
          7. *-commutativeN/A

            \[\leadsto \frac{{\color{blue}{\left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot n\right)\right)}}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \]
          8. lower-*.f64N/A

            \[\leadsto \frac{{\color{blue}{\left(\mathsf{PI}\left(\right) \cdot \left(2 \cdot n\right)\right)}}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \]
          9. lift-*.f64N/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(2 \cdot n\right)}\right)}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \]
          10. *-commutativeN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(n \cdot 2\right)}\right)}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \]
          11. lower-*.f64N/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \color{blue}{\left(n \cdot 2\right)}\right)}^{\left(\frac{1 - k}{2}\right)}}{\sqrt{k}} \]
          12. lift--.f64N/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{\color{blue}{1 - k}}{2}\right)}}{\sqrt{k}} \]
          13. clear-numN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\color{blue}{\left(\frac{1}{\frac{2}{1 - k}}\right)}}}{\sqrt{k}} \]
          14. associate-/r/N/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\color{blue}{\left(\frac{1}{2} \cdot \left(1 - k\right)\right)}}}{\sqrt{k}} \]
          15. metadata-evalN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\color{blue}{\frac{1}{2}} \cdot \left(1 - k\right)\right)}}{\sqrt{k}} \]
          16. lift--.f64N/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} \cdot \color{blue}{\left(1 - k\right)}\right)}}{\sqrt{k}} \]
          17. sub-negN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} \cdot \color{blue}{\left(1 + \left(\mathsf{neg}\left(k\right)\right)\right)}\right)}}{\sqrt{k}} \]
          18. +-commutativeN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} \cdot \color{blue}{\left(\left(\mathsf{neg}\left(k\right)\right) + 1\right)}\right)}}{\sqrt{k}} \]
          19. distribute-rgt-inN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\color{blue}{\left(\left(\mathsf{neg}\left(k\right)\right) \cdot \frac{1}{2} + 1 \cdot \frac{1}{2}\right)}}}{\sqrt{k}} \]
          20. distribute-lft-neg-inN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\color{blue}{\left(\mathsf{neg}\left(k \cdot \frac{1}{2}\right)\right)} + 1 \cdot \frac{1}{2}\right)}}{\sqrt{k}} \]
          21. distribute-rgt-neg-inN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\color{blue}{k \cdot \left(\mathsf{neg}\left(\frac{1}{2}\right)\right)} + 1 \cdot \frac{1}{2}\right)}}{\sqrt{k}} \]
          22. metadata-evalN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(k \cdot \color{blue}{\frac{-1}{2}} + 1 \cdot \frac{1}{2}\right)}}{\sqrt{k}} \]
          23. metadata-evalN/A

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(k \cdot \frac{-1}{2} + \color{blue}{\frac{1}{2}}\right)}}{\sqrt{k}} \]
          24. lower-fma.f6499.6

            \[\leadsto \frac{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\color{blue}{\left(\mathsf{fma}\left(k, -0.5, 0.5\right)\right)}}}{\sqrt{k}} \]
        8. Applied rewrites99.6%

          \[\leadsto \frac{\color{blue}{{\left(\mathsf{PI}\left(\right) \cdot \left(n \cdot 2\right)\right)}^{\left(\mathsf{fma}\left(k, -0.5, 0.5\right)\right)}}}{\sqrt{k}} \]
        9. Final simplification99.6%

          \[\leadsto \frac{{\left(\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)\right)}^{\left(\mathsf{fma}\left(k, -0.5, 0.5\right)\right)}}{\sqrt{k}} \]
        10. Add Preprocessing

        Alternative 5: 50.0% accurate, 3.6× speedup?

        \[\begin{array}{l} \\ \frac{\sqrt{\mathsf{PI}\left(\right) \cdot n}}{\sqrt{0.5 \cdot k}} \end{array} \]
        (FPCore (k n) :precision binary64 (/ (sqrt (* (PI) n)) (sqrt (* 0.5 k))))
        \begin{array}{l}
        
        \\
        \frac{\sqrt{\mathsf{PI}\left(\right) \cdot n}}{\sqrt{0.5 \cdot k}}
        \end{array}
        
        Derivation
        1. Initial program 99.5%

          \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
        2. Add Preprocessing
        3. Taylor expanded in k around 0

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

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

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

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

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

            \[\leadsto \sqrt{2} \cdot \sqrt{\color{blue}{\frac{n \cdot \mathsf{PI}\left(\right)}{k}}} \]
          6. *-commutativeN/A

            \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
          7. lower-*.f64N/A

            \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
          8. lower-PI.f6438.4

            \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right)} \cdot n}{k}} \]
        5. Applied rewrites38.4%

          \[\leadsto \color{blue}{\sqrt{2} \cdot \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k}}} \]
        6. Step-by-step derivation
          1. Applied rewrites38.5%

            \[\leadsto \color{blue}{\sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2}} \]
          2. Step-by-step derivation
            1. Applied rewrites38.5%

              \[\leadsto \sqrt{\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot \frac{n}{k}} \]
            2. Step-by-step derivation
              1. Applied rewrites49.1%

                \[\leadsto \frac{\sqrt{n \cdot \mathsf{PI}\left(\right)}}{\color{blue}{\sqrt{0.5 \cdot k}}} \]
              2. Final simplification49.1%

                \[\leadsto \frac{\sqrt{\mathsf{PI}\left(\right) \cdot n}}{\sqrt{0.5 \cdot k}} \]
              3. Add Preprocessing

              Alternative 6: 50.0% accurate, 3.6× speedup?

              \[\begin{array}{l} \\ \sqrt{\frac{2 \cdot \mathsf{PI}\left(\right)}{k}} \cdot \sqrt{n} \end{array} \]
              (FPCore (k n) :precision binary64 (* (sqrt (/ (* 2.0 (PI)) k)) (sqrt n)))
              \begin{array}{l}
              
              \\
              \sqrt{\frac{2 \cdot \mathsf{PI}\left(\right)}{k}} \cdot \sqrt{n}
              \end{array}
              
              Derivation
              1. Initial program 99.5%

                \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
              2. Add Preprocessing
              3. Taylor expanded in k around 0

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

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

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

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

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

                  \[\leadsto \sqrt{2} \cdot \sqrt{\color{blue}{\frac{n \cdot \mathsf{PI}\left(\right)}{k}}} \]
                6. *-commutativeN/A

                  \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                7. lower-*.f64N/A

                  \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                8. lower-PI.f6438.4

                  \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right)} \cdot n}{k}} \]
              5. Applied rewrites38.4%

                \[\leadsto \color{blue}{\sqrt{2} \cdot \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k}}} \]
              6. Step-by-step derivation
                1. Applied rewrites48.7%

                  \[\leadsto \sqrt{n} \cdot \color{blue}{\sqrt{\frac{2 \cdot \mathsf{PI}\left(\right)}{k}}} \]
                2. Final simplification48.7%

                  \[\leadsto \sqrt{\frac{2 \cdot \mathsf{PI}\left(\right)}{k}} \cdot \sqrt{n} \]
                3. Add Preprocessing

                Alternative 7: 49.9% accurate, 3.6× speedup?

                \[\begin{array}{l} \\ \sqrt{\frac{2}{k} \cdot \mathsf{PI}\left(\right)} \cdot \sqrt{n} \end{array} \]
                (FPCore (k n) :precision binary64 (* (sqrt (* (/ 2.0 k) (PI))) (sqrt n)))
                \begin{array}{l}
                
                \\
                \sqrt{\frac{2}{k} \cdot \mathsf{PI}\left(\right)} \cdot \sqrt{n}
                \end{array}
                
                Derivation
                1. Initial program 99.5%

                  \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
                2. Add Preprocessing
                3. Taylor expanded in k around 0

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

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

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

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

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

                    \[\leadsto \sqrt{2} \cdot \sqrt{\color{blue}{\frac{n \cdot \mathsf{PI}\left(\right)}{k}}} \]
                  6. *-commutativeN/A

                    \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                  7. lower-*.f64N/A

                    \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                  8. lower-PI.f6438.4

                    \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right)} \cdot n}{k}} \]
                5. Applied rewrites38.4%

                  \[\leadsto \color{blue}{\sqrt{2} \cdot \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k}}} \]
                6. Step-by-step derivation
                  1. Applied rewrites38.5%

                    \[\leadsto \color{blue}{\sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2}} \]
                  2. Step-by-step derivation
                    1. Applied rewrites38.5%

                      \[\leadsto \sqrt{\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot \frac{n}{k}} \]
                    2. Step-by-step derivation
                      1. Applied rewrites48.7%

                        \[\leadsto \sqrt{n} \cdot \color{blue}{\sqrt{\frac{2}{k} \cdot \mathsf{PI}\left(\right)}} \]
                      2. Final simplification48.7%

                        \[\leadsto \sqrt{\frac{2}{k} \cdot \mathsf{PI}\left(\right)} \cdot \sqrt{n} \]
                      3. Add Preprocessing

                      Alternative 8: 38.3% accurate, 4.8× speedup?

                      \[\begin{array}{l} \\ \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2} \end{array} \]
                      (FPCore (k n) :precision binary64 (sqrt (* (/ (* (PI) n) k) 2.0)))
                      \begin{array}{l}
                      
                      \\
                      \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2}
                      \end{array}
                      
                      Derivation
                      1. Initial program 99.5%

                        \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
                      2. Add Preprocessing
                      3. Taylor expanded in k around 0

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

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

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

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

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

                          \[\leadsto \sqrt{2} \cdot \sqrt{\color{blue}{\frac{n \cdot \mathsf{PI}\left(\right)}{k}}} \]
                        6. *-commutativeN/A

                          \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                        7. lower-*.f64N/A

                          \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                        8. lower-PI.f6438.4

                          \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right)} \cdot n}{k}} \]
                      5. Applied rewrites38.4%

                        \[\leadsto \color{blue}{\sqrt{2} \cdot \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k}}} \]
                      6. Step-by-step derivation
                        1. Applied rewrites38.5%

                          \[\leadsto \color{blue}{\sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2}} \]
                        2. Add Preprocessing

                        Alternative 9: 38.3% accurate, 4.8× speedup?

                        \[\begin{array}{l} \\ \sqrt{\frac{n}{k} \cdot \left(2 \cdot \mathsf{PI}\left(\right)\right)} \end{array} \]
                        (FPCore (k n) :precision binary64 (sqrt (* (/ n k) (* 2.0 (PI)))))
                        \begin{array}{l}
                        
                        \\
                        \sqrt{\frac{n}{k} \cdot \left(2 \cdot \mathsf{PI}\left(\right)\right)}
                        \end{array}
                        
                        Derivation
                        1. Initial program 99.5%

                          \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)} \]
                        2. Add Preprocessing
                        3. Taylor expanded in k around 0

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

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

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

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

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

                            \[\leadsto \sqrt{2} \cdot \sqrt{\color{blue}{\frac{n \cdot \mathsf{PI}\left(\right)}{k}}} \]
                          6. *-commutativeN/A

                            \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                          7. lower-*.f64N/A

                            \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right) \cdot n}}{k}} \]
                          8. lower-PI.f6438.4

                            \[\leadsto \sqrt{2} \cdot \sqrt{\frac{\color{blue}{\mathsf{PI}\left(\right)} \cdot n}{k}} \]
                        5. Applied rewrites38.4%

                          \[\leadsto \color{blue}{\sqrt{2} \cdot \sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k}}} \]
                        6. Step-by-step derivation
                          1. Applied rewrites38.5%

                            \[\leadsto \color{blue}{\sqrt{\frac{\mathsf{PI}\left(\right) \cdot n}{k} \cdot 2}} \]
                          2. Step-by-step derivation
                            1. Applied rewrites38.5%

                              \[\leadsto \sqrt{\left(\mathsf{PI}\left(\right) \cdot 2\right) \cdot \frac{n}{k}} \]
                            2. Final simplification38.5%

                              \[\leadsto \sqrt{\frac{n}{k} \cdot \left(2 \cdot \mathsf{PI}\left(\right)\right)} \]
                            3. Add Preprocessing

                            Reproduce

                            ?
                            herbie shell --seed 2024241 
                            (FPCore (k n)
                              :name "Migdal et al, Equation (51)"
                              :precision binary64
                              (* (/ 1.0 (sqrt k)) (pow (* (* 2.0 (PI)) n) (/ (- 1.0 k) 2.0))))