_Bool isJFC(double q, double e, double i) {
  if (q < 1.3)
    return 0;
  double t =
      5.2 * (1 - e) / q + 2 * sqrt(q * (1 + e) / 5.2) *
                              cos(i * 3.14159265358979323846264338327 / 180.);
  return t < 3 && t > 2;
}


x86 assembly (clang, -O0):
<isJFC>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  0xec1(%rip),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    1153 <isJFC+0x33>
movb   $0x0,-0x1(%rbp)
jmp    1220 <isJFC+0x100>
movsd  0xed5(%rip),%xmm1
subsd  -0x18(%rbp),%xmm1
movsd  0xec0(%rip),%xmm0
mulsd  %xmm1,%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  0xead(%rip),%xmm1
addsd  -0x18(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe94(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xe78(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe5f(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xe4a(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <cos@plt>
movsd  -0x40(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x38(%rbp),%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
xor    %eax,%eax
movsd  0xe15(%rip),%xmm0
ucomisd -0x28(%rbp),%xmm0
mov    %al,-0x29(%rbp)
jbe    1218 <isJFC+0xf8>
movsd  -0x28(%rbp),%xmm0
movsd  0xe12(%rip),%xmm1
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x29(%rbp)
mov    -0x29(%rbp),%al
and    $0x1,%al
mov    %al,-0x1(%rbp)
mov    -0x1(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<isJFC>:
sub    $0x18,%rsp
movsd  0xed4(%rip),%xmm3
ucomisd %xmm0,%xmm3
jbe    1139 <isJFC+0x19>
xor    %eax,%eax
add    $0x18,%rsp
ret
movsd  0xec7(%rip),%xmm4
movapd %xmm4,%xmm5
subsd  %xmm1,%xmm5
movsd  0xebf(%rip),%xmm3
mulsd  %xmm3,%xmm5
divsd  %xmm0,%xmm5
movsd  %xmm5,0x10(%rsp)
addsd  %xmm4,%xmm1
mulsd  %xmm0,%xmm1
divsd  %xmm3,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jb     117e <isJFC+0x5e>
xorps  %xmm0,%xmm0
sqrtsd %xmm1,%xmm0
jmp    1193 <isJFC+0x73>
movapd %xmm1,%xmm0
movsd  %xmm2,0x8(%rsp)
call   1040 <sqrt@plt>
movsd  0x8(%rsp),%xmm2
addsd  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
mulsd  0xe73(%rip),%xmm2
divsd  0xe73(%rip),%xmm2
movapd %xmm2,%xmm0
call   1030 <cos@plt>
mulsd  0x8(%rsp),%xmm0
addsd  0x10(%rsp),%xmm0
movsd  0xe5e(%rip),%xmm1
ucomisd %xmm0,%xmm1
seta   %cl
ucomisd 0xe57(%rip),%xmm0
seta   %al
and    %cl,%al
add    $0x18,%rsp
ret

LLVM IR (-O0):
define zeroext i1 @isJFC(double %0, double %1, double %2) {
%4 = alloca i1
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca double
store double %0, double* %5
store double %1, double* %6
store double %2, double* %7
%9 = load double, double* %5
%10 = fcmp olt double %9, 1.300000e+00
br i1 %10, label %11, label %12
11: ; preds = %3
store i1 false, i1* %4
br label %37
12: ; preds = %3
%13 = load double, double* %6
%14 = fsub double 1.000000e+00, %13
%15 = fmul double 5.200000e+00, %14
%16 = load double, double* %5
%17 = fdiv double %15, %16
%18 = load double, double* %5
%19 = load double, double* %6
%20 = fadd double 1.000000e+00, %19
%21 = fmul double %18, %20
%22 = fdiv double %21, 5.200000e+00
%23 = call double @sqrt(double %22)
%24 = fmul double 2.000000e+00, %23
%25 = load double, double* %7
%26 = fmul double %25, 0x400921FB54442D18
%27 = fdiv double %26, 1.800000e+02
%28 = call double @cos(double %27)
%29 = call double @llvm.fmuladd.f64(double %24, double %28, double %17)
store double %29, double* %8
%30 = load double, double* %8
%31 = fcmp olt double %30, 3.000000e+00
br i1 %31, label %32, label %35
32: ; preds = %12
%33 = load double, double* %8
%34 = fcmp ogt double %33, 2.000000e+00
br label %35
35: ; preds = %32, %12
%36 = phi i1 [ false, %12 ], [ %34, %32 ]
store i1 %36, i1* %4
br label %37
37: ; preds = %35, %11
%38 = load i1, i1* %4
ret i1 %38
}

LLVM IR (-O2):
define zeroext i1 @isJFC(double %0, double %1, double %2) {
%4 = fcmp olt double %0, 1.300000e+00
br i1 %4, label %21, label %5
5: ; preds = %3
%6 = fsub double 1.000000e+00, %1
%7 = fmul double %6, 5.200000e+00
%8 = fdiv double %7, %0
%9 = fadd double %1, 1.000000e+00
%10 = fmul double %9, %0
%11 = fdiv double %10, 5.200000e+00
%12 = tail call double @sqrt(double %11)
%13 = fmul double %12, 2.000000e+00
%14 = fmul double %2, 0x400921FB54442D18
%15 = fdiv double %14, 1.800000e+02
%16 = tail call double @cos(double %15)
%17 = tail call double @llvm.fmuladd.f64(double %13, double %16, double %8)
%18 = fcmp olt double %17, 3.000000e+00
%19 = fcmp ogt double %17, 2.000000e+00
%20 = and i1 %18, %19
br label %21
21: ; preds = %3, %5
%22 = phi i1 [ %20, %5 ], [ false, %3 ]
ret i1 %22
}