#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
float incbetf(float, float, double);
int /*<<< orphan*/ mtherrf(char *, int);
float powf(double, float);
float bdtrf(int k, int n, float pp) {
  float p, dk, dn;
  p = pp;
  if ((p < 0.0) || (p > 1.0))
    goto domerr;
  if ((k < 0) || (n < k)) {
  domerr:
    mtherrf("bdtrf", 1);
    return (0.0);
  }
  if (k == n)
    return (1.0);
  dn = n - k;
  if (k == 0) {
    dk = powf(1.0 - p, dn);
  } else {
    dk = k + 1;
    dk = incbetf(dn, dk, 1.0 - p);
  }
  return (dk);
}


x86 assembly (clang, -O0):
<bdtrf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
ja     117e <bdtrf+0x4e>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe8c(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1183 <bdtrf+0x53>
jmp    119e <bdtrf+0x6e>
cmpl   $0x0,-0x8(%rbp)
jl     1199 <bdtrf+0x69>
mov    -0xc(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11bc <bdtrf+0x8c>
jmp    119e <bdtrf+0x6e>
lea    0xe67(%rip),%rdi
mov    $0x1,%esi
call   1050 <mtherrf@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    125e <bdtrf+0x12e>
mov    -0x8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jne    11da <bdtrf+0xaa>
movss  0xe38(%rip),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    125e <bdtrf+0x12e>
mov    -0xc(%rbp),%eax
sub    -0x8(%rbp),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x1c(%rbp)
cmpl   $0x0,-0x8(%rbp)
jne    121c <bdtrf+0xec>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdfc(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x1c(%rbp),%xmm1
call   1030 <powf@plt>
movss  %xmm0,-0x18(%rbp)
jmp    1254 <bdtrf+0x124>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x1c(%rbp),%xmm0
movss  -0x18(%rbp),%xmm1
movss  -0x14(%rbp),%xmm2
cvtss2sd %xmm2,%xmm3
movsd  0xdba(%rip),%xmm2
subsd  %xmm3,%xmm2
call   1040 <incbetf@plt>
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bdtrf>:
push   %rax
xorps  %xmm1,%xmm1
ucomiss %xmm0,%xmm1
jbe    114f <bdtrf+0x1f>
lea    0xed0(%rip),%rdi
mov    $0x1,%esi
call   1050 <mtherrf@plt>
xorps  %xmm0,%xmm0
pop    %rax
ret
cmp    %edi,%esi
jl     1139 <bdtrf+0x9>
test   %edi,%edi
js     1139 <bdtrf+0x9>
ucomiss 0xea2(%rip),%xmm0
ja     1139 <bdtrf+0x9>
sub    %edi,%esi
jne    116e <bdtrf+0x3e>
movss  0xe94(%rip),%xmm0
pop    %rax
ret
cvtss2sd %xmm0,%xmm4
cvtsi2ss %esi,%xmm3
test   %edi,%edi
je     1199 <bdtrf+0x69>
add    $0x1,%edi
xorps  %xmm1,%xmm1
cvtsi2ss %edi,%xmm1
movsd  0xe7c(%rip),%xmm2
subsd  %xmm4,%xmm2
movaps %xmm3,%xmm0
pop    %rax
jmp    1040 <incbetf@plt>
movsd  0xe67(%rip),%xmm0
subsd  %xmm4,%xmm0
movaps %xmm3,%xmm1
pop    %rax
jmp    1030 <powf@plt>

LLVM IR (-O0):
@.str = [6 x i8] c"bdtrf\00"
define float @bdtrf(i32 %0, i32 %1, float %2) {
%4 = alloca float
%5 = alloca i32
%6 = alloca i32
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
store i32 %0, i32* %5
store i32 %1, i32* %6
store float %2, float* %7
%11 = load float, float* %7
store float %11, float* %8
%12 = load float, float* %8
%13 = fpext float %12 to double
%14 = fcmp olt double %13, 0.000000e+00
br i1 %14, label %19, label %15
15: ; preds = %3
%16 = load float, float* %8
%17 = fpext float %16 to double
%18 = fcmp ogt double %17, 1.000000e+00
br i1 %18, label %19, label %20
19: ; preds = %15, %3
br label %28
20: ; preds = %15
%21 = load i32, i32* %5
%22 = icmp slt i32 %21, 0
br i1 %22, label %27, label %23
23: ; preds = %20
%24 = load i32, i32* %6
%25 = load i32, i32* %5
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %30
27: ; preds = %23, %20
br label %28
28: ; preds = %27, %19
%29 = call i32 @mtherrf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str, i64 0, i64 0), i32 1)
store float 0.000000e+00, float* %4
br label %60
30: ; preds = %23
%31 = load i32, i32* %5
%32 = load i32, i32* %6
%33 = icmp eq i32 %31, %32
br i1 %33, label %34, label %35
34: ; preds = %30
store float 1.000000e+00, float* %4
br label %60
35: ; preds = %30
%36 = load i32, i32* %6
%37 = load i32, i32* %5
%38 = sub nsw i32 %36, %37
%39 = sitofp i32 %38 to float
store float %39, float* %10
%40 = load i32, i32* %5
%41 = icmp eq i32 %40, 0
br i1 %41, label %42, label %48
42: ; preds = %35
%43 = load float, float* %8
%44 = fpext float %43 to double
%45 = fsub double 1.000000e+00, %44
%46 = load float, float* %10
%47 = call float @powf(double %45, float %46)
store float %47, float* %9
br label %58
48: ; preds = %35
%49 = load i32, i32* %5
%50 = add nsw i32 %49, 1
%51 = sitofp i32 %50 to float
store float %51, float* %9
%52 = load float, float* %10
%53 = load float, float* %9
%54 = load float, float* %8
%55 = fpext float %54 to double
%56 = fsub double 1.000000e+00, %55
%57 = call float @incbetf(float %52, float %53, double %56)
store float %57, float* %9
br label %58
58: ; preds = %48, %42
%59 = load float, float* %9
store float %59, float* %4
br label %60
60: ; preds = %58, %34, %28
%61 = load float, float* %4
ret float %61
}

LLVM IR (-O2):
@.str = [6 x i8] c"bdtrf\00"
define float @bdtrf(i32 %0, i32 %1, float %2) {
%4 = fpext float %2 to double
%5 = fcmp olt float %2, 0.000000e+00
br i1 %5, label %12, label %6
6: ; preds = %3
%7 = fcmp ogt float %2, 1.000000e+00
%8 = icmp slt i32 %0, 0
%9 = or i1 %8, %7
%10 = icmp slt i32 %1, %0
%11 = or i1 %10, %9
br i1 %11, label %12, label %14
12: ; preds = %3, %6
%13 = tail call i32 @mtherrf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str, i64 0, i64 0), i32 1)
br label %28
14: ; preds = %6
%15 = icmp eq i32 %1, %0
br i1 %15, label %28, label %16
16: ; preds = %14
%17 = sub nsw i32 %1, %0
%18 = sitofp i32 %17 to float
%19 = icmp eq i32 %0, 0
br i1 %19, label %20, label %23
20: ; preds = %16
%21 = fsub double 1.000000e+00, %4
%22 = tail call float @powf(double %21, float %18)
br label %28
23: ; preds = %16
%24 = add nuw nsw i32 %0, 1
%25 = sitofp i32 %24 to float
%26 = fsub double 1.000000e+00, %4
%27 = tail call float @incbetf(float %18, float %25, double %26)
br label %28
28: ; preds = %20, %23, %14, %12
%29 = phi float [ 0.000000e+00, %12 ], [ 1.000000e+00, %14 ], [ %22, %20 ], [ %27, %23 ]
ret float %29
}