
(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 22 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}
cosTheta_O\_m = (fabs.f32 cosTheta_O)
cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
(FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v)
:precision binary32
(*
cosTheta_O_s
(/
(*
(* cosTheta_O_m (* (/ 1.0 v) cosTheta_i))
(fma sinTheta_O (/ (* -0.5 sinTheta_i) (* v v)) (/ 0.5 v)))
(sinh (/ 1.0 v)))))cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (((cosTheta_O_m * ((1.0f / v) * cosTheta_i)) * fmaf(sinTheta_O, ((-0.5f * sinTheta_i) / (v * v)), (0.5f / v))) / sinhf((1.0f / v)));
}
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(Float32(cosTheta_O_m * Float32(Float32(Float32(1.0) / v) * cosTheta_i)) * fma(sinTheta_O, Float32(Float32(Float32(-0.5) * sinTheta_i) / Float32(v * v)), Float32(Float32(0.5) / v))) / sinh(Float32(Float32(1.0) / v)))) end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{\left(cosTheta\_O\_m \cdot \left(\frac{1}{v} \cdot cosTheta\_i\right)\right) \cdot \mathsf{fma}\left(sinTheta\_O, \frac{-0.5 \cdot sinTheta\_i}{v \cdot v}, \frac{0.5}{v}\right)}{\sinh \left(\frac{1}{v}\right)}
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-*l*N/A
times-fracN/A
Applied egg-rr99.0%
clear-numN/A
associate-/r/N/A
lift-/.f32N/A
lower-*.f3299.1
Applied egg-rr99.1%
Taylor expanded in sinTheta_i around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f32N/A
associate-*r/N/A
lower-/.f32N/A
lower-*.f32N/A
unpow2N/A
lower-*.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f3299.1
Simplified99.1%
cosTheta_O\_m = (fabs.f32 cosTheta_O)
cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
(FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v)
:precision binary32
(*
cosTheta_O_s
(/
(*
(* cosTheta_O_m (* (/ 1.0 v) cosTheta_i))
(/ (fma (* sinTheta_O sinTheta_i) (/ -0.5 v) 0.5) v))
(sinh (/ 1.0 v)))))cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (((cosTheta_O_m * ((1.0f / v) * cosTheta_i)) * (fmaf((sinTheta_O * sinTheta_i), (-0.5f / v), 0.5f) / v)) / sinhf((1.0f / v)));
}
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(Float32(cosTheta_O_m * Float32(Float32(Float32(1.0) / v) * cosTheta_i)) * Float32(fma(Float32(sinTheta_O * sinTheta_i), Float32(Float32(-0.5) / v), Float32(0.5)) / v)) / sinh(Float32(Float32(1.0) / v)))) end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{\left(cosTheta\_O\_m \cdot \left(\frac{1}{v} \cdot cosTheta\_i\right)\right) \cdot \frac{\mathsf{fma}\left(sinTheta\_O \cdot sinTheta\_i, \frac{-0.5}{v}, 0.5\right)}{v}}{\sinh \left(\frac{1}{v}\right)}
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-*l*N/A
times-fracN/A
Applied egg-rr99.0%
clear-numN/A
associate-/r/N/A
lift-/.f32N/A
lower-*.f3299.1
Applied egg-rr99.1%
Taylor expanded in v around inf
lower-/.f32N/A
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
lower-fma.f32N/A
lower-*.f32N/A
lower-/.f3299.1
Simplified99.1%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ (* (* cosTheta_O_m (* (/ 1.0 v) cosTheta_i)) (/ 0.5 v)) (sinh (/ 1.0 v)))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (((cosTheta_O_m * ((1.0f / v) * cosTheta_i)) * (0.5f / v)) / sinhf((1.0f / v)));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (((costheta_o_m * ((1.0e0 / v) * costheta_i)) * (0.5e0 / v)) / sinh((1.0e0 / v)))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(Float32(cosTheta_O_m * Float32(Float32(Float32(1.0) / v) * cosTheta_i)) * Float32(Float32(0.5) / v)) / sinh(Float32(Float32(1.0) / v)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (((cosTheta_O_m * ((single(1.0) / v) * cosTheta_i)) * (single(0.5) / v)) / sinh((single(1.0) / v)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{\left(cosTheta\_O\_m \cdot \left(\frac{1}{v} \cdot cosTheta\_i\right)\right) \cdot \frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)}
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-*l*N/A
times-fracN/A
Applied egg-rr99.0%
clear-numN/A
associate-/r/N/A
lift-/.f32N/A
lower-*.f3299.1
Applied egg-rr99.1%
Taylor expanded in sinTheta_i around 0
lower-/.f3298.9
Simplified98.9%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ (* (/ cosTheta_i v) (/ cosTheta_O_m v)) (* (sinh (/ 1.0 v)) 2.0))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (((cosTheta_i / v) * (cosTheta_O_m / v)) / (sinhf((1.0f / v)) * 2.0f));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (((costheta_i / v) * (costheta_o_m / v)) / (sinh((1.0e0 / v)) * 2.0e0))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(Float32(cosTheta_i / v) * Float32(cosTheta_O_m / v)) / Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (((cosTheta_i / v) * (cosTheta_O_m / v)) / (sinh((single(1.0) / v)) * single(2.0)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{\frac{cosTheta\_i}{v} \cdot \frac{cosTheta\_O\_m}{v}}{\sinh \left(\frac{1}{v}\right) \cdot 2}
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
neg-mul-1N/A
lift-exp.f32N/A
neg-mul-1N/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
lift-exp.f32N/A
lift--.f32N/A
lift-*.f32N/A
times-fracN/A
associate-*l/N/A
lower-/.f32N/A
Applied egg-rr98.4%
lift-*.f32N/A
associate-*r/N/A
lift-*.f32N/A
times-fracN/A
lower-*.f32N/A
lower-/.f32N/A
lower-/.f3298.7
Applied egg-rr98.7%
Final simplification98.7%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* (/ cosTheta_i v) (/ cosTheta_O_m (* v (* (sinh (/ 1.0 v)) 2.0))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * ((cosTheta_i / v) * (cosTheta_O_m / (v * (sinhf((1.0f / v)) * 2.0f))));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * ((costheta_i / v) * (costheta_o_m / (v * (sinh((1.0e0 / v)) * 2.0e0))))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(cosTheta_i / v) * Float32(cosTheta_O_m / Float32(v * Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0)))))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * ((cosTheta_i / v) * (cosTheta_O_m / (v * (sinh((single(1.0) / v)) * single(2.0)))));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(\frac{cosTheta\_i}{v} \cdot \frac{cosTheta\_O\_m}{v \cdot \left(\sinh \left(\frac{1}{v}\right) \cdot 2\right)}\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
neg-mul-1N/A
lift-exp.f32N/A
neg-mul-1N/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
lift-exp.f32N/A
lift--.f32N/A
associate-*r*N/A
times-fracN/A
un-div-invN/A
lift-/.f32N/A
*-commutativeN/A
lift-*.f32N/A
Applied egg-rr98.6%
Final simplification98.6%
cosTheta_O\_m = (fabs.f32 cosTheta_O)
cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
(FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v)
:precision binary32
(*
cosTheta_O_s
(if (<= v 0.33500000834465027)
(/ (* cosTheta_O_m cosTheta_i) (* (* v v) (+ (exp (/ 1.0 v)) -1.0)))
(/
(* cosTheta_O_m cosTheta_i)
(*
v
(-
(/ (+ 0.3333333333333333 (/ 0.016666666666666666 (* v v))) (* v v))
-2.0))))))cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (v <= 0.33500000834465027f) {
tmp = (cosTheta_O_m * cosTheta_i) / ((v * v) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = (cosTheta_O_m * cosTheta_i) / (v * (((0.3333333333333333f + (0.016666666666666666f / (v * v))) / (v * v)) - -2.0f));
}
return cosTheta_O_s * tmp;
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: tmp
if (v <= 0.33500000834465027e0) then
tmp = (costheta_o_m * costheta_i) / ((v * v) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = (costheta_o_m * costheta_i) / (v * (((0.3333333333333333e0 + (0.016666666666666666e0 / (v * v))) / (v * v)) - (-2.0e0)))
end if
code = costheta_o_s * tmp
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (v <= Float32(0.33500000834465027)) tmp = Float32(Float32(cosTheta_O_m * cosTheta_i) / Float32(Float32(v * v) * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0)))); else tmp = Float32(Float32(cosTheta_O_m * cosTheta_i) / Float32(v * Float32(Float32(Float32(Float32(0.3333333333333333) + Float32(Float32(0.016666666666666666) / Float32(v * v))) / Float32(v * v)) - Float32(-2.0)))); end return Float32(cosTheta_O_s * tmp) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp_2 = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = single(0.0);
if (v <= single(0.33500000834465027))
tmp = (cosTheta_O_m * cosTheta_i) / ((v * v) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = (cosTheta_O_m * cosTheta_i) / (v * (((single(0.3333333333333333) + (single(0.016666666666666666) / (v * v))) / (v * v)) - single(-2.0)));
end
tmp_2 = cosTheta_O_s * tmp;
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \begin{array}{l}
\mathbf{if}\;v \leq 0.33500000834465027:\\
\;\;\;\;\frac{cosTheta\_O\_m \cdot cosTheta\_i}{\left(v \cdot v\right) \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta\_O\_m \cdot cosTheta\_i}{v \cdot \left(\frac{0.3333333333333333 + \frac{0.016666666666666666}{v \cdot v}}{v \cdot v} - -2\right)}\\
\end{array}
\end{array}
if v < 0.335000008Initial program 98.1%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.2%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3297.9
Simplified97.9%
Taylor expanded in v around inf
Simplified75.1%
if 0.335000008 < v Initial program 99.0%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr99.0%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.9
Simplified98.9%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
metadata-evalN/A
lower-/.f32N/A
Simplified79.3%
Final simplification77.4%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* cosTheta_i (/ cosTheta_O_m (* v (* v (* (sinh (/ 1.0 v)) 2.0)))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (cosTheta_i * (cosTheta_O_m / (v * (v * (sinhf((1.0f / v)) * 2.0f)))));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (costheta_i * (costheta_o_m / (v * (v * (sinh((1.0e0 / v)) * 2.0e0)))))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(cosTheta_i * Float32(cosTheta_O_m / Float32(v * Float32(v * Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0))))))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (cosTheta_i * (cosTheta_O_m / (v * (v * (sinh((single(1.0) / v)) * single(2.0))))));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(cosTheta\_i \cdot \frac{cosTheta\_O\_m}{v \cdot \left(v \cdot \left(\sinh \left(\frac{1}{v}\right) \cdot 2\right)\right)}\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
neg-mul-1N/A
lift-exp.f32N/A
neg-mul-1N/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
lift-exp.f32N/A
lift--.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied egg-rr98.5%
Final simplification98.5%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* cosTheta_O_m (/ cosTheta_i (* v (* v (* (sinh (/ 1.0 v)) 2.0)))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (cosTheta_O_m * (cosTheta_i / (v * (v * (sinhf((1.0f / v)) * 2.0f)))));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (costheta_o_m * (costheta_i / (v * (v * (sinh((1.0e0 / v)) * 2.0e0)))))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(cosTheta_O_m * Float32(cosTheta_i / Float32(v * Float32(v * Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0))))))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (cosTheta_O_m * (cosTheta_i / (v * (v * (sinh((single(1.0) / v)) * single(2.0))))));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(cosTheta\_O\_m \cdot \frac{cosTheta\_i}{v \cdot \left(v \cdot \left(\sinh \left(\frac{1}{v}\right) \cdot 2\right)\right)}\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
*-commutativeN/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
neg-mul-1N/A
lift-exp.f32N/A
neg-mul-1N/A
lift-/.f32N/A
div-invN/A
lift-/.f32N/A
lift-exp.f32N/A
lift--.f32N/A
lift-*.f32N/A
lift-*.f32N/A
associate-/l*N/A
lower-*.f32N/A
Applied egg-rr98.5%
Final simplification98.5%
cosTheta_O\_m = (fabs.f32 cosTheta_O)
cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
(FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v)
:precision binary32
(*
cosTheta_O_s
(/
(* cosTheta_O_m cosTheta_i)
(*
v
(-
(/ (+ 0.3333333333333333 (/ 0.016666666666666666 (* v v))) (* v v))
-2.0)))))cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * ((cosTheta_O_m * cosTheta_i) / (v * (((0.3333333333333333f + (0.016666666666666666f / (v * v))) / (v * v)) - -2.0f)));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * ((costheta_o_m * costheta_i) / (v * (((0.3333333333333333e0 + (0.016666666666666666e0 / (v * v))) / (v * v)) - (-2.0e0))))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(cosTheta_O_m * cosTheta_i) / Float32(v * Float32(Float32(Float32(Float32(0.3333333333333333) + Float32(Float32(0.016666666666666666) / Float32(v * v))) / Float32(v * v)) - Float32(-2.0))))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * ((cosTheta_O_m * cosTheta_i) / (v * (((single(0.3333333333333333) + (single(0.016666666666666666) / (v * v))) / (v * v)) - single(-2.0))));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{cosTheta\_O\_m \cdot cosTheta\_i}{v \cdot \left(\frac{0.3333333333333333 + \frac{0.016666666666666666}{v \cdot v}}{v \cdot v} - -2\right)}
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f32N/A
mul-1-negN/A
lower-neg.f32N/A
sub-negN/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
lower--.f32N/A
metadata-evalN/A
lower-/.f32N/A
Simplified70.8%
Final simplification70.8%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ (* cosTheta_i (/ cosTheta_O_m (+ 2.0 (/ 0.3333333333333333 (* v v))))) v)))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * ((cosTheta_i * (cosTheta_O_m / (2.0f + (0.3333333333333333f / (v * v))))) / v);
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * ((costheta_i * (costheta_o_m / (2.0e0 + (0.3333333333333333e0 / (v * v))))) / v)
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(cosTheta_i * Float32(cosTheta_O_m / Float32(Float32(2.0) + Float32(Float32(0.3333333333333333) / Float32(v * v))))) / v)) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * ((cosTheta_i * (cosTheta_O_m / (single(2.0) + (single(0.3333333333333333) / (v * v))))) / v);
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{cosTheta\_i \cdot \frac{cosTheta\_O\_m}{2 + \frac{0.3333333333333333}{v \cdot v}}}{v}
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
Taylor expanded in v around inf
lower-*.f32N/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
unpow2N/A
lower-*.f3265.0
Simplified65.0%
lift-*.f32N/A
lift-/.f32N/A
lift-+.f32N/A
times-fracN/A
associate-*l/N/A
lower-/.f32N/A
lower-*.f32N/A
lower-/.f3265.0
Applied egg-rr65.0%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* (/ cosTheta_i v) (/ cosTheta_O_m (+ 2.0 (/ 0.3333333333333333 (* v v)))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * ((cosTheta_i / v) * (cosTheta_O_m / (2.0f + (0.3333333333333333f / (v * v)))));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * ((costheta_i / v) * (costheta_o_m / (2.0e0 + (0.3333333333333333e0 / (v * v)))))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(cosTheta_i / v) * Float32(cosTheta_O_m / Float32(Float32(2.0) + Float32(Float32(0.3333333333333333) / Float32(v * v)))))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * ((cosTheta_i / v) * (cosTheta_O_m / (single(2.0) + (single(0.3333333333333333) / (v * v)))));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(\frac{cosTheta\_i}{v} \cdot \frac{cosTheta\_O\_m}{2 + \frac{0.3333333333333333}{v \cdot v}}\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
Taylor expanded in v around inf
lower-*.f32N/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
unpow2N/A
lower-*.f3265.0
Simplified65.0%
lift-*.f32N/A
lift-/.f32N/A
lift-+.f32N/A
times-fracN/A
lower-*.f32N/A
lower-/.f32N/A
lower-/.f3265.0
Applied egg-rr65.0%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* cosTheta_i (/ cosTheta_O_m (fma v 2.0 (/ 0.3333333333333333 v))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (cosTheta_i * (cosTheta_O_m / fmaf(v, 2.0f, (0.3333333333333333f / v))));
}
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(cosTheta_i * Float32(cosTheta_O_m / fma(v, Float32(2.0), Float32(Float32(0.3333333333333333) / v))))) end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(cosTheta\_i \cdot \frac{cosTheta\_O\_m}{\mathsf{fma}\left(v, 2, \frac{0.3333333333333333}{v}\right)}\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
Taylor expanded in v around inf
lower-*.f32N/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
unpow2N/A
lower-*.f3265.0
Simplified65.0%
lift-*.f32N/A
lift-/.f32N/A
lift-+.f32N/A
lift-*.f32N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f32N/A
Applied egg-rr65.0%
Final simplification65.0%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* cosTheta_O_m (/ cosTheta_i (fma v 2.0 (/ 0.3333333333333333 v))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (cosTheta_O_m * (cosTheta_i / fmaf(v, 2.0f, (0.3333333333333333f / v))));
}
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(cosTheta_O_m * Float32(cosTheta_i / fma(v, Float32(2.0), Float32(Float32(0.3333333333333333) / v))))) end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(cosTheta\_O\_m \cdot \frac{cosTheta\_i}{\mathsf{fma}\left(v, 2, \frac{0.3333333333333333}{v}\right)}\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
Taylor expanded in v around inf
lower-*.f32N/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
unpow2N/A
lower-*.f3265.0
Simplified65.0%
lift-*.f32N/A
lift-/.f32N/A
lift-+.f32N/A
*-commutativeN/A
lift-*.f32N/A
associate-/l*N/A
lower-*.f32N/A
lower-/.f3265.0
lift-*.f32N/A
lift-+.f32N/A
distribute-rgt-inN/A
*-commutativeN/A
lift-/.f32N/A
div-invN/A
associate-*l*N/A
lift-*.f32N/A
pow2N/A
pow-flipN/A
pow-plusN/A
metadata-evalN/A
metadata-evalN/A
inv-powN/A
Applied egg-rr65.0%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ 1.0 (/ (* v 2.0) (* cosTheta_O_m cosTheta_i)))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (1.0f / ((v * 2.0f) / (cosTheta_O_m * cosTheta_i)));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (1.0e0 / ((v * 2.0e0) / (costheta_o_m * costheta_i)))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(1.0) / Float32(Float32(v * Float32(2.0)) / Float32(cosTheta_O_m * cosTheta_i)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (single(1.0) / ((v * single(2.0)) / (cosTheta_O_m * cosTheta_i)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{1}{\frac{v \cdot 2}{cosTheta\_O\_m \cdot cosTheta\_i}}
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
associate-/l*N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f3259.5
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.5
Applied egg-rr59.5%
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
un-div-invN/A
clear-numN/A
lower-/.f32N/A
lower-/.f3259.8
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.8
Applied egg-rr59.8%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ 1.0 (/ v (* 0.5 (* cosTheta_O_m cosTheta_i))))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (1.0f / (v / (0.5f * (cosTheta_O_m * cosTheta_i))));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (1.0e0 / (v / (0.5e0 * (costheta_o_m * costheta_i))))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(1.0) / Float32(v / Float32(Float32(0.5) * Float32(cosTheta_O_m * cosTheta_i))))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (single(1.0) / (v / (single(0.5) * (cosTheta_O_m * cosTheta_i))));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{1}{\frac{v}{0.5 \cdot \left(cosTheta\_O\_m \cdot cosTheta\_i\right)}}
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
lift-*.f32N/A
clear-numN/A
lower-/.f32N/A
lower-/.f3259.8
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-*.f3259.8
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.8
Applied egg-rr59.8%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ 0.5 (/ v (* cosTheta_O_m cosTheta_i)))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (0.5f / (v / (cosTheta_O_m * cosTheta_i)));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (0.5e0 / (v / (costheta_o_m * costheta_i)))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(0.5) / Float32(v / Float32(cosTheta_O_m * cosTheta_i)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (single(0.5) / (v / (cosTheta_O_m * cosTheta_i)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{0.5}{\frac{v}{cosTheta\_O\_m \cdot cosTheta\_i}}
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f32N/A
lower-/.f3259.8
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.8
Applied egg-rr59.8%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (/ (* 0.5 (* cosTheta_O_m cosTheta_i)) v)))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * ((0.5f * (cosTheta_O_m * cosTheta_i)) / v);
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * ((0.5e0 * (costheta_o_m * costheta_i)) / v)
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(Float32(0.5) * Float32(cosTheta_O_m * cosTheta_i)) / v)) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * ((single(0.5) * (cosTheta_O_m * cosTheta_i)) / v);
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \frac{0.5 \cdot \left(cosTheta\_O\_m \cdot cosTheta\_i\right)}{v}
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
Final simplification59.5%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* 0.5 (/ (* cosTheta_O_m cosTheta_i) v))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (0.5f * ((cosTheta_O_m * cosTheta_i) / v));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (0.5e0 * ((costheta_o_m * costheta_i) / v))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(0.5) * Float32(Float32(cosTheta_O_m * cosTheta_i) / v))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (single(0.5) * ((cosTheta_O_m * cosTheta_i) / v));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(0.5 \cdot \frac{cosTheta\_O\_m \cdot cosTheta\_i}{v}\right)
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
associate-/l*N/A
lift-*.f32N/A
*-commutativeN/A
associate-*r/N/A
lift-/.f32N/A
lift-*.f32N/A
*-commutativeN/A
lower-*.f3259.5
lift-*.f32N/A
lift-/.f32N/A
associate-*r/N/A
*-commutativeN/A
lift-*.f32N/A
lower-/.f3259.5
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.5
Applied egg-rr59.5%
Final simplification59.5%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* (/ 0.5 v) (* cosTheta_O_m cosTheta_i))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * ((0.5f / v) * (cosTheta_O_m * cosTheta_i));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * ((0.5e0 / v) * (costheta_o_m * costheta_i))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(Float32(Float32(0.5) / v) * Float32(cosTheta_O_m * cosTheta_i))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * ((single(0.5) / v) * (cosTheta_O_m * cosTheta_i));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(\frac{0.5}{v} \cdot \left(cosTheta\_O\_m \cdot cosTheta\_i\right)\right)
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
associate-/l*N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f3259.5
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.5
Applied egg-rr59.5%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* cosTheta_O_m (* cosTheta_i (/ 0.5 v)))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (cosTheta_O_m * (cosTheta_i * (0.5f / v)));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (costheta_o_m * (costheta_i * (0.5e0 / v)))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(cosTheta_O_m * Float32(cosTheta_i * Float32(Float32(0.5) / v)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (cosTheta_O_m * (cosTheta_i * (single(0.5) / v)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(cosTheta\_O\_m \cdot \left(cosTheta\_i \cdot \frac{0.5}{v}\right)\right)
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
associate-/l*N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f3259.5
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.5
Applied egg-rr59.5%
lift-/.f32N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
lift-/.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
lower-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
lift-/.f3259.4
Applied egg-rr59.4%
Final simplification59.4%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* cosTheta_i (* cosTheta_O_m (/ 0.5 v)))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (cosTheta_i * (cosTheta_O_m * (0.5f / v)));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (costheta_i * (costheta_o_m * (0.5e0 / v)))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(cosTheta_i * Float32(cosTheta_O_m * Float32(Float32(0.5) / v)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (cosTheta_i * (cosTheta_O_m * (single(0.5) / v)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(cosTheta\_i \cdot \left(cosTheta\_O\_m \cdot \frac{0.5}{v}\right)\right)
\end{array}
Initial program 98.6%
Taylor expanded in v around inf
associate-*r/N/A
lower-/.f32N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f3259.5
Simplified59.5%
lift-*.f32N/A
associate-/l*N/A
lift-*.f32N/A
*-commutativeN/A
lift-*.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
*-commutativeN/A
lower-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
lower-/.f3259.5
lift-*.f32N/A
*-commutativeN/A
lift-*.f3259.5
Applied egg-rr59.5%
lift-/.f32N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
lift-/.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
lift-*.f32N/A
lower-*.f32N/A
lift-*.f32N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
lift-/.f3259.4
Applied egg-rr59.4%
lift-/.f32N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f32N/A
lower-*.f3259.4
Applied egg-rr59.4%
Final simplification59.4%
cosTheta_O\_m = (fabs.f32 cosTheta_O) cosTheta_O\_s = (copysign.f32 #s(literal 1 binary32) cosTheta_O) NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function. (FPCore (cosTheta_O_s cosTheta_i cosTheta_O_m sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O_s (* v (* (* cosTheta_O_m cosTheta_i) 3.0))))
cosTheta_O\_m = fabs(cosTheta_O);
cosTheta_O\_s = copysign(1.0, cosTheta_O);
assert(cosTheta_i < cosTheta_O_m && cosTheta_O_m < sinTheta_i && sinTheta_i < sinTheta_O && sinTheta_O < v);
float code(float cosTheta_O_s, float cosTheta_i, float cosTheta_O_m, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O_s * (v * ((cosTheta_O_m * cosTheta_i) * 3.0f));
}
cosTheta_O\_m = abs(costheta_o)
cosTheta_O\_s = copysign(1.0d0, costheta_o)
NOTE: cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, and v should be sorted in increasing order before calling this function.
real(4) function code(costheta_o_s, costheta_i, costheta_o_m, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_o_s
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o_m
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_o_s * (v * ((costheta_o_m * costheta_i) * 3.0e0))
end function
cosTheta_O\_m = abs(cosTheta_O) cosTheta_O\_s = copysign(1.0, cosTheta_O) cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v]) function code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O_s * Float32(v * Float32(Float32(cosTheta_O_m * cosTheta_i) * Float32(3.0)))) end
cosTheta_O\_m = abs(cosTheta_O);
cosTheta_O\_s = sign(cosTheta_O) * abs(1.0);
cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v = num2cell(sort([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])){:}
function tmp = code(cosTheta_O_s, cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_O_s * (v * ((cosTheta_O_m * cosTheta_i) * single(3.0)));
end
\begin{array}{l}
cosTheta_O\_m = \left|cosTheta\_O\right|
\\
cosTheta_O\_s = \mathsf{copysign}\left(1, cosTheta\_O\right)
\\
[cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v] = \mathsf{sort}([cosTheta_i, cosTheta_O_m, sinTheta_i, sinTheta_O, v])\\
\\
cosTheta\_O\_s \cdot \left(v \cdot \left(\left(cosTheta\_O\_m \cdot cosTheta\_i\right) \cdot 3\right)\right)
\end{array}
Initial program 98.6%
lift-*.f32N/A
lift-/.f32N/A
lift-neg.f32N/A
lift-exp.f32N/A
lift-*.f32N/A
lift-/.f32N/A
lift-/.f32N/A
lift-sinh.f32N/A
lift-*.f32N/A
times-fracN/A
Applied egg-rr98.6%
Taylor expanded in sinTheta_i around 0
lower--.f32N/A
lower-exp.f32N/A
lower-/.f32N/A
rec-expN/A
distribute-neg-fracN/A
metadata-evalN/A
lower-exp.f32N/A
lower-/.f3298.4
Simplified98.4%
Taylor expanded in v around inf
lower-*.f32N/A
lower-+.f32N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f32N/A
unpow2N/A
lower-*.f3265.0
Simplified65.0%
Taylor expanded in v around 0
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-*.f32N/A
lower-*.f32N/A
lower-*.f3257.9
Simplified57.9%
Final simplification57.9%
herbie shell --seed 2024207
(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)))