
(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 6 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 (* (log1p (* u -4.0)) (- s)))
float code(float s, float u) {
return log1pf((u * -4.0f)) * -s;
}
function code(s, u) return Float32(log1p(Float32(u * Float32(-4.0))) * Float32(-s)) end
\begin{array}{l}
\\
\mathsf{log1p}\left(u \cdot -4\right) \cdot \left(-s\right)
\end{array}
Initial program 59.5%
*-commutative59.5%
log-rec62.3%
distribute-lft-neg-out62.3%
distribute-rgt-neg-in62.3%
sub-neg62.3%
log1p-def99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
metadata-eval99.4%
Simplified99.4%
Final simplification99.4%
(FPCore (s u) :precision binary32 (* s (+ (* u (* u 8.0)) (* u 4.0))))
float code(float s, float u) {
return s * ((u * (u * 8.0f)) + (u * 4.0f));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * ((u * (u * 8.0e0)) + (u * 4.0e0))
end function
function code(s, u) return Float32(s * Float32(Float32(u * Float32(u * Float32(8.0))) + Float32(u * Float32(4.0)))) end
function tmp = code(s, u) tmp = s * ((u * (u * single(8.0))) + (u * single(4.0))); end
\begin{array}{l}
\\
s \cdot \left(u \cdot \left(u \cdot 8\right) + u \cdot 4\right)
\end{array}
Initial program 59.5%
Taylor expanded in u around 0 92.0%
associate-+r+92.1%
distribute-rgt-in92.0%
fma-def92.0%
unpow292.0%
*-commutative92.0%
*-commutative92.0%
Applied egg-rr92.0%
Taylor expanded in u around 0 87.9%
*-commutative87.9%
unpow287.9%
associate-*l*87.9%
Simplified87.9%
distribute-rgt-out88.0%
Applied egg-rr88.0%
Final simplification88.0%
(FPCore (s u) :precision binary32 (* s (* u (+ (* u 8.0) 4.0))))
float code(float s, float u) {
return s * (u * ((u * 8.0f) + 4.0f));
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * (u * ((u * 8.0e0) + 4.0e0))
end function
function code(s, u) return Float32(s * Float32(u * Float32(Float32(u * Float32(8.0)) + Float32(4.0)))) end
function tmp = code(s, u) tmp = s * (u * ((u * single(8.0)) + single(4.0))); end
\begin{array}{l}
\\
s \cdot \left(u \cdot \left(u \cdot 8 + 4\right)\right)
\end{array}
Initial program 59.5%
flip3--59.0%
associate-/r/58.8%
log-prod58.9%
metadata-eval58.9%
*-commutative58.9%
unpow-prod-down58.9%
metadata-eval58.9%
metadata-eval58.9%
log1p-udef96.1%
*-un-lft-identity96.1%
distribute-lft1-in96.0%
fma-def96.0%
Applied egg-rr96.0%
+-commutative96.0%
log-rec96.7%
sub-neg96.7%
log1p-def98.8%
unsub-neg98.8%
*-commutative98.8%
*-commutative98.8%
associate-*l*98.8%
distribute-rgt-neg-in98.8%
metadata-eval98.8%
Simplified98.8%
Taylor expanded in u around 0 87.6%
*-commutative87.6%
associate-*r*87.9%
*-commutative87.9%
*-commutative87.9%
associate-*l*87.9%
unpow287.9%
distribute-rgt-out88.0%
*-commutative88.0%
associate-*r*88.0%
distribute-rgt-out87.7%
*-commutative87.7%
Simplified87.7%
Final simplification87.7%
(FPCore (s u) :precision binary32 (* 4.0 (* u s)))
float code(float s, float u) {
return 4.0f * (u * s);
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = 4.0e0 * (u * s)
end function
function code(s, u) return Float32(Float32(4.0) * Float32(u * s)) end
function tmp = code(s, u) tmp = single(4.0) * (u * s); end
\begin{array}{l}
\\
4 \cdot \left(u \cdot s\right)
\end{array}
Initial program 59.5%
Taylor expanded in u around 0 75.2%
*-commutative75.2%
Simplified75.2%
Final simplification75.2%
(FPCore (s u) :precision binary32 (* s (* u 4.0)))
float code(float s, float u) {
return s * (u * 4.0f);
}
real(4) function code(s, u)
real(4), intent (in) :: s
real(4), intent (in) :: u
code = s * (u * 4.0e0)
end function
function code(s, u) return Float32(s * Float32(u * Float32(4.0))) end
function tmp = code(s, u) tmp = s * (u * single(4.0)); end
\begin{array}{l}
\\
s \cdot \left(u \cdot 4\right)
\end{array}
Initial program 59.5%
Taylor expanded in u around 0 75.5%
Final simplification75.5%
(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.5%
Applied egg-rr17.0%
+-inverses17.0%
Simplified17.0%
Final simplification17.0%
herbie shell --seed 2023258
(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))))))