
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1.0 v)) 2.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinhf((1.0f / v)) * 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(-((sintheta_i * sintheta_o) / v)) * ((costheta_i * costheta_o) / v)) / ((sinh((1.0e0 / v)) * 2.0e0) * v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(-Float32(Float32(sinTheta_i * sinTheta_O) / v))) * Float32(Float32(cosTheta_i * cosTheta_O) / v)) / Float32(Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0)) * v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinh((single(1.0) / v)) * single(2.0)) * v); end
\begin{array}{l}
\\
\frac{e^{-\frac{sinTheta\_i \cdot sinTheta\_O}{v}} \cdot \frac{cosTheta\_i \cdot cosTheta\_O}{v}}{\left(\sinh \left(\frac{1}{v}\right) \cdot 2\right) \cdot v}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1.0 v)) 2.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinhf((1.0f / v)) * 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(-((sintheta_i * sintheta_o) / v)) * ((costheta_i * costheta_o) / v)) / ((sinh((1.0e0 / v)) * 2.0e0) * v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(-Float32(Float32(sinTheta_i * sinTheta_O) / v))) * Float32(Float32(cosTheta_i * cosTheta_O) / v)) / Float32(Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0)) * v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinh((single(1.0) / v)) * single(2.0)) * v); end
\begin{array}{l}
\\
\frac{e^{-\frac{sinTheta\_i \cdot sinTheta\_O}{v}} \cdot \frac{cosTheta\_i \cdot cosTheta\_O}{v}}{\left(\sinh \left(\frac{1}{v}\right) \cdot 2\right) \cdot v}
\end{array}
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* (/ (pow (exp sinTheta_i) (/ sinTheta_O v)) v) (/ cosTheta_i (sinh (/ 1.0 v)))) (/ cosTheta_O (* v 2.0))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((powf(expf(sinTheta_i), (sinTheta_O / v)) / v) * (cosTheta_i / sinhf((1.0f / v)))) * (cosTheta_O / (v * 2.0f));
}
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(sintheta_i) ** (sintheta_o / v)) / v) * (costheta_i / sinh((1.0e0 / v)))) * (costheta_o / (v * 2.0e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32((exp(sinTheta_i) ^ Float32(sinTheta_O / v)) / v) * Float32(cosTheta_i / sinh(Float32(Float32(1.0) / v)))) * Float32(cosTheta_O / Float32(v * Float32(2.0)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (((exp(sinTheta_i) ^ (sinTheta_O / v)) / v) * (cosTheta_i / sinh((single(1.0) / v)))) * (cosTheta_O / (v * single(2.0))); end
\begin{array}{l}
\\
\left(\frac{{\left(e^{sinTheta\_i}\right)}^{\left(\frac{sinTheta\_O}{v}\right)}}{v} \cdot \frac{cosTheta\_i}{\sinh \left(\frac{1}{v}\right)}\right) \cdot \frac{cosTheta\_O}{v \cdot 2}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ cosTheta_O (* v 2.0)) (/ (* cosTheta_i (/ (exp (* sinTheta_O (/ sinTheta_i v))) v)) (sinh (/ 1.0 v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O / (v * 2.0f)) * ((cosTheta_i * (expf((sinTheta_O * (sinTheta_i / v))) / v)) / sinhf((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 = (costheta_o / (v * 2.0e0)) * ((costheta_i * (exp((sintheta_o * (sintheta_i / v))) / v)) / sinh((1.0e0 / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O / Float32(v * Float32(2.0))) * Float32(Float32(cosTheta_i * Float32(exp(Float32(sinTheta_O * Float32(sinTheta_i / v))) / v)) / sinh(Float32(Float32(1.0) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O / (v * single(2.0))) * ((cosTheta_i * (exp((sinTheta_O * (sinTheta_i / v))) / v)) / sinh((single(1.0) / v))); end
\begin{array}{l}
\\
\frac{cosTheta\_O}{v \cdot 2} \cdot \frac{cosTheta\_i \cdot \frac{e^{sinTheta\_O \cdot \frac{sinTheta\_i}{v}}}{v}}{\sinh \left(\frac{1}{v}\right)}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
associate-*r/98.7%
Applied egg-rr98.7%
Taylor expanded in sinTheta_i around inf 98.5%
associate-/l*98.7%
associate-/l*98.7%
Simplified98.7%
Final simplification98.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ cosTheta_O (* v 2.0)) (* 2.0 (* cosTheta_i (/ 1.0 (* (sinh (/ 1.0 v)) (* v 2.0)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O / (v * 2.0f)) * (2.0f * (cosTheta_i * (1.0f / (sinhf((1.0f / v)) * (v * 2.0f)))));
}
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 = (costheta_o / (v * 2.0e0)) * (2.0e0 * (costheta_i * (1.0e0 / (sinh((1.0e0 / v)) * (v * 2.0e0)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O / Float32(v * Float32(2.0))) * Float32(Float32(2.0) * Float32(cosTheta_i * Float32(Float32(1.0) / Float32(sinh(Float32(Float32(1.0) / v)) * Float32(v * Float32(2.0))))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O / (v * single(2.0))) * (single(2.0) * (cosTheta_i * (single(1.0) / (sinh((single(1.0) / v)) * (v * single(2.0)))))); end
\begin{array}{l}
\\
\frac{cosTheta\_O}{v \cdot 2} \cdot \left(2 \cdot \left(cosTheta\_i \cdot \frac{1}{\sinh \left(\frac{1}{v}\right) \cdot \left(v \cdot 2\right)}\right)\right)
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
associate-*r/98.7%
Applied egg-rr98.7%
Taylor expanded in sinTheta_i around 0 98.5%
rec-exp98.6%
Simplified98.6%
div-inv98.6%
sinh-undef98.6%
Applied egg-rr98.6%
associate-*r*98.6%
*-commutative98.6%
*-commutative98.6%
Simplified98.6%
Final simplification98.6%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ cosTheta_O (* v 2.0)) (/ (/ cosTheta_i v) (sinh (/ 1.0 v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O / (v * 2.0f)) * ((cosTheta_i / v) / sinhf((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 = (costheta_o / (v * 2.0e0)) * ((costheta_i / v) / sinh((1.0e0 / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O / Float32(v * Float32(2.0))) * Float32(Float32(cosTheta_i / v) / sinh(Float32(Float32(1.0) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O / (v * single(2.0))) * ((cosTheta_i / v) / sinh((single(1.0) / v))); end
\begin{array}{l}
\\
\frac{cosTheta\_O}{v \cdot 2} \cdot \frac{\frac{cosTheta\_i}{v}}{\sinh \left(\frac{1}{v}\right)}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
associate-*r/98.7%
Applied egg-rr98.7%
Taylor expanded in sinTheta_i around 0 98.5%
Final simplification98.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (pow (/ v (* cosTheta_i cosTheta_O)) -1.0)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * powf((v / (cosTheta_i * cosTheta_O)), -1.0f);
}
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 / (costheta_i * costheta_o)) ** (-1.0e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * (Float32(v / Float32(cosTheta_i * cosTheta_O)) ^ Float32(-1.0))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * ((v / (cosTheta_i * cosTheta_O)) ^ single(-1.0)); end
\begin{array}{l}
\\
0.5 \cdot {\left(\frac{v}{cosTheta\_i \cdot cosTheta\_O}\right)}^{-1}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
Taylor expanded in v around inf 54.1%
*-commutative54.1%
Simplified54.1%
clear-num54.8%
inv-pow54.8%
Applied egg-rr54.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(/
(*
0.5
(+
(* cosTheta_i cosTheta_O)
(* cosTheta_O (* cosTheta_i (* sinTheta_O (/ sinTheta_i v))))))
v))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f * ((cosTheta_i * cosTheta_O) + (cosTheta_O * (cosTheta_i * (sinTheta_O * (sinTheta_i / 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 * ((costheta_i * costheta_o) + (costheta_o * (costheta_i * (sintheta_o * (sintheta_i / v)))))) / v
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) * Float32(Float32(cosTheta_i * cosTheta_O) + Float32(cosTheta_O * Float32(cosTheta_i * Float32(sinTheta_O * Float32(sinTheta_i / v)))))) / v) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) * ((cosTheta_i * cosTheta_O) + (cosTheta_O * (cosTheta_i * (sinTheta_O * (sinTheta_i / v)))))) / v; end
\begin{array}{l}
\\
\frac{0.5 \cdot \left(cosTheta\_i \cdot cosTheta\_O + cosTheta\_O \cdot \left(cosTheta\_i \cdot \left(sinTheta\_O \cdot \frac{sinTheta\_i}{v}\right)\right)\right)}{v}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
Taylor expanded in v around inf 54.1%
distribute-lft-out54.1%
*-commutative54.1%
associate-/l*54.1%
associate-/l*54.1%
associate-*r/54.1%
Simplified54.1%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* 0.5 (* cosTheta_i cosTheta_O)) v))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f * (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 = (0.5e0 * (costheta_i * costheta_o)) / v
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) * Float32(cosTheta_i * cosTheta_O)) / v) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) * (cosTheta_i * cosTheta_O)) / v; end
\begin{array}{l}
\\
\frac{0.5 \cdot \left(cosTheta\_i \cdot cosTheta\_O\right)}{v}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
Taylor expanded in v around inf 54.1%
*-commutative54.1%
Simplified54.1%
associate-*r/54.1%
Applied egg-rr54.1%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_i (/ cosTheta_O (* v 2.0))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i * (cosTheta_O / (v * 2.0f));
}
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 = costheta_i * (costheta_o / (v * 2.0e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i * Float32(cosTheta_O / Float32(v * Float32(2.0)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = cosTheta_i * (cosTheta_O / (v * single(2.0))); end
\begin{array}{l}
\\
cosTheta\_i \cdot \frac{cosTheta\_O}{v \cdot 2}
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
Taylor expanded in v around inf 54.1%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (* cosTheta_i (/ cosTheta_O v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (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 = 0.5e0 * (costheta_i * (costheta_o / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_i * Float32(cosTheta_O / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * (cosTheta_i * (cosTheta_O / v)); end
\begin{array}{l}
\\
0.5 \cdot \left(cosTheta\_i \cdot \frac{cosTheta\_O}{v}\right)
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
Taylor expanded in v around inf 54.1%
*-commutative54.1%
Simplified54.1%
associate-*r/54.1%
*-commutative54.1%
Applied egg-rr54.1%
Final simplification54.1%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (* cosTheta_O (/ cosTheta_i v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (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 = 0.5e0 * (costheta_o * (costheta_i / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_O * Float32(cosTheta_i / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * (cosTheta_O * (cosTheta_i / v)); end
\begin{array}{l}
\\
0.5 \cdot \left(cosTheta\_O \cdot \frac{cosTheta\_i}{v}\right)
\end{array}
Initial program 98.5%
associate-*r/98.5%
associate-/l/98.5%
remove-double-neg98.5%
distribute-rgt-neg-out98.5%
distribute-rgt-neg-out98.5%
distribute-lft-neg-in98.5%
associate-*r/98.5%
associate-/l/98.5%
associate-*r/98.5%
Simplified98.5%
associate-*r*98.5%
associate-*r*98.5%
times-frac98.6%
Applied egg-rr98.5%
times-frac98.8%
associate-/r*98.8%
Simplified98.8%
associate-*r/98.7%
Applied egg-rr98.7%
Taylor expanded in sinTheta_i around inf 98.5%
associate-/l*98.7%
associate-/l*98.7%
Simplified98.7%
Taylor expanded in v around inf 54.1%
associate-/l*54.1%
Simplified54.1%
herbie shell --seed 2024098
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:name "HairBSDF, Mp, upper"
:precision binary32
:pre (and (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))) (< 0.1 v)) (<= v 1.5707964))
(/ (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1.0 v)) 2.0) v)))