Midpoint on a great circle

Percentage Accurate: 98.5% → 99.6%

Time: 13.6s

Alternatives: 19

Speedup: 0.7×

Specification

Math

\[\begin{array}{l} \\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  lambda1
  (atan2
   (* (cos phi2) (sin (- lambda1 lambda2)))
   (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), (Math.cos(phi1) + (Math.cos(phi2) * Math.cos((lambda1 - lambda2)))));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.cos(phi1) + (math.cos(phi2) * math.cos((lambda1 - lambda2)))))

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(cos(phi1) + Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Initial Program: 98.5% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  lambda1
  (atan2
   (* (cos phi2) (sin (- lambda1 lambda2)))
   (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), (Math.cos(phi1) + (Math.cos(phi2) * Math.cos((lambda1 - lambda2)))));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (math.cos(phi1) + (math.cos(phi2) * math.cos((lambda1 - lambda2)))))

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(cos(phi1) + Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Alternative 1: 99.6% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \cos \lambda_1 \cdot \sin \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  lambda1
  (atan2
   (*
    (cos phi2)
    (- (* (cos lambda2) (sin lambda1)) (* (cos lambda1) (sin lambda2))))
   (+
    (cos phi1)
    (*
     (cos phi2)
     (fma (sin lambda2) (sin lambda1) (* (cos lambda2) (cos lambda1))))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + atan2((cos(phi2) * ((cos(lambda2) * sin(lambda1)) - (cos(lambda1) * sin(lambda2)))), (cos(phi1) + (cos(phi2) * fma(sin(lambda2), sin(lambda1), (cos(lambda2) * cos(lambda1))))));
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(lambda1 + atan(Float64(cos(phi2) * Float64(Float64(cos(lambda2) * sin(lambda1)) - Float64(cos(lambda1) * sin(lambda2)))), Float64(cos(phi1) + Float64(cos(phi2) * fma(sin(lambda2), sin(lambda1), Float64(cos(lambda2) * cos(lambda1)))))))
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[lambda1], $MachinePrecision] * N[Sin[lambda2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[(N[Sin[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision] + N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
2. Step-by-step derivation

   1. lift-sin.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\sin \left(\lambda_1 - \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \color{blue}{\left(\lambda_1 - \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. sin-diff N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. sub-to-mult N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   5. lower-special-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   6. lower-special--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\left(1 - \frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}\right)} \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. lower-special-/.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \color{blue}{\frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   8. *-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\color{blue}{\sin \lambda_2 \cdot \cos \lambda_1}}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   9. lower-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\color{blue}{\sin \lambda_2 \cdot \cos \lambda_1}}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   10. lower-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\color{blue}{\sin \lambda_2} \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   11. lower-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \color{blue}{\cos \lambda_1}}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   12. lower-*.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\color{blue}{\sin \lambda_1 \cdot \cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   13. lower-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\color{blue}{\sin \lambda_1} \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   14. lower-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \color{blue}{\cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   15. lower-*.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2\right)}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   16. lower-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\color{blue}{\sin \lambda_1} \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   17. lower-cos.f64 98.6

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \color{blue}{\cos \lambda_2}\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

3. Applied rewrites 98.6%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
4. Step-by-step derivation

   1. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

   3. cos-diff N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \color{blue}{\left(\cos \lambda_1 \cdot \cos \lambda_2 + \sin \lambda_1 \cdot \sin \lambda_2\right)}} \]

   4. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\color{blue}{\cos \lambda_1} \cdot \cos \lambda_2 + \sin \lambda_1 \cdot \sin \lambda_2\right)} \]

   5. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\cos \lambda_1 \cdot \color{blue}{\cos \lambda_2} + \sin \lambda_1 \cdot \sin \lambda_2\right)} \]

   6. *-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\color{blue}{\cos \lambda_2 \cdot \cos \lambda_1} + \sin \lambda_1 \cdot \sin \lambda_2\right)} \]

   7. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\color{blue}{\cos \lambda_2} \cdot \cos \lambda_1 + \sin \lambda_1 \cdot \sin \lambda_2\right)} \]

   8. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \color{blue}{\cos \lambda_1} + \sin \lambda_1 \cdot \sin \lambda_2\right)} \]

   9. +-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \color{blue}{\left(\sin \lambda_1 \cdot \sin \lambda_2 + \cos \lambda_2 \cdot \cos \lambda_1\right)}} \]

   10. lift-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\color{blue}{\sin \lambda_1} \cdot \sin \lambda_2 + \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   11. lift-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \color{blue}{\sin \lambda_2} + \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   12. *-commutative N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \left(\color{blue}{\sin \lambda_2 \cdot \sin \lambda_1} + \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   13. lower-fma.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \color{blue}{\mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)}} \]

   14. lift-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \color{blue}{\cos \lambda_2} \cdot \cos \lambda_1\right)} \]

   15. lift-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \color{blue}{\cos \lambda_1}\right)} \]

   16. lower-*.f64 99.6

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \color{blue}{\cos \lambda_2 \cdot \cos \lambda_1}\right)} \]

5. Applied rewrites 99.6%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \color{blue}{\mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)}} \]
6. Step-by-step derivation

   1. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   2. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right)} \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   3. lift-/.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \color{blue}{\frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   4. sub-to-mult-rev N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2 - \sin \lambda_2 \cdot \cos \lambda_1\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   5. lower--.f64 99.6

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2 - \sin \lambda_2 \cdot \cos \lambda_1\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   6. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\sin \lambda_1 \cdot \cos \lambda_2} - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   7. *-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\cos \lambda_2 \cdot \sin \lambda_1} - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   8. lower-*.f64 99.6

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\cos \lambda_2 \cdot \sin \lambda_1} - \sin \lambda_2 \cdot \cos \lambda_1\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   9. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \color{blue}{\sin \lambda_2 \cdot \cos \lambda_1}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   10. *-commutative N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \color{blue}{\cos \lambda_1 \cdot \sin \lambda_2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

   11. lower-*.f64 99.6

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\cos \lambda_2 \cdot \sin \lambda_1 - \color{blue}{\cos \lambda_1 \cdot \sin \lambda_2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]

7. Applied rewrites 99.6%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\cos \lambda_2 \cdot \sin \lambda_1 - \cos \lambda_1 \cdot \sin \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \mathsf{fma}\left(\sin \lambda_2, \sin \lambda_1, \cos \lambda_2 \cdot \cos \lambda_1\right)} \]
8. Add Preprocessing

Alternative 2: 98.6% accurate, 0.6× speedup

Math

