Timeout for inverse-phi

Specification

Math

\[\tan^{-1} \left(\frac{\tan lamdp \cdot \cos lamt - ca \cdot \sin lamt}{one\_es \cdot sa}\right) \]

FPCore

(FPCore (lamdp lamt ca one_es sa)
  :precision binary64
  :pre TRUE
  (atan
 (/ (- (* (tan lamdp) (cos lamt)) (* ca (sin lamt))) (* one_es sa))))

C

double code(double lamdp, double lamt, double ca, double one_es, double sa) {
	return atan((((tan(lamdp) * cos(lamt)) - (ca * sin(lamt))) / (one_es * sa)));
}

Fortran

real(8) function code(lamdp, lamt, ca, one_es, sa)
use fmin_fmax_functions
    real(8), intent (in) :: lamdp
    real(8), intent (in) :: lamt
    real(8), intent (in) :: ca
    real(8), intent (in) :: one_es
    real(8), intent (in) :: sa
    code = atan((((tan(lamdp) * cos(lamt)) - (ca * sin(lamt))) / (one_es * sa)))
end function

Java

public static double code(double lamdp, double lamt, double ca, double one_es, double sa) {
	return Math.atan((((Math.tan(lamdp) * Math.cos(lamt)) - (ca * Math.sin(lamt))) / (one_es * sa)));
}

Python

def code(lamdp, lamt, ca, one_es, sa):
	return math.atan((((math.tan(lamdp) * math.cos(lamt)) - (ca * math.sin(lamt))) / (one_es * sa)))

Julia

function code(lamdp, lamt, ca, one_es, sa)
	return atan(Float64(Float64(Float64(tan(lamdp) * cos(lamt)) - Float64(ca * sin(lamt))) / Float64(one_es * sa)))
end

MATLAB

function tmp = code(lamdp, lamt, ca, one_es, sa)
	tmp = atan((((tan(lamdp) * cos(lamt)) - (ca * sin(lamt))) / (one_es * sa)));
end

Wolfram

code[lamdp_, lamt_, ca_, one$95$es_, sa_] := N[ArcTan[N[(N[(N[(N[Tan[lamdp], $MachinePrecision] * N[Cos[lamt], $MachinePrecision]), $MachinePrecision] - N[(ca * N[Sin[lamt], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(one$95$es * sa), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]

ReFlow

f(lamdp, lamt, ca, one_es, sa):
	lamdp in [-inf, +inf],
	lamt in [-inf, +inf],
	ca in [-inf, +inf],
	one_es in [-inf, +inf],
	sa in [-inf, +inf]
code: THEORY
BEGIN
f(lamdp, lamt, ca, one_es, sa: real): real =
	atan(((((tan(lamdp)) * (cos(lamt))) - (ca * (sin(lamt)))) / (one_es * sa)))
END code