
(FPCore (d1) :precision binary64 (* (* d1 (* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1)) d1))
double code(double d1) {
return (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1;
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1
end function
public static double code(double d1) {
return (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1;
}
def code(d1): return (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1
function code(d1) return Float64(Float64(d1 * Float64(Float64(Float64(Float64(Float64(d1 * Float64(d1 * d1)) * d1) * d1) * Float64(d1 * d1)) * d1)) * d1) end
function tmp = code(d1) tmp = (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1; end
code[d1_] := N[(N[(d1 * N[(N[(N[(N[(N[(d1 * N[(d1 * d1), $MachinePrecision]), $MachinePrecision] * d1), $MachinePrecision] * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision] * d1), $MachinePrecision]), $MachinePrecision] * d1), $MachinePrecision]
\begin{array}{l}
\\
\left(d1 \cdot \left(\left(\left(\left(\left(d1 \cdot \left(d1 \cdot d1\right)\right) \cdot d1\right) \cdot d1\right) \cdot \left(d1 \cdot d1\right)\right) \cdot d1\right)\right) \cdot d1
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (d1) :precision binary64 (* (* d1 (* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1)) d1))
double code(double d1) {
return (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1;
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1
end function
public static double code(double d1) {
return (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1;
}
def code(d1): return (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1
function code(d1) return Float64(Float64(d1 * Float64(Float64(Float64(Float64(Float64(d1 * Float64(d1 * d1)) * d1) * d1) * Float64(d1 * d1)) * d1)) * d1) end
function tmp = code(d1) tmp = (d1 * (((((d1 * (d1 * d1)) * d1) * d1) * (d1 * d1)) * d1)) * d1; end
code[d1_] := N[(N[(d1 * N[(N[(N[(N[(N[(d1 * N[(d1 * d1), $MachinePrecision]), $MachinePrecision] * d1), $MachinePrecision] * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision] * d1), $MachinePrecision]), $MachinePrecision] * d1), $MachinePrecision]
\begin{array}{l}
\\
\left(d1 \cdot \left(\left(\left(\left(\left(d1 \cdot \left(d1 \cdot d1\right)\right) \cdot d1\right) \cdot d1\right) \cdot \left(d1 \cdot d1\right)\right) \cdot d1\right)\right) \cdot d1
\end{array}
(FPCore (d1) :precision binary64 (pow d1 10.0))
double code(double d1) {
return pow(d1, 10.0);
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 ** 10.0d0
end function
public static double code(double d1) {
return Math.pow(d1, 10.0);
}
def code(d1): return math.pow(d1, 10.0)
function code(d1) return d1 ^ 10.0 end
function tmp = code(d1) tmp = d1 ^ 10.0; end
code[d1_] := N[Power[d1, 10.0], $MachinePrecision]
\begin{array}{l}
\\
{d1}^{10}
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Taylor expanded in d1 around 0
pow-lowering-pow.f64100.0%
Simplified100.0%
(FPCore (d1) :precision binary64 (* d1 (* d1 (* d1 (* (* d1 d1) (* d1 (* d1 (* d1 (* d1 d1)))))))))
double code(double d1) {
return d1 * (d1 * (d1 * ((d1 * d1) * (d1 * (d1 * (d1 * (d1 * d1)))))));
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 * (d1 * (d1 * ((d1 * d1) * (d1 * (d1 * (d1 * (d1 * d1)))))))
end function
public static double code(double d1) {
return d1 * (d1 * (d1 * ((d1 * d1) * (d1 * (d1 * (d1 * (d1 * d1)))))));
}
def code(d1): return d1 * (d1 * (d1 * ((d1 * d1) * (d1 * (d1 * (d1 * (d1 * d1)))))))
function code(d1) return Float64(d1 * Float64(d1 * Float64(d1 * Float64(Float64(d1 * d1) * Float64(d1 * Float64(d1 * Float64(d1 * Float64(d1 * d1)))))))) end
function tmp = code(d1) tmp = d1 * (d1 * (d1 * ((d1 * d1) * (d1 * (d1 * (d1 * (d1 * d1))))))); end
code[d1_] := N[(d1 * N[(d1 * N[(d1 * N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * N[(d1 * N[(d1 * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d1 \cdot \left(d1 \cdot \left(d1 \cdot \left(\left(d1 \cdot d1\right) \cdot \left(d1 \cdot \left(d1 \cdot \left(d1 \cdot \left(d1 \cdot d1\right)\right)\right)\right)\right)\right)\right)
\end{array}
Initial program 99.9%
Final simplification99.9%
(FPCore (d1) :precision binary64 (let* ((t_0 (* (* d1 d1) (* d1 d1)))) (* (* d1 d1) (* t_0 t_0))))
double code(double d1) {
double t_0 = (d1 * d1) * (d1 * d1);
return (d1 * d1) * (t_0 * t_0);
}
real(8) function code(d1)
real(8), intent (in) :: d1
real(8) :: t_0
t_0 = (d1 * d1) * (d1 * d1)
code = (d1 * d1) * (t_0 * t_0)
end function
public static double code(double d1) {
double t_0 = (d1 * d1) * (d1 * d1);
return (d1 * d1) * (t_0 * t_0);
}
def code(d1): t_0 = (d1 * d1) * (d1 * d1) return (d1 * d1) * (t_0 * t_0)
function code(d1) t_0 = Float64(Float64(d1 * d1) * Float64(d1 * d1)) return Float64(Float64(d1 * d1) * Float64(t_0 * t_0)) end
function tmp = code(d1) t_0 = (d1 * d1) * (d1 * d1); tmp = (d1 * d1) * (t_0 * t_0); end
code[d1_] := Block[{t$95$0 = N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]}, N[(N[(d1 * d1), $MachinePrecision] * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(d1 \cdot d1\right) \cdot \left(d1 \cdot d1\right)\\
\left(d1 \cdot d1\right) \cdot \left(t\_0 \cdot t\_0\right)
\end{array}
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
(FPCore (d1) :precision binary64 (let* ((t_0 (* (* d1 d1) (* d1 d1)))) (* d1 (* t_0 t_0))))
double code(double d1) {
double t_0 = (d1 * d1) * (d1 * d1);
return d1 * (t_0 * t_0);
}
real(8) function code(d1)
real(8), intent (in) :: d1
real(8) :: t_0
t_0 = (d1 * d1) * (d1 * d1)
code = d1 * (t_0 * t_0)
end function
public static double code(double d1) {
double t_0 = (d1 * d1) * (d1 * d1);
return d1 * (t_0 * t_0);
}
def code(d1): t_0 = (d1 * d1) * (d1 * d1) return d1 * (t_0 * t_0)
function code(d1) t_0 = Float64(Float64(d1 * d1) * Float64(d1 * d1)) return Float64(d1 * Float64(t_0 * t_0)) end
function tmp = code(d1) t_0 = (d1 * d1) * (d1 * d1); tmp = d1 * (t_0 * t_0); end
code[d1_] := Block[{t$95$0 = N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]}, N[(d1 * N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(d1 \cdot d1\right) \cdot \left(d1 \cdot d1\right)\\
d1 \cdot \left(t\_0 \cdot t\_0\right)
\end{array}
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr71.1%
(FPCore (d1) :precision binary64 (* d1 (* (* (* d1 d1) (* d1 d1)) (* d1 (* d1 d1)))))
double code(double d1) {
return d1 * (((d1 * d1) * (d1 * d1)) * (d1 * (d1 * d1)));
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 * (((d1 * d1) * (d1 * d1)) * (d1 * (d1 * d1)))
end function
public static double code(double d1) {
return d1 * (((d1 * d1) * (d1 * d1)) * (d1 * (d1 * d1)));
}
def code(d1): return d1 * (((d1 * d1) * (d1 * d1)) * (d1 * (d1 * d1)))
function code(d1) return Float64(d1 * Float64(Float64(Float64(d1 * d1) * Float64(d1 * d1)) * Float64(d1 * Float64(d1 * d1)))) end
function tmp = code(d1) tmp = d1 * (((d1 * d1) * (d1 * d1)) * (d1 * (d1 * d1))); end
code[d1_] := N[(d1 * N[(N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision] * N[(d1 * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d1 \cdot \left(\left(\left(d1 \cdot d1\right) \cdot \left(d1 \cdot d1\right)\right) \cdot \left(d1 \cdot \left(d1 \cdot d1\right)\right)\right)
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr71.1%
Applied egg-rr89.0%
Final simplification89.0%
(FPCore (d1) :precision binary64 (* (* d1 d1) (* d1 (* (* d1 d1) (* d1 d1)))))
double code(double d1) {
return (d1 * d1) * (d1 * ((d1 * d1) * (d1 * d1)));
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = (d1 * d1) * (d1 * ((d1 * d1) * (d1 * d1)))
end function
public static double code(double d1) {
return (d1 * d1) * (d1 * ((d1 * d1) * (d1 * d1)));
}
def code(d1): return (d1 * d1) * (d1 * ((d1 * d1) * (d1 * d1)))
function code(d1) return Float64(Float64(d1 * d1) * Float64(d1 * Float64(Float64(d1 * d1) * Float64(d1 * d1)))) end
function tmp = code(d1) tmp = (d1 * d1) * (d1 * ((d1 * d1) * (d1 * d1))); end
code[d1_] := N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(d1 \cdot d1\right) \cdot \left(d1 \cdot \left(\left(d1 \cdot d1\right) \cdot \left(d1 \cdot d1\right)\right)\right)
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr68.1%
Final simplification68.1%
(FPCore (d1) :precision binary64 (* d1 (* d1 (* (* d1 d1) (* d1 d1)))))
double code(double d1) {
return d1 * (d1 * ((d1 * d1) * (d1 * d1)));
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 * (d1 * ((d1 * d1) * (d1 * d1)))
end function
public static double code(double d1) {
return d1 * (d1 * ((d1 * d1) * (d1 * d1)));
}
def code(d1): return d1 * (d1 * ((d1 * d1) * (d1 * d1)))
function code(d1) return Float64(d1 * Float64(d1 * Float64(Float64(d1 * d1) * Float64(d1 * d1)))) end
function tmp = code(d1) tmp = d1 * (d1 * ((d1 * d1) * (d1 * d1))); end
code[d1_] := N[(d1 * N[(d1 * N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d1 \cdot \left(d1 \cdot \left(\left(d1 \cdot d1\right) \cdot \left(d1 \cdot d1\right)\right)\right)
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr71.1%
Applied egg-rr86.9%
Final simplification86.9%
(FPCore (d1) :precision binary64 (* d1 (* d1 (* d1 (* d1 d1)))))
double code(double d1) {
return d1 * (d1 * (d1 * (d1 * d1)));
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 * (d1 * (d1 * (d1 * d1)))
end function
public static double code(double d1) {
return d1 * (d1 * (d1 * (d1 * d1)));
}
def code(d1): return d1 * (d1 * (d1 * (d1 * d1)))
function code(d1) return Float64(d1 * Float64(d1 * Float64(d1 * Float64(d1 * d1)))) end
function tmp = code(d1) tmp = d1 * (d1 * (d1 * (d1 * d1))); end
code[d1_] := N[(d1 * N[(d1 * N[(d1 * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d1 \cdot \left(d1 \cdot \left(d1 \cdot \left(d1 \cdot d1\right)\right)\right)
\end{array}
Initial program 99.9%
Applied egg-rr86.9%
rem-exp-logN/A
log-prodN/A
flip-+N/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
metadata-evalN/A
+-inversesN/A
+-inversesN/A
associate-*l/N/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
+-inversesN/A
flip-+N/A
log-prodN/A
pow-to-expN/A
unpow-prod-downN/A
sqr-powN/A
unpow163.7%
Applied egg-rr63.7%
Final simplification63.7%
(FPCore (d1) :precision binary64 (* (* d1 d1) (* d1 d1)))
double code(double d1) {
return (d1 * d1) * (d1 * d1);
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = (d1 * d1) * (d1 * d1)
end function
public static double code(double d1) {
return (d1 * d1) * (d1 * d1);
}
def code(d1): return (d1 * d1) * (d1 * d1)
function code(d1) return Float64(Float64(d1 * d1) * Float64(d1 * d1)) end
function tmp = code(d1) tmp = (d1 * d1) * (d1 * d1); end
code[d1_] := N[(N[(d1 * d1), $MachinePrecision] * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(d1 \cdot d1\right) \cdot \left(d1 \cdot d1\right)
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr68.1%
Applied egg-rr75.8%
(FPCore (d1) :precision binary64 (* d1 (* d1 d1)))
double code(double d1) {
return d1 * (d1 * d1);
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 * (d1 * d1)
end function
public static double code(double d1) {
return d1 * (d1 * d1);
}
def code(d1): return d1 * (d1 * d1)
function code(d1) return Float64(d1 * Float64(d1 * d1)) end
function tmp = code(d1) tmp = d1 * (d1 * d1); end
code[d1_] := N[(d1 * N[(d1 * d1), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
d1 \cdot \left(d1 \cdot d1\right)
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr53.0%
Final simplification53.0%
(FPCore (d1) :precision binary64 (* d1 d1))
double code(double d1) {
return d1 * d1;
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 * d1
end function
public static double code(double d1) {
return d1 * d1;
}
def code(d1): return d1 * d1
function code(d1) return Float64(d1 * d1) end
function tmp = code(d1) tmp = d1 * d1; end
code[d1_] := N[(d1 * d1), $MachinePrecision]
\begin{array}{l}
\\
d1 \cdot d1
\end{array}
Initial program 99.9%
Applied egg-rr51.4%
(FPCore (d1) :precision binary64 d1)
double code(double d1) {
return d1;
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1
end function
public static double code(double d1) {
return d1;
}
def code(d1): return d1
function code(d1) return d1 end
function tmp = code(d1) tmp = d1; end
code[d1_] := d1
\begin{array}{l}
\\
d1
\end{array}
Initial program 99.9%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.8%
Simplified99.8%
Applied egg-rr4.5%
(FPCore (d1) :precision binary64 (pow d1 10.0))
double code(double d1) {
return pow(d1, 10.0);
}
real(8) function code(d1)
real(8), intent (in) :: d1
code = d1 ** 10.0d0
end function
public static double code(double d1) {
return Math.pow(d1, 10.0);
}
def code(d1): return math.pow(d1, 10.0)
function code(d1) return d1 ^ 10.0 end
function tmp = code(d1) tmp = d1 ^ 10.0; end
code[d1_] := N[Power[d1, 10.0], $MachinePrecision]
\begin{array}{l}
\\
{d1}^{10}
\end{array}
herbie shell --seed 2024184
(FPCore (d1)
:name "FastMath test5"
:precision binary64
:alt
(! :herbie-platform default (pow d1 10))
(* (* d1 (* (* (* (* (* d1 (* d1 d1)) d1) d1) (* d1 d1)) d1)) d1))