
(FPCore (u s)
:precision binary32
(let* ((t_0 (/ 1.0 (+ 1.0 (exp (/ (PI) s))))))
(*
(- s)
(log
(-
(/ 1.0 (+ (* u (- (/ 1.0 (+ 1.0 (exp (/ (- (PI)) s)))) t_0)) t_0))
1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{1 + e^{\frac{\mathsf{PI}\left(\right)}{s}}}\\
\left(-s\right) \cdot \log \left(\frac{1}{u \cdot \left(\frac{1}{1 + e^{\frac{-\mathsf{PI}\left(\right)}{s}}} - t\_0\right) + t\_0} - 1\right)
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (u s)
:precision binary32
(let* ((t_0 (/ 1.0 (+ 1.0 (exp (/ (PI) s))))))
(*
(- s)
(log
(-
(/ 1.0 (+ (* u (- (/ 1.0 (+ 1.0 (exp (/ (- (PI)) s)))) t_0)) t_0))
1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{1 + e^{\frac{\mathsf{PI}\left(\right)}{s}}}\\
\left(-s\right) \cdot \log \left(\frac{1}{u \cdot \left(\frac{1}{1 + e^{\frac{-\mathsf{PI}\left(\right)}{s}}} - t\_0\right) + t\_0} - 1\right)
\end{array}
\end{array}
(FPCore (u s)
:precision binary32
(let* ((t_0 (exp (- (log1p (exp (/ (PI) s))))))
(t_1 (fma (- (exp (- (log1p (exp (/ (- (PI)) s))))) t_0) u t_0)))
(* (- s) (log (/ (- (/ 1.0 (pow t_1 2.0)) 1.0) (+ (/ 1.0 t_1) 1.0))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{-\mathsf{log1p}\left(e^{\frac{\mathsf{PI}\left(\right)}{s}}\right)}\\
t_1 := \mathsf{fma}\left(e^{-\mathsf{log1p}\left(e^{\frac{-\mathsf{PI}\left(\right)}{s}}\right)} - t\_0, u, t\_0\right)\\
\left(-s\right) \cdot \log \left(\frac{\frac{1}{{t\_1}^{2}} - 1}{\frac{1}{t\_1} + 1}\right)
\end{array}
\end{array}
Initial program 99.1%
lift--.f32N/A
flip--N/A
lower-/.f32N/A
Applied rewrites99.1%
Taylor expanded in s around 0
lower-/.f32N/A
Applied rewrites99.1%
(FPCore (u s)
:precision binary32
(let* ((t_0 (exp (- (log1p (exp (/ (PI) s))))))
(t_1 (fma (- (exp (- (log1p (exp (/ (- (PI)) s))))) t_0) u t_0)))
(* (- s) (log (/ (- (pow t_1 -2.0) 1.0) (+ (/ 1.0 t_1) 1.0))))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{-\mathsf{log1p}\left(e^{\frac{\mathsf{PI}\left(\right)}{s}}\right)}\\
t_1 := \mathsf{fma}\left(e^{-\mathsf{log1p}\left(e^{\frac{-\mathsf{PI}\left(\right)}{s}}\right)} - t\_0, u, t\_0\right)\\
\left(-s\right) \cdot \log \left(\frac{{t\_1}^{-2} - 1}{\frac{1}{t\_1} + 1}\right)
\end{array}
\end{array}
Initial program 99.1%
lift--.f32N/A
flip--N/A
lower-/.f32N/A
Applied rewrites99.1%
(FPCore (u s)
:precision binary32
(*
(- s)
(log
(-
(/
1.0
(fma
(- (/ 1.0 (+ (exp (/ (- (PI)) s)) 1.0)) (/ 1.0 (+ (exp (/ (PI) s)) 1.0)))
u
(/ 1.0 (+ (exp (/ (pow (cbrt (PI)) 3.0) s)) 1.0))))
1.0))))\begin{array}{l}
\\
\left(-s\right) \cdot \log \left(\frac{1}{\mathsf{fma}\left(\frac{1}{e^{\frac{-\mathsf{PI}\left(\right)}{s}} + 1} - \frac{1}{e^{\frac{\mathsf{PI}\left(\right)}{s}} + 1}, u, \frac{1}{e^{\frac{{\left(\sqrt[3]{\mathsf{PI}\left(\right)}\right)}^{3}}{s}} + 1}\right)} - 1\right)
\end{array}
Initial program 99.1%
Taylor expanded in u around 0
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites99.1%
Applied rewrites99.1%
(FPCore (u s)
:precision binary32
(let* ((t_0 (/ 1.0 (+ (exp (/ (PI) s)) 1.0))))
(*
(- s)
(log
(-
(/ 1.0 (fma (- (/ 1.0 (+ (exp (/ (- (PI)) s)) 1.0)) t_0) u t_0))
1.0)))))\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{e^{\frac{\mathsf{PI}\left(\right)}{s}} + 1}\\
\left(-s\right) \cdot \log \left(\frac{1}{\mathsf{fma}\left(\frac{1}{e^{\frac{-\mathsf{PI}\left(\right)}{s}} + 1} - t\_0, u, t\_0\right)} - 1\right)
\end{array}
\end{array}
Initial program 99.1%
Taylor expanded in u around 0
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites99.1%
(FPCore (u s)
:precision binary32
(*
(- s)
(log
(-
(/
1.0
(*
(- (/ 1.0 (+ (exp (/ (- (PI)) s)) 1.0)) (/ 1.0 (+ (exp (/ (PI) s)) 1.0)))
u))
1.0))))\begin{array}{l}
\\
\left(-s\right) \cdot \log \left(\frac{1}{\left(\frac{1}{e^{\frac{-\mathsf{PI}\left(\right)}{s}} + 1} - \frac{1}{e^{\frac{\mathsf{PI}\left(\right)}{s}} + 1}\right) \cdot u} - 1\right)
\end{array}
Initial program 99.1%
Taylor expanded in u around inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites98.4%
(FPCore (u s) :precision binary32 (let* ((t_0 (/ 1.0 (+ (exp (/ (PI) s)) 1.0)))) (* (- s) (log (- (/ 1.0 (fma (- 0.5 t_0) u t_0)) 1.0)))))
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{e^{\frac{\mathsf{PI}\left(\right)}{s}} + 1}\\
\left(-s\right) \cdot \log \left(\frac{1}{\mathsf{fma}\left(0.5 - t\_0, u, t\_0\right)} - 1\right)
\end{array}
\end{array}
Initial program 99.1%
Taylor expanded in u around 0
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites99.1%
Taylor expanded in s around inf
Applied rewrites37.1%
(FPCore (u s)
:precision binary32
(*
(- s)
(log
(-
(/
1.0
(fma
(-
0.5
(/ 1.0 (+ (- 1.0 (/ (fma (/ (* (PI) (PI)) s) -0.5 (- (PI))) s)) 1.0)))
u
(/ 1.0 (+ (exp (/ (PI) s)) 1.0))))
1.0))))\begin{array}{l}
\\
\left(-s\right) \cdot \log \left(\frac{1}{\mathsf{fma}\left(0.5 - \frac{1}{\left(1 - \frac{\mathsf{fma}\left(\frac{\mathsf{PI}\left(\right) \cdot \mathsf{PI}\left(\right)}{s}, -0.5, -\mathsf{PI}\left(\right)\right)}{s}\right) + 1}, u, \frac{1}{e^{\frac{\mathsf{PI}\left(\right)}{s}} + 1}\right)} - 1\right)
\end{array}
Initial program 99.1%
Taylor expanded in u around 0
*-commutativeN/A
lower-fma.f32N/A
Applied rewrites99.1%
Taylor expanded in s around inf
Applied rewrites37.1%
Taylor expanded in s around -inf
Applied rewrites37.1%
(FPCore (u s) :precision binary32 (* (- s) (log (fma (* -0.25 (PI)) (/ -4.0 s) 1.0))))
\begin{array}{l}
\\
\left(-s\right) \cdot \log \left(\mathsf{fma}\left(-0.25 \cdot \mathsf{PI}\left(\right), \frac{-4}{s}, 1\right)\right)
\end{array}
Initial program 99.1%
Taylor expanded in s around inf
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
lower-fma.f32N/A
Applied rewrites25.1%
Taylor expanded in u around 0
Applied rewrites25.2%
(FPCore (u s) :precision binary32 (* (- s) (log (+ (/ (PI) s) 1.0))))
\begin{array}{l}
\\
\left(-s\right) \cdot \log \left(\frac{\mathsf{PI}\left(\right)}{s} + 1\right)
\end{array}
Initial program 99.1%
Taylor expanded in s around inf
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
lower-fma.f32N/A
Applied rewrites25.1%
Taylor expanded in u around 0
Applied rewrites25.2%
(FPCore (u s) :precision binary32 (* (fma (* u (PI)) 0.5 (* -0.25 (PI))) 4.0))
\begin{array}{l}
\\
\mathsf{fma}\left(u \cdot \mathsf{PI}\left(\right), 0.5, -0.25 \cdot \mathsf{PI}\left(\right)\right) \cdot 4
\end{array}
Initial program 99.1%
Taylor expanded in s around inf
*-commutativeN/A
lower-*.f32N/A
fp-cancel-sub-sign-invN/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f32N/A
lower-*.f32N/A
lower-PI.f32N/A
lower-*.f32N/A
lower-PI.f3210.9
Applied rewrites10.9%
(FPCore (u s) :precision binary32 (* (fma (* 0.5 (PI)) u (* -0.25 (PI))) 4.0))
\begin{array}{l}
\\
\mathsf{fma}\left(0.5 \cdot \mathsf{PI}\left(\right), u, -0.25 \cdot \mathsf{PI}\left(\right)\right) \cdot 4
\end{array}
Initial program 99.1%
Taylor expanded in u around 0
mul-1-negN/A
lower-neg.f32N/A
lower-PI.f3210.7
Applied rewrites10.7%
Taylor expanded in s around inf
*-commutativeN/A
lower-*.f32N/A
Applied rewrites10.9%
(FPCore (u s) :precision binary32 (- (PI)))
\begin{array}{l}
\\
-\mathsf{PI}\left(\right)
\end{array}
Initial program 99.1%
Taylor expanded in u around 0
mul-1-negN/A
lower-neg.f32N/A
lower-PI.f3210.7
Applied rewrites10.7%
herbie shell --seed 2024356
(FPCore (u s)
:name "Sample trimmed logistic on [-pi, pi]"
:precision binary32
:pre (and (and (<= 2.328306437e-10 u) (<= u 1.0)) (and (<= 0.0 s) (<= s 1.0651631)))
(* (- s) (log (- (/ 1.0 (+ (* u (- (/ 1.0 (+ 1.0 (exp (/ (- (PI)) s)))) (/ 1.0 (+ 1.0 (exp (/ (PI) s)))))) (/ 1.0 (+ 1.0 (exp (/ (PI) s)))))) 1.0))))