#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 */
double DF(double, double, double, double, double);
double genunf(double, double);
scalar_t__ is_positive(double);
scalar_t__ log(double);
double maximum_point(double, double, double, double);
scalar_t__ my_fmin(double, double);
double qacc(double, double, double, double, double, double, double, double);
double sgamma(double);
double sample_quad(double xkl, double xkk, double xll, double ckl, double clk,
                   double ckk, double cll, double bk, double bl) {
  double xlk, skl, slk, s2, s3, s, a1, a2, a3;
  double vbar, alpha, beta, v, w;
  double q, U;
  xlk = xkl;
  skl = xkk + xkl;
  slk = xll + xlk;
  if (skl <= slk) {
    s2 = skl;
    s3 = slk;
    s = s3 / s2;
    a1 = ckl + clk - 1.0;
    a2 = ckk + bk - 1.0;
    a3 = cll + bl - 1.0;
  } else {
    s2 = slk;
    s3 = skl;
    s = s3 / s2;
    a1 = ckl + clk - 1.0;
    a2 = cll + bl - 1.0;
    a3 = ckk + bk - 1.0;
  }
  if (is_positive(s - 1.0)) {
    vbar = maximum_point(s, a1, a2, a3);
    beta = -1.0 * DF(vbar, s, a1, a2, a3) * vbar;
    alpha = beta * vbar;
    if (is_positive(s2 - xkl)) {
      v = xkl / (s2 - xkl);
      if (is_positive(alpha) && is_positive(1.0 / beta)) {
        w = 1.0 / beta * sgamma(alpha);
        if (is_positive(w)) {
          if (!is_positive(v)) {
            return s2 * w / (1.0 + w);
          } else {
            q = qacc(w, v, s, a1, a2, a3, alpha, beta);
            U = genunf(0.0, 1.0);
            if (log(U) < my_fmin(0.0, q)) {
              return s2 * w / (1.0 + w);
            }
          }
        }
      }
    }
  }
  return xkl;
}