\[\begin{array}{l} \\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \frac{\left(\cos \lambda_2 \cdot \sin \lambda_1\right) \cdot 2 - \left(\sin \left(\lambda_2 + \lambda_1\right) + \sin \left(\lambda_2 - \lambda_1\right)\right)}{2}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  lambda1
  (atan2
   (*
    (cos phi2)
    (/
     (-
      (* (* (cos lambda2) (sin lambda1)) 2.0)
      (+ (sin (+ lambda2 lambda1)) (sin (- lambda2 lambda1))))
     2.0))
   (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + atan2((cos(phi2) * ((((cos(lambda2) * sin(lambda1)) * 2.0) - (sin((lambda2 + lambda1)) + sin((lambda2 - lambda1)))) / 2.0)), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    code = lambda1 + atan2((cos(phi2) * ((((cos(lambda2) * sin(lambda1)) * 2.0d0) - (sin((lambda2 + lambda1)) + sin((lambda2 - lambda1)))) / 2.0d0)), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + Math.atan2((Math.cos(phi2) * ((((Math.cos(lambda2) * Math.sin(lambda1)) * 2.0) - (Math.sin((lambda2 + lambda1)) + Math.sin((lambda2 - lambda1)))) / 2.0)), (Math.cos(phi1) + (Math.cos(phi2) * Math.cos((lambda1 - lambda2)))));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	return lambda1 + math.atan2((math.cos(phi2) * ((((math.cos(lambda2) * math.sin(lambda1)) * 2.0) - (math.sin((lambda2 + lambda1)) + math.sin((lambda2 - lambda1)))) / 2.0)), (math.cos(phi1) + (math.cos(phi2) * math.cos((lambda1 - lambda2)))))

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(lambda1 + atan(Float64(cos(phi2) * Float64(Float64(Float64(Float64(cos(lambda2) * sin(lambda1)) * 2.0) - Float64(sin(Float64(lambda2 + lambda1)) + sin(Float64(lambda2 - lambda1)))) / 2.0)), Float64(cos(phi1) + Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	tmp = lambda1 + atan2((cos(phi2) * ((((cos(lambda2) * sin(lambda1)) * 2.0) - (sin((lambda2 + lambda1)) + sin((lambda2 - lambda1)))) / 2.0)), (cos(phi1) + (cos(phi2) * cos((lambda1 - lambda2)))));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[(N[(N[(N[(N[Cos[lambda2], $MachinePrecision] * N[Sin[lambda1], $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] - N[(N[Sin[N[(lambda2 + lambda1), $MachinePrecision]], $MachinePrecision] + N[Sin[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
2. Step-by-step derivation

   1. lift-sin.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\sin \left(\lambda_1 - \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \color{blue}{\left(\lambda_1 - \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. sin-diff N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2 - \cos \lambda_1 \cdot \sin \lambda_2\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. sub-to-mult N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   5. lower-special-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   6. lower-special--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\left(1 - \frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}\right)} \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. lower-special-/.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \color{blue}{\frac{\cos \lambda_1 \cdot \sin \lambda_2}{\sin \lambda_1 \cdot \cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   8. *-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\color{blue}{\sin \lambda_2 \cdot \cos \lambda_1}}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   9. lower-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\color{blue}{\sin \lambda_2 \cdot \cos \lambda_1}}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   10. lower-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\color{blue}{\sin \lambda_2} \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   11. lower-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \color{blue}{\cos \lambda_1}}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   12. lower-*.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\color{blue}{\sin \lambda_1 \cdot \cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   13. lower-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\color{blue}{\sin \lambda_1} \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   14. lower-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \color{blue}{\cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   15. lower-*.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2\right)}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   16. lower-sin.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\color{blue}{\sin \lambda_1} \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   17. lower-cos.f64 98.6

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \color{blue}{\cos \lambda_2}\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

3. Applied rewrites 98.6%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
4. Step-by-step derivation

   1. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\color{blue}{\left(1 - \frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}\right)} \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lift-/.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\left(1 - \color{blue}{\frac{\sin \lambda_2 \cdot \cos \lambda_1}{\sin \lambda_1 \cdot \cos \lambda_2}}\right) \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2\right)\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. sub-to-mult-rev N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\left(\sin \lambda_1 \cdot \cos \lambda_2 - \sin \lambda_2 \cdot \cos \lambda_1\right)}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   5. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2 - \color{blue}{\sin \lambda_2 \cdot \cos \lambda_1}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   6. lift-sin.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2 - \color{blue}{\sin \lambda_2} \cdot \cos \lambda_1\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2 - \sin \lambda_2 \cdot \color{blue}{\cos \lambda_1}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   8. sin-cos-mult N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \left(\sin \lambda_1 \cdot \cos \lambda_2 - \color{blue}{\frac{\sin \left(\lambda_2 - \lambda_1\right) + \sin \left(\lambda_2 + \lambda_1\right)}{2}}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   9. sub-to-fraction N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\frac{\left(\sin \lambda_1 \cdot \cos \lambda_2\right) \cdot 2 - \left(\sin \left(\lambda_2 - \lambda_1\right) + \sin \left(\lambda_2 + \lambda_1\right)\right)}{2}}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   10. lower-/.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\frac{\left(\sin \lambda_1 \cdot \cos \lambda_2\right) \cdot 2 - \left(\sin \left(\lambda_2 - \lambda_1\right) + \sin \left(\lambda_2 + \lambda_1\right)\right)}{2}}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

5. Applied rewrites 98.6%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \color{blue}{\frac{\left(\cos \lambda_2 \cdot \sin \lambda_1\right) \cdot 2 - \left(\sin \left(\lambda_2 + \lambda_1\right) + \sin \left(\lambda_2 - \lambda_1\right)\right)}{2}}}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
6. Add Preprocessing

Alternative 3: 98.5% accurate, 0.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{t\_0}}{\cos \phi_2}\right) \cdot \left(t\_0 \cdot \cos \phi_2\right)} \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (cos (- lambda2 lambda1))))
   (+
    lambda1
    (atan2
     (* (cos phi2) (sin (- lambda1 lambda2)))
     (* (+ 1.0 (/ (/ (cos phi1) t_0) (cos phi2))) (* t_0 (cos phi2)))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = cos((lambda2 - lambda1));
	return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), ((1.0 + ((cos(phi1) / t_0) / cos(phi2))) * (t_0 * cos(phi2))));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    real(8) :: t_0
    t_0 = cos((lambda2 - lambda1))
    code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), ((1.0d0 + ((cos(phi1) / t_0) / cos(phi2))) * (t_0 * cos(phi2))))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = Math.cos((lambda2 - lambda1));
	return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), ((1.0 + ((Math.cos(phi1) / t_0) / Math.cos(phi2))) * (t_0 * Math.cos(phi2))));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	t_0 = math.cos((lambda2 - lambda1))
	return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), ((1.0 + ((math.cos(phi1) / t_0) / math.cos(phi2))) * (t_0 * math.cos(phi2))))

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = cos(Float64(lambda2 - lambda1))
	return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(1.0 + Float64(Float64(cos(phi1) / t_0) / cos(phi2))) * Float64(t_0 * cos(phi2)))))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	t_0 = cos((lambda2 - lambda1));
	tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), ((1.0 + ((cos(phi1) / t_0) / cos(phi2))) * (t_0 * cos(phi2))));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(N[(N[Cos[phi1], $MachinePrecision] / t$95$0), $MachinePrecision] / N[Cos[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
2. Step-by-step derivation

   1. lift-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. +-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right) + \cos \phi_1}} \]

   3. sum-to-mult N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}} \]

   4. lower-special-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}} \]

   5. lower-special-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\right)} \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   6. lower-special-/.f64 98.5

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \color{blue}{\frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   7. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   8. *-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   9. lower-*.f64 98.5

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   10. lift-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   11. cos-neg-rev N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\mathsf{neg}\left(\left(\lambda_1 - \lambda_2\right)\right)\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   12. lower-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\mathsf{neg}\left(\left(\lambda_1 - \lambda_2\right)\right)\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   13. lift--.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\mathsf{neg}\left(\color{blue}{\left(\lambda_1 - \lambda_2\right)}\right)\right) \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   14. sub-negate-rev N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   15. lower--.f64 98.5

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   16. lift-*.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \color{blue}{\left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}} \]

   17. *-commutative N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \color{blue}{\left(\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\right)}} \]

   18. lower-*.f64 98.5

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \color{blue}{\left(\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\right)}} \]

