
(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 9 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 (* (* (/ 0.5 v) (exp (+ (/ (* cosTheta_O cosTheta_i) v) (/ -1.0 v)))) (exp 0.6931)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((0.5f / v) * expf((((cosTheta_O * cosTheta_i) / v) + (-1.0f / v)))) * expf(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((((costheta_o * costheta_i) / v) + ((-1.0e0) / v)))) * exp(0.6931e0)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(Float32(0.5) / v) * exp(Float32(Float32(Float32(cosTheta_O * cosTheta_i) / v) + Float32(Float32(-1.0) / v)))) * exp(Float32(0.6931))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((single(0.5) / v) * exp((((cosTheta_O * cosTheta_i) / v) + (single(-1.0) / v)))) * exp(single(0.6931)); end
\begin{array}{l}
\\
\left(\frac{0.5}{v} \cdot e^{\frac{cosTheta\_O \cdot cosTheta\_i}{v} + \frac{-1}{v}}\right) \cdot e^{0.6931}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
associate-+r+99.6%
fma-define99.6%
prod-exp99.7%
add-exp-log99.7%
*-commutative99.7%
exp-sum99.7%
fma-define99.7%
associate-*r*99.7%
Applied egg-rr99.7%
Taylor expanded in cosTheta_i around inf 99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (exp 0.6931) (/ (* 0.5 (exp (/ -1.0 v))) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(0.6931f) * ((0.5f * expf((-1.0f / 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 = exp(0.6931e0) * ((0.5e0 * exp(((-1.0e0) / v))) / v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(exp(Float32(0.6931)) * Float32(Float32(Float32(0.5) * exp(Float32(Float32(-1.0) / v))) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(single(0.6931)) * ((single(0.5) * exp((single(-1.0) / v))) / v); end
\begin{array}{l}
\\
e^{0.6931} \cdot \frac{0.5 \cdot e^{\frac{-1}{v}}}{v}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
associate-+r+99.6%
fma-define99.6%
prod-exp99.7%
add-exp-log99.7%
*-commutative99.7%
exp-sum99.7%
fma-define99.7%
associate-*r*99.7%
Applied egg-rr99.7%
Taylor expanded in cosTheta_i around inf 99.7%
Taylor expanded in cosTheta_O around 0 99.7%
associate-*r/99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ 0.5 v) (exp (+ (+ (/ (* cosTheta_O cosTheta_i) v) 0.6931) (/ -1.0 v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f / v) * expf(((((cosTheta_O * cosTheta_i) / v) + 0.6931f) + (-1.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 = (0.5e0 / v) * exp(((((costheta_o * costheta_i) / v) + 0.6931e0) + ((-1.0e0) / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) / v) * exp(Float32(Float32(Float32(Float32(cosTheta_O * cosTheta_i) / v) + Float32(0.6931)) + Float32(Float32(-1.0) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) / v) * exp(((((cosTheta_O * cosTheta_i) / v) + single(0.6931)) + (single(-1.0) / v))); end
\begin{array}{l}
\\
\frac{0.5}{v} \cdot e^{\left(\frac{cosTheta\_O \cdot cosTheta\_i}{v} + 0.6931\right) + \frac{-1}{v}}
\end{array}
Initial program 99.6%
exp-sum99.7%
*-commutative99.7%
rem-exp-log99.7%
associate-/r*99.7%
metadata-eval99.7%
associate--l-99.7%
associate-/l*99.7%
associate-/l*99.7%
fma-define99.7%
Simplified99.7%
Taylor expanded in sinTheta_i around 0 99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (/ (exp (+ 0.6931 (/ -1.0 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 / 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) / 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(-1.0) / 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) / v))) / v); end
\begin{array}{l}
\\
0.5 \cdot \frac{e^{0.6931 + \frac{-1}{v}}}{v}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
associate-+r+99.6%
fma-define99.6%
prod-exp99.7%
add-exp-log99.7%
*-commutative99.7%
exp-sum99.7%
fma-define99.7%
associate-*r*99.7%
Applied egg-rr99.7%
Taylor expanded in cosTheta_i around inf 99.7%
Taylor expanded in cosTheta_O around 0 99.7%
associate-*r/99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in v around 0 99.7%
metadata-eval99.7%
distribute-neg-frac99.7%
*-commutative99.7%
prod-exp99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (/ (+ (* cosTheta_O cosTheta_i) -1.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf((((cosTheta_O * cosTheta_i) + -1.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_o * costheta_i) + (-1.0e0)) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(cosTheta_O * cosTheta_i) + Float32(-1.0)) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp((((cosTheta_O * cosTheta_i) + single(-1.0)) / v)); end
\begin{array}{l}
\\
e^{\frac{cosTheta\_O \cdot cosTheta\_i + -1}{v}}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in sinTheta_i around 0 99.6%
associate-+r+99.6%
associate-/l*99.6%
Simplified99.6%
Taylor expanded in v around 0 97.5%
Final simplification97.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (+ -1.0 (- 1.0 (* sinTheta_O (/ sinTheta_i v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return -1.0f + (1.0f - (sinTheta_O * (sinTheta_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 = (-1.0e0) + (1.0e0 - (sintheta_o * (sintheta_i / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(-1.0) + Float32(Float32(1.0) - Float32(sinTheta_O * Float32(sinTheta_i / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(-1.0) + (single(1.0) - (sinTheta_O * (sinTheta_i / v))); end
\begin{array}{l}
\\
-1 + \left(1 - sinTheta\_O \cdot \frac{sinTheta\_i}{v}\right)
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in sinTheta_i around inf 10.6%
*-commutative10.6%
associate-*r/10.6%
neg-mul-110.6%
distribute-lft-neg-in10.6%
Simplified10.6%
Taylor expanded in sinTheta_i around 0 6.2%
mul-1-neg6.2%
unsub-neg6.2%
associate-/l*6.2%
Simplified6.2%
Taylor expanded in sinTheta_O around inf 43.0%
mul-1-neg43.0%
associate-*r/22.6%
distribute-lft-neg-out22.6%
*-commutative22.6%
Simplified22.6%
expm1-log1p-u22.2%
expm1-undefine76.7%
Applied egg-rr76.7%
sub-neg76.7%
log1p-undefine76.7%
distribute-rgt-neg-in76.7%
neg-mul-176.7%
*-commutative76.7%
associate-/l*76.7%
rem-exp-log77.0%
mul-1-neg77.0%
unsub-neg77.0%
associate-/l*77.0%
metadata-eval77.0%
Simplified77.0%
Final simplification77.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* sinTheta_i (/ sinTheta_O (- v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return sinTheta_i * (sinTheta_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 = sintheta_i * (sintheta_o / -v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(sinTheta_i * Float32(sinTheta_O / Float32(-v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = sinTheta_i * (sinTheta_O / -v); end
\begin{array}{l}
\\
sinTheta\_i \cdot \frac{sinTheta\_O}{-v}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in sinTheta_i around inf 10.6%
*-commutative10.6%
associate-*r/10.6%
neg-mul-110.6%
distribute-lft-neg-in10.6%
Simplified10.6%
Taylor expanded in sinTheta_i around 0 6.2%
mul-1-neg6.2%
unsub-neg6.2%
associate-/l*6.2%
Simplified6.2%
Taylor expanded in sinTheta_O around inf 43.0%
mul-1-neg43.0%
associate-*r/22.6%
distribute-lft-neg-out22.6%
*-commutative22.6%
Simplified22.6%
Taylor expanded in sinTheta_i around 0 43.0%
mul-1-neg43.0%
distribute-frac-neg43.0%
*-commutative43.0%
distribute-rgt-neg-out43.0%
associate-/l*22.6%
Simplified22.6%
Final simplification22.6%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* sinTheta_O (/ sinTheta_i (- v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return sinTheta_O * (sinTheta_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 = sintheta_o * (sintheta_i / -v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(sinTheta_O * Float32(sinTheta_i / Float32(-v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = sinTheta_O * (sinTheta_i / -v); end
\begin{array}{l}
\\
sinTheta\_O \cdot \frac{sinTheta\_i}{-v}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in sinTheta_i around inf 10.6%
*-commutative10.6%
associate-*r/10.6%
neg-mul-110.6%
distribute-lft-neg-in10.6%
Simplified10.6%
Taylor expanded in sinTheta_i around 0 6.2%
mul-1-neg6.2%
unsub-neg6.2%
associate-/l*6.2%
Simplified6.2%
Taylor expanded in sinTheta_O around inf 43.0%
mul-1-neg43.0%
associate-*r/22.6%
distribute-lft-neg-out22.6%
*-commutative22.6%
Simplified22.6%
Final simplification22.6%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* sinTheta_O (- sinTheta_i)) v))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (sinTheta_O * -sinTheta_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 = (sintheta_o * -sintheta_i) / v
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(sinTheta_O * Float32(-sinTheta_i)) / v) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (sinTheta_O * -sinTheta_i) / v; end
\begin{array}{l}
\\
\frac{sinTheta\_O \cdot \left(-sinTheta\_i\right)}{v}
\end{array}
Initial program 99.6%
associate-+l+99.6%
associate--l-99.6%
associate-/l*99.6%
associate-/l*99.6%
associate-/r*99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in sinTheta_i around inf 10.6%
*-commutative10.6%
associate-*r/10.6%
neg-mul-110.6%
distribute-lft-neg-in10.6%
Simplified10.6%
Taylor expanded in sinTheta_i around 0 6.2%
mul-1-neg6.2%
unsub-neg6.2%
associate-/l*6.2%
Simplified6.2%
Taylor expanded in sinTheta_O around inf 43.0%
mul-1-neg43.0%
associate-*r/22.6%
distribute-lft-neg-out22.6%
*-commutative22.6%
Simplified22.6%
associate-*l/43.0%
Applied egg-rr43.0%
Final simplification43.0%
herbie shell --seed 2024044
(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))))))