Numeric.AD.Rank1.Halley:findZero from ad-4.2.4

Percentage Accurate: 82.4% → 82.4%

Time: 1.7s

Alternatives: 4

Speedup: 1.0×

Specification

\[\mathsf{TRUE}\left(\right)\]

Math

\[\begin{array}{l} \\ x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \end{array} \]

FPCore

(FPCore (x y z t)
 :precision binary64
 (- x (/ (* (* y 2.0) z) (- (* (* z 2.0) z) (* y t)))))

C

double code(double x, double y, double z, double t) {
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
}

Fortran

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 - (((y * 2.0d0) * z) / (((z * 2.0d0) * z) - (y * t)))
end function

Java

public static double code(double x, double y, double z, double t) {
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
}

Python

def code(x, y, z, t):
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)))

Julia

function code(x, y, z, t)
	return Float64(x - Float64(Float64(Float64(y * 2.0) * z) / Float64(Float64(Float64(z * 2.0) * z) - Float64(y * t))))
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
end

Wolfram

code[x_, y_, z_, t_] := N[(x - N[(N[(N[(y * 2.0), $MachinePrecision] * z), $MachinePrecision] / N[(N[(N[(z * 2.0), $MachinePrecision] * z), $MachinePrecision] - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Initial Program: 82.4% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \end{array} \]

FPCore

(FPCore (x y z t)
 :precision binary64
 (- x (/ (* (* y 2.0) z) (- (* (* z 2.0) z) (* y t)))))

C

double code(double x, double y, double z, double t) {
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
}

Fortran

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 - (((y * 2.0d0) * z) / (((z * 2.0d0) * z) - (y * t)))
end function

Java

public static double code(double x, double y, double z, double t) {
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
}

Python

def code(x, y, z, t):
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)))

Julia

function code(x, y, z, t)
	return Float64(x - Float64(Float64(Float64(y * 2.0) * z) / Float64(Float64(Float64(z * 2.0) * z) - Float64(y * t))))
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
end

Wolfram

code[x_, y_, z_, t_] := N[(x - N[(N[(N[(y * 2.0), $MachinePrecision] * z), $MachinePrecision] / N[(N[(N[(z * 2.0), $MachinePrecision] * z), $MachinePrecision] - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Alternative 1: 82.4% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \end{array} \]

FPCore

(FPCore (x y z t)
 :precision binary64
 (- x (/ (* (* y 2.0) z) (- (* (* z 2.0) z) (* y t)))))

C

double code(double x, double y, double z, double t) {
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
}

Fortran

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 - (((y * 2.0d0) * z) / (((z * 2.0d0) * z) - (y * t)))
end function

Java

public static double code(double x, double y, double z, double t) {
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
}

Python

def code(x, y, z, t):
	return x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)))

Julia

function code(x, y, z, t)
	return Float64(x - Float64(Float64(Float64(y * 2.0) * z) / Float64(Float64(Float64(z * 2.0) * z) - Float64(y * t))))
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = x - (((y * 2.0) * z) / (((z * 2.0) * z) - (y * t)));
end

Wolfram

code[x_, y_, z_, t_] := N[(x - N[(N[(N[(y * 2.0), $MachinePrecision] * z), $MachinePrecision] / N[(N[(N[(z * 2.0), $MachinePrecision] * z), $MachinePrecision] - N[(y * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 84.4%

   \[x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \]
2. Add Preprocessing
3. Add Preprocessing

Alternative 2: 11.7% accurate, 1.6× speedup

Math

\[\begin{array}{l} \\ \frac{\left(y \cdot 2\right) \cdot z}{y \cdot t} \end{array} \]

FPCore

(FPCore (x y z t) :precision binary64 (/ (* (* y 2.0) z) (* y t)))

C

double code(double x, double y, double z, double t) {
	return ((y * 2.0) * z) / (y * t);
}

Fortran

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 * 2.0d0) * z) / (y * t)
end function

Java

public static double code(double x, double y, double z, double t) {
	return ((y * 2.0) * z) / (y * t);
}

Python

def code(x, y, z, t):
	return ((y * 2.0) * z) / (y * t)

Julia

function code(x, y, z, t)
	return Float64(Float64(Float64(y * 2.0) * z) / Float64(y * t))
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = ((y * 2.0) * z) / (y * t);
end

Wolfram

code[x_, y_, z_, t_] := N[(N[(N[(y * 2.0), $MachinePrecision] * z), $MachinePrecision] / N[(y * t), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 84.4%

   \[x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \]
2. Add Preprocessing

3. Taylor expanded in x around 0

   \[\leadsto \color{blue}{x - 2 \cdot \frac{y \cdot z}{2 \cdot {z}^{2} - t \cdot y}} \]

5. Applied rewrites 2.9%

   \[\leadsto \color{blue}{\frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t}} \]

6. Taylor expanded in y around 0

   \[\leadsto \frac{\left(y \cdot 2\right) \cdot z}{2 \cdot \color{blue}{{z}^{2}}} \]

8. Applied rewrites 11.2%

   \[\leadsto \frac{\left(y \cdot 2\right) \cdot z}{y \cdot \color{blue}{t}} \]
9. Add Preprocessing

Alternative 3: 3.5% accurate, 7.2× speedup

Math

\[\begin{array}{l} \\ y \cdot 2 \end{array} \]

FPCore

(FPCore (x y z t) :precision binary64 (* y 2.0))

C

double code(double x, double y, double z, double t) {
	return y * 2.0;
}

Fortran

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 * 2.0d0
end function

Java

public static double code(double x, double y, double z, double t) {
	return y * 2.0;
}

Python

def code(x, y, z, t):
	return y * 2.0

Julia

function code(x, y, z, t)
	return Float64(y * 2.0)
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = y * 2.0;
end

Wolfram

code[x_, y_, z_, t_] := N[(y * 2.0), $MachinePrecision]

Derivation

1. Initial program 84.4%

   \[x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \]
2. Add Preprocessing

3. Taylor expanded in x around inf

   \[\leadsto \color{blue}{x \cdot \left(1 + -2 \cdot \frac{y \cdot z}{x \cdot \left(2 \cdot {z}^{2} - t \cdot y\right)}\right)} \]

5. Applied rewrites 3.3%

   \[\leadsto \color{blue}{y \cdot 2} \]
6. Add Preprocessing

Alternative 4: 3.2% accurate, 7.2× speedup

Math

\[\begin{array}{l} \\ y \cdot t \end{array} \]

FPCore

(FPCore (x y z t) :precision binary64 (* y t))

C

double code(double x, double y, double z, double t) {
	return y * t;
}

Fortran

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 * t
end function

Java

public static double code(double x, double y, double z, double t) {
	return y * t;
}

Python

def code(x, y, z, t):
	return y * t

Julia

function code(x, y, z, t)
	return Float64(y * t)
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = y * t;
end

Wolfram

code[x_, y_, z_, t_] := N[(y * t), $MachinePrecision]

Derivation

1. Initial program 84.4%

   \[x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t} \]
2. Add Preprocessing

3. Taylor expanded in x around -inf

   \[\leadsto \color{blue}{-1 \cdot \left(x \cdot \left(2 \cdot \frac{y \cdot z}{x \cdot \left(2 \cdot {z}^{2} - t \cdot y\right)} - 1\right)\right)} \]

5. Applied rewrites 2.8%

   \[\leadsto \color{blue}{\left(z \cdot 2\right) \cdot z - y \cdot t} \]

6. Taylor expanded in y around 0

   \[\leadsto 2 \cdot \color{blue}{{z}^{2}} \]

8. Applied rewrites 3.0%

   \[\leadsto y \cdot \color{blue}{t} \]
9. Add Preprocessing

Developer Target 1: 99.9% accurate, 0.8× speedup

Math

\[\begin{array}{l} \\ x - \frac{1}{\frac{z}{y} - \frac{\frac{t}{2}}{z}} \end{array} \]

FPCore

(FPCore (x y z t)
 :precision binary64
 (- x (/ 1.0 (- (/ z y) (/ (/ t 2.0) z)))))

C

double code(double x, double y, double z, double t) {
	return x - (1.0 / ((z / y) - ((t / 2.0) / z)));
}

Fortran

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 / ((z / y) - ((t / 2.0d0) / z)))
end function

Java

public static double code(double x, double y, double z, double t) {
	return x - (1.0 / ((z / y) - ((t / 2.0) / z)));
}

Python

def code(x, y, z, t):
	return x - (1.0 / ((z / y) - ((t / 2.0) / z)))

Julia

function code(x, y, z, t)
	return Float64(x - Float64(1.0 / Float64(Float64(z / y) - Float64(Float64(t / 2.0) / z))))
end

MATLAB

function tmp = code(x, y, z, t)
	tmp = x - (1.0 / ((z / y) - ((t / 2.0) / z)));
end

Wolfram

code[x_, y_, z_, t_] := N[(x - N[(1.0 / N[(N[(z / y), $MachinePrecision] - N[(N[(t / 2.0), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Reproduce

herbie shell --seed 2024321 
(FPCore (x y z t)
  :name "Numeric.AD.Rank1.Halley:findZero from ad-4.2.4"
  :precision binary64
  :pre (TRUE)

  :alt
  (! :herbie-platform default (- x (/ 1 (- (/ z y) (/ (/ t 2) z)))))

  (- x (/ (* (* y 2.0) z) (- (* (* z 2.0) z) (* y t)))))