3. Applied rewrites 98.5%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)}} \]
4. Step-by-step derivation

   1. lift-/.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \color{blue}{\frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   2. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   3. associate-/r* N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \color{blue}{\frac{\frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right)}}{\cos \phi_2}}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   4. lower-/.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \color{blue}{\frac{\frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right)}}{\cos \phi_2}}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   5. lift-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_2 - \lambda_1\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   6. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   7. sub-negate-rev N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\cos \color{blue}{\left(\mathsf{neg}\left(\left(\lambda_1 - \lambda_2\right)\right)\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   8. lift--.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\cos \left(\mathsf{neg}\left(\color{blue}{\left(\lambda_1 - \lambda_2\right)}\right)\right)}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   9. cos-neg-rev N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   10. lift-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   11. lower-/.f64 98.5

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\color{blue}{\frac{\cos \phi_1}{\cos \left(\lambda_1 - \lambda_2\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   12. lift-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   13. cos-neg-rev N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\color{blue}{\cos \left(\mathsf{neg}\left(\left(\lambda_1 - \lambda_2\right)\right)\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   14. lift--.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\cos \left(\mathsf{neg}\left(\color{blue}{\left(\lambda_1 - \lambda_2\right)}\right)\right)}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   15. sub-negate-rev N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   16. lift--.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

   17. lift-cos.f64 98.5

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_2 - \lambda_1\right)}}}{\cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]

5. Applied rewrites 98.5%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \color{blue}{\frac{\frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right)}}{\cos \phi_2}}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)} \]
6. Add Preprocessing

Alternative 4: 98.5% accurate, 0.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{t\_0}\right) \cdot t\_0} \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (* (cos (- lambda2 lambda1)) (cos phi2))))
   (+
    lambda1
    (atan2
     (* (cos phi2) (sin (- lambda1 lambda2)))
     (* (+ 1.0 (/ (cos phi1) t_0)) t_0)))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = cos((lambda2 - lambda1)) * cos(phi2);
	return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), ((1.0 + (cos(phi1) / t_0)) * t_0));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    real(8) :: t_0
    t_0 = cos((lambda2 - lambda1)) * cos(phi2)
    code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), ((1.0d0 + (cos(phi1) / t_0)) * t_0))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = Math.cos((lambda2 - lambda1)) * Math.cos(phi2);
	return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), ((1.0 + (Math.cos(phi1) / t_0)) * t_0));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	t_0 = math.cos((lambda2 - lambda1)) * math.cos(phi2)
	return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), ((1.0 + (math.cos(phi1) / t_0)) * t_0))

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = Float64(cos(Float64(lambda2 - lambda1)) * cos(phi2))
	return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(Float64(1.0 + Float64(cos(phi1) / t_0)) * t_0)))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	t_0 = cos((lambda2 - lambda1)) * cos(phi2);
	tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), ((1.0 + (cos(phi1) / t_0)) * t_0));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision]}, N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(N[Cos[phi1], $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
2. Step-by-step derivation

   1. lift-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. +-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right) + \cos \phi_1}} \]

   3. sum-to-mult N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}} \]

   4. lower-special-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}} \]

   5. lower-special-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\right)} \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   6. lower-special-/.f64 98.5

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \color{blue}{\frac{\cos \phi_1}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   7. lift-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   8. *-commutative N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   9. lower-*.f64 98.5

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   10. lift-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\lambda_1 - \lambda_2\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   11. cos-neg-rev N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\mathsf{neg}\left(\left(\lambda_1 - \lambda_2\right)\right)\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   12. lower-cos.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\color{blue}{\cos \left(\mathsf{neg}\left(\left(\lambda_1 - \lambda_2\right)\right)\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   13. lift--.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\mathsf{neg}\left(\color{blue}{\left(\lambda_1 - \lambda_2\right)}\right)\right) \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   14. sub-negate-rev N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   15. lower--.f64 98.5

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \color{blue}{\left(\lambda_2 - \lambda_1\right)} \cdot \cos \phi_2}\right) \cdot \left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)} \]

   16. lift-*.f64 N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \color{blue}{\left(\cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)\right)}} \]

   17. *-commutative N/A

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \color{blue}{\left(\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\right)}} \]

   18. lower-*.f64 98.5

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \color{blue}{\left(\cos \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2\right)}} \]

3. Applied rewrites 98.5%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\left(1 + \frac{\cos \phi_1}{\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2}\right) \cdot \left(\cos \left(\lambda_2 - \lambda_1\right) \cdot \cos \phi_2\right)}} \]
4. Add Preprocessing

Alternative 5: 98.5% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  (atan2
   (* (sin (- lambda1 lambda2)) (cos phi2))
   (fma (cos (- lambda2 lambda1)) (cos phi2) (cos phi1)))
  lambda1))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return atan2((sin((lambda1 - lambda2)) * cos(phi2)), fma(cos((lambda2 - lambda1)), cos(phi2), cos(phi1))) + lambda1;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), fma(cos(Float64(lambda2 - lambda1)), cos(phi2), cos(phi1))) + lambda1)
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
2. Step-by-step derivation

   1. lift-+.f64 N/A

      \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. +-commutative N/A

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. lower-+.f64 98.5

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

3. Applied rewrites 98.5%

   \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]
4. Add Preprocessing

Alternative 6: 97.8% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  (atan2
   (* (sin (- lambda1 lambda2)) (cos phi2))
   (fma (cos lambda2) (cos phi2) (cos phi1)))
  lambda1))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return atan2((sin((lambda1 - lambda2)) * cos(phi2)), fma(cos(lambda2), cos(phi2), cos(phi1))) + lambda1;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), fma(cos(lambda2), cos(phi2), cos(phi1))) + lambda1)
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
2. Step-by-step derivation

   1. lift-+.f64 N/A

      \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. +-commutative N/A

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. lower-+.f64 98.5

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

3. Applied rewrites 98.5%

   \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

4. Taylor expanded in lambda1 around 0

   \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \color{blue}{\lambda_2}, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
5. Step-by-step derivation

6. Applied rewrites 97.8%

   \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \color{blue}{\lambda_2}, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
7. Add Preprocessing

