
(FPCore (x y z t) :precision binary64 (- (+ (* (- x 1.0) (log y)) (* (- z 1.0) (log (- 1.0 y)))) t))
double code(double x, double y, double z, double t) {
return (((x - 1.0) * log(y)) + ((z - 1.0) * log((1.0 - y)))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x - 1.0d0) * log(y)) + ((z - 1.0d0) * log((1.0d0 - y)))) - t
end function
public static double code(double x, double y, double z, double t) {
return (((x - 1.0) * Math.log(y)) + ((z - 1.0) * Math.log((1.0 - y)))) - t;
}
def code(x, y, z, t): return (((x - 1.0) * math.log(y)) + ((z - 1.0) * math.log((1.0 - y)))) - t
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x - 1.0) * log(y)) + Float64(Float64(z - 1.0) * log(Float64(1.0 - y)))) - t) end
function tmp = code(x, y, z, t) tmp = (((x - 1.0) * log(y)) + ((z - 1.0) * log((1.0 - y)))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x - 1.0), $MachinePrecision] * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(N[(z - 1.0), $MachinePrecision] * N[Log[N[(1.0 - y), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - 1\right) \cdot \log y + \left(z - 1\right) \cdot \log \left(1 - y\right)\right) - t
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (- (+ (* (- x 1.0) (log y)) (* (- z 1.0) (log (- 1.0 y)))) t))
double code(double x, double y, double z, double t) {
return (((x - 1.0) * log(y)) + ((z - 1.0) * log((1.0 - y)))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x - 1.0d0) * log(y)) + ((z - 1.0d0) * log((1.0d0 - y)))) - t
end function
public static double code(double x, double y, double z, double t) {
return (((x - 1.0) * Math.log(y)) + ((z - 1.0) * Math.log((1.0 - y)))) - t;
}
def code(x, y, z, t): return (((x - 1.0) * math.log(y)) + ((z - 1.0) * math.log((1.0 - y)))) - t
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x - 1.0) * log(y)) + Float64(Float64(z - 1.0) * log(Float64(1.0 - y)))) - t) end
function tmp = code(x, y, z, t) tmp = (((x - 1.0) * log(y)) + ((z - 1.0) * log((1.0 - y)))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x - 1.0), $MachinePrecision] * N[Log[y], $MachinePrecision]), $MachinePrecision] + N[(N[(z - 1.0), $MachinePrecision] * N[Log[N[(1.0 - y), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - 1\right) \cdot \log y + \left(z - 1\right) \cdot \log \left(1 - y\right)\right) - t
\end{array}
(FPCore (x y z t) :precision binary64 (- (+ (* (log y) (+ x -1.0)) (* (log1p (- 0.0 y)) (+ z -1.0))) t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + (log1p((0.0 - y)) * (z + -1.0))) - t;
}
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + (Math.log1p((0.0 - y)) * (z + -1.0))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + (math.log1p((0.0 - y)) * (z + -1.0))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(log1p(Float64(0.0 - y)) * Float64(z + -1.0))) - t) end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(N[Log[1 + N[(0.0 - y), $MachinePrecision]], $MachinePrecision] * N[(z + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + \mathsf{log1p}\left(0 - y\right) \cdot \left(z + -1\right)\right) - t
\end{array}
Initial program 87.9%
sub-negN/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
neg-sub0N/A
--lowering--.f6499.8%
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (x y z t)
:precision binary64
(-
(+
(* (log y) (+ x -1.0))
(/
(*
(* y (+ z -1.0))
(+
-1.0
(*
(* y (* y y))
(+
(*
y
(+ (* y (+ (* y -0.28703703703703703) -0.3541666666666667)) -0.25))
-0.125))))
(+
1.0
(*
(* y (+ -0.5 (* y (+ -0.3333333333333333 (* y -0.25)))))
(- (* y (+ -0.5 (* y -0.3333333333333333))) -1.0)))))
t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * ((y * ((y * ((y * -0.28703703703703703) + -0.3541666666666667)) + -0.25)) + -0.125)))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + (((y * (z + (-1.0d0))) * ((-1.0d0) + ((y * (y * y)) * ((y * ((y * ((y * (-0.28703703703703703d0)) + (-0.3541666666666667d0))) + (-0.25d0))) + (-0.125d0))))) / (1.0d0 + ((y * ((-0.5d0) + (y * ((-0.3333333333333333d0) + (y * (-0.25d0)))))) * ((y * ((-0.5d0) + (y * (-0.3333333333333333d0)))) - (-1.0d0)))))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * ((y * ((y * ((y * -0.28703703703703703) + -0.3541666666666667)) + -0.25)) + -0.125)))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * ((y * ((y * ((y * -0.28703703703703703) + -0.3541666666666667)) + -0.25)) + -0.125)))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(Float64(Float64(y * Float64(z + -1.0)) * Float64(-1.0 + Float64(Float64(y * Float64(y * y)) * Float64(Float64(y * Float64(Float64(y * Float64(Float64(y * -0.28703703703703703) + -0.3541666666666667)) + -0.25)) + -0.125)))) / Float64(1.0 + Float64(Float64(y * Float64(-0.5 + Float64(y * Float64(-0.3333333333333333 + Float64(y * -0.25))))) * Float64(Float64(y * Float64(-0.5 + Float64(y * -0.3333333333333333))) - -1.0))))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * ((y * ((y * ((y * -0.28703703703703703) + -0.3541666666666667)) + -0.25)) + -0.125)))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * N[(z + -1.0), $MachinePrecision]), $MachinePrecision] * N[(-1.0 + N[(N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(N[(y * N[(N[(y * -0.28703703703703703), $MachinePrecision] + -0.3541666666666667), $MachinePrecision]), $MachinePrecision] + -0.25), $MachinePrecision]), $MachinePrecision] + -0.125), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(N[(y * N[(-0.5 + N[(y * N[(-0.3333333333333333 + N[(y * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(-0.5 + N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + \frac{\left(y \cdot \left(z + -1\right)\right) \cdot \left(-1 + \left(y \cdot \left(y \cdot y\right)\right) \cdot \left(y \cdot \left(y \cdot \left(y \cdot -0.28703703703703703 + -0.3541666666666667\right) + -0.25\right) + -0.125\right)\right)}{1 + \left(y \cdot \left(-0.5 + y \cdot \left(-0.3333333333333333 + y \cdot -0.25\right)\right)\right) \cdot \left(y \cdot \left(-0.5 + y \cdot -0.3333333333333333\right) - -1\right)}\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
associate-*r*N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
flip3-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr99.7%
Taylor expanded in y around 0
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Taylor expanded in y around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y z t)
:precision binary64
(-
(+
(* (log y) (+ x -1.0))
(/
(*
(* y (+ z -1.0))
(+
-1.0
(* (* y (* y y)) (+ -0.125 (* y (+ -0.25 (* y -0.3541666666666667)))))))
(+
1.0
(*
(* y (+ -0.5 (* y (+ -0.3333333333333333 (* y -0.25)))))
(- (* y (+ -0.5 (* y -0.3333333333333333))) -1.0)))))
t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * (-0.25 + (y * -0.3541666666666667))))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + (((y * (z + (-1.0d0))) * ((-1.0d0) + ((y * (y * y)) * ((-0.125d0) + (y * ((-0.25d0) + (y * (-0.3541666666666667d0)))))))) / (1.0d0 + ((y * ((-0.5d0) + (y * ((-0.3333333333333333d0) + (y * (-0.25d0)))))) * ((y * ((-0.5d0) + (y * (-0.3333333333333333d0)))) - (-1.0d0)))))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * (-0.25 + (y * -0.3541666666666667))))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * (-0.25 + (y * -0.3541666666666667))))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(Float64(Float64(y * Float64(z + -1.0)) * Float64(-1.0 + Float64(Float64(y * Float64(y * y)) * Float64(-0.125 + Float64(y * Float64(-0.25 + Float64(y * -0.3541666666666667))))))) / Float64(1.0 + Float64(Float64(y * Float64(-0.5 + Float64(y * Float64(-0.3333333333333333 + Float64(y * -0.25))))) * Float64(Float64(y * Float64(-0.5 + Float64(y * -0.3333333333333333))) - -1.0))))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * (-0.25 + (y * -0.3541666666666667))))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * N[(z + -1.0), $MachinePrecision]), $MachinePrecision] * N[(-1.0 + N[(N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision] * N[(-0.125 + N[(y * N[(-0.25 + N[(y * -0.3541666666666667), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(N[(y * N[(-0.5 + N[(y * N[(-0.3333333333333333 + N[(y * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(-0.5 + N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + \frac{\left(y \cdot \left(z + -1\right)\right) \cdot \left(-1 + \left(y \cdot \left(y \cdot y\right)\right) \cdot \left(-0.125 + y \cdot \left(-0.25 + y \cdot -0.3541666666666667\right)\right)\right)}{1 + \left(y \cdot \left(-0.5 + y \cdot \left(-0.3333333333333333 + y \cdot -0.25\right)\right)\right) \cdot \left(y \cdot \left(-0.5 + y \cdot -0.3333333333333333\right) - -1\right)}\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
associate-*r*N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
flip3-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr99.7%
Taylor expanded in y around 0
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Taylor expanded in y around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y z t)
:precision binary64
(-
(+
(* (log y) (+ x -1.0))
(/
(* (* y (+ z -1.0)) (+ -1.0 (* (* y (* y y)) (+ -0.125 (* y -0.25)))))
(+
1.0
(*
(* y (+ -0.5 (* y (+ -0.3333333333333333 (* y -0.25)))))
(- (* y (+ -0.5 (* y -0.3333333333333333))) -1.0)))))
t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * -0.25))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + (((y * (z + (-1.0d0))) * ((-1.0d0) + ((y * (y * y)) * ((-0.125d0) + (y * (-0.25d0)))))) / (1.0d0 + ((y * ((-0.5d0) + (y * ((-0.3333333333333333d0) + (y * (-0.25d0)))))) * ((y * ((-0.5d0) + (y * (-0.3333333333333333d0)))) - (-1.0d0)))))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * -0.25))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * -0.25))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(Float64(Float64(y * Float64(z + -1.0)) * Float64(-1.0 + Float64(Float64(y * Float64(y * y)) * Float64(-0.125 + Float64(y * -0.25))))) / Float64(1.0 + Float64(Float64(y * Float64(-0.5 + Float64(y * Float64(-0.3333333333333333 + Float64(y * -0.25))))) * Float64(Float64(y * Float64(-0.5 + Float64(y * -0.3333333333333333))) - -1.0))))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + (((y * (z + -1.0)) * (-1.0 + ((y * (y * y)) * (-0.125 + (y * -0.25))))) / (1.0 + ((y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))) * ((y * (-0.5 + (y * -0.3333333333333333))) - -1.0))))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * N[(z + -1.0), $MachinePrecision]), $MachinePrecision] * N[(-1.0 + N[(N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision] * N[(-0.125 + N[(y * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(N[(y * N[(-0.5 + N[(y * N[(-0.3333333333333333 + N[(y * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(y * N[(-0.5 + N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + \frac{\left(y \cdot \left(z + -1\right)\right) \cdot \left(-1 + \left(y \cdot \left(y \cdot y\right)\right) \cdot \left(-0.125 + y \cdot -0.25\right)\right)}{1 + \left(y \cdot \left(-0.5 + y \cdot \left(-0.3333333333333333 + y \cdot -0.25\right)\right)\right) \cdot \left(y \cdot \left(-0.5 + y \cdot -0.3333333333333333\right) - -1\right)}\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
associate-*r*N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
flip3-+N/A
associate-*r/N/A
/-lowering-/.f64N/A
Applied egg-rr99.7%
Taylor expanded in y around 0
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Taylor expanded in y around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y z t)
:precision binary64
(-
(+
(* (log y) (+ x -1.0))
(*
(* y (+ -1.0 (* y (+ -0.5 (* y (+ -0.3333333333333333 (* y -0.25)))))))
(+ z -1.0)))
t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + ((y * (-1.0 + (y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))))) * (z + -1.0))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + ((y * ((-1.0d0) + (y * ((-0.5d0) + (y * ((-0.3333333333333333d0) + (y * (-0.25d0)))))))) * (z + (-1.0d0)))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + ((y * (-1.0 + (y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))))) * (z + -1.0))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + ((y * (-1.0 + (y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))))) * (z + -1.0))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(Float64(y * Float64(-1.0 + Float64(y * Float64(-0.5 + Float64(y * Float64(-0.3333333333333333 + Float64(y * -0.25))))))) * Float64(z + -1.0))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + ((y * (-1.0 + (y * (-0.5 + (y * (-0.3333333333333333 + (y * -0.25))))))) * (z + -1.0))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(N[(y * N[(-1.0 + N[(y * N[(-0.5 + N[(y * N[(-0.3333333333333333 + N[(y * -0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(z + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + \left(y \cdot \left(-1 + y \cdot \left(-0.5 + y \cdot \left(-0.3333333333333333 + y \cdot -0.25\right)\right)\right)\right) \cdot \left(z + -1\right)\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x y z t)
:precision binary64
(-
(+
(* (log y) (+ x -1.0))
(*
y
(- (* (* y (+ -0.5 (* y -0.3333333333333333))) (+ z -1.0)) (+ z -1.0))))
t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + (y * (((y * (-0.5 + (y * -0.3333333333333333))) * (z + -1.0)) - (z + -1.0)))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + (y * (((y * ((-0.5d0) + (y * (-0.3333333333333333d0)))) * (z + (-1.0d0))) - (z + (-1.0d0))))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + (y * (((y * (-0.5 + (y * -0.3333333333333333))) * (z + -1.0)) - (z + -1.0)))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + (y * (((y * (-0.5 + (y * -0.3333333333333333))) * (z + -1.0)) - (z + -1.0)))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(y * Float64(Float64(Float64(y * Float64(-0.5 + Float64(y * -0.3333333333333333))) * Float64(z + -1.0)) - Float64(z + -1.0)))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + (y * (((y * (-0.5 + (y * -0.3333333333333333))) * (z + -1.0)) - (z + -1.0)))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(y * N[(N[(N[(y * N[(-0.5 + N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(z + -1.0), $MachinePrecision]), $MachinePrecision] - N[(z + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + y \cdot \left(\left(y \cdot \left(-0.5 + y \cdot -0.3333333333333333\right)\right) \cdot \left(z + -1\right) - \left(z + -1\right)\right)\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
Simplified99.6%
Final simplification99.6%
(FPCore (x y z t) :precision binary64 (- (+ (* (log y) (+ x -1.0)) (* (+ z -1.0) (* y (+ -1.0 (* y (+ -0.5 (* y -0.3333333333333333))))))) t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + ((z + -1.0) * (y * (-1.0 + (y * (-0.5 + (y * -0.3333333333333333))))))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + ((z + (-1.0d0)) * (y * ((-1.0d0) + (y * ((-0.5d0) + (y * (-0.3333333333333333d0)))))))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + ((z + -1.0) * (y * (-1.0 + (y * (-0.5 + (y * -0.3333333333333333))))))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + ((z + -1.0) * (y * (-1.0 + (y * (-0.5 + (y * -0.3333333333333333))))))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(Float64(z + -1.0) * Float64(y * Float64(-1.0 + Float64(y * Float64(-0.5 + Float64(y * -0.3333333333333333))))))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + ((z + -1.0) * (y * (-1.0 + (y * (-0.5 + (y * -0.3333333333333333))))))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(N[(z + -1.0), $MachinePrecision] * N[(y * N[(-1.0 + N[(y * N[(-0.5 + N[(y * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + \left(z + -1\right) \cdot \left(y \cdot \left(-1 + y \cdot \left(-0.5 + y \cdot -0.3333333333333333\right)\right)\right)\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.6%
Simplified99.6%
Final simplification99.6%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* (log y) (+ x -1.0)) t)))
(if (<= x -7.6e-7)
t_1
(if (<= x 1.7e-39) (- (- (log (/ 1.0 y)) (* y (+ z -1.0))) t) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (log(y) * (x + -1.0)) - t;
double tmp;
if (x <= -7.6e-7) {
tmp = t_1;
} else if (x <= 1.7e-39) {
tmp = (log((1.0 / y)) - (y * (z + -1.0))) - t;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (log(y) * (x + (-1.0d0))) - t
if (x <= (-7.6d-7)) then
tmp = t_1
else if (x <= 1.7d-39) then
tmp = (log((1.0d0 / y)) - (y * (z + (-1.0d0)))) - t
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (Math.log(y) * (x + -1.0)) - t;
double tmp;
if (x <= -7.6e-7) {
tmp = t_1;
} else if (x <= 1.7e-39) {
tmp = (Math.log((1.0 / y)) - (y * (z + -1.0))) - t;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (math.log(y) * (x + -1.0)) - t tmp = 0 if x <= -7.6e-7: tmp = t_1 elif x <= 1.7e-39: tmp = (math.log((1.0 / y)) - (y * (z + -1.0))) - t else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(log(y) * Float64(x + -1.0)) - t) tmp = 0.0 if (x <= -7.6e-7) tmp = t_1; elseif (x <= 1.7e-39) tmp = Float64(Float64(log(Float64(1.0 / y)) - Float64(y * Float64(z + -1.0))) - t); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (log(y) * (x + -1.0)) - t; tmp = 0.0; if (x <= -7.6e-7) tmp = t_1; elseif (x <= 1.7e-39) tmp = (log((1.0 / y)) - (y * (z + -1.0))) - t; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]}, If[LessEqual[x, -7.6e-7], t$95$1, If[LessEqual[x, 1.7e-39], N[(N[(N[Log[N[(1.0 / y), $MachinePrecision]], $MachinePrecision] - N[(y * N[(z + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log y \cdot \left(x + -1\right) - t\\
\mathbf{if}\;x \leq -7.6 \cdot 10^{-7}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 1.7 \cdot 10^{-39}:\\
\;\;\;\;\left(\log \left(\frac{1}{y}\right) - y \cdot \left(z + -1\right)\right) - t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -7.60000000000000029e-7 or 1.7e-39 < x Initial program 94.4%
Taylor expanded in y around 0
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6493.2%
Simplified93.2%
if -7.60000000000000029e-7 < x < 1.7e-39Initial program 80.8%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.1%
Simplified98.1%
Taylor expanded in x around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
log-lowering-log.f6498.1%
Simplified98.1%
sub0-negN/A
neg-logN/A
log-lowering-log.f64N/A
/-lowering-/.f6498.1%
Applied egg-rr98.1%
Final simplification95.6%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (log y) (+ x -1.0))))
(if (<= x -1.9e-5)
(- t_1 t)
(if (<= x 5.5e-44) (- (- 0.0 t) (+ (log y) (* y z))) (- (+ y t_1) t)))))
double code(double x, double y, double z, double t) {
double t_1 = log(y) * (x + -1.0);
double tmp;
if (x <= -1.9e-5) {
tmp = t_1 - t;
} else if (x <= 5.5e-44) {
tmp = (0.0 - t) - (log(y) + (y * z));
} else {
tmp = (y + t_1) - t;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = log(y) * (x + (-1.0d0))
if (x <= (-1.9d-5)) then
tmp = t_1 - t
else if (x <= 5.5d-44) then
tmp = (0.0d0 - t) - (log(y) + (y * z))
else
tmp = (y + t_1) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = Math.log(y) * (x + -1.0);
double tmp;
if (x <= -1.9e-5) {
tmp = t_1 - t;
} else if (x <= 5.5e-44) {
tmp = (0.0 - t) - (Math.log(y) + (y * z));
} else {
tmp = (y + t_1) - t;
}
return tmp;
}
def code(x, y, z, t): t_1 = math.log(y) * (x + -1.0) tmp = 0 if x <= -1.9e-5: tmp = t_1 - t elif x <= 5.5e-44: tmp = (0.0 - t) - (math.log(y) + (y * z)) else: tmp = (y + t_1) - t return tmp
function code(x, y, z, t) t_1 = Float64(log(y) * Float64(x + -1.0)) tmp = 0.0 if (x <= -1.9e-5) tmp = Float64(t_1 - t); elseif (x <= 5.5e-44) tmp = Float64(Float64(0.0 - t) - Float64(log(y) + Float64(y * z))); else tmp = Float64(Float64(y + t_1) - t); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = log(y) * (x + -1.0); tmp = 0.0; if (x <= -1.9e-5) tmp = t_1 - t; elseif (x <= 5.5e-44) tmp = (0.0 - t) - (log(y) + (y * z)); else tmp = (y + t_1) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.9e-5], N[(t$95$1 - t), $MachinePrecision], If[LessEqual[x, 5.5e-44], N[(N[(0.0 - t), $MachinePrecision] - N[(N[Log[y], $MachinePrecision] + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y + t$95$1), $MachinePrecision] - t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \log y \cdot \left(x + -1\right)\\
\mathbf{if}\;x \leq -1.9 \cdot 10^{-5}:\\
\;\;\;\;t\_1 - t\\
\mathbf{elif}\;x \leq 5.5 \cdot 10^{-44}:\\
\;\;\;\;\left(0 - t\right) - \left(\log y + y \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y + t\_1\right) - t\\
\end{array}
\end{array}
if x < -1.9000000000000001e-5Initial program 95.0%
Taylor expanded in y around 0
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6494.3%
Simplified94.3%
if -1.9000000000000001e-5 < x < 5.49999999999999993e-44Initial program 80.3%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.1%
Simplified98.1%
Taylor expanded in x around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
log-lowering-log.f6498.1%
Simplified98.1%
Taylor expanded in z around inf
*-lowering-*.f6497.0%
Simplified97.0%
if 5.49999999999999993e-44 < x Initial program 94.0%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.8%
Simplified98.8%
Taylor expanded in z around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6492.3%
Simplified92.3%
Final simplification95.1%
(FPCore (x y z t) :precision binary64 (- (+ (* (log y) (+ x -1.0)) (* y (* (+ z -1.0) (+ -1.0 (* y -0.5))))) t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) + (y * ((z + -1.0) * (-1.0 + (y * -0.5))))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) + (y * ((z + (-1.0d0)) * ((-1.0d0) + (y * (-0.5d0)))))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) + (y * ((z + -1.0) * (-1.0 + (y * -0.5))))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) + (y * ((z + -1.0) * (-1.0 + (y * -0.5))))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) + Float64(y * Float64(Float64(z + -1.0) * Float64(-1.0 + Float64(y * -0.5))))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) + (y * ((z + -1.0) * (-1.0 + (y * -0.5))))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] + N[(y * N[(N[(z + -1.0), $MachinePrecision] * N[(-1.0 + N[(y * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) + y \cdot \left(\left(z + -1\right) \cdot \left(-1 + y \cdot -0.5\right)\right)\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6499.4%
Simplified99.4%
Final simplification99.4%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (- (* x (log y)) t))) (if (<= x -1450.0) t_1 (if (<= x 7400000000.0) (- (log (/ 1.0 y)) t) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (x * log(y)) - t;
double tmp;
if (x <= -1450.0) {
tmp = t_1;
} else if (x <= 7400000000.0) {
tmp = log((1.0 / y)) - t;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * log(y)) - t
if (x <= (-1450.0d0)) then
tmp = t_1
else if (x <= 7400000000.0d0) then
tmp = log((1.0d0 / y)) - t
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * Math.log(y)) - t;
double tmp;
if (x <= -1450.0) {
tmp = t_1;
} else if (x <= 7400000000.0) {
tmp = Math.log((1.0 / y)) - t;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * math.log(y)) - t tmp = 0 if x <= -1450.0: tmp = t_1 elif x <= 7400000000.0: tmp = math.log((1.0 / y)) - t else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * log(y)) - t) tmp = 0.0 if (x <= -1450.0) tmp = t_1; elseif (x <= 7400000000.0) tmp = Float64(log(Float64(1.0 / y)) - t); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * log(y)) - t; tmp = 0.0; if (x <= -1450.0) tmp = t_1; elseif (x <= 7400000000.0) tmp = log((1.0 / y)) - t; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]}, If[LessEqual[x, -1450.0], t$95$1, If[LessEqual[x, 7400000000.0], N[(N[Log[N[(1.0 / y), $MachinePrecision]], $MachinePrecision] - t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \log y - t\\
\mathbf{if}\;x \leq -1450:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 7400000000:\\
\;\;\;\;\log \left(\frac{1}{y}\right) - t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -1450 or 7.4e9 < x Initial program 93.8%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6492.5%
Simplified92.5%
if -1450 < x < 7.4e9Initial program 83.0%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.1%
Simplified98.1%
Taylor expanded in x around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
log-lowering-log.f6497.1%
Simplified97.1%
Taylor expanded in y around 0
mul-1-negN/A
log-recN/A
log-lowering-log.f64N/A
/-lowering-/.f6478.3%
Simplified78.3%
Final simplification84.7%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (* x (log y)))) (if (<= x -6.9e+78) t_1 (if (<= x 1.95e+127) (- (log (/ 1.0 y)) t) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x * log(y);
double tmp;
if (x <= -6.9e+78) {
tmp = t_1;
} else if (x <= 1.95e+127) {
tmp = log((1.0 / y)) - t;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x * log(y)
if (x <= (-6.9d+78)) then
tmp = t_1
else if (x <= 1.95d+127) then
tmp = log((1.0d0 / y)) - t
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x * Math.log(y);
double tmp;
if (x <= -6.9e+78) {
tmp = t_1;
} else if (x <= 1.95e+127) {
tmp = Math.log((1.0 / y)) - t;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x * math.log(y) tmp = 0 if x <= -6.9e+78: tmp = t_1 elif x <= 1.95e+127: tmp = math.log((1.0 / y)) - t else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x * log(y)) tmp = 0.0 if (x <= -6.9e+78) tmp = t_1; elseif (x <= 1.95e+127) tmp = Float64(log(Float64(1.0 / y)) - t); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x * log(y); tmp = 0.0; if (x <= -6.9e+78) tmp = t_1; elseif (x <= 1.95e+127) tmp = log((1.0 / y)) - t; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -6.9e+78], t$95$1, If[LessEqual[x, 1.95e+127], N[(N[Log[N[(1.0 / y), $MachinePrecision]], $MachinePrecision] - t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \log y\\
\mathbf{if}\;x \leq -6.9 \cdot 10^{+78}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 1.95 \cdot 10^{+127}:\\
\;\;\;\;\log \left(\frac{1}{y}\right) - t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -6.8999999999999998e78 or 1.94999999999999991e127 < x Initial program 96.5%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6474.4%
Simplified74.4%
if -6.8999999999999998e78 < x < 1.94999999999999991e127Initial program 84.3%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.3%
Simplified98.3%
Taylor expanded in x around 0
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
log-lowering-log.f6490.7%
Simplified90.7%
Taylor expanded in y around 0
mul-1-negN/A
log-recN/A
log-lowering-log.f64N/A
/-lowering-/.f6473.5%
Simplified73.5%
Final simplification73.8%
(FPCore (x y z t) :precision binary64 (if (<= z -3.4e+247) (- (* y (- 1.0 z)) t) (- (+ y (* (log y) (+ x -1.0))) t)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -3.4e+247) {
tmp = (y * (1.0 - z)) - t;
} else {
tmp = (y + (log(y) * (x + -1.0))) - t;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-3.4d+247)) then
tmp = (y * (1.0d0 - z)) - t
else
tmp = (y + (log(y) * (x + (-1.0d0)))) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -3.4e+247) {
tmp = (y * (1.0 - z)) - t;
} else {
tmp = (y + (Math.log(y) * (x + -1.0))) - t;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -3.4e+247: tmp = (y * (1.0 - z)) - t else: tmp = (y + (math.log(y) * (x + -1.0))) - t return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -3.4e+247) tmp = Float64(Float64(y * Float64(1.0 - z)) - t); else tmp = Float64(Float64(y + Float64(log(y) * Float64(x + -1.0))) - t); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -3.4e+247) tmp = (y * (1.0 - z)) - t; else tmp = (y + (log(y) * (x + -1.0))) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -3.4e+247], N[(N[(y * N[(1.0 - z), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[(y + N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.4 \cdot 10^{+247}:\\
\;\;\;\;y \cdot \left(1 - z\right) - t\\
\mathbf{else}:\\
\;\;\;\;\left(y + \log y \cdot \left(x + -1\right)\right) - t\\
\end{array}
\end{array}
if z < -3.39999999999999981e247Initial program 43.3%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6497.1%
Simplified97.1%
Taylor expanded in y around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6492.1%
Simplified92.1%
if -3.39999999999999981e247 < z Initial program 91.6%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.8%
Simplified98.8%
Taylor expanded in z around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6490.3%
Simplified90.3%
Final simplification90.4%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (* x (log y)))) (if (<= x -2.05e+79) t_1 (if (<= x 1e+175) (- (* y (- 1.0 z)) t) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = x * log(y);
double tmp;
if (x <= -2.05e+79) {
tmp = t_1;
} else if (x <= 1e+175) {
tmp = (y * (1.0 - z)) - t;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = x * log(y)
if (x <= (-2.05d+79)) then
tmp = t_1
else if (x <= 1d+175) then
tmp = (y * (1.0d0 - z)) - t
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = x * Math.log(y);
double tmp;
if (x <= -2.05e+79) {
tmp = t_1;
} else if (x <= 1e+175) {
tmp = (y * (1.0 - z)) - t;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = x * math.log(y) tmp = 0 if x <= -2.05e+79: tmp = t_1 elif x <= 1e+175: tmp = (y * (1.0 - z)) - t else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(x * log(y)) tmp = 0.0 if (x <= -2.05e+79) tmp = t_1; elseif (x <= 1e+175) tmp = Float64(Float64(y * Float64(1.0 - z)) - t); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = x * log(y); tmp = 0.0; if (x <= -2.05e+79) tmp = t_1; elseif (x <= 1e+175) tmp = (y * (1.0 - z)) - t; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(x * N[Log[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -2.05e+79], t$95$1, If[LessEqual[x, 1e+175], N[(N[(y * N[(1.0 - z), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot \log y\\
\mathbf{if}\;x \leq -2.05 \cdot 10^{+79}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 10^{+175}:\\
\;\;\;\;y \cdot \left(1 - z\right) - t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -2.05e79 or 9.9999999999999994e174 < x Initial program 99.0%
Taylor expanded in x around inf
*-commutativeN/A
*-lowering-*.f64N/A
log-lowering-log.f6478.3%
Simplified78.3%
if -2.05e79 < x < 9.9999999999999994e174Initial program 84.0%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.4%
Simplified98.4%
Taylor expanded in y around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6464.5%
Simplified64.5%
Final simplification68.1%
(FPCore (x y z t) :precision binary64 (- (- (* (log y) (+ x -1.0)) (* y (+ z -1.0))) t))
double code(double x, double y, double z, double t) {
return ((log(y) * (x + -1.0)) - (y * (z + -1.0))) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((log(y) * (x + (-1.0d0))) - (y * (z + (-1.0d0)))) - t
end function
public static double code(double x, double y, double z, double t) {
return ((Math.log(y) * (x + -1.0)) - (y * (z + -1.0))) - t;
}
def code(x, y, z, t): return ((math.log(y) * (x + -1.0)) - (y * (z + -1.0))) - t
function code(x, y, z, t) return Float64(Float64(Float64(log(y) * Float64(x + -1.0)) - Float64(y * Float64(z + -1.0))) - t) end
function tmp = code(x, y, z, t) tmp = ((log(y) * (x + -1.0)) - (y * (z + -1.0))) - t; end
code[x_, y_, z_, t_] := N[(N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] - N[(y * N[(z + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
\left(\log y \cdot \left(x + -1\right) - y \cdot \left(z + -1\right)\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.6%
Simplified98.6%
Final simplification98.6%
(FPCore (x y z t) :precision binary64 (if (<= z -2.75e+247) (- (* y (- 1.0 z)) t) (- (* (log y) (+ x -1.0)) t)))
double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.75e+247) {
tmp = (y * (1.0 - z)) - t;
} else {
tmp = (log(y) * (x + -1.0)) - t;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (z <= (-2.75d+247)) then
tmp = (y * (1.0d0 - z)) - t
else
tmp = (log(y) * (x + (-1.0d0))) - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (z <= -2.75e+247) {
tmp = (y * (1.0 - z)) - t;
} else {
tmp = (Math.log(y) * (x + -1.0)) - t;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if z <= -2.75e+247: tmp = (y * (1.0 - z)) - t else: tmp = (math.log(y) * (x + -1.0)) - t return tmp
function code(x, y, z, t) tmp = 0.0 if (z <= -2.75e+247) tmp = Float64(Float64(y * Float64(1.0 - z)) - t); else tmp = Float64(Float64(log(y) * Float64(x + -1.0)) - t); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (z <= -2.75e+247) tmp = (y * (1.0 - z)) - t; else tmp = (log(y) * (x + -1.0)) - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[z, -2.75e+247], N[(N[(y * N[(1.0 - z), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision], N[(N[(N[Log[y], $MachinePrecision] * N[(x + -1.0), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.75 \cdot 10^{+247}:\\
\;\;\;\;y \cdot \left(1 - z\right) - t\\
\mathbf{else}:\\
\;\;\;\;\log y \cdot \left(x + -1\right) - t\\
\end{array}
\end{array}
if z < -2.7499999999999999e247Initial program 43.3%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6497.1%
Simplified97.1%
Taylor expanded in y around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6492.1%
Simplified92.1%
if -2.7499999999999999e247 < z Initial program 91.6%
Taylor expanded in y around 0
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6489.7%
Simplified89.7%
Final simplification89.9%
(FPCore (x y z t) :precision binary64 (if (<= t -65000000000.0) (- 0.0 t) (if (<= t 200.0) (* y (- 1.0 z)) (- 0.0 t))))
double code(double x, double y, double z, double t) {
double tmp;
if (t <= -65000000000.0) {
tmp = 0.0 - t;
} else if (t <= 200.0) {
tmp = y * (1.0 - z);
} else {
tmp = 0.0 - t;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-65000000000.0d0)) then
tmp = 0.0d0 - t
else if (t <= 200.0d0) then
tmp = y * (1.0d0 - z)
else
tmp = 0.0d0 - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (t <= -65000000000.0) {
tmp = 0.0 - t;
} else if (t <= 200.0) {
tmp = y * (1.0 - z);
} else {
tmp = 0.0 - t;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if t <= -65000000000.0: tmp = 0.0 - t elif t <= 200.0: tmp = y * (1.0 - z) else: tmp = 0.0 - t return tmp
function code(x, y, z, t) tmp = 0.0 if (t <= -65000000000.0) tmp = Float64(0.0 - t); elseif (t <= 200.0) tmp = Float64(y * Float64(1.0 - z)); else tmp = Float64(0.0 - t); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (t <= -65000000000.0) tmp = 0.0 - t; elseif (t <= 200.0) tmp = y * (1.0 - z); else tmp = 0.0 - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[t, -65000000000.0], N[(0.0 - t), $MachinePrecision], If[LessEqual[t, 200.0], N[(y * N[(1.0 - z), $MachinePrecision]), $MachinePrecision], N[(0.0 - t), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -65000000000:\\
\;\;\;\;0 - t\\
\mathbf{elif}\;t \leq 200:\\
\;\;\;\;y \cdot \left(1 - z\right)\\
\mathbf{else}:\\
\;\;\;\;0 - t\\
\end{array}
\end{array}
if t < -6.5e10 or 200 < t Initial program 95.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6475.8%
Simplified75.8%
sub0-negN/A
neg-lowering-neg.f6475.8%
Applied egg-rr75.8%
if -6.5e10 < t < 200Initial program 79.0%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6497.7%
Simplified97.7%
Taylor expanded in y around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6424.2%
Simplified24.2%
Final simplification51.2%
(FPCore (x y z t) :precision binary64 (if (<= t -72000000000.0) (- 0.0 t) (if (<= t 200.0) (- 0.0 (* y z)) (- 0.0 t))))
double code(double x, double y, double z, double t) {
double tmp;
if (t <= -72000000000.0) {
tmp = 0.0 - t;
} else if (t <= 200.0) {
tmp = 0.0 - (y * z);
} else {
tmp = 0.0 - t;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (t <= (-72000000000.0d0)) then
tmp = 0.0d0 - t
else if (t <= 200.0d0) then
tmp = 0.0d0 - (y * z)
else
tmp = 0.0d0 - t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (t <= -72000000000.0) {
tmp = 0.0 - t;
} else if (t <= 200.0) {
tmp = 0.0 - (y * z);
} else {
tmp = 0.0 - t;
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if t <= -72000000000.0: tmp = 0.0 - t elif t <= 200.0: tmp = 0.0 - (y * z) else: tmp = 0.0 - t return tmp
function code(x, y, z, t) tmp = 0.0 if (t <= -72000000000.0) tmp = Float64(0.0 - t); elseif (t <= 200.0) tmp = Float64(0.0 - Float64(y * z)); else tmp = Float64(0.0 - t); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (t <= -72000000000.0) tmp = 0.0 - t; elseif (t <= 200.0) tmp = 0.0 - (y * z); else tmp = 0.0 - t; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[t, -72000000000.0], N[(0.0 - t), $MachinePrecision], If[LessEqual[t, 200.0], N[(0.0 - N[(y * z), $MachinePrecision]), $MachinePrecision], N[(0.0 - t), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -72000000000:\\
\;\;\;\;0 - t\\
\mathbf{elif}\;t \leq 200:\\
\;\;\;\;0 - y \cdot z\\
\mathbf{else}:\\
\;\;\;\;0 - t\\
\end{array}
\end{array}
if t < -7.2e10 or 200 < t Initial program 95.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6475.8%
Simplified75.8%
sub0-negN/A
neg-lowering-neg.f6475.8%
Applied egg-rr75.8%
if -7.2e10 < t < 200Initial program 79.0%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6497.7%
Simplified97.7%
Taylor expanded in z around inf
mul-1-negN/A
distribute-rgt-neg-inN/A
neg-mul-1N/A
*-lowering-*.f64N/A
neg-mul-1N/A
neg-sub0N/A
--lowering--.f6423.4%
Simplified23.4%
sub0-negN/A
neg-lowering-neg.f6423.4%
Applied egg-rr23.4%
Final simplification50.9%
(FPCore (x y z t) :precision binary64 (- (* y (- 1.0 z)) t))
double code(double x, double y, double z, double t) {
return (y * (1.0 - z)) - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (y * (1.0d0 - z)) - t
end function
public static double code(double x, double y, double z, double t) {
return (y * (1.0 - z)) - t;
}
def code(x, y, z, t): return (y * (1.0 - z)) - t
function code(x, y, z, t) return Float64(Float64(y * Float64(1.0 - z)) - t) end
function tmp = code(x, y, z, t) tmp = (y * (1.0 - z)) - t; end
code[x_, y_, z_, t_] := N[(N[(y * N[(1.0 - z), $MachinePrecision]), $MachinePrecision] - t), $MachinePrecision]
\begin{array}{l}
\\
y \cdot \left(1 - z\right) - t
\end{array}
Initial program 87.9%
Taylor expanded in y around 0
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
log-lowering-log.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f6498.6%
Simplified98.6%
Taylor expanded in y around inf
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
--lowering--.f6453.4%
Simplified53.4%
(FPCore (x y z t) :precision binary64 (- 0.0 t))
double code(double x, double y, double z, double t) {
return 0.0 - t;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 0.0d0 - t
end function
public static double code(double x, double y, double z, double t) {
return 0.0 - t;
}
def code(x, y, z, t): return 0.0 - t
function code(x, y, z, t) return Float64(0.0 - t) end
function tmp = code(x, y, z, t) tmp = 0.0 - t; end
code[x_, y_, z_, t_] := N[(0.0 - t), $MachinePrecision]
\begin{array}{l}
\\
0 - t
\end{array}
Initial program 87.9%
Taylor expanded in t around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f6441.1%
Simplified41.1%
sub0-negN/A
neg-lowering-neg.f6441.1%
Applied egg-rr41.1%
Final simplification41.1%
herbie shell --seed 2024138
(FPCore (x y z t)
:name "Statistics.Distribution.Beta:$cdensity from math-functions-0.1.5.2"
:precision binary64
(- (+ (* (- x 1.0) (log y)) (* (- z 1.0) (log (- 1.0 y)))) t))