
(FPCore (s u) :precision binary32 (* s (log (/ 1.0 (- 1.0 (* 4.0 u))))))
float code(float s, float u) {
return s * logf((1.0f / (1.0f - (4.0f * u))));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * log((1.0e0 / (1.0e0 - (4.0e0 * u))))
end function
function code(s, u) return Float32(s * log(Float32(Float32(1.0) / Float32(Float32(1.0) - Float32(Float32(4.0) * u))))) end
function tmp = code(s, u) tmp = s * log((single(1.0) / (single(1.0) - (single(4.0) * u)))); end
\begin{array}{l}
\\
s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (s u) :precision binary32 (* s (log (/ 1.0 (- 1.0 (* 4.0 u))))))
float code(float s, float u) {
return s * logf((1.0f / (1.0f - (4.0f * u))));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * log((1.0e0 / (1.0e0 - (4.0e0 * u))))
end function
function code(s, u) return Float32(s * log(Float32(Float32(1.0) / Float32(Float32(1.0) - Float32(Float32(4.0) * u))))) end
function tmp = code(s, u) tmp = s * log((single(1.0) / (single(1.0) - (single(4.0) * u)))); end
\begin{array}{l}
\\
s \cdot \log \left(\frac{1}{1 - 4 \cdot u}\right)
\end{array}
(FPCore (s u) :precision binary32 (* (- s) (log1p (* u -4.0))))
float code(float s, float u) {
return -s * log1pf((u * -4.0f));
}
function code(s, u) return Float32(Float32(-s) * log1p(Float32(u * Float32(-4.0)))) end
\begin{array}{l}
\\
\left(-s\right) \cdot \mathsf{log1p}\left(u \cdot -4\right)
\end{array}
Initial program 59.9%
log-rec62.7%
distribute-rgt-neg-out62.7%
distribute-lft-neg-out62.7%
cancel-sign-sub-inv62.7%
log1p-def99.2%
*-commutative99.2%
metadata-eval99.2%
Simplified99.2%
Final simplification99.2%
(FPCore (s u) :precision binary32 (* s (* u (- (- -4.0) (+ (* u -8.0) (* u (* u -21.333333333333332)))))))
float code(float s, float u) {
return s * (u * (-(-4.0f) - ((u * -8.0f) + (u * (u * -21.333333333333332f)))));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * (u * (-(-4.0e0) - ((u * (-8.0e0)) + (u * (u * (-21.333333333333332e0))))))
end function
function code(s, u) return Float32(s * Float32(u * Float32(Float32(-Float32(-4.0)) - Float32(Float32(u * Float32(-8.0)) + Float32(u * Float32(u * Float32(-21.333333333333332))))))) end
function tmp = code(s, u) tmp = s * (u * (-single(-4.0) - ((u * single(-8.0)) + (u * (u * single(-21.333333333333332)))))); end
\begin{array}{l}
\\
s \cdot \left(u \cdot \left(\left(--4\right) - \left(u \cdot -8 + u \cdot \left(u \cdot -21.333333333333332\right)\right)\right)\right)
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Taylor expanded in u around 0 91.1%
+-commutative91.1%
+-commutative91.1%
associate-+l+91.1%
*-commutative91.1%
unpow291.1%
associate-*r*91.1%
unpow391.1%
unpow291.1%
associate-*r*91.1%
distribute-rgt-out91.1%
distribute-lft-out90.9%
unpow290.9%
associate-*r*90.9%
distribute-rgt-out90.9%
*-commutative90.9%
Simplified90.9%
distribute-lft-in90.9%
Applied egg-rr90.9%
Final simplification90.9%
(FPCore (s u) :precision binary32 (* s (* (+ -4.0 (* u (+ -8.0 (* u -21.333333333333332)))) (- u))))
float code(float s, float u) {
return s * ((-4.0f + (u * (-8.0f + (u * -21.333333333333332f)))) * -u);
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * (((-4.0e0) + (u * ((-8.0e0) + (u * (-21.333333333333332e0))))) * -u)
end function
function code(s, u) return Float32(s * Float32(Float32(Float32(-4.0) + Float32(u * Float32(Float32(-8.0) + Float32(u * Float32(-21.333333333333332))))) * Float32(-u))) end
function tmp = code(s, u) tmp = s * ((single(-4.0) + (u * (single(-8.0) + (u * single(-21.333333333333332))))) * -u); end
\begin{array}{l}
\\
s \cdot \left(\left(-4 + u \cdot \left(-8 + u \cdot -21.333333333333332\right)\right) \cdot \left(-u\right)\right)
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Taylor expanded in u around 0 91.1%
+-commutative91.1%
+-commutative91.1%
associate-+l+91.1%
*-commutative91.1%
unpow291.1%
associate-*r*91.1%
unpow391.1%
unpow291.1%
associate-*r*91.1%
distribute-rgt-out91.1%
distribute-lft-out90.9%
unpow290.9%
associate-*r*90.9%
distribute-rgt-out90.9%
*-commutative90.9%
Simplified90.9%
Final simplification90.9%
(FPCore (s u) :precision binary32 (* s (+ (* u 4.0) (* u (* u 8.0)))))
float code(float s, float u) {
return s * ((u * 4.0f) + (u * (u * 8.0f)));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * ((u * 4.0e0) + (u * (u * 8.0e0)))
end function
function code(s, u) return Float32(s * Float32(Float32(u * Float32(4.0)) + Float32(u * Float32(u * Float32(8.0))))) end
function tmp = code(s, u) tmp = s * ((u * single(4.0)) + (u * (u * single(8.0)))); end
\begin{array}{l}
\\
s \cdot \left(u \cdot 4 + u \cdot \left(u \cdot 8\right)\right)
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Taylor expanded in u around 0 86.3%
*-commutative86.3%
associate-*l*86.6%
*-commutative86.6%
associate-*l*86.7%
distribute-lft-out86.7%
unpow286.7%
associate-*l*86.7%
distribute-lft-out86.5%
Simplified86.5%
distribute-rgt-in86.7%
*-commutative86.7%
Applied egg-rr86.7%
Final simplification86.7%
(FPCore (s u) :precision binary32 (* s (* u (+ 4.0 (* u 8.0)))))
float code(float s, float u) {
return s * (u * (4.0f + (u * 8.0f)));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * (u * (4.0e0 + (u * 8.0e0)))
end function
function code(s, u) return Float32(s * Float32(u * Float32(Float32(4.0) + Float32(u * Float32(8.0))))) end
function tmp = code(s, u) tmp = s * (u * (single(4.0) + (u * single(8.0)))); end
\begin{array}{l}
\\
s \cdot \left(u \cdot \left(4 + u \cdot 8\right)\right)
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Taylor expanded in u around 0 86.3%
*-commutative86.3%
associate-*l*86.6%
*-commutative86.6%
associate-*l*86.7%
distribute-lft-out86.7%
unpow286.7%
associate-*l*86.7%
distribute-lft-out86.5%
Simplified86.5%
Final simplification86.5%
(FPCore (s u) :precision binary32 (* 4.0 (* s u)))
float code(float s, float u) {
return 4.0f * (s * u);
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = 4.0e0 * (s * u)
end function
function code(s, u) return Float32(Float32(4.0) * Float32(s * u)) end
function tmp = code(s, u) tmp = single(4.0) * (s * u); end
\begin{array}{l}
\\
4 \cdot \left(s \cdot u\right)
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Taylor expanded in u around 0 74.0%
Final simplification74.0%
(FPCore (s u) :precision binary32 (* u (* s 4.0)))
float code(float s, float u) {
return u * (s * 4.0f);
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = u * (s * 4.0e0)
end function
function code(s, u) return Float32(u * Float32(s * Float32(4.0))) end
function tmp = code(s, u) tmp = u * (s * single(4.0)); end
\begin{array}{l}
\\
u \cdot \left(s \cdot 4\right)
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Taylor expanded in u around 0 74.0%
associate-*r*74.3%
Simplified74.3%
Final simplification74.3%
(FPCore (s u) :precision binary32 (* s 0.0))
float code(float s, float u) {
return s * 0.0f;
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * 0.0e0
end function
function code(s, u) return Float32(s * Float32(0.0)) end
function tmp = code(s, u) tmp = s * single(0.0); end
\begin{array}{l}
\\
s \cdot 0
\end{array}
Initial program 59.9%
log-rec62.7%
sub-neg62.7%
distribute-rgt-neg-out62.7%
+-commutative62.7%
distribute-rgt-neg-out62.7%
neg-mul-162.7%
metadata-eval62.7%
distribute-lft-in62.7%
metadata-eval62.7%
sub-neg62.7%
log-prod-0.0%
distribute-neg-in-0.0%
unsub-neg-0.0%
unsub-neg-0.0%
distribute-neg-in-0.0%
log-prod62.7%
sub-neg62.7%
metadata-eval62.7%
Simplified62.7%
Applied egg-rr16.6%
Final simplification16.6%
herbie shell --seed 2024020
(FPCore (s u)
:name "Disney BSSRDF, sample scattering profile, lower"
:precision binary32
:pre (and (and (<= 0.0 s) (<= s 256.0)) (and (<= 2.328306437e-10 u) (<= u 0.25)))
(* s (log (/ 1.0 (- 1.0 (* 4.0 u))))))