Alternative 7: 93.1% accurate, 0.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;\cos \phi_1 \leq 0.999996:\\ \;\;\;\;\tan^{-1}_* \frac{\sin \left(-\lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \cos \phi_1\right)} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{2}\right)} + \lambda_1\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (if (<= (cos phi1) 0.999996)
   (+
    (atan2
     (* (sin (- lambda2)) (cos phi2))
     (fma (cos lambda2) (cos phi2) (cos phi1)))
    lambda1)
   (+
    (atan2
     (* (sin (- lambda1 lambda2)) (cos phi2))
     (fma
      (cos (- lambda2 lambda1))
      (cos phi2)
      (+ 1.0 (* -0.5 (pow phi1 2.0)))))
    lambda1)))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double tmp;
	if (cos(phi1) <= 0.999996) {
		tmp = atan2((sin(-lambda2) * cos(phi2)), fma(cos(lambda2), cos(phi2), cos(phi1))) + lambda1;
	} else {
		tmp = atan2((sin((lambda1 - lambda2)) * cos(phi2)), fma(cos((lambda2 - lambda1)), cos(phi2), (1.0 + (-0.5 * pow(phi1, 2.0))))) + lambda1;
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	tmp = 0.0
	if (cos(phi1) <= 0.999996)
		tmp = Float64(atan(Float64(sin(Float64(-lambda2)) * cos(phi2)), fma(cos(lambda2), cos(phi2), cos(phi1))) + lambda1);
	else
		tmp = Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), fma(cos(Float64(lambda2 - lambda1)), cos(phi2), Float64(1.0 + Float64(-0.5 * (phi1 ^ 2.0))))) + lambda1);
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.999996], N[(N[ArcTan[N[(N[Sin[(-lambda2)], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(1.0 + N[(-0.5 * N[Power[phi1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if (cos.f64 phi1) < 0.999995999999999996

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in lambda1 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \color{blue}{\lambda_2}, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

   6. Applied rewrites 97.8%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \color{blue}{\lambda_2}, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in lambda1 around 0

      \[\leadsto \tan^{-1}_* \frac{\color{blue}{\sin \left(\mathsf{neg}\left(\lambda_2\right)\right)} \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-sin.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\mathsf{neg}\left(\lambda_2\right)\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-neg.f64 87.9

         \[\leadsto \tan^{-1}_* \frac{\sin \left(-\lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 87.9%

      \[\leadsto \tan^{-1}_* \frac{\color{blue}{\sin \left(-\lambda_2\right)} \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   if 0.999995999999999996 < (cos.f64 phi1)

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi1 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \color{blue}{1 + \frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_1}^{2}}\right)} + \lambda_1 \]

      3. lower-pow.f64 79.8

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{\color{blue}{2}}\right)} + \lambda_1 \]

   6. Applied rewrites 79.8%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \color{blue}{1 + -0.5 \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 8: 88.3% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\ t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\ \mathbf{if}\;\phi_2 \leq 44000:\\ \;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot t\_0} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot \cos \phi_2}{\mathsf{fma}\left(t\_0, \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{2}\right)} + \lambda_1\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (cos (- lambda2 lambda1))) (t_1 (sin (- lambda1 lambda2))))
   (if (<= phi2 44000.0)
     (+
      (atan2
       (* t_1 (+ 1.0 (* -0.5 (pow phi2 2.0))))
       (+ (cos phi1) (* (cos phi2) t_0)))
      lambda1)
     (+
      (atan2
       (* t_1 (cos phi2))
       (fma t_0 (cos phi2) (+ 1.0 (* -0.5 (pow phi1 2.0)))))
      lambda1))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = cos((lambda2 - lambda1));
	double t_1 = sin((lambda1 - lambda2));
	double tmp;
	if (phi2 <= 44000.0) {
		tmp = atan2((t_1 * (1.0 + (-0.5 * pow(phi2, 2.0)))), (cos(phi1) + (cos(phi2) * t_0))) + lambda1;
	} else {
		tmp = atan2((t_1 * cos(phi2)), fma(t_0, cos(phi2), (1.0 + (-0.5 * pow(phi1, 2.0))))) + lambda1;
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = cos(Float64(lambda2 - lambda1))
	t_1 = sin(Float64(lambda1 - lambda2))
	tmp = 0.0
	if (phi2 <= 44000.0)
		tmp = Float64(atan(Float64(t_1 * Float64(1.0 + Float64(-0.5 * (phi2 ^ 2.0)))), Float64(cos(phi1) + Float64(cos(phi2) * t_0))) + lambda1);
	else
		tmp = Float64(atan(Float64(t_1 * cos(phi2)), fma(t_0, cos(phi2), Float64(1.0 + Float64(-0.5 * (phi1 ^ 2.0))))) + lambda1);
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, 44000.0], N[(N[ArcTan[N[(t$95$1 * N[(1.0 + N[(-0.5 * N[Power[phi2, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[phi1], $MachinePrecision] + N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Cos[phi2], $MachinePrecision] + N[(1.0 + N[(-0.5 * N[Power[phi1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]]

Derivation

1. Split input into 2 regimes

2. if phi2 < 44000

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 76.3

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   6. Applied rewrites 76.3%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + \frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 77.0

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 77.0%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + -0.5 \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

   10. Taylor expanded in lambda1 around inf

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\color{blue}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_2 - \lambda_1\right)}} + \lambda_1 \]
   11. Step-by-step derivation

       1. lower-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \color{blue}{\cos \phi_2 \cdot \cos \left(\lambda_2 - \lambda_1\right)}} + \lambda_1 \]

       2. lower-cos.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \color{blue}{\cos \phi_2} \cdot \cos \left(\lambda_2 - \lambda_1\right)} + \lambda_1 \]

       3. lower-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \color{blue}{\cos \left(\lambda_2 - \lambda_1\right)}} + \lambda_1 \]

       4. lower-cos.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \color{blue}{\left(\lambda_2 - \lambda_1\right)}} + \lambda_1 \]

       5. lower-cos.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_2 - \lambda_1\right)} + \lambda_1 \]

       6. lower--.f64 76.3

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_2 - \lambda_1\right)} + \lambda_1 \]

   12. Applied rewrites 76.3%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\color{blue}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_2 - \lambda_1\right)}} + \lambda_1 \]

   if 44000 < phi2

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi1 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \color{blue}{1 + \frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_1}^{2}}\right)} + \lambda_1 \]

      3. lower-pow.f64 79.8

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{\color{blue}{2}}\right)} + \lambda_1 \]

   6. Applied rewrites 79.8%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \color{blue}{1 + -0.5 \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 9: 83.8% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\ t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\ t_2 := \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)\\ \mathbf{if}\;\phi_2 \leq 44000:\\ \;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot t\_2}{\mathsf{fma}\left(t\_0, t\_2, \cos \phi_1\right)} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot \cos \phi_2}{\mathsf{fma}\left(t\_0, \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{2}\right)} + \lambda_1\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (cos (- lambda2 lambda1)))
        (t_1 (sin (- lambda1 lambda2)))
        (t_2 (fma (* phi2 phi2) -0.5 1.0)))
   (if (<= phi2 44000.0)
     (+ (atan2 (* t_1 t_2) (fma t_0 t_2 (cos phi1))) lambda1)
     (+
      (atan2
       (* t_1 (cos phi2))
       (fma t_0 (cos phi2) (+ 1.0 (* -0.5 (pow phi1 2.0)))))
      lambda1))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = cos((lambda2 - lambda1));
	double t_1 = sin((lambda1 - lambda2));
	double t_2 = fma((phi2 * phi2), -0.5, 1.0);
	double tmp;
	if (phi2 <= 44000.0) {
		tmp = atan2((t_1 * t_2), fma(t_0, t_2, cos(phi1))) + lambda1;
	} else {
		tmp = atan2((t_1 * cos(phi2)), fma(t_0, cos(phi2), (1.0 + (-0.5 * pow(phi1, 2.0))))) + lambda1;
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = cos(Float64(lambda2 - lambda1))
	t_1 = sin(Float64(lambda1 - lambda2))
	t_2 = fma(Float64(phi2 * phi2), -0.5, 1.0)
	tmp = 0.0
	if (phi2 <= 44000.0)
		tmp = Float64(atan(Float64(t_1 * t_2), fma(t_0, t_2, cos(phi1))) + lambda1);
	else
		tmp = Float64(atan(Float64(t_1 * cos(phi2)), fma(t_0, cos(phi2), Float64(1.0 + Float64(-0.5 * (phi1 ^ 2.0))))) + lambda1);
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(phi2 * phi2), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, If[LessEqual[phi2, 44000.0], N[(N[ArcTan[N[(t$95$1 * t$95$2), $MachinePrecision] / N[(t$95$0 * t$95$2 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * N[Cos[phi2], $MachinePrecision] + N[(1.0 + N[(-0.5 * N[Power[phi1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]]]

Derivation

1. Split input into 2 regimes

2. if phi2 < 44000

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 76.3

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   6. Applied rewrites 76.3%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + \frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 77.0

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 77.0%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + -0.5 \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   10. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left({\phi_2}^{2} \cdot \frac{-1}{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

   11. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]
   12. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), {\phi_2}^{2} \cdot \frac{-1}{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right), \cos \phi_1\right)} + \lambda_1 \]

   13. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}, \cos \phi_1\right)} + \lambda_1 \]

   if 44000 < phi2

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi1 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \color{blue}{1 + \frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_1}^{2}}\right)} + \lambda_1 \]

      3. lower-pow.f64 79.8

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{\color{blue}{2}}\right)} + \lambda_1 \]

   6. Applied rewrites 79.8%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \color{blue}{1 + -0.5 \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 10: 83.7% accurate, 1.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\ t_1 := \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)\\ \mathbf{if}\;\phi_2 \leq 44000:\\ \;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot t\_1}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), t\_1, \cos \phi_1\right)} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{2}\right)} + \lambda_1\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (sin (- lambda1 lambda2))) (t_1 (fma (* phi2 phi2) -0.5 1.0)))
   (if (<= phi2 44000.0)
     (+
      (atan2 (* t_0 t_1) (fma (cos (- lambda2 lambda1)) t_1 (cos phi1)))
      lambda1)
     (+
      (atan2
       (* t_0 (cos phi2))
       (fma (cos lambda2) (cos phi2) (+ 1.0 (* -0.5 (pow phi1 2.0)))))
      lambda1))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = sin((lambda1 - lambda2));
	double t_1 = fma((phi2 * phi2), -0.5, 1.0);
	double tmp;
	if (phi2 <= 44000.0) {
		tmp = atan2((t_0 * t_1), fma(cos((lambda2 - lambda1)), t_1, cos(phi1))) + lambda1;
	} else {
		tmp = atan2((t_0 * cos(phi2)), fma(cos(lambda2), cos(phi2), (1.0 + (-0.5 * pow(phi1, 2.0))))) + lambda1;
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = sin(Float64(lambda1 - lambda2))
	t_1 = fma(Float64(phi2 * phi2), -0.5, 1.0)
	tmp = 0.0
	if (phi2 <= 44000.0)
		tmp = Float64(atan(Float64(t_0 * t_1), fma(cos(Float64(lambda2 - lambda1)), t_1, cos(phi1))) + lambda1);
	else
		tmp = Float64(atan(Float64(t_0 * cos(phi2)), fma(cos(lambda2), cos(phi2), Float64(1.0 + Float64(-0.5 * (phi1 ^ 2.0))))) + lambda1);
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(phi2 * phi2), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, If[LessEqual[phi2, 44000.0], N[(N[ArcTan[N[(t$95$0 * t$95$1), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * t$95$1 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[ArcTan[N[(t$95$0 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[lambda2], $MachinePrecision] * N[Cos[phi2], $MachinePrecision] + N[(1.0 + N[(-0.5 * N[Power[phi1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]]]]

Derivation

1. Split input into 2 regimes

2. if phi2 < 44000

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 76.3

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   6. Applied rewrites 76.3%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + \frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 77.0

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 77.0%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + -0.5 \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   10. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left({\phi_2}^{2} \cdot \frac{-1}{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

   11. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]
   12. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), {\phi_2}^{2} \cdot \frac{-1}{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right), \cos \phi_1\right)} + \lambda_1 \]

   13. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}, \cos \phi_1\right)} + \lambda_1 \]

   if 44000 < phi2

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in lambda1 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \color{blue}{\lambda_2}, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

   6. Applied rewrites 97.8%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \color{blue}{\lambda_2}, \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi1 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \color{blue}{1 + \frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_1}^{2}}\right)} + \lambda_1 \]

      3. lower-pow.f64 79.5

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, 1 + -0.5 \cdot {\phi_1}^{\color{blue}{2}}\right)} + \lambda_1 \]

   9. Applied rewrites 79.5%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \lambda_2, \cos \phi_2, \color{blue}{1 + -0.5 \cdot {\phi_1}^{2}}\right)} + \lambda_1 \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 11: 83.3% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)\\ t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\ \mathbf{if}\;\cos \phi_1 \leq 0.9999:\\ \;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot t\_0}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), t\_0, \cos \phi_1\right)} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot t\_1}{1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (fma (* phi2 phi2) -0.5 1.0)) (t_1 (sin (- lambda1 lambda2))))
   (if (<= (cos phi1) 0.9999)
     (+
      (atan2 (* t_1 t_0) (fma (cos (- lambda2 lambda1)) t_0 (cos phi1)))
      lambda1)
     (+
      lambda1
      (atan2
       (* (cos phi2) t_1)
       (+ 1.0 (* (cos phi2) (cos (- lambda1 lambda2)))))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = fma((phi2 * phi2), -0.5, 1.0);
	double t_1 = sin((lambda1 - lambda2));
	double tmp;
	if (cos(phi1) <= 0.9999) {
		tmp = atan2((t_1 * t_0), fma(cos((lambda2 - lambda1)), t_0, cos(phi1))) + lambda1;
	} else {
		tmp = lambda1 + atan2((cos(phi2) * t_1), (1.0 + (cos(phi2) * cos((lambda1 - lambda2)))));
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = fma(Float64(phi2 * phi2), -0.5, 1.0)
	t_1 = sin(Float64(lambda1 - lambda2))
	tmp = 0.0
	if (cos(phi1) <= 0.9999)
		tmp = Float64(atan(Float64(t_1 * t_0), fma(cos(Float64(lambda2 - lambda1)), t_0, cos(phi1))) + lambda1);
	else
		tmp = Float64(lambda1 + atan(Float64(cos(phi2) * t_1), Float64(1.0 + Float64(cos(phi2) * cos(Float64(lambda1 - lambda2))))));
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[(N[(phi2 * phi2), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[Cos[phi1], $MachinePrecision], 0.9999], N[(N[ArcTan[N[(t$95$1 * t$95$0), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * t$95$0 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$1), $MachinePrecision] / N[(1.0 + N[(N[Cos[phi2], $MachinePrecision] * N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]

Derivation

1. Split input into 2 regimes

2. if (cos.f64 phi1) < 0.99990000000000001

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 76.3

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   6. Applied rewrites 76.3%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + \frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 77.0

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 77.0%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + -0.5 \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   10. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left({\phi_2}^{2} \cdot \frac{-1}{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

   11. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]
   12. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), {\phi_2}^{2} \cdot \frac{-1}{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right), \cos \phi_1\right)} + \lambda_1 \]

   13. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}, \cos \phi_1\right)} + \lambda_1 \]

   if 0.99990000000000001 < (cos.f64 phi1)

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. Taylor expanded in phi1 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{1} + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   3. Step-by-step derivation

   4. Applied rewrites 78.2%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{1} + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 12: 79.0% accurate, 1.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\ t_1 := \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)\\ \mathbf{if}\;\phi_2 \leq 45000:\\ \;\;\;\;\tan^{-1}_* \frac{t\_0 \cdot t\_1}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), t\_1, \cos \phi_1\right)} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot t\_0}{\left(1 + -0.5 \cdot {\phi_1}^{2}\right) + \cos \left(\lambda_1 - \lambda_2\right)}\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (sin (- lambda1 lambda2))) (t_1 (fma (* phi2 phi2) -0.5 1.0)))
   (if (<= phi2 45000.0)
     (+
      (atan2 (* t_0 t_1) (fma (cos (- lambda2 lambda1)) t_1 (cos phi1)))
      lambda1)
     (+
      lambda1
      (atan2
       (* (cos phi2) t_0)
       (+ (+ 1.0 (* -0.5 (pow phi1 2.0))) (cos (- lambda1 lambda2))))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = sin((lambda1 - lambda2));
	double t_1 = fma((phi2 * phi2), -0.5, 1.0);
	double tmp;
	if (phi2 <= 45000.0) {
		tmp = atan2((t_0 * t_1), fma(cos((lambda2 - lambda1)), t_1, cos(phi1))) + lambda1;
	} else {
		tmp = lambda1 + atan2((cos(phi2) * t_0), ((1.0 + (-0.5 * pow(phi1, 2.0))) + cos((lambda1 - lambda2))));
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = sin(Float64(lambda1 - lambda2))
	t_1 = fma(Float64(phi2 * phi2), -0.5, 1.0)
	tmp = 0.0
	if (phi2 <= 45000.0)
		tmp = Float64(atan(Float64(t_0 * t_1), fma(cos(Float64(lambda2 - lambda1)), t_1, cos(phi1))) + lambda1);
	else
		tmp = Float64(lambda1 + atan(Float64(cos(phi2) * t_0), Float64(Float64(1.0 + Float64(-0.5 * (phi1 ^ 2.0))) + cos(Float64(lambda1 - lambda2)))));
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(phi2 * phi2), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, If[LessEqual[phi2, 45000.0], N[(N[ArcTan[N[(t$95$0 * t$95$1), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] * t$95$1 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision] / N[(N[(1.0 + N[(-0.5 * N[Power[phi1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]

Derivation

1. Split input into 2 regimes

2. if phi2 < 45000

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 76.3

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   6. Applied rewrites 76.3%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + \frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 77.0

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 77.0%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + -0.5 \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   10. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left({\phi_2}^{2} \cdot \frac{-1}{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

   11. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]
   12. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), {\phi_2}^{2} \cdot \frac{-1}{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right), \cos \phi_1\right)} + \lambda_1 \]

   13. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}, \cos \phi_1\right)} + \lambda_1 \]

   if 45000 < phi2

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. Taylor expanded in phi2 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
   3. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

      3. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      4. lower--.f64 78.4

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. Applied rewrites 78.4%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

   5. Taylor expanded in phi1 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{-1}{2} \cdot {\phi_1}^{2}\right) + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]
   6. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{-1}{2} \cdot {\phi_1}^{2}\right) + \cos \left(\lambda_1 - \color{blue}{\lambda_2}\right)} \]

      2. lower-*.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + \frac{-1}{2} \cdot {\phi_1}^{2}\right) + \cos \left(\lambda_1 - \lambda_2\right)} \]

      3. lower-pow.f64 69.9

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + -0.5 \cdot {\phi_1}^{2}\right) + \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. Applied rewrites 69.9%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\left(1 + -0.5 \cdot {\phi_1}^{2}\right) + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 13: 77.7% accurate, 1.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \cos \left(\lambda_2 - \lambda_1\right)\\ t_1 := \sin \left(\lambda_1 - \lambda_2\right)\\ t_2 := \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)\\ \mathbf{if}\;\phi_2 \leq 5.2 \cdot 10^{+23}:\\ \;\;\;\;\tan^{-1}_* \frac{t\_1 \cdot t\_2}{\mathsf{fma}\left(t\_0, t\_2, \cos \phi_1\right)} + \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\left(1 + \frac{\tan^{-1}_* \frac{t\_1 \cdot \cos \phi_2}{t\_0 - -1}}{\lambda_1}\right) \cdot \lambda_1\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (cos (- lambda2 lambda1)))
        (t_1 (sin (- lambda1 lambda2)))
        (t_2 (fma (* phi2 phi2) -0.5 1.0)))
   (if (<= phi2 5.2e+23)
     (+ (atan2 (* t_1 t_2) (fma t_0 t_2 (cos phi1))) lambda1)
     (* (+ 1.0 (/ (atan2 (* t_1 (cos phi2)) (- t_0 -1.0)) lambda1)) lambda1))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = cos((lambda2 - lambda1));
	double t_1 = sin((lambda1 - lambda2));
	double t_2 = fma((phi2 * phi2), -0.5, 1.0);
	double tmp;
	if (phi2 <= 5.2e+23) {
		tmp = atan2((t_1 * t_2), fma(t_0, t_2, cos(phi1))) + lambda1;
	} else {
		tmp = (1.0 + (atan2((t_1 * cos(phi2)), (t_0 - -1.0)) / lambda1)) * lambda1;
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = cos(Float64(lambda2 - lambda1))
	t_1 = sin(Float64(lambda1 - lambda2))
	t_2 = fma(Float64(phi2 * phi2), -0.5, 1.0)
	tmp = 0.0
	if (phi2 <= 5.2e+23)
		tmp = Float64(atan(Float64(t_1 * t_2), fma(t_0, t_2, cos(phi1))) + lambda1);
	else
		tmp = Float64(Float64(1.0 + Float64(atan(Float64(t_1 * cos(phi2)), Float64(t_0 - -1.0)) / lambda1)) * lambda1);
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(phi2 * phi2), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]}, If[LessEqual[phi2, 5.2e+23], N[(N[ArcTan[N[(t$95$1 * t$95$2), $MachinePrecision] / N[(t$95$0 * t$95$2 + N[Cos[phi1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision], N[(N[(1.0 + N[(N[ArcTan[N[(t$95$1 * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(t$95$0 - -1.0), $MachinePrecision]], $MachinePrecision] / lambda1), $MachinePrecision]), $MachinePrecision] * lambda1), $MachinePrecision]]]]]

Derivation

1. Split input into 2 regimes

2. if phi2 < 5.19999999999999983e23

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]
   2. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

      3. lower-+.f64 98.5

         \[\leadsto \color{blue}{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

   3. Applied rewrites 98.5%

      \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1} \]

   4. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 76.3

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   6. Applied rewrites 76.3%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \cos \phi_2, \cos \phi_1\right)} + \lambda_1 \]

   7. Taylor expanded in phi2 around 0

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + \frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   8. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      2. lower-*.f64 N/A

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

      3. lower-pow.f64 77.0

         \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}, \cos \phi_1\right)} + \lambda_1 \]

   9. Applied rewrites 77.0%

      \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + -0.5 \cdot {\phi_2}^{2}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{1 + -0.5 \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]
   10. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left(\frac{-1}{2} \cdot {\phi_2}^{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \left({\phi_2}^{2} \cdot \frac{-1}{2} + 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \frac{-1}{2} \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]

   11. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + -0.5 \cdot {\phi_2}^{2}, \cos \phi_1\right)} + \lambda_1 \]
   12. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), 1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}, \cos \phi_1\right)} + \lambda_1 \]

       2. +-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + \color{blue}{1}, \cos \phi_1\right)} + \lambda_1 \]

       3. lift-*.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \frac{-1}{2} \cdot {\phi_2}^{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       4. *-commutative N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), {\phi_2}^{2} \cdot \frac{-1}{2} + 1, \cos \phi_1\right)} + \lambda_1 \]

       5. lower-fma.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \color{blue}{-0.5}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       6. lift-pow.f64 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left({\phi_2}^{2}, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       7. unpow2 N/A

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, \frac{-1}{2}, 1\right), \cos \phi_1\right)} + \lambda_1 \]

       8. lower-*.f64 77.0

          \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right), \cos \phi_1\right)} + \lambda_1 \]

   13. Applied rewrites 77.0%

       \[\leadsto \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}{\mathsf{fma}\left(\cos \left(\lambda_2 - \lambda_1\right), \color{blue}{\mathsf{fma}\left(\phi_2 \cdot \phi_2, -0.5, 1\right)}, \cos \phi_1\right)} + \lambda_1 \]

   if 5.19999999999999983e23 < phi2

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. Taylor expanded in phi2 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
   3. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

      3. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      4. lower--.f64 78.4

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. Applied rewrites 78.4%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

   5. Taylor expanded in phi1 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
   6. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      3. lower--.f64 68.2

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. Applied rewrites 68.2%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
   8. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. sum-to-mult N/A

         \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

      3. lower-special-*.f64 N/A

         \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

   9. Applied rewrites 68.2%

      \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 14: 68.2% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \left(1 + \frac{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1 \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (*
  (+
   1.0
   (/
    (atan2
     (* (sin (- lambda1 lambda2)) (cos phi2))
     (- (cos (- lambda2 lambda1)) -1.0))
    lambda1))
  lambda1))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return (1.0 + (atan2((sin((lambda1 - lambda2)) * cos(phi2)), (cos((lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1;
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    code = (1.0d0 + (atan2((sin((lambda1 - lambda2)) * cos(phi2)), (cos((lambda2 - lambda1)) - (-1.0d0))) / lambda1)) * lambda1
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	return (1.0 + (Math.atan2((Math.sin((lambda1 - lambda2)) * Math.cos(phi2)), (Math.cos((lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1;
}

Python

def code(lambda1, lambda2, phi1, phi2):
	return (1.0 + (math.atan2((math.sin((lambda1 - lambda2)) * math.cos(phi2)), (math.cos((lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(Float64(1.0 + Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * cos(phi2)), Float64(cos(Float64(lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1)
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	tmp = (1.0 + (atan2((sin((lambda1 - lambda2)) * cos(phi2)), (cos((lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1;
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[(1.0 + N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[Cos[phi2], $MachinePrecision]), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision]], $MachinePrecision] / lambda1), $MachinePrecision]), $MachinePrecision] * lambda1), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

2. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
3. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

   3. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. lower--.f64 78.4

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

4. Applied rewrites 78.4%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

5. Taylor expanded in phi1 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
6. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower--.f64 68.2

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

7. Applied rewrites 68.2%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
8. Step-by-step derivation

   1. lift-+.f64 N/A

      \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. sum-to-mult N/A

      \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

   3. lower-special-*.f64 N/A

      \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

9. Applied rewrites 68.2%

   \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \cos \phi_2}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1} \]
10. Add Preprocessing

Alternative 15: 68.2% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  lambda1
  (atan2
   (* (cos phi2) (sin (- lambda1 lambda2)))
   (+ 1.0 (cos (- lambda1 lambda2))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (1.0 + cos((lambda1 - lambda2))));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (1.0d0 + cos((lambda1 - lambda2))))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), (1.0 + Math.cos((lambda1 - lambda2))));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (1.0 + math.cos((lambda1 - lambda2))))

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(1.0 + cos(Float64(lambda1 - lambda2)))))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (1.0 + cos((lambda1 - lambda2))));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[Cos[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

2. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
3. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

   3. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. lower--.f64 78.4

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

4. Applied rewrites 78.4%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

5. Taylor expanded in phi1 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
6. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower--.f64 68.2

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

7. Applied rewrites 68.2%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
8. Add Preprocessing

Alternative 16: 68.1% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \sin \left(\lambda_1 - \lambda_2\right)\\ \mathbf{if}\;\phi_2 \leq 195000:\\ \;\;\;\;\left(1 + \frac{\tan^{-1}_* \frac{t\_0 \cdot \mathsf{fma}\left(-0.5 \cdot \phi_2, \phi_2, 1\right)}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1\\ \mathbf{else}:\\ \;\;\;\;\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot t\_0}{1 + \cos \lambda_1}\\ \end{array} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (let* ((t_0 (sin (- lambda1 lambda2))))
   (if (<= phi2 195000.0)
     (*
      (+
       1.0
       (/
        (atan2
         (* t_0 (fma (* -0.5 phi2) phi2 1.0))
         (- (cos (- lambda2 lambda1)) -1.0))
        lambda1))
      lambda1)
     (+ lambda1 (atan2 (* (cos phi2) t_0) (+ 1.0 (cos lambda1)))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	double t_0 = sin((lambda1 - lambda2));
	double tmp;
	if (phi2 <= 195000.0) {
		tmp = (1.0 + (atan2((t_0 * fma((-0.5 * phi2), phi2, 1.0)), (cos((lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1;
	} else {
		tmp = lambda1 + atan2((cos(phi2) * t_0), (1.0 + cos(lambda1)));
	}
	return tmp;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	t_0 = sin(Float64(lambda1 - lambda2))
	tmp = 0.0
	if (phi2 <= 195000.0)
		tmp = Float64(Float64(1.0 + Float64(atan(Float64(t_0 * fma(Float64(-0.5 * phi2), phi2, 1.0)), Float64(cos(Float64(lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1);
	else
		tmp = Float64(lambda1 + atan(Float64(cos(phi2) * t_0), Float64(1.0 + cos(lambda1))));
	end
	return tmp
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := Block[{t$95$0 = N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[phi2, 195000.0], N[(N[(1.0 + N[(N[ArcTan[N[(t$95$0 * N[(N[(-0.5 * phi2), $MachinePrecision] * phi2 + 1.0), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision]], $MachinePrecision] / lambda1), $MachinePrecision]), $MachinePrecision] * lambda1), $MachinePrecision], N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * t$95$0), $MachinePrecision] / N[(1.0 + N[Cos[lambda1], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if phi2 < 195000

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. Taylor expanded in phi2 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
   3. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

      3. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      4. lower--.f64 78.4

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. Applied rewrites 78.4%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

   5. Taylor expanded in phi1 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
   6. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      3. lower--.f64 68.2

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. Applied rewrites 68.2%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   8. Taylor expanded in phi2 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)} \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]
   9. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      2. lower-*.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      3. lower-pow.f64 65.9

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   10. Applied rewrites 65.9%

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)} \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]
   11. Step-by-step derivation

       1. lift-+.f64 N/A

          \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

       2. sum-to-mult N/A

          \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

       3. lower-special-*.f64 N/A

          \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

   12. Applied rewrites 65.9%

       \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(-0.5 \cdot \phi_2, \phi_2, 1\right)}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1} \]

   if 195000 < phi2

   1. Initial program 98.5%

      \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. Taylor expanded in phi2 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
   3. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

      3. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      4. lower--.f64 78.4

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. Applied rewrites 78.4%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

   5. Taylor expanded in phi1 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
   6. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      2. lower-cos.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

      3. lower--.f64 68.2

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   7. Applied rewrites 68.2%

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   8. Taylor expanded in lambda2 around 0

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \lambda_1} \]
   9. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \lambda_1} \]

      2. lower-cos.f64 62.8

         \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \lambda_1} \]

   10. Applied rewrites 62.8%

       \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \lambda_1} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 17: 67.0% accurate, 1.6× speedup

Math

\[\begin{array}{l} \\ \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(-\lambda_2\right)} \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  lambda1
  (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ 1.0 (cos (- lambda2))))))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (1.0 + cos(-lambda2)));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(lambda1, lambda2, phi1, phi2)
use fmin_fmax_functions
    real(8), intent (in) :: lambda1
    real(8), intent (in) :: lambda2
    real(8), intent (in) :: phi1
    real(8), intent (in) :: phi2
    code = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (1.0d0 + cos(-lambda2)))
end function

Java

public static double code(double lambda1, double lambda2, double phi1, double phi2) {
	return lambda1 + Math.atan2((Math.cos(phi2) * Math.sin((lambda1 - lambda2))), (1.0 + Math.cos(-lambda2)));
}

Python

def code(lambda1, lambda2, phi1, phi2):
	return lambda1 + math.atan2((math.cos(phi2) * math.sin((lambda1 - lambda2))), (1.0 + math.cos(-lambda2)))

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(lambda1 + atan(Float64(cos(phi2) * sin(Float64(lambda1 - lambda2))), Float64(1.0 + cos(Float64(-lambda2)))))
end

MATLAB

function tmp = code(lambda1, lambda2, phi1, phi2)
	tmp = lambda1 + atan2((cos(phi2) * sin((lambda1 - lambda2))), (1.0 + cos(-lambda2)));
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(lambda1 + N[ArcTan[N[(N[Cos[phi2], $MachinePrecision] * N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[Cos[(-lambda2)], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

2. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
3. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

   3. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. lower--.f64 78.4

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

4. Applied rewrites 78.4%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

5. Taylor expanded in phi1 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
6. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower--.f64 68.2

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

7. Applied rewrites 68.2%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

8. Taylor expanded in lambda1 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\mathsf{neg}\left(\lambda_2\right)\right)} \]
9. Step-by-step derivation

   1. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\mathsf{neg}\left(\lambda_2\right)\right)} \]

   2. lower-neg.f64 68.1

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(-\lambda_2\right)} \]

10. Applied rewrites 68.1%

    \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(-\lambda_2\right)} \]
11. Add Preprocessing

Alternative 18: 65.9% accurate, 1.8× speedup

Math

\[\begin{array}{l} \\ \left(1 + \frac{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(-0.5 \cdot \phi_2, \phi_2, 1\right)}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1 \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (*
  (+
   1.0
   (/
    (atan2
     (* (sin (- lambda1 lambda2)) (fma (* -0.5 phi2) phi2 1.0))
     (- (cos (- lambda2 lambda1)) -1.0))
    lambda1))
  lambda1))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return (1.0 + (atan2((sin((lambda1 - lambda2)) * fma((-0.5 * phi2), phi2, 1.0)), (cos((lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(Float64(1.0 + Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * fma(Float64(-0.5 * phi2), phi2, 1.0)), Float64(cos(Float64(lambda2 - lambda1)) - -1.0)) / lambda1)) * lambda1)
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[(1.0 + N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[(N[(-0.5 * phi2), $MachinePrecision] * phi2 + 1.0), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision]], $MachinePrecision] / lambda1), $MachinePrecision]), $MachinePrecision] * lambda1), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

2. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
3. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

   3. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. lower--.f64 78.4

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

4. Applied rewrites 78.4%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

5. Taylor expanded in phi1 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
6. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower--.f64 68.2

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

7. Applied rewrites 68.2%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

8. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)} \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]
9. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower-pow.f64 65.9

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

10. Applied rewrites 65.9%

    \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)} \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]
11. Step-by-step derivation

    1. lift-+.f64 N/A

       \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

    2. sum-to-mult N/A

       \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

    3. lower-special-*.f64 N/A

       \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}}{\lambda_1}\right) \cdot \lambda_1} \]

