Migdal et al, Equation (51)

Percentage Accurate: 99.5% → 99.5%
Time: 9.4s
Alternatives: 12
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 12 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.5% 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.5% accurate, 0.6× speedup?

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

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

    \[\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.6%

    \[\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. Add Preprocessing

Alternative 2: 99.5% accurate, 0.6× speedup?

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

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

    \[\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(\sqrt{k}\right)}^{-1} \cdot \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) \]
  6. Add Preprocessing

Alternative 3: 99.5% accurate, 1.0× speedup?

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

\\
\frac{{\left(\sqrt{\left(2 \cdot n\right) \cdot \mathsf{PI}\left(\right)}\right)}^{\left(1 - k\right)}}{\sqrt{k}}
\end{array}
Derivation
  1. Initial program 99.3%

    \[\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.6%

    \[\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. Step-by-step derivation
    1. lift-/.f64N/A

      \[\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)}^{\frac{1}{2}}}} \]
    2. lift-*.f64N/A

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

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

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

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

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

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

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

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

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

Alternative 4: 98.1% accurate, 1.1× speedup?

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;k \leq 1:\\ \;\;\;\;\frac{\sqrt{n \cdot \mathsf{PI}\left(\right)}}{\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 (* n (PI))) (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{n \cdot \mathsf{PI}\left(\right)}}{\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 98.8%

      \[\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.f6472.0

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

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

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

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

            \[\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. lower-*.f6499.2

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

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

              \[\leadsto \color{blue}{{\left(\left(2 \cdot \mathsf{PI}\left(\right)\right) \cdot n\right)}^{\left(\frac{-1}{2} \cdot k\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}{2} \cdot k\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}{2} \cdot k\right)}}{\sqrt{k}}} \]
            5. lower-/.f6499.2

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

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

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

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

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

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

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

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

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

            \[\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. Add Preprocessing

        Alternative 5: 99.5% accurate, 1.1× speedup?

        \[\begin{array}{l} \\ \frac{{\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}} \end{array} \]
        (FPCore (k n)
         :precision binary64
         (/ (pow (* (PI) (* n 2.0)) (fma k -0.5 0.5)) (sqrt k)))
        \begin{array}{l}
        
        \\
        \frac{{\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}}
        \end{array}
        
        Derivation
        1. Initial program 99.3%

          \[\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.6%

          \[\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. Step-by-step derivation
          1. lift-/.f64N/A

            \[\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)}^{\frac{1}{2}}}} \]
          2. lift-*.f64N/A

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

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

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

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

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

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

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

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

          \[\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.4

            \[\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.4%

          \[\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. Add Preprocessing

        Alternative 6: 49.8% accurate, 3.6× speedup?

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

          \[\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.f6440.6

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

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

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

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

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

              Alternative 7: 49.7% 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.3%

                \[\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.f6440.6

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

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

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

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

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

                    Alternative 8: 49.7% accurate, 3.6× speedup?

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

                      \[\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.f6440.6

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

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

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

                      Alternative 9: 38.4% 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.3%

                        \[\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.f6440.6

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

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

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

                        Alternative 10: 38.4% accurate, 4.8× speedup?

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

                          \[\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.f6440.6

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

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

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

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

                            Alternative 11: 38.4% accurate, 4.8× speedup?

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

                              \[\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.f6440.6

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

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

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

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

                                Alternative 12: 38.3% accurate, 4.8× speedup?

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

                                  \[\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.f6440.6

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

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

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

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

                                    Reproduce

                                    ?
                                    herbie shell --seed 2024340 
                                    (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))))