
(FPCore (x tau) :precision binary32 (let* ((t_1 (* (* x PI) tau))) (* (/ (sin t_1) t_1) (/ (sin (* x PI)) (* x PI)))))
float code(float x, float tau) {
float t_1 = (x * ((float) M_PI)) * tau;
return (sinf(t_1) / t_1) * (sinf((x * ((float) M_PI))) / (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(Float32(x * Float32(pi)) * tau) return Float32(Float32(sin(t_1) / t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = (x * single(pi)) * tau; tmp = (sin(t_1) / t_1) * (sin((x * single(pi))) / (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x \cdot \pi\right) \cdot tau\\
\frac{\sin t\_1}{t\_1} \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \pi}
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x tau) :precision binary32 (let* ((t_1 (* (* x PI) tau))) (* (/ (sin t_1) t_1) (/ (sin (* x PI)) (* x PI)))))
float code(float x, float tau) {
float t_1 = (x * ((float) M_PI)) * tau;
return (sinf(t_1) / t_1) * (sinf((x * ((float) M_PI))) / (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(Float32(x * Float32(pi)) * tau) return Float32(Float32(sin(t_1) / t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = (x * single(pi)) * tau; tmp = (sin(t_1) / t_1) * (sin((x * single(pi))) / (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x \cdot \pi\right) \cdot tau\\
\frac{\sin t\_1}{t\_1} \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \pi}
\end{array}
\end{array}
(FPCore (x tau) :precision binary32 (let* ((t_1 (* x (* PI tau)))) (/ (/ (sin t_1) (* x PI)) (/ t_1 (sin (* x PI))))))
float code(float x, float tau) {
float t_1 = x * (((float) M_PI) * tau);
return (sinf(t_1) / (x * ((float) M_PI))) / (t_1 / sinf((x * ((float) M_PI))));
}
function code(x, tau) t_1 = Float32(x * Float32(Float32(pi) * tau)) return Float32(Float32(sin(t_1) / Float32(x * Float32(pi))) / Float32(t_1 / sin(Float32(x * Float32(pi))))) end
function tmp = code(x, tau) t_1 = x * (single(pi) * tau); tmp = (sin(t_1) / (x * single(pi))) / (t_1 / sin((x * single(pi)))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(\pi \cdot tau\right)\\
\frac{\frac{\sin t\_1}{x \cdot \pi}}{\frac{t\_1}{\sin \left(x \cdot \pi\right)}}
\end{array}
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
associate-/l/N/A
associate-*r/N/A
associate-/r/N/A
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
Applied egg-rr97.4%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* tau (* x PI)))) (* (/ (sin t_1) t_1) (/ (sin (* x PI)) (* x PI)))))
float code(float x, float tau) {
float t_1 = tau * (x * ((float) M_PI));
return (sinf(t_1) / t_1) * (sinf((x * ((float) M_PI))) / (x * ((float) M_PI)));
}
function code(x, tau) t_1 = Float32(tau * Float32(x * Float32(pi))) return Float32(Float32(sin(t_1) / t_1) * Float32(sin(Float32(x * Float32(pi))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) t_1 = tau * (x * single(pi)); tmp = (sin(t_1) / t_1) * (sin((x * single(pi))) / (x * single(pi))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := tau \cdot \left(x \cdot \pi\right)\\
\frac{\sin t\_1}{t\_1} \cdot \frac{\sin \left(x \cdot \pi\right)}{x \cdot \pi}
\end{array}
\end{array}
Initial program 97.3%
Final simplification97.3%
(FPCore (x tau) :precision binary32 (let* ((t_1 (* x (* PI tau)))) (* (/ (sin t_1) x) (/ (/ (sin (* x PI)) t_1) PI))))
float code(float x, float tau) {
float t_1 = x * (((float) M_PI) * tau);
return (sinf(t_1) / x) * ((sinf((x * ((float) M_PI))) / t_1) / ((float) M_PI));
}
function code(x, tau) t_1 = Float32(x * Float32(Float32(pi) * tau)) return Float32(Float32(sin(t_1) / x) * Float32(Float32(sin(Float32(x * Float32(pi))) / t_1) / Float32(pi))) end
function tmp = code(x, tau) t_1 = x * (single(pi) * tau); tmp = (sin(t_1) / x) * ((sin((x * single(pi))) / t_1) / single(pi)); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \left(\pi \cdot tau\right)\\
\frac{\sin t\_1}{x} \cdot \frac{\frac{\sin \left(x \cdot \pi\right)}{t\_1}}{\pi}
\end{array}
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
associate-/l/N/A
associate-*r/N/A
associate-/r/N/A
div-invN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr97.1%
(FPCore (x tau)
:precision binary32
(let* ((t_1 (* tau (* x PI))))
(*
(/ (sin t_1) t_1)
(+
1.0
(*
x
(*
x
(+
(* -0.16666666666666666 (* PI PI))
(* (* (* PI PI) (* PI PI)) (* (* x x) 0.008333333333333333)))))))))
float code(float x, float tau) {
float t_1 = tau * (x * ((float) M_PI));
return (sinf(t_1) / t_1) * (1.0f + (x * (x * ((-0.16666666666666666f * (((float) M_PI) * ((float) M_PI))) + (((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))) * ((x * x) * 0.008333333333333333f))))));
}
function code(x, tau) t_1 = Float32(tau * Float32(x * Float32(pi))) return Float32(Float32(sin(t_1) / t_1) * Float32(Float32(1.0) + Float32(x * Float32(x * Float32(Float32(Float32(-0.16666666666666666) * Float32(Float32(pi) * Float32(pi))) + Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))) * Float32(Float32(x * x) * Float32(0.008333333333333333)))))))) end
function tmp = code(x, tau) t_1 = tau * (x * single(pi)); tmp = (sin(t_1) / t_1) * (single(1.0) + (x * (x * ((single(-0.16666666666666666) * (single(pi) * single(pi))) + (((single(pi) * single(pi)) * (single(pi) * single(pi))) * ((x * x) * single(0.008333333333333333))))))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := tau \cdot \left(x \cdot \pi\right)\\
\frac{\sin t\_1}{t\_1} \cdot \left(1 + x \cdot \left(x \cdot \left(-0.16666666666666666 \cdot \left(\pi \cdot \pi\right) + \left(\left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(\left(x \cdot x\right) \cdot 0.008333333333333333\right)\right)\right)\right)
\end{array}
\end{array}
Initial program 97.3%
associate-/l/N/A
div-invN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
/-lowering-/.f3296.9%
Applied egg-rr96.9%
Taylor expanded in x around 0
+-lowering-+.f32N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
+-lowering-+.f32N/A
Simplified89.6%
Final simplification89.6%
(FPCore (x tau)
:precision binary32
(let* ((t_1 (* tau (* x PI))))
(*
(/ (sin t_1) t_1)
(+ 1.0 (* x (* x (* PI (* PI -0.16666666666666666))))))))
float code(float x, float tau) {
float t_1 = tau * (x * ((float) M_PI));
return (sinf(t_1) / t_1) * (1.0f + (x * (x * (((float) M_PI) * (((float) M_PI) * -0.16666666666666666f)))));
}
function code(x, tau) t_1 = Float32(tau * Float32(x * Float32(pi))) return Float32(Float32(sin(t_1) / t_1) * Float32(Float32(1.0) + Float32(x * Float32(x * Float32(Float32(pi) * Float32(Float32(pi) * Float32(-0.16666666666666666))))))) end
function tmp = code(x, tau) t_1 = tau * (x * single(pi)); tmp = (sin(t_1) / t_1) * (single(1.0) + (x * (x * (single(pi) * (single(pi) * single(-0.16666666666666666)))))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := tau \cdot \left(x \cdot \pi\right)\\
\frac{\sin t\_1}{t\_1} \cdot \left(1 + x \cdot \left(x \cdot \left(\pi \cdot \left(\pi \cdot -0.16666666666666666\right)\right)\right)\right)
\end{array}
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.7%
Simplified82.7%
Final simplification82.7%
(FPCore (x tau) :precision binary32 (* (/ (sin (* x (* PI tau))) x) (/ (+ 1.0 (* x (* x (* -0.16666666666666666 (* PI PI))))) (* PI tau))))
float code(float x, float tau) {
return (sinf((x * (((float) M_PI) * tau))) / x) * ((1.0f + (x * (x * (-0.16666666666666666f * (((float) M_PI) * ((float) M_PI)))))) / (((float) M_PI) * tau));
}
function code(x, tau) return Float32(Float32(sin(Float32(x * Float32(Float32(pi) * tau))) / x) * Float32(Float32(Float32(1.0) + Float32(x * Float32(x * Float32(Float32(-0.16666666666666666) * Float32(Float32(pi) * Float32(pi)))))) / Float32(Float32(pi) * tau))) end
function tmp = code(x, tau) tmp = (sin((x * (single(pi) * tau))) / x) * ((single(1.0) + (x * (x * (single(-0.16666666666666666) * (single(pi) * single(pi)))))) / (single(pi) * tau)); end
\begin{array}{l}
\\
\frac{\sin \left(x \cdot \left(\pi \cdot tau\right)\right)}{x} \cdot \frac{1 + x \cdot \left(x \cdot \left(-0.16666666666666666 \cdot \left(\pi \cdot \pi\right)\right)\right)}{\pi \cdot tau}
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.7%
Simplified82.7%
associate-*r*N/A
associate-*l/N/A
associate-*r*N/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr82.5%
Final simplification82.5%
(FPCore (x tau) :precision binary32 (* (/ (sin (* x (* PI tau))) tau) (/ (+ 1.0 (* x (* x (* -0.16666666666666666 (* PI PI))))) (* x PI))))
float code(float x, float tau) {
return (sinf((x * (((float) M_PI) * tau))) / tau) * ((1.0f + (x * (x * (-0.16666666666666666f * (((float) M_PI) * ((float) M_PI)))))) / (x * ((float) M_PI)));
}
function code(x, tau) return Float32(Float32(sin(Float32(x * Float32(Float32(pi) * tau))) / tau) * Float32(Float32(Float32(1.0) + Float32(x * Float32(x * Float32(Float32(-0.16666666666666666) * Float32(Float32(pi) * Float32(pi)))))) / Float32(x * Float32(pi)))) end
function tmp = code(x, tau) tmp = (sin((x * (single(pi) * tau))) / tau) * ((single(1.0) + (x * (x * (single(-0.16666666666666666) * (single(pi) * single(pi)))))) / (x * single(pi))); end
\begin{array}{l}
\\
\frac{\sin \left(x \cdot \left(\pi \cdot tau\right)\right)}{tau} \cdot \frac{1 + x \cdot \left(x \cdot \left(-0.16666666666666666 \cdot \left(\pi \cdot \pi\right)\right)\right)}{x \cdot \pi}
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.7%
Simplified82.7%
associate-*r*N/A
associate-*l/N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
times-fracN/A
*-lowering-*.f32N/A
Applied egg-rr82.3%
Final simplification82.3%
(FPCore (x tau) :precision binary32 (* x (* (/ (sin (* x (* PI tau))) tau) (- (/ -1.0 (* x (- 0.0 (* x PI)))) (* PI 0.16666666666666666)))))
float code(float x, float tau) {
return x * ((sinf((x * (((float) M_PI) * tau))) / tau) * ((-1.0f / (x * (0.0f - (x * ((float) M_PI))))) - (((float) M_PI) * 0.16666666666666666f)));
}
function code(x, tau) return Float32(x * Float32(Float32(sin(Float32(x * Float32(Float32(pi) * tau))) / tau) * Float32(Float32(Float32(-1.0) / Float32(x * Float32(Float32(0.0) - Float32(x * Float32(pi))))) - Float32(Float32(pi) * Float32(0.16666666666666666))))) end
function tmp = code(x, tau) tmp = x * ((sin((x * (single(pi) * tau))) / tau) * ((single(-1.0) / (x * (single(0.0) - (x * single(pi))))) - (single(pi) * single(0.16666666666666666)))); end
\begin{array}{l}
\\
x \cdot \left(\frac{\sin \left(x \cdot \left(\pi \cdot tau\right)\right)}{tau} \cdot \left(\frac{-1}{x \cdot \left(0 - x \cdot \pi\right)} - \pi \cdot 0.16666666666666666\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.6%
Simplified82.6%
Taylor expanded in x around inf
Simplified82.0%
Final simplification82.0%
(FPCore (x tau)
:precision binary32
(+
1.0
(*
(* x x)
(+
(* (* PI PI) (+ -0.16666666666666666 (* -0.16666666666666666 (* tau tau))))
(*
(* x x)
(*
(* (* PI PI) (* PI PI))
(+
(* 0.008333333333333333 (* (* tau tau) (* tau tau)))
(+ 0.008333333333333333 (* (* tau tau) 0.027777777777777776)))))))))
float code(float x, float tau) {
return 1.0f + ((x * x) * (((((float) M_PI) * ((float) M_PI)) * (-0.16666666666666666f + (-0.16666666666666666f * (tau * tau)))) + ((x * x) * (((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))) * ((0.008333333333333333f * ((tau * tau) * (tau * tau))) + (0.008333333333333333f + ((tau * tau) * 0.027777777777777776f)))))));
}
function code(x, tau) return Float32(Float32(1.0) + Float32(Float32(x * x) * Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-0.16666666666666666) + Float32(Float32(-0.16666666666666666) * Float32(tau * tau)))) + Float32(Float32(x * x) * Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))) * Float32(Float32(Float32(0.008333333333333333) * Float32(Float32(tau * tau) * Float32(tau * tau))) + Float32(Float32(0.008333333333333333) + Float32(Float32(tau * tau) * Float32(0.027777777777777776))))))))) end
function tmp = code(x, tau) tmp = single(1.0) + ((x * x) * (((single(pi) * single(pi)) * (single(-0.16666666666666666) + (single(-0.16666666666666666) * (tau * tau)))) + ((x * x) * (((single(pi) * single(pi)) * (single(pi) * single(pi))) * ((single(0.008333333333333333) * ((tau * tau) * (tau * tau))) + (single(0.008333333333333333) + ((tau * tau) * single(0.027777777777777776)))))))); end
\begin{array}{l}
\\
1 + \left(x \cdot x\right) \cdot \left(\left(\pi \cdot \pi\right) \cdot \left(-0.16666666666666666 + -0.16666666666666666 \cdot \left(tau \cdot tau\right)\right) + \left(x \cdot x\right) \cdot \left(\left(\left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(0.008333333333333333 \cdot \left(\left(tau \cdot tau\right) \cdot \left(tau \cdot tau\right)\right) + \left(0.008333333333333333 + \left(tau \cdot tau\right) \cdot 0.027777777777777776\right)\right)\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
associate-/l/N/A
associate-*r/N/A
associate-/r/N/A
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
Applied egg-rr97.4%
associate-/l/N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
associate-*l/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
Applied egg-rr97.2%
Taylor expanded in x around 0
Simplified81.9%
(FPCore (x tau)
:precision binary32
(+
1.0
(*
(* x x)
(+
(* (* PI PI) (+ -0.16666666666666666 (* -0.16666666666666666 (* tau tau))))
(*
(* x x)
(*
(* (* PI PI) (* PI PI))
(+
(* 0.008333333333333333 (* (* tau tau) (* tau tau)))
(* (* tau tau) 0.027777777777777776))))))))
float code(float x, float tau) {
return 1.0f + ((x * x) * (((((float) M_PI) * ((float) M_PI)) * (-0.16666666666666666f + (-0.16666666666666666f * (tau * tau)))) + ((x * x) * (((((float) M_PI) * ((float) M_PI)) * (((float) M_PI) * ((float) M_PI))) * ((0.008333333333333333f * ((tau * tau) * (tau * tau))) + ((tau * tau) * 0.027777777777777776f))))));
}
function code(x, tau) return Float32(Float32(1.0) + Float32(Float32(x * x) * Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-0.16666666666666666) + Float32(Float32(-0.16666666666666666) * Float32(tau * tau)))) + Float32(Float32(x * x) * Float32(Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(pi) * Float32(pi))) * Float32(Float32(Float32(0.008333333333333333) * Float32(Float32(tau * tau) * Float32(tau * tau))) + Float32(Float32(tau * tau) * Float32(0.027777777777777776)))))))) end
function tmp = code(x, tau) tmp = single(1.0) + ((x * x) * (((single(pi) * single(pi)) * (single(-0.16666666666666666) + (single(-0.16666666666666666) * (tau * tau)))) + ((x * x) * (((single(pi) * single(pi)) * (single(pi) * single(pi))) * ((single(0.008333333333333333) * ((tau * tau) * (tau * tau))) + ((tau * tau) * single(0.027777777777777776))))))); end
\begin{array}{l}
\\
1 + \left(x \cdot x\right) \cdot \left(\left(\pi \cdot \pi\right) \cdot \left(-0.16666666666666666 + -0.16666666666666666 \cdot \left(tau \cdot tau\right)\right) + \left(x \cdot x\right) \cdot \left(\left(\left(\pi \cdot \pi\right) \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(0.008333333333333333 \cdot \left(\left(tau \cdot tau\right) \cdot \left(tau \cdot tau\right)\right) + \left(tau \cdot tau\right) \cdot 0.027777777777777776\right)\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.6%
Simplified82.6%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-+r+N/A
+-lowering-+.f32N/A
Simplified81.2%
(FPCore (x tau)
:precision binary32
(let* ((t_1 (* (* PI PI) (* x x))))
(+
1.0
(*
-0.16666666666666666
(+
t_1
(*
(* tau tau)
(* t_1 (+ 1.0 (* (* PI PI) (* -0.16666666666666666 (* x x)))))))))))
float code(float x, float tau) {
float t_1 = (((float) M_PI) * ((float) M_PI)) * (x * x);
return 1.0f + (-0.16666666666666666f * (t_1 + ((tau * tau) * (t_1 * (1.0f + ((((float) M_PI) * ((float) M_PI)) * (-0.16666666666666666f * (x * x))))))));
}
function code(x, tau) t_1 = Float32(Float32(Float32(pi) * Float32(pi)) * Float32(x * x)) return Float32(Float32(1.0) + Float32(Float32(-0.16666666666666666) * Float32(t_1 + Float32(Float32(tau * tau) * Float32(t_1 * Float32(Float32(1.0) + Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-0.16666666666666666) * Float32(x * x))))))))) end
function tmp = code(x, tau) t_1 = (single(pi) * single(pi)) * (x * x); tmp = single(1.0) + (single(-0.16666666666666666) * (t_1 + ((tau * tau) * (t_1 * (single(1.0) + ((single(pi) * single(pi)) * (single(-0.16666666666666666) * (x * x)))))))); end
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\pi \cdot \pi\right) \cdot \left(x \cdot x\right)\\
1 + -0.16666666666666666 \cdot \left(t\_1 + \left(tau \cdot tau\right) \cdot \left(t\_1 \cdot \left(1 + \left(\pi \cdot \pi\right) \cdot \left(-0.16666666666666666 \cdot \left(x \cdot x\right)\right)\right)\right)\right)
\end{array}
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.7%
Simplified82.7%
Taylor expanded in tau around 0
+-lowering-+.f32N/A
distribute-lft-outN/A
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
Simplified77.0%
Final simplification77.0%
(FPCore (x tau) :precision binary32 (* (+ 1.0 (* (* tau (* tau (* PI PI))) (* -0.16666666666666666 (* x x)))) (+ 1.0 (* (* -0.16666666666666666 (* PI PI)) (* x x)))))
float code(float x, float tau) {
return (1.0f + ((tau * (tau * (((float) M_PI) * ((float) M_PI)))) * (-0.16666666666666666f * (x * x)))) * (1.0f + ((-0.16666666666666666f * (((float) M_PI) * ((float) M_PI))) * (x * x)));
}
function code(x, tau) return Float32(Float32(Float32(1.0) + Float32(Float32(tau * Float32(tau * Float32(Float32(pi) * Float32(pi)))) * Float32(Float32(-0.16666666666666666) * Float32(x * x)))) * Float32(Float32(1.0) + Float32(Float32(Float32(-0.16666666666666666) * Float32(Float32(pi) * Float32(pi))) * Float32(x * x)))) end
function tmp = code(x, tau) tmp = (single(1.0) + ((tau * (tau * (single(pi) * single(pi)))) * (single(-0.16666666666666666) * (x * x)))) * (single(1.0) + ((single(-0.16666666666666666) * (single(pi) * single(pi))) * (x * x))); end
\begin{array}{l}
\\
\left(1 + \left(tau \cdot \left(tau \cdot \left(\pi \cdot \pi\right)\right)\right) \cdot \left(-0.16666666666666666 \cdot \left(x \cdot x\right)\right)\right) \cdot \left(1 + \left(-0.16666666666666666 \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(x \cdot x\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.6%
Simplified82.6%
Taylor expanded in tau around 0
+-commutativeN/A
associate-+r+N/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
Simplified76.6%
Final simplification76.6%
(FPCore (x tau) :precision binary32 (* (+ 1.0 (* x (* x (* PI (* PI -0.16666666666666666))))) (+ 1.0 (* (* -0.16666666666666666 (* tau tau)) (* (* PI PI) (* x x))))))
float code(float x, float tau) {
return (1.0f + (x * (x * (((float) M_PI) * (((float) M_PI) * -0.16666666666666666f))))) * (1.0f + ((-0.16666666666666666f * (tau * tau)) * ((((float) M_PI) * ((float) M_PI)) * (x * x))));
}
function code(x, tau) return Float32(Float32(Float32(1.0) + Float32(x * Float32(x * Float32(Float32(pi) * Float32(Float32(pi) * Float32(-0.16666666666666666)))))) * Float32(Float32(1.0) + Float32(Float32(Float32(-0.16666666666666666) * Float32(tau * tau)) * Float32(Float32(Float32(pi) * Float32(pi)) * Float32(x * x))))) end
function tmp = code(x, tau) tmp = (single(1.0) + (x * (x * (single(pi) * (single(pi) * single(-0.16666666666666666)))))) * (single(1.0) + ((single(-0.16666666666666666) * (tau * tau)) * ((single(pi) * single(pi)) * (x * x)))); end
\begin{array}{l}
\\
\left(1 + x \cdot \left(x \cdot \left(\pi \cdot \left(\pi \cdot -0.16666666666666666\right)\right)\right)\right) \cdot \left(1 + \left(-0.16666666666666666 \cdot \left(tau \cdot tau\right)\right) \cdot \left(\left(\pi \cdot \pi\right) \cdot \left(x \cdot x\right)\right)\right)
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3282.7%
Simplified82.7%
Taylor expanded in x around 0
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3276.6%
Simplified76.6%
Final simplification76.6%
(FPCore (x tau) :precision binary32 (+ 1.0 (* (* x x) (* PI (* PI (* -0.16666666666666666 (+ 1.0 (* tau tau))))))))
float code(float x, float tau) {
return 1.0f + ((x * x) * (((float) M_PI) * (((float) M_PI) * (-0.16666666666666666f * (1.0f + (tau * tau))))));
}
function code(x, tau) return Float32(Float32(1.0) + Float32(Float32(x * x) * Float32(Float32(pi) * Float32(Float32(pi) * Float32(Float32(-0.16666666666666666) * Float32(Float32(1.0) + Float32(tau * tau))))))) end
function tmp = code(x, tau) tmp = single(1.0) + ((x * x) * (single(pi) * (single(pi) * (single(-0.16666666666666666) * (single(1.0) + (tau * tau)))))); end
\begin{array}{l}
\\
1 + \left(x \cdot x\right) \cdot \left(\pi \cdot \left(\pi \cdot \left(-0.16666666666666666 \cdot \left(1 + tau \cdot tau\right)\right)\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3275.6%
Simplified75.6%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
*-commutativeN/A
distribute-lft1-inN/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f3275.6%
Applied egg-rr75.6%
Final simplification75.6%
(FPCore (x tau)
:precision binary32
(+
1.0
(*
(* x x)
(*
(* PI PI)
(+ -0.16666666666666666 (* -0.16666666666666666 (* tau tau)))))))
float code(float x, float tau) {
return 1.0f + ((x * x) * ((((float) M_PI) * ((float) M_PI)) * (-0.16666666666666666f + (-0.16666666666666666f * (tau * tau)))));
}
function code(x, tau) return Float32(Float32(1.0) + Float32(Float32(x * x) * Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-0.16666666666666666) + Float32(Float32(-0.16666666666666666) * Float32(tau * tau)))))) end
function tmp = code(x, tau) tmp = single(1.0) + ((x * x) * ((single(pi) * single(pi)) * (single(-0.16666666666666666) + (single(-0.16666666666666666) * (tau * tau))))); end
\begin{array}{l}
\\
1 + \left(x \cdot x\right) \cdot \left(\left(\pi \cdot \pi\right) \cdot \left(-0.16666666666666666 + -0.16666666666666666 \cdot \left(tau \cdot tau\right)\right)\right)
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3275.6%
Simplified75.6%
(FPCore (x tau) :precision binary32 (+ 1.0 (* (* -0.16666666666666666 (* tau tau)) (* (* PI PI) (* x x)))))
float code(float x, float tau) {
return 1.0f + ((-0.16666666666666666f * (tau * tau)) * ((((float) M_PI) * ((float) M_PI)) * (x * x)));
}
function code(x, tau) return Float32(Float32(1.0) + Float32(Float32(Float32(-0.16666666666666666) * Float32(tau * tau)) * Float32(Float32(Float32(pi) * Float32(pi)) * Float32(x * x)))) end
function tmp = code(x, tau) tmp = single(1.0) + ((single(-0.16666666666666666) * (tau * tau)) * ((single(pi) * single(pi)) * (x * x))); end
\begin{array}{l}
\\
1 + \left(-0.16666666666666666 \cdot \left(tau \cdot tau\right)\right) \cdot \left(\left(\pi \cdot \pi\right) \cdot \left(x \cdot x\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3275.6%
Simplified75.6%
Taylor expanded in tau around inf
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3266.4%
Simplified66.4%
Final simplification66.4%
(FPCore (x tau) :precision binary32 (* x (/ (+ 1.0 (* (* -0.16666666666666666 (* PI PI)) (* x x))) x)))
float code(float x, float tau) {
return x * ((1.0f + ((-0.16666666666666666f * (((float) M_PI) * ((float) M_PI))) * (x * x))) / x);
}
function code(x, tau) return Float32(x * Float32(Float32(Float32(1.0) + Float32(Float32(Float32(-0.16666666666666666) * Float32(Float32(pi) * Float32(pi))) * Float32(x * x))) / x)) end
function tmp = code(x, tau) tmp = x * ((single(1.0) + ((single(-0.16666666666666666) * (single(pi) * single(pi))) * (x * x))) / x); end
\begin{array}{l}
\\
x \cdot \frac{1 + \left(-0.16666666666666666 \cdot \left(\pi \cdot \pi\right)\right) \cdot \left(x \cdot x\right)}{x}
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3275.6%
Simplified75.6%
Taylor expanded in tau around 0
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3261.8%
Simplified61.8%
flip-+N/A
clear-numN/A
clear-numN/A
flip-+N/A
*-inversesN/A
associate-/r*N/A
clear-numN/A
associate-/l*N/A
Applied egg-rr61.8%
Final simplification61.8%
(FPCore (x tau) :precision binary32 (- (* PI (* (* (* x x) 0.16666666666666666) (- PI))) -1.0))
float code(float x, float tau) {
return (((float) M_PI) * (((x * x) * 0.16666666666666666f) * -((float) M_PI))) - -1.0f;
}
function code(x, tau) return Float32(Float32(Float32(pi) * Float32(Float32(Float32(x * x) * Float32(0.16666666666666666)) * Float32(-Float32(pi)))) - Float32(-1.0)) end
function tmp = code(x, tau) tmp = (single(pi) * (((x * x) * single(0.16666666666666666)) * -single(pi))) - single(-1.0); end
\begin{array}{l}
\\
\pi \cdot \left(\left(\left(x \cdot x\right) \cdot 0.16666666666666666\right) \cdot \left(-\pi\right)\right) - -1
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3275.6%
Simplified75.6%
Taylor expanded in tau around 0
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3261.8%
Simplified61.8%
flip-+N/A
clear-numN/A
clear-numN/A
flip-+N/A
*-inversesN/A
associate-/r*N/A
clear-numN/A
frac-2negN/A
Applied egg-rr61.8%
Final simplification61.8%
(FPCore (x tau) :precision binary32 (+ 1.0 (* (* PI PI) (* -0.16666666666666666 (* x x)))))
float code(float x, float tau) {
return 1.0f + ((((float) M_PI) * ((float) M_PI)) * (-0.16666666666666666f * (x * x)));
}
function code(x, tau) return Float32(Float32(1.0) + Float32(Float32(Float32(pi) * Float32(pi)) * Float32(Float32(-0.16666666666666666) * Float32(x * x)))) end
function tmp = code(x, tau) tmp = single(1.0) + ((single(pi) * single(pi)) * (single(-0.16666666666666666) * (x * x))); end
\begin{array}{l}
\\
1 + \left(\pi \cdot \pi\right) \cdot \left(-0.16666666666666666 \cdot \left(x \cdot x\right)\right)
\end{array}
Initial program 97.3%
associate-/l/N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sin-lowering-sin.f32N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3297.2%
Applied egg-rr97.2%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3275.6%
Simplified75.6%
Taylor expanded in tau around 0
+-lowering-+.f32N/A
associate-*r*N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
PI-lowering-PI.f32N/A
PI-lowering-PI.f3261.8%
Simplified61.8%
Final simplification61.8%
(FPCore (x tau) :precision binary32 1.0)
float code(float x, float tau) {
return 1.0f;
}
real(4) function code(x, tau)
real(4), intent (in) :: x
real(4), intent (in) :: tau
code = 1.0e0
end function
function code(x, tau) return Float32(1.0) end
function tmp = code(x, tau) tmp = single(1.0); end
\begin{array}{l}
\\
1
\end{array}
Initial program 97.3%
Taylor expanded in x around 0
Simplified60.9%
herbie shell --seed 2024138
(FPCore (x tau)
:name "Lanczos kernel"
:precision binary32
:pre (and (and (<= 1e-5 x) (<= x 1.0)) (and (<= 1.0 tau) (<= tau 5.0)))
(* (/ (sin (* (* x PI) tau)) (* (* x PI) tau)) (/ (sin (* x PI)) (* x PI))))