double get_precision(double FPR, double TPR, long P, long N) {
  if (2 - FPR - TPR > 0)
    return (1 - FPR) * N / ((1 - FPR) * N + (1 - TPR) * P);
  else
    return 1;
}


x86 assembly (clang, -O0):
<get_precision>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %rdi,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
movsd  0xee2(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    118c <get_precision+0x8c>
movsd  0xecb(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
cvtsi2sdq -0x28(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xeb4(%rip),%xmm1
subsd  -0x10(%rbp),%xmm1
cvtsi2sdq -0x28(%rbp),%xmm3
movsd  0xea1(%rip),%xmm2
subsd  -0x18(%rbp),%xmm2
cvtsi2sdq -0x20(%rbp),%xmm4
mulsd  %xmm4,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1199 <get_precision+0x99>
movsd  0xe74(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_precision>:
movapd %xmm0,%xmm2
movsd  0xef4(%rip),%xmm0
subsd  %xmm2,%xmm0
subsd  %xmm1,%xmm0
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm0
jbe    1153 <get_precision+0x53>
movsd  0xee2(%rip),%xmm3
movapd %xmm3,%xmm4
xorps  %xmm0,%xmm0
cvtsi2sd %rsi,%xmm0
subsd  %xmm2,%xmm4
mulsd  %xmm4,%xmm0
xorps  %xmm2,%xmm2
cvtsi2sd %rdi,%xmm2
subsd  %xmm1,%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
divsd  %xmm2,%xmm0
ret
movsd  0xead(%rip),%xmm0
ret

LLVM IR (-O0):
define double @get_precision(double %0, double %1, i64 %2, i64 %3) {
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca i64
%9 = alloca i64
store double %0, double* %6
store double %1, double* %7
store i64 %2, i64* %8
store i64 %3, i64* %9
%10 = load double, double* %6
%11 = fsub double 2.000000e+00, %10
%12 = load double, double* %7
%13 = fsub double %11, %12
%14 = fcmp ogt double %13, 0.000000e+00
br i1 %14, label %15, label %32
15: ; preds = %4
%16 = load double, double* %6
%17 = fsub double 1.000000e+00, %16
%18 = load i64, i64* %9
%19 = sitofp i64 %18 to double
%20 = fmul double %17, %19
%21 = load double, double* %6
%22 = fsub double 1.000000e+00, %21
%23 = load i64, i64* %9
%24 = sitofp i64 %23 to double
%25 = load double, double* %7
%26 = fsub double 1.000000e+00, %25
%27 = load i64, i64* %8
%28 = sitofp i64 %27 to double
%29 = fmul double %26, %28
%30 = call double @llvm.fmuladd.f64(double %22, double %24, double %29)
%31 = fdiv double %20, %30
store double %31, double* %5
br label %33
32: ; preds = %4
store double 1.000000e+00, double* %5
br label %33
33: ; preds = %32, %15
%34 = load double, double* %5
ret double %34
}

LLVM IR (-O2):
define double @get_precision(double %0, double %1, i64 %2, i64 %3) {
%5 = fsub double 2.000000e+00, %0
%6 = fsub double %5, %1
%7 = fcmp ogt double %6, 0.000000e+00
br i1 %7, label %8, label %17
8: ; preds = %4
%9 = fsub double 1.000000e+00, %0
%10 = sitofp i64 %3 to double
%11 = fmul double %9, %10
%12 = fsub double 1.000000e+00, %1
%13 = sitofp i64 %2 to double
%14 = fmul double %12, %13
%15 = tail call double @llvm.fmuladd.f64(double %9, double %10, double %14)
%16 = fdiv double %11, %15
br label %17
17: ; preds = %4, %8
%18 = phi double [ %16, %8 ], [ 1.000000e+00, %4 ]
ret double %18
}