
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(exp
(+
(+
(-
(- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_i sinTheta_O) v))
(/ 1.0 v))
0.6931)
(log (/ 1.0 (* 2.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (1.0f / v)) + 0.6931f) + logf((1.0f / (2.0f * v)))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = exp(((((((costheta_i * costheta_o) / v) - ((sintheta_i * sintheta_o) / v)) - (1.0e0 / v)) + 0.6931e0) + log((1.0e0 / (2.0e0 * v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) / v) - Float32(Float32(sinTheta_i * sinTheta_O) / v)) - Float32(Float32(1.0) / v)) + Float32(0.6931)) + log(Float32(Float32(1.0) / Float32(Float32(2.0) * v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (single(1.0) / v)) + single(0.6931)) + log((single(1.0) / (single(2.0) * v))))); end
\begin{array}{l}
\\
e^{\left(\left(\left(\frac{cosTheta\_i \cdot cosTheta\_O}{v} - \frac{sinTheta\_i \cdot sinTheta\_O}{v}\right) - \frac{1}{v}\right) + 0.6931\right) + \log \left(\frac{1}{2 \cdot v}\right)}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(exp
(+
(+
(-
(- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_i sinTheta_O) v))
(/ 1.0 v))
0.6931)
(log (/ 1.0 (* 2.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (1.0f / v)) + 0.6931f) + logf((1.0f / (2.0f * v)))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = exp(((((((costheta_i * costheta_o) / v) - ((sintheta_i * sintheta_o) / v)) - (1.0e0 / v)) + 0.6931e0) + log((1.0e0 / (2.0e0 * v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) / v) - Float32(Float32(sinTheta_i * sinTheta_O) / v)) - Float32(Float32(1.0) / v)) + Float32(0.6931)) + log(Float32(Float32(1.0) / Float32(Float32(2.0) * v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (single(1.0) / v)) + single(0.6931)) + log((single(1.0) / (single(2.0) * v))))); end
\begin{array}{l}
\\
e^{\left(\left(\left(\frac{cosTheta\_i \cdot cosTheta\_O}{v} - \frac{sinTheta\_i \cdot sinTheta\_O}{v}\right) - \frac{1}{v}\right) + 0.6931\right) + \log \left(\frac{1}{2 \cdot v}\right)}
\end{array}
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 1.0 (* (exp (* (log (/ 0.5 v)) -1.0)) (exp (- (/ (- -1.0 (* sinTheta_i sinTheta_O)) (- v)) 0.6931)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 1.0f / (expf((logf((0.5f / v)) * -1.0f)) * expf((((-1.0f - (sinTheta_i * sinTheta_O)) / -v) - 0.6931f)));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 1.0e0 / (exp((log((0.5e0 / v)) * (-1.0e0))) * exp(((((-1.0e0) - (sintheta_i * sintheta_o)) / -v) - 0.6931e0)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(1.0) / Float32(exp(Float32(log(Float32(Float32(0.5) / v)) * Float32(-1.0))) * exp(Float32(Float32(Float32(Float32(-1.0) - Float32(sinTheta_i * sinTheta_O)) / Float32(-v)) - Float32(0.6931))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(1.0) / (exp((log((single(0.5) / v)) * single(-1.0))) * exp((((single(-1.0) - (sinTheta_i * sinTheta_O)) / -v) - single(0.6931)))); end
\begin{array}{l}
\\
\frac{1}{e^{\log \left(\frac{0.5}{v}\right) \cdot -1} \cdot e^{\frac{-1 - sinTheta\_i \cdot sinTheta\_O}{-v} - 0.6931}}
\end{array}
Initial program 99.7%
lift-exp.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-log.f32N/A
lift-/.f32N/A
log-divN/A
associate-+l-N/A
exp-diffN/A
metadata-evalN/A
1-expN/A
lower-/.f32N/A
Applied rewrites99.7%
Taylor expanded in cosTheta_i around 0
distribute-lft-inN/A
metadata-evalN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
*-commutativeN/A
lower-*.f3299.7
Applied rewrites99.7%
lift-exp.f32N/A
lift--.f32N/A
exp-diffN/A
lift-log.f32N/A
rem-exp-logN/A
div-invN/A
lower-*.f32N/A
rec-expN/A
lower-exp.f32N/A
lower-neg.f3299.8
Applied rewrites99.8%
lift-*.f32N/A
metadata-evalN/A
associate-/r/N/A
lift-/.f32N/A
inv-powN/A
pow-to-expN/A
lower-exp.f32N/A
lift-/.f32N/A
lower-*.f32N/A
lift-/.f32N/A
lower-log.f3299.8
Applied rewrites99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 0.5 (* (exp (- -0.6931 (/ (- -1.0 (* sinTheta_i sinTheta_O)) v))) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f / (expf((-0.6931f - ((-1.0f - (sinTheta_i * sinTheta_O)) / v))) * v);
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 0.5e0 / (exp(((-0.6931e0) - (((-1.0e0) - (sintheta_i * sintheta_o)) / v))) * v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) / Float32(exp(Float32(Float32(-0.6931) - Float32(Float32(Float32(-1.0) - Float32(sinTheta_i * sinTheta_O)) / v))) * v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) / (exp((single(-0.6931) - ((single(-1.0) - (sinTheta_i * sinTheta_O)) / v))) * v); end
\begin{array}{l}
\\
\frac{0.5}{e^{-0.6931 - \frac{-1 - sinTheta\_i \cdot sinTheta\_O}{v}} \cdot v}
\end{array}
Initial program 99.7%
lift-exp.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-log.f32N/A
lift-/.f32N/A
log-divN/A
associate-+l-N/A
exp-diffN/A
metadata-evalN/A
1-expN/A
lower-/.f32N/A
Applied rewrites99.7%
Taylor expanded in cosTheta_i around 0
distribute-lft-inN/A
metadata-evalN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
*-commutativeN/A
lower-*.f3299.7
Applied rewrites99.7%
lift-exp.f32N/A
lift--.f32N/A
exp-diffN/A
lift-log.f32N/A
rem-exp-logN/A
div-invN/A
lower-*.f32N/A
rec-expN/A
lower-exp.f32N/A
lower-neg.f3299.8
Applied rewrites99.8%
lift-/.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-*l*N/A
associate-/r*N/A
metadata-evalN/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f3299.8
Applied rewrites99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (/ 0.5 v) (exp (- (/ 1.0 v) 0.6931))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f / v) / expf(((1.0f / v) - 0.6931f));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (0.5e0 / v) / exp(((1.0e0 / v) - 0.6931e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) / v) / exp(Float32(Float32(Float32(1.0) / v) - Float32(0.6931)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) / v) / exp(((single(1.0) / v) - single(0.6931))); end
\begin{array}{l}
\\
\frac{\frac{0.5}{v}}{e^{\frac{1}{v} - 0.6931}}
\end{array}
Initial program 99.7%
lift-exp.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-+.f32N/A
lift--.f32N/A
associate-+l-N/A
associate-+r-N/A
exp-diffN/A
lower-/.f32N/A
Applied rewrites86.9%
Taylor expanded in v around -inf
+-commutativeN/A
exp-sumN/A
rem-exp-logN/A
lower-*.f32N/A
rem-exp-logN/A
lower-/.f3299.8
Applied rewrites99.8%
Applied rewrites99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 1.0 (exp (/ (- 1.0 (* cosTheta_i cosTheta_O)) v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 1.0f / expf(((1.0f - (cosTheta_i * cosTheta_O)) / v));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 1.0e0 / exp(((1.0e0 - (costheta_i * costheta_o)) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(1.0) / exp(Float32(Float32(Float32(1.0) - Float32(cosTheta_i * cosTheta_O)) / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(1.0) / exp(((single(1.0) - (cosTheta_i * cosTheta_O)) / v)); end
\begin{array}{l}
\\
\frac{1}{e^{\frac{1 - cosTheta\_i \cdot cosTheta\_O}{v}}}
\end{array}
Initial program 99.7%
lift-exp.f32N/A
lift-+.f32N/A
+-commutativeN/A
lift-log.f32N/A
lift-/.f32N/A
log-divN/A
associate-+l-N/A
exp-diffN/A
metadata-evalN/A
1-expN/A
lower-/.f32N/A
Applied rewrites99.7%
Taylor expanded in v around 0
lower-/.f32N/A
lower--.f32N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f32N/A
*-commutativeN/A
lower-*.f3297.4
Applied rewrites97.4%
Taylor expanded in sinTheta_i around 0
Applied rewrites98.1%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (/ (- (fma cosTheta_i cosTheta_O -1.0) (* sinTheta_i sinTheta_O)) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((fmaf(cosTheta_i, cosTheta_O, -1.0f) - (sinTheta_i * sinTheta_O)) / v));
}
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(fma(cosTheta_i, cosTheta_O, Float32(-1.0)) - Float32(sinTheta_i * sinTheta_O)) / v)) end
\begin{array}{l}
\\
e^{\frac{\mathsf{fma}\left(cosTheta\_i, cosTheta\_O, -1\right) - sinTheta\_i \cdot sinTheta\_O}{v}}
\end{array}
Initial program 99.7%
lift-+.f32N/A
lift-log.f32N/A
lift-/.f32N/A
log-recN/A
unsub-negN/A
lower--.f32N/A
Applied rewrites99.7%
Taylor expanded in v around 0
lower-/.f32N/A
associate--r+N/A
lower--.f32N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f32N/A
*-commutativeN/A
lower-*.f3298.1
Applied rewrites97.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (* cosTheta_O (/ cosTheta_i v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf((cosTheta_O * (cosTheta_i / v)));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = exp((costheta_o * (costheta_i / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(cosTheta_O * Float32(cosTheta_i / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp((cosTheta_O * (cosTheta_i / v))); end
\begin{array}{l}
\\
e^{cosTheta\_O \cdot \frac{cosTheta\_i}{v}}
\end{array}
Initial program 99.7%
lift-+.f32N/A
lift-log.f32N/A
lift-/.f32N/A
log-recN/A
unsub-negN/A
lower--.f32N/A
Applied rewrites99.7%
Taylor expanded in cosTheta_i around inf
lower-/.f32N/A
*-commutativeN/A
lower-*.f3213.5
Applied rewrites13.5%
Applied rewrites13.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 0.5 (* (exp -0.6931) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f / (expf(-0.6931f) * v);
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 0.5e0 / (exp((-0.6931e0)) * v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) / Float32(exp(Float32(-0.6931)) * v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) / (exp(single(-0.6931)) * v); end
\begin{array}{l}
\\
\frac{0.5}{e^{-0.6931} \cdot v}
\end{array}
Initial program 99.7%
lift-exp.f32N/A
lift-+.f32N/A
+-commutativeN/A
exp-sumN/A
lift-log.f32N/A
rem-exp-logN/A
lower-*.f32N/A
lift-/.f32N/A
lift-*.f32N/A
associate-/r*N/A
lower-/.f32N/A
metadata-evalN/A
lower-exp.f3299.8
lift-+.f32N/A
Applied rewrites99.8%
Applied rewrites87.3%
Taylor expanded in v around inf
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-exp.f324.6
Applied rewrites4.6%
herbie shell --seed 2024318
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:name "HairBSDF, Mp, lower"
:precision binary32
:pre (and (and (and (and (and (<= -1.0 cosTheta_i) (<= cosTheta_i 1.0)) (and (<= -1.0 cosTheta_O) (<= cosTheta_O 1.0))) (and (<= -1.0 sinTheta_i) (<= sinTheta_i 1.0))) (and (<= -1.0 sinTheta_O) (<= sinTheta_O 1.0))) (and (<= -1.5707964 v) (<= v 0.1)))
(exp (+ (+ (- (- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_i sinTheta_O) v)) (/ 1.0 v)) 0.6931) (log (/ 1.0 (* 2.0 v))))))