12. Applied rewrites 65.9%

    \[\leadsto \color{blue}{\left(1 + \frac{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(-0.5 \cdot \phi_2, \phi_2, 1\right)}{\cos \left(\lambda_2 - \lambda_1\right) - -1}}{\lambda_1}\right) \cdot \lambda_1} \]
13. Add Preprocessing

Alternative 19: 65.9% accurate, 1.9× speedup

Math

\[\begin{array}{l} \\ \tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(-0.5 \cdot \phi_2, \phi_2, 1\right)}{\cos \left(\lambda_2 - \lambda_1\right) - -1} + \lambda_1 \end{array} \]

FPCore

(FPCore (lambda1 lambda2 phi1 phi2)
 :precision binary64
 (+
  (atan2
   (* (sin (- lambda1 lambda2)) (fma (* -0.5 phi2) phi2 1.0))
   (- (cos (- lambda2 lambda1)) -1.0))
  lambda1))

C

double code(double lambda1, double lambda2, double phi1, double phi2) {
	return atan2((sin((lambda1 - lambda2)) * fma((-0.5 * phi2), phi2, 1.0)), (cos((lambda2 - lambda1)) - -1.0)) + lambda1;
}

Julia

function code(lambda1, lambda2, phi1, phi2)
	return Float64(atan(Float64(sin(Float64(lambda1 - lambda2)) * fma(Float64(-0.5 * phi2), phi2, 1.0)), Float64(cos(Float64(lambda2 - lambda1)) - -1.0)) + lambda1)
