| Alternative 1 | |
|---|---|
| Error | 0.5 |
| Cost | 10272 |
\[\frac{1}{\left(1 + c\right) + e^{cosTheta \cdot \left(-cosTheta\right)} \cdot \frac{1}{cosTheta \cdot \sqrt{\frac{\pi}{1 - \left(cosTheta + cosTheta\right)}}}}
\]
(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(*
(* (/ 1.0 (sqrt PI)) (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta))
(exp (* (- cosTheta) cosTheta))))))(FPCore (cosTheta c)
:precision binary32
(/
1.0
(+
(+ 1.0 c)
(*
(*
(/ (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta) (cbrt PI))
(/ 1.0 (cbrt (sqrt PI))))
(exp (* cosTheta (- cosTheta)))))))float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + (((1.0f / sqrtf(((float) M_PI))) * (sqrtf(((1.0f - cosTheta) - cosTheta)) / cosTheta)) * expf((-cosTheta * cosTheta))));
}
float code(float cosTheta, float c) {
return 1.0f / ((1.0f + c) + ((((sqrtf(((1.0f - cosTheta) - cosTheta)) / cosTheta) / cbrtf(((float) M_PI))) * (1.0f / cbrtf(sqrtf(((float) M_PI))))) * expf((cosTheta * -cosTheta))));
}
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(Float32(1.0) / sqrt(Float32(pi))) * Float32(sqrt(Float32(Float32(Float32(1.0) - cosTheta) - cosTheta)) / cosTheta)) * exp(Float32(Float32(-cosTheta) * cosTheta))))) end
function code(cosTheta, c) return Float32(Float32(1.0) / Float32(Float32(Float32(1.0) + c) + Float32(Float32(Float32(Float32(sqrt(Float32(Float32(Float32(1.0) - cosTheta) - cosTheta)) / cosTheta) / cbrt(Float32(pi))) * Float32(Float32(1.0) / cbrt(sqrt(Float32(pi))))) * exp(Float32(cosTheta * Float32(-cosTheta)))))) end
\frac{1}{\left(1 + c\right) + \left(\frac{1}{\sqrt{\pi}} \cdot \frac{\sqrt{\left(1 - cosTheta\right) - cosTheta}}{cosTheta}\right) \cdot e^{\left(-cosTheta\right) \cdot cosTheta}}
\frac{1}{\left(1 + c\right) + \left(\frac{\frac{\sqrt{\left(1 - cosTheta\right) - cosTheta}}{cosTheta}}{\sqrt[3]{\pi}} \cdot \frac{1}{\sqrt[3]{\sqrt{\pi}}}\right) \cdot e^{cosTheta \cdot \left(-cosTheta\right)}}
Results
Initial program 0.7
Applied egg-rr0.5
Applied egg-rr0.5
Final simplification0.5
| Alternative 1 | |
|---|---|
| Error | 0.5 |
| Cost | 10272 |
| Alternative 2 | |
|---|---|
| Error | 0.7 |
| Cost | 10240 |
| Alternative 3 | |
|---|---|
| Error | 1.1 |
| Cost | 6976 |
| Alternative 4 | |
|---|---|
| Error | 1.5 |
| Cost | 6848 |
| Alternative 5 | |
|---|---|
| Error | 2.2 |
| Cost | 6464 |
| Alternative 6 | |
|---|---|
| Error | 28.5 |
| Cost | 3424 |
| Alternative 7 | |
|---|---|
| Error | 28.6 |
| Cost | 32 |
herbie shell --seed 2023187
(FPCore (cosTheta c)
:name "Beckmann Sample, normalization factor"
:precision binary32
:pre (and (and (< 0.0 cosTheta) (< cosTheta 0.9999)) (and (< -1.0 c) (< c 1.0)))
(/ 1.0 (+ (+ 1.0 c) (* (* (/ 1.0 (sqrt PI)) (/ (sqrt (- (- 1.0 cosTheta) cosTheta)) cosTheta)) (exp (* (- cosTheta) cosTheta))))))