Timeout for init-u0

Specification

Math

\[\begin{array}{l} t_0 := \left|ArB \cdot \tan^{-1} \left(\frac{\sqrt{D \cdot D - 1}}{\cos alpha\_c}\right)\right|\\ \mathbf{if}\;phi0 < 0:\\ \;\;\;\;-t\_0\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \]

FPCore

(FPCore (ArB D alpha_c phi0)
  :precision binary64
  :pre TRUE
  (let* ((t_0
        (fabs
         (* ArB (atan (/ (sqrt (- (* D D) 1.0)) (cos alpha_c)))))))
  (if (< phi0 0.0) (- t_0) t_0)))

C

double code(double ArB, double D, double alpha_c, double phi0) {
	double t_0 = fabs((ArB * atan((sqrt(((D * D) - 1.0)) / cos(alpha_c)))));
	double tmp;
	if (phi0 < 0.0) {
		tmp = -t_0;
	} else {
		tmp = t_0;
	}
	return tmp;
}

Fortran

real(8) function code(arb, d, alpha_c, phi0)
use fmin_fmax_functions
    real(8), intent (in) :: arb
    real(8), intent (in) :: d
    real(8), intent (in) :: alpha_c
    real(8), intent (in) :: phi0
    real(8) :: t_0
    real(8) :: tmp
    t_0 = abs((arb * atan((sqrt(((d * d) - 1.0d0)) / cos(alpha_c)))))
    if (phi0 < 0.0d0) then
        tmp = -t_0
    else
        tmp = t_0
    end if
    code = tmp
end function

Java

public static double code(double ArB, double D, double alpha_c, double phi0) {
	double t_0 = Math.abs((ArB * Math.atan((Math.sqrt(((D * D) - 1.0)) / Math.cos(alpha_c)))));
	double tmp;
	if (phi0 < 0.0) {
		tmp = -t_0;
	} else {
		tmp = t_0;
	}
	return tmp;
}

Python

def code(ArB, D, alpha_c, phi0):
	t_0 = math.fabs((ArB * math.atan((math.sqrt(((D * D) - 1.0)) / math.cos(alpha_c)))))
	tmp = 0
	if phi0 < 0.0:
		tmp = -t_0
	else:
		tmp = t_0
	return tmp

Julia

function code(ArB, D, alpha_c, phi0)
	t_0 = abs(Float64(ArB * atan(Float64(sqrt(Float64(Float64(D * D) - 1.0)) / cos(alpha_c)))))
	tmp = 0.0
	if (phi0 < 0.0)
		tmp = Float64(-t_0);
	else
		tmp = t_0;
	end
	return tmp
end

MATLAB

function tmp_2 = code(ArB, D, alpha_c, phi0)
	t_0 = abs((ArB * atan((sqrt(((D * D) - 1.0)) / cos(alpha_c)))));
	tmp = 0.0;
	if (phi0 < 0.0)
		tmp = -t_0;
	else
		tmp = t_0;
	end
	tmp_2 = tmp;
end

Wolfram

code[ArB_, D_, alpha$95$c_, phi0_] := Block[{t$95$0 = N[Abs[N[(ArB * N[ArcTan[N[(N[Sqrt[N[(N[(D * D), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision] / N[Cos[alpha$95$c], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[Less[phi0, 0.0], (-t$95$0), t$95$0]]

ReFlow

f(ArB, D, alpha_c, phi0):
	ArB in [-inf, +inf],
	D in [-inf, +inf],
	alpha_c in [-inf, +inf],
	phi0 in [-inf, +inf]
code: THEORY
BEGIN
f(ArB, D, alpha_c, phi0: real): real =
	LET t_0 = (abs((ArB * (atan(((sqrt(((D * D) - (1)))) / (cos(alpha_c)))))))) IN
		LET tmp = IF (phi0 < (0)) THEN (- t_0) ELSE t_0 ENDIF IN
	tmp
END code