end

Wolfram

code[lambda1_, lambda2_, phi1_, phi2_] := N[(N[ArcTan[N[(N[Sin[N[(lambda1 - lambda2), $MachinePrecision]], $MachinePrecision] * N[(N[(-0.5 * phi2), $MachinePrecision] * phi2 + 1.0), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[N[(lambda2 - lambda1), $MachinePrecision]], $MachinePrecision] - -1.0), $MachinePrecision]], $MachinePrecision] + lambda1), $MachinePrecision]

Derivation

1. Initial program 98.5%

   \[\lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \phi_2 \cdot \cos \left(\lambda_1 - \lambda_2\right)} \]

2. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]
3. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \color{blue}{\left(\lambda_1 - \lambda_2\right)}} \]

   3. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   4. lower--.f64 78.4

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

4. Applied rewrites 78.4%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{\color{blue}{\cos \phi_1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

5. Taylor expanded in phi1 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]
6. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-cos.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower--.f64 68.2

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

7. Applied rewrites 68.2%

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\cos \phi_2 \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \color{blue}{\cos \left(\lambda_1 - \lambda_2\right)}} \]

8. Taylor expanded in phi2 around 0

   \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\color{blue}{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right)} \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]
9. Step-by-step derivation

   1. lower-+.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + \color{blue}{\frac{-1}{2} \cdot {\phi_2}^{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   2. lower-*.f64 N/A

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot \color{blue}{{\phi_2}^{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

   3. lower-pow.f64 65.9

      \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\left(1 + -0.5 \cdot {\phi_2}^{\color{blue}{2}}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]

10. Applied rewrites 65.9%

    \[\leadsto \lambda_1 + \tan^{-1}_* \frac{\color{blue}{\left(1 + -0.5 \cdot {\phi_2}^{2}\right)} \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} \]
11. Step-by-step derivation

    1. lift-+.f64 N/A

       \[\leadsto \color{blue}{\lambda_1 + \tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)}} \]

    2. +-commutative N/A

       \[\leadsto \color{blue}{\tan^{-1}_* \frac{\left(1 + \frac{-1}{2} \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

    3. lower-+.f64 65.9

       \[\leadsto \color{blue}{\tan^{-1}_* \frac{\left(1 + -0.5 \cdot {\phi_2}^{2}\right) \cdot \sin \left(\lambda_1 - \lambda_2\right)}{1 + \cos \left(\lambda_1 - \lambda_2\right)} + \lambda_1} \]

12. Applied rewrites 65.9%

    \[\leadsto \color{blue}{\tan^{-1}_* \frac{\sin \left(\lambda_1 - \lambda_2\right) \cdot \mathsf{fma}\left(-0.5 \cdot \phi_2, \phi_2, 1\right)}{\cos \left(\lambda_2 - \lambda_1\right) - -1} + \lambda_1} \]
13. Add Preprocessing

Reproduce

herbie shell --seed 2025150 
(FPCore (lambda1 lambda2 phi1 phi2)
  :name "Midpoint on a great circle"
  :precision binary64
  (+ lambda1 (atan2 (* (cos phi2) (sin (- lambda1 lambda2))) (+ (cos phi1) (* (cos phi2) (cos (- lambda1 lambda2)))))))