
(FPCore (v) :precision binary64 (/ 4 (* (* (* 3 PI) (- 1 (* v v))) (sqrt (- 2 (* 6 (* v v)))))))
double code(double v) {
return 4.0 / (((3.0 * ((double) M_PI)) * (1.0 - (v * v))) * sqrt((2.0 - (6.0 * (v * v)))));
}
public static double code(double v) {
return 4.0 / (((3.0 * Math.PI) * (1.0 - (v * v))) * Math.sqrt((2.0 - (6.0 * (v * v)))));
}
def code(v): return 4.0 / (((3.0 * math.pi) * (1.0 - (v * v))) * math.sqrt((2.0 - (6.0 * (v * v)))))
function code(v) return Float64(4.0 / Float64(Float64(Float64(3.0 * pi) * Float64(1.0 - Float64(v * v))) * sqrt(Float64(2.0 - Float64(6.0 * Float64(v * v)))))) end
function tmp = code(v) tmp = 4.0 / (((3.0 * pi) * (1.0 - (v * v))) * sqrt((2.0 - (6.0 * (v * v))))); end
code[v_] := N[(4 / N[(N[(N[(3 * Pi), $MachinePrecision] * N[(1 - N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(2 - N[(6 * N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (v) :precision binary64 (/ 4 (* (* (* 3 PI) (- 1 (* v v))) (sqrt (- 2 (* 6 (* v v)))))))
double code(double v) {
return 4.0 / (((3.0 * ((double) M_PI)) * (1.0 - (v * v))) * sqrt((2.0 - (6.0 * (v * v)))));
}
public static double code(double v) {
return 4.0 / (((3.0 * Math.PI) * (1.0 - (v * v))) * Math.sqrt((2.0 - (6.0 * (v * v)))));
}
def code(v): return 4.0 / (((3.0 * math.pi) * (1.0 - (v * v))) * math.sqrt((2.0 - (6.0 * (v * v)))))
function code(v) return Float64(4.0 / Float64(Float64(Float64(3.0 * pi) * Float64(1.0 - Float64(v * v))) * sqrt(Float64(2.0 - Float64(6.0 * Float64(v * v)))))) end
function tmp = code(v) tmp = 4.0 / (((3.0 * pi) * (1.0 - (v * v))) * sqrt((2.0 - (6.0 * (v * v))))); end
code[v_] := N[(4 / N[(N[(N[(3 * Pi), $MachinePrecision] * N[(1 - N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(2 - N[(6 * N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{4}{\left(\left(3 \cdot \pi\right) \cdot \left(1 - v \cdot v\right)\right) \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}
(FPCore (v) :precision binary64 (/ 4/3 (* (* (- 1 (* v v)) PI) (sqrt (- 2 (* 6 (* v v)))))))
double code(double v) {
return 1.3333333333333333 / (((1.0 - (v * v)) * ((double) M_PI)) * sqrt((2.0 - (6.0 * (v * v)))));
}
public static double code(double v) {
return 1.3333333333333333 / (((1.0 - (v * v)) * Math.PI) * Math.sqrt((2.0 - (6.0 * (v * v)))));
}
def code(v): return 1.3333333333333333 / (((1.0 - (v * v)) * math.pi) * math.sqrt((2.0 - (6.0 * (v * v)))))
function code(v) return Float64(1.3333333333333333 / Float64(Float64(Float64(1.0 - Float64(v * v)) * pi) * sqrt(Float64(2.0 - Float64(6.0 * Float64(v * v)))))) end
function tmp = code(v) tmp = 1.3333333333333333 / (((1.0 - (v * v)) * pi) * sqrt((2.0 - (6.0 * (v * v))))); end
code[v_] := N[(4/3 / N[(N[(N[(1 - N[(v * v), $MachinePrecision]), $MachinePrecision] * Pi), $MachinePrecision] * N[Sqrt[N[(2 - N[(6 * N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{\frac{4}{3}}{\left(\left(1 - v \cdot v\right) \cdot \pi\right) \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}
Initial program 98.5%
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
associate-/r*N/A
lower-/.f64N/A
metadata-evalN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64100.0%
Applied rewrites100.0%
(FPCore (v) :precision binary64 (/ 4/3 (* PI (sqrt (- (- 1 (* (* v v) 6)) -1)))))
double code(double v) {
return 1.3333333333333333 / (((double) M_PI) * sqrt(((1.0 - ((v * v) * 6.0)) - -1.0)));
}
public static double code(double v) {
return 1.3333333333333333 / (Math.PI * Math.sqrt(((1.0 - ((v * v) * 6.0)) - -1.0)));
}
def code(v): return 1.3333333333333333 / (math.pi * math.sqrt(((1.0 - ((v * v) * 6.0)) - -1.0)))
function code(v) return Float64(1.3333333333333333 / Float64(pi * sqrt(Float64(Float64(1.0 - Float64(Float64(v * v) * 6.0)) - -1.0)))) end
function tmp = code(v) tmp = 1.3333333333333333 / (pi * sqrt(((1.0 - ((v * v) * 6.0)) - -1.0))); end
code[v_] := N[(4/3 / N[(Pi * N[Sqrt[N[(N[(1 - N[(N[(v * v), $MachinePrecision] * 6), $MachinePrecision]), $MachinePrecision] - -1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{\frac{4}{3}}{\pi \cdot \sqrt{\left(1 - \left(v \cdot v\right) \cdot 6\right) - -1}}
Initial program 98.5%
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
associate-/r*N/A
lower-/.f64N/A
metadata-evalN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64100.0%
Applied rewrites100.0%
Taylor expanded in v around 0
lower-PI.f6499.0%
Applied rewrites99.0%
lift--.f64N/A
metadata-evalN/A
associate--l-N/A
+-commutativeN/A
associate--l-N/A
lift--.f64N/A
lift--.f6499.0%
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.0%
Applied rewrites99.0%
(FPCore (v) :precision binary64 (/ 4/3 (* PI (sqrt (- 2 (* 6 (* v v)))))))
double code(double v) {
return 1.3333333333333333 / (((double) M_PI) * sqrt((2.0 - (6.0 * (v * v)))));
}
public static double code(double v) {
return 1.3333333333333333 / (Math.PI * Math.sqrt((2.0 - (6.0 * (v * v)))));
}
def code(v): return 1.3333333333333333 / (math.pi * math.sqrt((2.0 - (6.0 * (v * v)))))
function code(v) return Float64(1.3333333333333333 / Float64(pi * sqrt(Float64(2.0 - Float64(6.0 * Float64(v * v)))))) end
function tmp = code(v) tmp = 1.3333333333333333 / (pi * sqrt((2.0 - (6.0 * (v * v))))); end
code[v_] := N[(4/3 / N[(Pi * N[Sqrt[N[(2 - N[(6 * N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{\frac{4}{3}}{\pi \cdot \sqrt{2 - 6 \cdot \left(v \cdot v\right)}}
Initial program 98.5%
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
associate-*l*N/A
associate-/r*N/A
lower-/.f64N/A
metadata-evalN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64100.0%
Applied rewrites100.0%
Taylor expanded in v around 0
lower-PI.f6499.0%
Applied rewrites99.0%
(FPCore (v) :precision binary64 (/ 4/3 (* PI (sqrt 2))))
double code(double v) {
return 1.3333333333333333 / (((double) M_PI) * sqrt(2.0));
}
public static double code(double v) {
return 1.3333333333333333 / (Math.PI * Math.sqrt(2.0));
}
def code(v): return 1.3333333333333333 / (math.pi * math.sqrt(2.0))
function code(v) return Float64(1.3333333333333333 / Float64(pi * sqrt(2.0))) end
function tmp = code(v) tmp = 1.3333333333333333 / (pi * sqrt(2.0)); end
code[v_] := N[(4/3 / N[(Pi * N[Sqrt[2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{\frac{4}{3}}{\pi \cdot \sqrt{2}}
Initial program 98.5%
Taylor expanded in v around 0
lower-/.f64N/A
lower-*.f64N/A
lower-PI.f64N/A
lower-sqrt.f6498.9%
Applied rewrites98.9%
herbie shell --seed 2025271 -o generate:evaluate
(FPCore (v)
:name "Falkner and Boettcher, Equation (22+)"
:precision binary64
(/ 4 (* (* (* 3 PI) (- 1 (* v v))) (sqrt (- 2 (* 6 (* v v)))))))