x86 assembly (clang, -O0):
<sample_quad>:
push   %rbp
mov    %rsp,%rbp
sub    $0xe0,%rsp
movsd  0x10(%rbp),%xmm8
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
movsd  %xmm4,-0x30(%rbp)
movsd  %xmm5,-0x38(%rbp)
movsd  %xmm6,-0x40(%rbp)
movsd  %xmm7,-0x48(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x20(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     1274 <sample_quad+0xf4>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  0xdd6(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  0xdbb(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
addsd  0x10(%rbp),%xmm0
movsd  0xd9d(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
jmp    12ee <sample_quad+0x16e>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  0xd57(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  0x10(%rbp),%xmm0
movsd  0xd3c(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  0xd1e(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  0xd05(%rip),%xmm1
subsd  %xmm1,%xmm0
call   10a0 <is_positive@plt>
cmp    $0x0,%rax
je     156c <sample_quad+0x3ec>
movsd  -0x78(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm1
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x90(%rbp),%xmm3
ff
call   1040 <maximum_point@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x78(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm3
ff
movsd  -0x90(%rbp),%xmm4
ff
call   1050 <DF@plt>
movaps %xmm0,%xmm1
movsd  0xca1(%rip),%xmm0
mulsd  %xmm1,%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
call   10a0 <is_positive@plt>
cmp    $0x0,%rax
je     1567 <sample_quad+0x3e7>
movsd  -0x10(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
call   10a0 <is_positive@plt>
cmp    $0x0,%rax
je     1562 <sample_quad+0x3e2>
movsd  0xc1a(%rip),%xmm0
divsd  -0xa8(%rbp),%xmm0
ff
call   10a0 <is_positive@plt>
cmp    $0x0,%rax
je     1562 <sample_quad+0x3e2>
movsd  0xbfb(%rip),%xmm0
divsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
call   1090 <sgamma@plt>
movaps %xmm0,%xmm1
movsd  -0xd0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
call   10a0 <is_positive@plt>
cmp    $0x0,%rax
je     155d <sample_quad+0x3dd>
movsd  -0xb0(%rbp),%xmm0
ff
call   10a0 <is_positive@plt>
cmp    $0x0,%rax
jne    1492 <sample_quad+0x312>
movsd  -0x68(%rbp),%xmm0
mulsd  -0xb8(%rbp),%xmm0
ff
movsd  0xb84(%rip),%xmm1
addsd  -0xb8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1576 <sample_quad+0x3f6>
movsd  -0xb8(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm1
ff
movsd  -0x78(%rbp),%xmm2
movsd  -0x80(%rbp),%xmm3
movsd  -0x88(%rbp),%xmm4
ff
movsd  -0x90(%rbp),%xmm5
ff
movsd  -0xa0(%rbp),%xmm6
ff
movsd  -0xa8(%rbp),%xmm7
ff
call   1060 <qacc@plt>
movsd  %xmm0,-0xc0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  0xb1c(%rip),%xmm1
call   1080 <genunf@plt>
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1070 <log@plt>
mov    %rax,-0xd8(%rbp)
movsd  -0xc0(%rbp),%xmm1
ff
xorps  %xmm0,%xmm0
call   1030 <my_fmin@plt>
mov    %rax,%rcx
mov    -0xd8(%rbp),%rax
cmp    %rcx,%rax
jge    1553 <sample_quad+0x3d3>
movsd  -0x68(%rbp),%xmm0
mulsd  -0xb8(%rbp),%xmm0
ff
movsd  0xac3(%rip),%xmm1
addsd  -0xb8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1576 <sample_quad+0x3f6>
jmp    1558 <sample_quad+0x3d8>
jmp    155d <sample_quad+0x3dd>
jmp    1562 <sample_quad+0x3e2>
jmp    1567 <sample_quad+0x3e7>
jmp    156c <sample_quad+0x3ec>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0xe0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sample_quad>:
push   %rbx
sub    $0x60,%rsp
movsd  %xmm4,0x50(%rsp)
movsd  %xmm3,0x38(%rsp)
movsd  0x70(%rsp),%xmm3
addsd  %xmm0,%xmm1
movsd  %xmm0,0x18(%rsp)
addsd  %xmm0,%xmm2
ucomisd %xmm1,%xmm2
jae    11cf <sample_quad+0x4f>
divsd  %xmm2,%xmm1
movapd %xmm6,%xmm0
addsd  %xmm3,%xmm0
addsd  %xmm7,%xmm5
movsd  %xmm5,0x20(%rsp)
movsd  %xmm0,0x8(%rsp)
movsd  %xmm2,0x30(%rsp)
jmp    11f1 <sample_quad+0x71>
movsd  %xmm1,0x30(%rsp)
divsd  %xmm1,%xmm2
addsd  %xmm7,%xmm5
movsd  %xmm5,0x8(%rsp)
addsd  %xmm3,%xmm6
movsd  %xmm6,0x20(%rsp)
movapd %xmm2,%xmm1
movapd %xmm1,%xmm0
movsd  %xmm1,0x48(%rsp)
addsd  0xdfd(%rip),%xmm0
call   10a0 <is_positive@plt>
test   %rax,%rax
je     13d1 <sample_quad+0x251>
movsd  0x20(%rsp),%xmm0
movsd  0xde1(%rip),%xmm4
addsd  %xmm4,%xmm0
movapd %xmm0,%xmm3
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,0x8(%rsp)
movsd  0x38(%rsp),%xmm1
addsd  0x50(%rsp),%xmm1
addsd  %xmm4,%xmm1
movsd  %xmm1,0x38(%rsp)
movsd  0x48(%rsp),%xmm0
call   1040 <maximum_point@plt>
movsd  %xmm0,0x28(%rsp)
movsd  0x48(%rsp),%xmm1
movsd  0x38(%rsp),%xmm2
movsd  0x8(%rsp),%xmm3
movsd  0x20(%rsp),%xmm4
call   1050 <DF@plt>
movaps %xmm0,0x50(%rsp)
movsd  0x30(%rsp),%xmm0
subsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x40(%rsp)
call   10a0 <is_positive@plt>
test   %rax,%rax
je     13d1 <sample_quad+0x251>
movapd 0x50(%rsp),%xmm1
xorpd  0xd5c(%rip),%xmm1
movsd  0x28(%rsp),%xmm0
mulsd  %xmm0,%xmm1
movapd %xmm1,0x50(%rsp)
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x28(%rsp)
call   10a0 <is_positive@plt>
test   %rax,%rax
je     13d1 <sample_quad+0x251>
movsd  0xd24(%rip),%xmm0
divsd  0x50(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
call   10a0 <is_positive@plt>
test   %rax,%rax
je     13d1 <sample_quad+0x251>
movsd  0x28(%rsp),%xmm0
call   1090 <sgamma@plt>
movsd  0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
movapd %xmm1,%xmm0
call   10a0 <is_positive@plt>
test   %rax,%rax
je     13d1 <sample_quad+0x251>
movsd  0x18(%rsp),%xmm0
divsd  0x40(%rsp),%xmm0
movsd  %xmm0,0x40(%rsp)
call   10a0 <is_positive@plt>
test   %rax,%rax
je     13af <sample_quad+0x22f>
movsd  0x10(%rsp),%xmm0
movsd  0x40(%rsp),%xmm1
movsd  0x48(%rsp),%xmm2
movsd  0x38(%rsp),%xmm3
movsd  0x8(%rsp),%xmm4
movsd  0x20(%rsp),%xmm5
movsd  0x28(%rsp),%xmm6
movapd 0x50(%rsp),%xmm7
call   1060 <qacc@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0xc7e(%rip),%xmm1
xorpd  %xmm0,%xmm0
call   1080 <genunf@plt>
call   1070 <log@plt>
mov    %rax,%rbx
xorpd  %xmm0,%xmm0
movsd  0x8(%rsp),%xmm1
call   1030 <my_fmin@plt>
cmp    %rax,%rbx
jge    13d1 <sample_quad+0x251>
movsd  0x30(%rsp),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm1,%xmm0
addsd  0xc41(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0x18(%rsp),%xmm0
add    $0x60,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define double @sample_quad(double %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8) {
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
store double %0, double* %11
store double %1, double* %12
store double %2, double* %13
store double %3, double* %14
store double %4, double* %15
store double %5, double* %16
store double %6, double* %17
store double %7, double* %18
store double %8, double* %19
%36 = load double, double* %11
store double %36, double* %20
%37 = load double, double* %12
%38 = load double, double* %11
%39 = fadd double %37, %38
store double %39, double* %21
%40 = load double, double* %13
%41 = load double, double* %20
%42 = fadd double %40, %41
store double %42, double* %22
%43 = load double, double* %21
%44 = load double, double* %22
%45 = fcmp ole double %43, %44
br i1 %45, label %46, label %64
46: ; preds = %9
%47 = load double, double* %21
store double %47, double* %23
%48 = load double, double* %22
store double %48, double* %24
%49 = load double, double* %24
%50 = load double, double* %23
%51 = fdiv double %49, %50
store double %51, double* %25
%52 = load double, double* %14
%53 = load double, double* %15
%54 = fadd double %52, %53
%55 = fsub double %54, 1.000000e+00
store double %55, double* %26
%56 = load double, double* %16
%57 = load double, double* %18
%58 = fadd double %56, %57
%59 = fsub double %58, 1.000000e+00
store double %59, double* %27
%60 = load double, double* %17
%61 = load double, double* %19
%62 = fadd double %60, %61
%63 = fsub double %62, 1.000000e+00
store double %63, double* %28
br label %82
64: ; preds = %9
%65 = load double, double* %22
store double %65, double* %23
%66 = load double, double* %21
store double %66, double* %24
%67 = load double, double* %24
%68 = load double, double* %23
%69 = fdiv double %67, %68
store double %69, double* %25
%70 = load double, double* %14
%71 = load double, double* %15
%72 = fadd double %70, %71
%73 = fsub double %72, 1.000000e+00
store double %73, double* %26
%74 = load double, double* %17
%75 = load double, double* %19
%76 = fadd double %74, %75
%77 = fsub double %76, 1.000000e+00
store double %77, double* %27
%78 = load double, double* %16
%79 = load double, double* %18
%80 = fadd double %78, %79
%81 = fsub double %80, 1.000000e+00
store double %81, double* %28
br label %82
82: ; preds = %64, %46
%83 = load double, double* %25
%84 = fsub double %83, 1.000000e+00
%85 = call i64 @is_positive(double %84)
%86 = icmp ne i64 %85, 0
br i1 %86, label %87, label %172
87: ; preds = %82
%88 = load double, double* %25
%89 = load double, double* %26
%90 = load double, double* %27
%91 = load double, double* %28
%92 = call double @maximum_point(double %88, double %89, double %90, double %91)
store double %92, double* %29
%93 = load double, double* %29
%94 = load double, double* %25
%95 = load double, double* %26
%96 = load double, double* %27
%97 = load double, double* %28
%98 = call double @DF(double %93, double %94, double %95, double %96, double %97)
%99 = fmul double -1.000000e+00, %98
%100 = load double, double* %29
%101 = fmul double %99, %100
store double %101, double* %31
%102 = load double, double* %31
%103 = load double, double* %29
%104 = fmul double %102, %103
store double %104, double* %30
%105 = load double, double* %23
%106 = load double, double* %11
%107 = fsub double %105, %106
%108 = call i64 @is_positive(double %107)
%109 = icmp ne i64 %108, 0
br i1 %109, label %110, label %171
110: ; preds = %87
%111 = load double, double* %11
%112 = load double, double* %23
%113 = load double, double* %11
%114 = fsub double %112, %113
%115 = fdiv double %111, %114
store double %115, double* %32
%116 = load double, double* %30
%117 = call i64 @is_positive(double %116)
%118 = icmp ne i64 %117, 0
br i1 %118, label %119, label %170
119: ; preds = %110
%120 = load double, double* %31
%121 = fdiv double 1.000000e+00, %120
%122 = call i64 @is_positive(double %121)
%123 = icmp ne i64 %122, 0
br i1 %123, label %124, label %170
124: ; preds = %119
%125 = load double, double* %31
%126 = fdiv double 1.000000e+00, %125
%127 = load double, double* %30
%128 = call double @sgamma(double %127)
%129 = fmul double %126, %128
store double %129, double* %33
%130 = load double, double* %33
%131 = call i64 @is_positive(double %130)
%132 = icmp ne i64 %131, 0
br i1 %132, label %133, label %169
133: ; preds = %124
%134 = load double, double* %32
%135 = call i64 @is_positive(double %134)
%136 = icmp ne i64 %135, 0
br i1 %136, label %144, label %137
137: ; preds = %133
%138 = load double, double* %23
%139 = load double, double* %33
%140 = fmul double %138, %139
%141 = load double, double* %33
%142 = fadd double 1.000000e+00, %141
%143 = fdiv double %140, %142
store double %143, double* %10
br label %174
144: ; preds = %133
%145 = load double, double* %33
%146 = load double, double* %32
%147 = load double, double* %25
%148 = load double, double* %26
%149 = load double, double* %27
%150 = load double, double* %28
%151 = load double, double* %30
%152 = load double, double* %31
%153 = call double @qacc(double %145, double %146, double %147, double %148, double %149, double %150, double %151, double %152)
store double %153, double* %34
%154 = call double @genunf(double 0.000000e+00, double 1.000000e+00)
store double %154, double* %35
%155 = load double, double* %35
%156 = call i64 @log(double %155)
%157 = load double, double* %34
%158 = call i64 @my_fmin(double 0.000000e+00, double %157)
%159 = icmp slt i64 %156, %158
br i1 %159, label %160, label %167
160: ; preds = %144
%161 = load double, double* %23
%162 = load double, double* %33
%163 = fmul double %161, %162
%164 = load double, double* %33
%165 = fadd double 1.000000e+00, %164
%166 = fdiv double %163, %165
store double %166, double* %10
br label %174
167: ; preds = %144
br label %168
168: ; preds = %167
br label %169
169: ; preds = %168, %124
br label %170
170: ; preds = %169, %119, %110
br label %171
171: ; preds = %170, %87
br label %172
172: ; preds = %171, %82
%173 = load double, double* %11
store double %173, double* %10
br label %174
174: ; preds = %172, %160, %137
%175 = load double, double* %10
ret double %175
}

LLVM IR (-O2):
define double @sample_quad(double %0, double %1, double %2, double %3, double %4, double %5, double %6, double %7, double %8) {
%10 = fadd double %0, %1
%11 = fadd double %0, %2
%12 = fcmp ugt double %10, %11
br i1 %12, label %17, label %13
13: ; preds = %9
%14 = fdiv double %11, %10
%15 = fadd double %5, %7
%16 = fadd double %6, %8
br label %21
17: ; preds = %9
%18 = fdiv double %10, %11
%19 = fadd double %6, %8
%20 = fadd double %5, %7
br label %21
21: ; preds = %17, %13
%22 = phi double [ %16, %13 ], [ %20, %17 ]
%23 = phi double [ %15, %13 ], [ %19, %17 ]
%24 = phi double [ %14, %13 ], [ %18, %17 ]
%25 = phi double [ %10, %13 ], [ %11, %17 ]
%26 = fadd double %22, -1.000000e+00
%27 = fadd double %23, -1.000000e+00
%28 = fadd double %3, %4
%29 = fadd double %28, -1.000000e+00
%30 = fadd double %24, -1.000000e+00
%31 = tail call i64 @is_positive(double %30)
%32 = icmp eq i64 %31, 0
br i1 %32, label %72, label %33
33: ; preds = %21
%34 = tail call double @maximum_point(double %24, double %29, double %27, double %26)
%35 = tail call double @DF(double %34, double %24, double %29, double %27, double %26)
%36 = fneg double %35
%37 = fmul double %34, %36
%38 = fmul double %34, %37
%39 = fsub double %25, %0
%40 = tail call i64 @is_positive(double %39)
%41 = icmp eq i64 %40, 0
br i1 %41, label %72, label %42
42: ; preds = %33
%43 = fdiv double %0, %39
%44 = tail call i64 @is_positive(double %38)
%45 = icmp eq i64 %44, 0
br i1 %45, label %72, label %46
46: ; preds = %42
%47 = fdiv double 1.000000e+00, %37
%48 = tail call i64 @is_positive(double %47)
%49 = icmp eq i64 %48, 0
br i1 %49, label %72, label %50
50: ; preds = %46
%51 = tail call double @sgamma(double %38)
%52 = fmul double %47, %51
%53 = tail call i64 @is_positive(double %52)
%54 = icmp eq i64 %53, 0
br i1 %54, label %72, label %55
55: ; preds = %50
%56 = tail call i64 @is_positive(double %43)
%57 = icmp eq i64 %56, 0
br i1 %57, label %58, label %62
58: ; preds = %55
%59 = fmul double %25, %52
%60 = fadd double %52, 1.000000e+00
%61 = fdiv double %59, %60
br label %72
62: ; preds = %55
%63 = tail call double @qacc(double %52, double %43, double %24, double %29, double %27, double %26, double %38, double %37)
%64 = tail call double @genunf(double 0.000000e+00, double 1.000000e+00)
%65 = tail call i64 @log(double %64)
%66 = tail call i64 @my_fmin(double 0.000000e+00, double %63)
%67 = icmp slt i64 %65, %66
br i1 %67, label %68, label %72
68: ; preds = %62
%69 = fmul double %25, %52
%70 = fadd double %52, 1.000000e+00
%71 = fdiv double %69, %70
br label %72
72: ; preds = %21, %42, %46, %62, %50, %33, %68, %58
%73 = phi double [ %71, %68 ], [ %61, %58 ], [ %0, %33 ], [ %0, %50 ], [ %0, %62 ], [ %0, %46 ], [ %0, %42 ], [ %0, %21 ]
ret double %73
}