
(FPCore (s r) :precision binary32 (+ (/ (* 0.25 (exp (/ (- r) s))) (* (* (* 2.0 PI) s) r)) (/ (* 0.75 (exp (/ (- r) (* 3.0 s)))) (* (* (* 6.0 PI) s) r))))
float code(float s, float r) {
return ((0.25f * expf((-r / s))) / (((2.0f * ((float) M_PI)) * s) * r)) + ((0.75f * expf((-r / (3.0f * s)))) / (((6.0f * ((float) M_PI)) * s) * r));
}
function code(s, r) return Float32(Float32(Float32(Float32(0.25) * exp(Float32(Float32(-r) / s))) / Float32(Float32(Float32(Float32(2.0) * Float32(pi)) * s) * r)) + Float32(Float32(Float32(0.75) * exp(Float32(Float32(-r) / Float32(Float32(3.0) * s)))) / Float32(Float32(Float32(Float32(6.0) * Float32(pi)) * s) * r))) end
function tmp = code(s, r) tmp = ((single(0.25) * exp((-r / s))) / (((single(2.0) * single(pi)) * s) * r)) + ((single(0.75) * exp((-r / (single(3.0) * s)))) / (((single(6.0) * single(pi)) * s) * r)); end
\begin{array}{l}
\\
\frac{0.25 \cdot e^{\frac{-r}{s}}}{\left(\left(2 \cdot \pi\right) \cdot s\right) \cdot r} + \frac{0.75 \cdot e^{\frac{-r}{3 \cdot s}}}{\left(\left(6 \cdot \pi\right) \cdot s\right) \cdot r}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (s r) :precision binary32 (+ (/ (* 0.25 (exp (/ (- r) s))) (* (* (* 2.0 PI) s) r)) (/ (* 0.75 (exp (/ (- r) (* 3.0 s)))) (* (* (* 6.0 PI) s) r))))
float code(float s, float r) {
return ((0.25f * expf((-r / s))) / (((2.0f * ((float) M_PI)) * s) * r)) + ((0.75f * expf((-r / (3.0f * s)))) / (((6.0f * ((float) M_PI)) * s) * r));
}
function code(s, r) return Float32(Float32(Float32(Float32(0.25) * exp(Float32(Float32(-r) / s))) / Float32(Float32(Float32(Float32(2.0) * Float32(pi)) * s) * r)) + Float32(Float32(Float32(0.75) * exp(Float32(Float32(-r) / Float32(Float32(3.0) * s)))) / Float32(Float32(Float32(Float32(6.0) * Float32(pi)) * s) * r))) end
function tmp = code(s, r) tmp = ((single(0.25) * exp((-r / s))) / (((single(2.0) * single(pi)) * s) * r)) + ((single(0.75) * exp((-r / (single(3.0) * s)))) / (((single(6.0) * single(pi)) * s) * r)); end
\begin{array}{l}
\\
\frac{0.25 \cdot e^{\frac{-r}{s}}}{\left(\left(2 \cdot \pi\right) \cdot s\right) \cdot r} + \frac{0.75 \cdot e^{\frac{-r}{3 \cdot s}}}{\left(\left(6 \cdot \pi\right) \cdot s\right) \cdot r}
\end{array}
(FPCore (s r) :precision binary32 (+ (/ (* 0.25 (exp (- 0.0 (/ r s)))) (* r (* s (* 2.0 PI)))) (* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r))))
float code(float s, float r) {
return ((0.25f * expf((0.0f - (r / s)))) / (r * (s * (2.0f * ((float) M_PI))))) + ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r));
}
function code(s, r) return Float32(Float32(Float32(Float32(0.25) * exp(Float32(Float32(0.0) - Float32(r / s)))) / Float32(r * Float32(s * Float32(Float32(2.0) * Float32(pi))))) + Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r))) end
function tmp = code(s, r) tmp = ((single(0.25) * exp((single(0.0) - (r / s)))) / (r * (s * (single(2.0) * single(pi))))) + ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)); end
\begin{array}{l}
\\
\frac{0.25 \cdot e^{0 - \frac{r}{s}}}{r \cdot \left(s \cdot \left(2 \cdot \pi\right)\right)} + \frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (s r) :precision binary32 (+ (* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r)) (/ (/ (/ (/ 0.125 (exp (/ r s))) PI) s) r)))
float code(float s, float r) {
return ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r)) + ((((0.125f / expf((r / s))) / ((float) M_PI)) / s) / r);
}
function code(s, r) return Float32(Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r)) + Float32(Float32(Float32(Float32(Float32(0.125) / exp(Float32(r / s))) / Float32(pi)) / s) / r)) end
function tmp = code(s, r) tmp = ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)) + ((((single(0.125) / exp((r / s))) / single(pi)) / s) / r); end
\begin{array}{l}
\\
\frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r} + \frac{\frac{\frac{\frac{0.125}{e^{\frac{r}{s}}}}{\pi}}{s}}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
exp-lowering-exp.f32N/A
distribute-frac-negN/A
neg-sub0N/A
--lowering--.f32N/A
/-lowering-/.f3299.6%
Applied egg-rr99.6%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3299.6%
Applied egg-rr99.6%
Final simplification99.6%
(FPCore (s r) :precision binary32 (/ (* (/ 0.125 PI) (+ (/ (exp (- 0.0 (/ r s))) s) (/ (exp (* (/ r s) -0.3333333333333333)) s))) r))
float code(float s, float r) {
return ((0.125f / ((float) M_PI)) * ((expf((0.0f - (r / s))) / s) + (expf(((r / s) * -0.3333333333333333f)) / s))) / r;
}
function code(s, r) return Float32(Float32(Float32(Float32(0.125) / Float32(pi)) * Float32(Float32(exp(Float32(Float32(0.0) - Float32(r / s))) / s) + Float32(exp(Float32(Float32(r / s) * Float32(-0.3333333333333333))) / s))) / r) end
function tmp = code(s, r) tmp = ((single(0.125) / single(pi)) * ((exp((single(0.0) - (r / s))) / s) + (exp(((r / s) * single(-0.3333333333333333))) / s))) / r; end
\begin{array}{l}
\\
\frac{\frac{0.125}{\pi} \cdot \left(\frac{e^{0 - \frac{r}{s}}}{s} + \frac{e^{\frac{r}{s} \cdot -0.3333333333333333}}{s}\right)}{r}
\end{array}
Initial program 99.6%
Taylor expanded in r around inf
Simplified99.6%
(FPCore (s r) :precision binary32 (/ (* 0.125 (+ (exp (- 0.0 (/ r s))) (exp (/ (/ r -3.0) s)))) (* PI (* r s))))
float code(float s, float r) {
return (0.125f * (expf((0.0f - (r / s))) + expf(((r / -3.0f) / s)))) / (((float) M_PI) * (r * s));
}
function code(s, r) return Float32(Float32(Float32(0.125) * Float32(exp(Float32(Float32(0.0) - Float32(r / s))) + exp(Float32(Float32(r / Float32(-3.0)) / s)))) / Float32(Float32(pi) * Float32(r * s))) end
function tmp = code(s, r) tmp = (single(0.125) * (exp((single(0.0) - (r / s))) + exp(((r / single(-3.0)) / s)))) / (single(pi) * (r * s)); end
\begin{array}{l}
\\
\frac{0.125 \cdot \left(e^{0 - \frac{r}{s}} + e^{\frac{\frac{r}{-3}}{s}}\right)}{\pi \cdot \left(r \cdot s\right)}
\end{array}
Initial program 99.6%
Simplified96.5%
associate-/l*N/A
frac-timesN/A
/-lowering-/.f32N/A
Applied egg-rr99.5%
*-commutativeN/A
associate-/r*N/A
remove-double-negN/A
metadata-evalN/A
frac-2negN/A
/-lowering-/.f32N/A
frac-2negN/A
remove-double-negN/A
metadata-evalN/A
/-lowering-/.f3299.5%
Applied egg-rr99.5%
Final simplification99.5%
(FPCore (s r) :precision binary32 (/ (* 0.125 (+ (exp (- 0.0 (/ r s))) (exp (/ r (* s -3.0))))) (* PI (* r s))))
float code(float s, float r) {
return (0.125f * (expf((0.0f - (r / s))) + expf((r / (s * -3.0f))))) / (((float) M_PI) * (r * s));
}
function code(s, r) return Float32(Float32(Float32(0.125) * Float32(exp(Float32(Float32(0.0) - Float32(r / s))) + exp(Float32(r / Float32(s * Float32(-3.0)))))) / Float32(Float32(pi) * Float32(r * s))) end
function tmp = code(s, r) tmp = (single(0.125) * (exp((single(0.0) - (r / s))) + exp((r / (s * single(-3.0)))))) / (single(pi) * (r * s)); end
\begin{array}{l}
\\
\frac{0.125 \cdot \left(e^{0 - \frac{r}{s}} + e^{\frac{r}{s \cdot -3}}\right)}{\pi \cdot \left(r \cdot s\right)}
\end{array}
Initial program 99.6%
Simplified96.5%
associate-/l*N/A
frac-timesN/A
/-lowering-/.f32N/A
Applied egg-rr99.5%
Final simplification99.5%
(FPCore (s r)
:precision binary32
(+
(* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r))
(/
(/
(/
(/
0.125
(+
(*
r
(+
(/ 1.0 s)
(*
r
(+ (* 0.16666666666666666 (/ r (* s (* s s)))) (/ 0.5 (* s s))))))
1.0))
PI)
s)
r)))
float code(float s, float r) {
return ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r)) + ((((0.125f / ((r * ((1.0f / s) + (r * ((0.16666666666666666f * (r / (s * (s * s)))) + (0.5f / (s * s)))))) + 1.0f)) / ((float) M_PI)) / s) / r);
}
function code(s, r) return Float32(Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r)) + Float32(Float32(Float32(Float32(Float32(0.125) / Float32(Float32(r * Float32(Float32(Float32(1.0) / s) + Float32(r * Float32(Float32(Float32(0.16666666666666666) * Float32(r / Float32(s * Float32(s * s)))) + Float32(Float32(0.5) / Float32(s * s)))))) + Float32(1.0))) / Float32(pi)) / s) / r)) end
function tmp = code(s, r) tmp = ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)) + ((((single(0.125) / ((r * ((single(1.0) / s) + (r * ((single(0.16666666666666666) * (r / (s * (s * s)))) + (single(0.5) / (s * s)))))) + single(1.0))) / single(pi)) / s) / r); end
\begin{array}{l}
\\
\frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r} + \frac{\frac{\frac{\frac{0.125}{r \cdot \left(\frac{1}{s} + r \cdot \left(0.16666666666666666 \cdot \frac{r}{s \cdot \left(s \cdot s\right)} + \frac{0.5}{s \cdot s}\right)\right) + 1}}{\pi}}{s}}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
exp-lowering-exp.f32N/A
distribute-frac-negN/A
neg-sub0N/A
--lowering--.f32N/A
/-lowering-/.f3299.6%
Applied egg-rr99.6%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3299.6%
Applied egg-rr99.6%
Taylor expanded in r around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
cube-multN/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3272.3%
Simplified72.3%
Final simplification72.3%
(FPCore (s r)
:precision binary32
(+
(* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r))
(/
(/
(/
(/
0.125
(+
(/
(+
r
(/ (+ (/ (* 0.16666666666666666 (* r (* r r))) s) (* r (* r 0.5))) s))
s)
1.0))
PI)
s)
r)))
float code(float s, float r) {
return ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r)) + ((((0.125f / (((r + ((((0.16666666666666666f * (r * (r * r))) / s) + (r * (r * 0.5f))) / s)) / s) + 1.0f)) / ((float) M_PI)) / s) / r);
}
function code(s, r) return Float32(Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r)) + Float32(Float32(Float32(Float32(Float32(0.125) / Float32(Float32(Float32(r + Float32(Float32(Float32(Float32(Float32(0.16666666666666666) * Float32(r * Float32(r * r))) / s) + Float32(r * Float32(r * Float32(0.5)))) / s)) / s) + Float32(1.0))) / Float32(pi)) / s) / r)) end
function tmp = code(s, r) tmp = ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)) + ((((single(0.125) / (((r + ((((single(0.16666666666666666) * (r * (r * r))) / s) + (r * (r * single(0.5)))) / s)) / s) + single(1.0))) / single(pi)) / s) / r); end
\begin{array}{l}
\\
\frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r} + \frac{\frac{\frac{\frac{0.125}{\frac{r + \frac{\frac{0.16666666666666666 \cdot \left(r \cdot \left(r \cdot r\right)\right)}{s} + r \cdot \left(r \cdot 0.5\right)}{s}}{s} + 1}}{\pi}}{s}}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
exp-lowering-exp.f32N/A
distribute-frac-negN/A
neg-sub0N/A
--lowering--.f32N/A
/-lowering-/.f3299.6%
Applied egg-rr99.6%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3299.6%
Applied egg-rr99.6%
Taylor expanded in s around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
/-lowering-/.f32N/A
Simplified69.8%
Final simplification69.8%
(FPCore (s r) :precision binary32 (+ (* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r)) (/ (/ (/ (/ 0.125 (+ (* r (+ (/ 1.0 s) (/ (/ (* r 0.5) s) s))) 1.0)) PI) s) r)))
float code(float s, float r) {
return ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r)) + ((((0.125f / ((r * ((1.0f / s) + (((r * 0.5f) / s) / s))) + 1.0f)) / ((float) M_PI)) / s) / r);
}
function code(s, r) return Float32(Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r)) + Float32(Float32(Float32(Float32(Float32(0.125) / Float32(Float32(r * Float32(Float32(Float32(1.0) / s) + Float32(Float32(Float32(r * Float32(0.5)) / s) / s))) + Float32(1.0))) / Float32(pi)) / s) / r)) end
function tmp = code(s, r) tmp = ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)) + ((((single(0.125) / ((r * ((single(1.0) / s) + (((r * single(0.5)) / s) / s))) + single(1.0))) / single(pi)) / s) / r); end
\begin{array}{l}
\\
\frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r} + \frac{\frac{\frac{\frac{0.125}{r \cdot \left(\frac{1}{s} + \frac{\frac{r \cdot 0.5}{s}}{s}\right) + 1}}{\pi}}{s}}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
exp-lowering-exp.f32N/A
distribute-frac-negN/A
neg-sub0N/A
--lowering--.f32N/A
/-lowering-/.f3299.6%
Applied egg-rr99.6%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3299.6%
Applied egg-rr99.6%
Taylor expanded in r around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
associate-*r/N/A
associate-*l/N/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
unpow2N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3255.2%
Simplified55.2%
Final simplification55.2%
(FPCore (s r) :precision binary32 (+ (* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r)) (/ (/ (/ (/ 0.125 (+ (/ (- r (/ (* (* r r) -0.5) s)) s) 1.0)) PI) s) r)))
float code(float s, float r) {
return ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r)) + ((((0.125f / (((r - (((r * r) * -0.5f) / s)) / s) + 1.0f)) / ((float) M_PI)) / s) / r);
}
function code(s, r) return Float32(Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r)) + Float32(Float32(Float32(Float32(Float32(0.125) / Float32(Float32(Float32(r - Float32(Float32(Float32(r * r) * Float32(-0.5)) / s)) / s) + Float32(1.0))) / Float32(pi)) / s) / r)) end
function tmp = code(s, r) tmp = ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)) + ((((single(0.125) / (((r - (((r * r) * single(-0.5)) / s)) / s) + single(1.0))) / single(pi)) / s) / r); end
\begin{array}{l}
\\
\frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r} + \frac{\frac{\frac{\frac{0.125}{\frac{r - \frac{\left(r \cdot r\right) \cdot -0.5}{s}}{s} + 1}}{\pi}}{s}}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
exp-lowering-exp.f32N/A
distribute-frac-negN/A
neg-sub0N/A
--lowering--.f32N/A
/-lowering-/.f3299.6%
Applied egg-rr99.6%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3299.6%
Applied egg-rr99.6%
Taylor expanded in s around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
/-lowering-/.f32N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3253.7%
Simplified53.7%
Final simplification53.7%
(FPCore (s r)
:precision binary32
(let* ((t_0 (* r (* (* s PI) -6.0))))
(if (<= s 1.4459999874816276e-5)
(/
(- (/ (/ 0.015625 (* (* s PI) (* s PI))) (* r r)) (/ 0.5625 (* t_0 t_0)))
(- (/ (/ 0.125 (* s PI)) r) (/ 0.75 (* 6.0 (* r (* s PI))))))
(/
(*
0.125
(+
(exp (- 0.0 (/ r s)))
(-
1.0
(/
(+ (* (/ (* r r) s) -0.05555555555555555) (* r 0.3333333333333333))
s))))
(* PI (* r s))))))
float code(float s, float r) {
float t_0 = r * ((s * ((float) M_PI)) * -6.0f);
float tmp;
if (s <= 1.4459999874816276e-5f) {
tmp = (((0.015625f / ((s * ((float) M_PI)) * (s * ((float) M_PI)))) / (r * r)) - (0.5625f / (t_0 * t_0))) / (((0.125f / (s * ((float) M_PI))) / r) - (0.75f / (6.0f * (r * (s * ((float) M_PI))))));
} else {
tmp = (0.125f * (expf((0.0f - (r / s))) + (1.0f - (((((r * r) / s) * -0.05555555555555555f) + (r * 0.3333333333333333f)) / s)))) / (((float) M_PI) * (r * s));
}
return tmp;
}
function code(s, r) t_0 = Float32(r * Float32(Float32(s * Float32(pi)) * Float32(-6.0))) tmp = Float32(0.0) if (s <= Float32(1.4459999874816276e-5)) tmp = Float32(Float32(Float32(Float32(Float32(0.015625) / Float32(Float32(s * Float32(pi)) * Float32(s * Float32(pi)))) / Float32(r * r)) - Float32(Float32(0.5625) / Float32(t_0 * t_0))) / Float32(Float32(Float32(Float32(0.125) / Float32(s * Float32(pi))) / r) - Float32(Float32(0.75) / Float32(Float32(6.0) * Float32(r * Float32(s * Float32(pi))))))); else tmp = Float32(Float32(Float32(0.125) * Float32(exp(Float32(Float32(0.0) - Float32(r / s))) + Float32(Float32(1.0) - Float32(Float32(Float32(Float32(Float32(r * r) / s) * Float32(-0.05555555555555555)) + Float32(r * Float32(0.3333333333333333))) / s)))) / Float32(Float32(pi) * Float32(r * s))); end return tmp end
function tmp_2 = code(s, r) t_0 = r * ((s * single(pi)) * single(-6.0)); tmp = single(0.0); if (s <= single(1.4459999874816276e-5)) tmp = (((single(0.015625) / ((s * single(pi)) * (s * single(pi)))) / (r * r)) - (single(0.5625) / (t_0 * t_0))) / (((single(0.125) / (s * single(pi))) / r) - (single(0.75) / (single(6.0) * (r * (s * single(pi)))))); else tmp = (single(0.125) * (exp((single(0.0) - (r / s))) + (single(1.0) - (((((r * r) / s) * single(-0.05555555555555555)) + (r * single(0.3333333333333333))) / s)))) / (single(pi) * (r * s)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := r \cdot \left(\left(s \cdot \pi\right) \cdot -6\right)\\
\mathbf{if}\;s \leq 1.4459999874816276 \cdot 10^{-5}:\\
\;\;\;\;\frac{\frac{\frac{0.015625}{\left(s \cdot \pi\right) \cdot \left(s \cdot \pi\right)}}{r \cdot r} - \frac{0.5625}{t\_0 \cdot t\_0}}{\frac{\frac{0.125}{s \cdot \pi}}{r} - \frac{0.75}{6 \cdot \left(r \cdot \left(s \cdot \pi\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.125 \cdot \left(e^{0 - \frac{r}{s}} + \left(1 - \frac{\frac{r \cdot r}{s} \cdot -0.05555555555555555 + r \cdot 0.3333333333333333}{s}\right)\right)}{\pi \cdot \left(r \cdot s\right)}\\
\end{array}
\end{array}
if s < 1.44599999e-5Initial program 99.9%
Taylor expanded in r around 0
Simplified4.6%
Taylor expanded in r around 0
Simplified4.6%
frac-2negN/A
div-invN/A
*-lowering-*.f32N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-neg-inN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
metadata-eval4.6%
Applied egg-rr4.6%
Applied egg-rr8.4%
if 1.44599999e-5 < s Initial program 98.2%
Simplified98.0%
associate-/l*N/A
frac-timesN/A
/-lowering-/.f32N/A
Applied egg-rr97.3%
Taylor expanded in s around -inf
mul-1-negN/A
unsub-negN/A
--lowering--.f32N/A
/-lowering-/.f32N/A
+-lowering-+.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f3245.1%
Simplified45.1%
Final simplification14.7%
(FPCore (s r) :precision binary32 (+ (* (/ (exp (/ r (* s -3.0))) (* s (* PI 6.0))) (/ 0.75 r)) (/ (/ (/ (/ 0.125 (+ (/ r s) 1.0)) PI) s) r)))
float code(float s, float r) {
return ((expf((r / (s * -3.0f))) / (s * (((float) M_PI) * 6.0f))) * (0.75f / r)) + ((((0.125f / ((r / s) + 1.0f)) / ((float) M_PI)) / s) / r);
}
function code(s, r) return Float32(Float32(Float32(exp(Float32(r / Float32(s * Float32(-3.0)))) / Float32(s * Float32(Float32(pi) * Float32(6.0)))) * Float32(Float32(0.75) / r)) + Float32(Float32(Float32(Float32(Float32(0.125) / Float32(Float32(r / s) + Float32(1.0))) / Float32(pi)) / s) / r)) end
function tmp = code(s, r) tmp = ((exp((r / (s * single(-3.0)))) / (s * (single(pi) * single(6.0)))) * (single(0.75) / r)) + ((((single(0.125) / ((r / s) + single(1.0))) / single(pi)) / s) / r); end
\begin{array}{l}
\\
\frac{e^{\frac{r}{s \cdot -3}}}{s \cdot \left(\pi \cdot 6\right)} \cdot \frac{0.75}{r} + \frac{\frac{\frac{\frac{0.125}{\frac{r}{s} + 1}}{\pi}}{s}}{r}
\end{array}
Initial program 99.6%
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr99.7%
times-fracN/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
exp-lowering-exp.f32N/A
distribute-frac-negN/A
neg-sub0N/A
--lowering--.f32N/A
/-lowering-/.f3299.6%
Applied egg-rr99.6%
associate-*l/N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-diffN/A
1-expN/A
un-div-invN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
PI-lowering-PI.f3299.6%
Applied egg-rr99.6%
Taylor expanded in r around 0
+-lowering-+.f32N/A
/-lowering-/.f3215.5%
Simplified15.5%
Final simplification15.5%
(FPCore (s r)
:precision binary32
(let* ((t_0 (* r (* (* s PI) -6.0))))
(if (<= s 1.4459999874816276e-5)
(/
(- (/ (/ 0.015625 (* (* s PI) (* s PI))) (* r r)) (/ 0.5625 (* t_0 t_0)))
(- (/ (/ 0.125 (* s PI)) r) (/ 0.75 (* 6.0 (* r (* s PI))))))
(/
(/
(+
(/ 0.25 r)
(/
(+ -0.16666666666666666 (* 0.125 (/ (* r 0.5555555555555556) s)))
s))
s)
PI))))
float code(float s, float r) {
float t_0 = r * ((s * ((float) M_PI)) * -6.0f);
float tmp;
if (s <= 1.4459999874816276e-5f) {
tmp = (((0.015625f / ((s * ((float) M_PI)) * (s * ((float) M_PI)))) / (r * r)) - (0.5625f / (t_0 * t_0))) / (((0.125f / (s * ((float) M_PI))) / r) - (0.75f / (6.0f * (r * (s * ((float) M_PI))))));
} else {
tmp = (((0.25f / r) + ((-0.16666666666666666f + (0.125f * ((r * 0.5555555555555556f) / s))) / s)) / s) / ((float) M_PI);
}
return tmp;
}
function code(s, r) t_0 = Float32(r * Float32(Float32(s * Float32(pi)) * Float32(-6.0))) tmp = Float32(0.0) if (s <= Float32(1.4459999874816276e-5)) tmp = Float32(Float32(Float32(Float32(Float32(0.015625) / Float32(Float32(s * Float32(pi)) * Float32(s * Float32(pi)))) / Float32(r * r)) - Float32(Float32(0.5625) / Float32(t_0 * t_0))) / Float32(Float32(Float32(Float32(0.125) / Float32(s * Float32(pi))) / r) - Float32(Float32(0.75) / Float32(Float32(6.0) * Float32(r * Float32(s * Float32(pi))))))); else tmp = Float32(Float32(Float32(Float32(Float32(0.25) / r) + Float32(Float32(Float32(-0.16666666666666666) + Float32(Float32(0.125) * Float32(Float32(r * Float32(0.5555555555555556)) / s))) / s)) / s) / Float32(pi)); end return tmp end
function tmp_2 = code(s, r) t_0 = r * ((s * single(pi)) * single(-6.0)); tmp = single(0.0); if (s <= single(1.4459999874816276e-5)) tmp = (((single(0.015625) / ((s * single(pi)) * (s * single(pi)))) / (r * r)) - (single(0.5625) / (t_0 * t_0))) / (((single(0.125) / (s * single(pi))) / r) - (single(0.75) / (single(6.0) * (r * (s * single(pi)))))); else tmp = (((single(0.25) / r) + ((single(-0.16666666666666666) + (single(0.125) * ((r * single(0.5555555555555556)) / s))) / s)) / s) / single(pi); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := r \cdot \left(\left(s \cdot \pi\right) \cdot -6\right)\\
\mathbf{if}\;s \leq 1.4459999874816276 \cdot 10^{-5}:\\
\;\;\;\;\frac{\frac{\frac{0.015625}{\left(s \cdot \pi\right) \cdot \left(s \cdot \pi\right)}}{r \cdot r} - \frac{0.5625}{t\_0 \cdot t\_0}}{\frac{\frac{0.125}{s \cdot \pi}}{r} - \frac{0.75}{6 \cdot \left(r \cdot \left(s \cdot \pi\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{0.25}{r} + \frac{-0.16666666666666666 + 0.125 \cdot \frac{r \cdot 0.5555555555555556}{s}}{s}}{s}}{\pi}\\
\end{array}
\end{array}
if s < 1.44599999e-5Initial program 99.9%
Taylor expanded in r around 0
Simplified4.6%
Taylor expanded in r around 0
Simplified4.6%
frac-2negN/A
div-invN/A
*-lowering-*.f32N/A
metadata-evalN/A
/-lowering-/.f32N/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-neg-inN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
metadata-eval4.6%
Applied egg-rr4.6%
Applied egg-rr8.4%
if 1.44599999e-5 < s Initial program 98.2%
Simplified98.0%
Taylor expanded in s around inf
Simplified41.8%
Final simplification14.1%
(FPCore (s r)
:precision binary32
(/
(/
(+
(/ 0.25 r)
(/ (+ -0.16666666666666666 (* 0.125 (/ (* r 0.5555555555555556) s))) s))
s)
PI))
float code(float s, float r) {
return (((0.25f / r) + ((-0.16666666666666666f + (0.125f * ((r * 0.5555555555555556f) / s))) / s)) / s) / ((float) M_PI);
}
function code(s, r) return Float32(Float32(Float32(Float32(Float32(0.25) / r) + Float32(Float32(Float32(-0.16666666666666666) + Float32(Float32(0.125) * Float32(Float32(r * Float32(0.5555555555555556)) / s))) / s)) / s) / Float32(pi)) end
function tmp = code(s, r) tmp = (((single(0.25) / r) + ((single(-0.16666666666666666) + (single(0.125) * ((r * single(0.5555555555555556)) / s))) / s)) / s) / single(pi); end
\begin{array}{l}
\\
\frac{\frac{\frac{0.25}{r} + \frac{-0.16666666666666666 + 0.125 \cdot \frac{r \cdot 0.5555555555555556}{s}}{s}}{s}}{\pi}
\end{array}
Initial program 99.6%
Simplified96.5%
Taylor expanded in s around inf
Simplified10.0%
(FPCore (s r) :precision binary32 (/ (/ (+ (/ 0.25 r) (/ -0.16666666666666666 s)) s) PI))
float code(float s, float r) {
return (((0.25f / r) + (-0.16666666666666666f / s)) / s) / ((float) M_PI);
}
function code(s, r) return Float32(Float32(Float32(Float32(Float32(0.25) / r) + Float32(Float32(-0.16666666666666666) / s)) / s) / Float32(pi)) end
function tmp = code(s, r) tmp = (((single(0.25) / r) + (single(-0.16666666666666666) / s)) / s) / single(pi); end
\begin{array}{l}
\\
\frac{\frac{\frac{0.25}{r} + \frac{-0.16666666666666666}{s}}{s}}{\pi}
\end{array}
Initial program 99.6%
Simplified96.5%
Taylor expanded in s around inf
/-lowering-/.f32N/A
sub-negN/A
+-lowering-+.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f328.9%
Simplified8.9%
(FPCore (s r) :precision binary32 (/ 1.0 (/ (* s PI) (/ 0.25 r))))
float code(float s, float r) {
return 1.0f / ((s * ((float) M_PI)) / (0.25f / r));
}
function code(s, r) return Float32(Float32(1.0) / Float32(Float32(s * Float32(pi)) / Float32(Float32(0.25) / r))) end
function tmp = code(s, r) tmp = single(1.0) / ((s * single(pi)) / (single(0.25) / r)); end
\begin{array}{l}
\\
\frac{1}{\frac{s \cdot \pi}{\frac{0.25}{r}}}
\end{array}
Initial program 99.6%
Simplified96.5%
Taylor expanded in r around 0
/-lowering-/.f32N/A
*-lowering-*.f328.7%
Simplified8.7%
associate-/r*N/A
/-lowering-/.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f328.7%
Applied egg-rr8.7%
associate-/r*N/A
clear-numN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
/-lowering-/.f328.7%
Applied egg-rr8.7%
(FPCore (s r) :precision binary32 (/ 0.25 (* s (* r PI))))
float code(float s, float r) {
return 0.25f / (s * (r * ((float) M_PI)));
}
function code(s, r) return Float32(Float32(0.25) / Float32(s * Float32(r * Float32(pi)))) end
function tmp = code(s, r) tmp = single(0.25) / (s * (r * single(pi))); end
\begin{array}{l}
\\
\frac{0.25}{s \cdot \left(r \cdot \pi\right)}
\end{array}
Initial program 99.6%
Simplified96.5%
Taylor expanded in r around 0
/-lowering-/.f32N/A
*-lowering-*.f328.7%
Simplified8.7%
associate-/r*N/A
/-lowering-/.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f328.7%
Applied egg-rr8.7%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
PI-lowering-PI.f328.7%
Applied egg-rr8.7%
Final simplification8.7%
(FPCore (s r) :precision binary32 (/ 0.25 (* r (* s PI))))
float code(float s, float r) {
return 0.25f / (r * (s * ((float) M_PI)));
}
function code(s, r) return Float32(Float32(0.25) / Float32(r * Float32(s * Float32(pi)))) end
function tmp = code(s, r) tmp = single(0.25) / (r * (s * single(pi))); end
\begin{array}{l}
\\
\frac{0.25}{r \cdot \left(s \cdot \pi\right)}
\end{array}
Initial program 99.6%
Simplified96.5%
Taylor expanded in r around 0
/-lowering-/.f32N/A
*-lowering-*.f328.7%
Simplified8.7%
associate-/r*N/A
/-lowering-/.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f328.7%
Applied egg-rr8.7%
herbie shell --seed 2024145
(FPCore (s r)
:name "Disney BSSRDF, PDF of scattering profile"
:precision binary32
:pre (and (and (<= 0.0 s) (<= s 256.0)) (and (< 1e-6 r) (< r 1000000.0)))
(+ (/ (* 0.25 (exp (/ (- r) s))) (* (* (* 2.0 PI) s) r)) (/ (* 0.75 (exp (/ (- r) (* 3.0 s)))) (* (* (* 6.0 PI) s) r))))