Alternative 1 | |
---|---|
Error | 0.3 |
Cost | 576 |
\[\frac{deltaP}{1 + \left(deltaP \cdot deltaP\right) \cdot -0.5}
\]
(FPCore (deltaP) :precision binary64 (/ deltaP (sqrt (- 1.0 (* deltaP deltaP)))))
(FPCore (deltaP) :precision binary64 (/ deltaP (sqrt (- 1.0 (* deltaP deltaP)))))
double code(double deltaP) { return deltaP / sqrt((1.0 - (deltaP * deltaP))); }
double code(double deltaP) { return deltaP / sqrt((1.0 - (deltaP * deltaP))); }
real(8) function code(deltap) real(8), intent (in) :: deltap code = deltap / sqrt((1.0d0 - (deltap * deltap))) end function
real(8) function code(deltap) real(8), intent (in) :: deltap code = deltap / sqrt((1.0d0 - (deltap * deltap))) end function
public static double code(double deltaP) { return deltaP / Math.sqrt((1.0 - (deltaP * deltaP))); }
public static double code(double deltaP) { return deltaP / Math.sqrt((1.0 - (deltaP * deltaP))); }
def code(deltaP): return deltaP / math.sqrt((1.0 - (deltaP * deltaP)))
def code(deltaP): return deltaP / math.sqrt((1.0 - (deltaP * deltaP)))
function code(deltaP) return Float64(deltaP / sqrt(Float64(1.0 - Float64(deltaP * deltaP)))) end
function code(deltaP) return Float64(deltaP / sqrt(Float64(1.0 - Float64(deltaP * deltaP)))) end
function tmp = code(deltaP) tmp = deltaP / sqrt((1.0 - (deltaP * deltaP))); end
function tmp = code(deltaP) tmp = deltaP / sqrt((1.0 - (deltaP * deltaP))); end
code[deltaP_] := N[(deltaP / N[Sqrt[N[(1.0 - N[(deltaP * deltaP), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
code[deltaP_] := N[(deltaP / N[Sqrt[N[(1.0 - N[(deltaP * deltaP), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\frac{deltaP}{\sqrt{1 - deltaP \cdot deltaP}}
\frac{deltaP}{\sqrt{1 - deltaP \cdot deltaP}}
Results
Initial program 0.0
Final simplification0.0
Alternative 1 | |
---|---|
Error | 0.3 |
Cost | 576 |
Alternative 2 | |
---|---|
Error | 0.7 |
Cost | 64 |
herbie shell --seed 1
(FPCore (deltaP)
:name "deltaP / sqrt(1. - deltaP * deltaP)"
:precision binary64
:pre (and (<= -1000000000.0 deltaP) (<= deltaP 1000000000.0))
(/ deltaP (sqrt (- 1.0 (* deltaP deltaP)))))