#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 M_PI;
int exp(double);
int pow(double, double);
int sqrt(int);
double sho_wf(double r, double mw, int n, int l) {
  int i, npl;
  double x, mwr2, cur, last, last2, factor;
  if ((r < 0) || (mw < 0) || (n < 0) || (l < 0))
    return 0.;
  if (r == 0) {
    if (l != 0)
      return 0.;
    factor = mw;
    for (i = 1; i <= n; i++)
      factor *= (i + 0.5) / i;
    return 2 * sqrt(sqrt(mw / M_PI) * factor);
  }
  x = sqrt(mw) * r;
  mwr2 = x * x;
  npl = n + l;
  factor = mw;
  for (i = 3; i <= 2 * npl + 1; i += 2)
    factor *= 2 * mwr2 / i;
  last2 = sqrt(factor);
  if (l == npl)
    cur = last2;
  last = ((npl + 0.5) / x - x) * last2;
  if (l == npl - 1)
    cur = last;
  for (i = npl - 2; i >= l; i--) {
    cur = (((i + 1.5) / x - x) * last - sqrt(npl - i - 1) * last2) /
          sqrt(npl - i);
    last2 = last;
    last = cur;
  }
  return 2 * pow(mw / M_PI, 0.25) * exp(-0.5 * mwr2) * cur;
}


x86 assembly (clang, -O0):
<sho_wf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %edi,-0x1c(%rbp)
mov    %esi,-0x20(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
ja     117b <sho_wf+0x4b>
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
ja     117b <sho_wf+0x4b>
cmpl   $0x0,-0x1c(%rbp)
jl     117b <sho_wf+0x4b>
cmpl   $0x0,-0x20(%rbp)
jge    1188 <sho_wf+0x58>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    143c <sho_wf+0x30c>
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1242 <sho_wf+0x112>
jp     1242 <sho_wf+0x112>
cmpl   $0x0,-0x20(%rbp)
je     11b7 <sho_wf+0x87>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    143c <sho_wf+0x30c>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1206 <sho_wf+0xd6>
cvtsi2sdl -0x24(%rbp),%xmm0
movsd  0xe1f(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11c8 <sho_wf+0x98>
movsd  -0x18(%rbp),%xmm0
mov    0x2dd6(%rip),%rax
divsd  (%rax),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
cvtsi2sd %eax,%xmm0
mulsd  -0x58(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
shl    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    143c <sho_wf+0x30c>
cvttsd2si -0x18(%rbp),%edi
call   1050 <sqrt@plt>
cvtsi2sd %eax,%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x3,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     12c5 <sho_wf+0x195>
movsd  0xd81(%rip),%xmm0
mulsd  -0x38(%rbp),%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    1283 <sho_wf+0x153>
cvttsd2si -0x58(%rbp),%edi
call   1050 <sqrt@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x28(%rbp),%eax
jne    12ee <sho_wf+0x1be>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x28(%rbp),%xmm0
movsd  0xd05(%rip),%xmm1
addsd  %xmm1,%xmm0
divsd  -0x30(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    132e <sho_wf+0x1fe>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jl     13f1 <sho_wf+0x2c1>
cvtsi2sdl -0x24(%rbp),%xmm0
movsd  0xcc8(%rip),%xmm1
addsd  %xmm1,%xmm0
divsd  -0x30(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x28(%rbp),%edi
sub    -0x24(%rbp),%edi
sub    $0x1,%edi
call   1050 <sqrt@plt>
movsd  -0x70(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
cvtsi2sd %eax,%xmm1
mulsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%edi
sub    -0x24(%rbp),%edi
call   1050 <sqrt@plt>
movsd  -0x60(%rbp),%xmm0
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1337 <sho_wf+0x207>
movsd  -0x18(%rbp),%xmm0
mov    0x2beb(%rip),%rax
divsd  (%rax),%xmm0
movsd  0xc07(%rip),%xmm1
call   1040 <pow@plt>
shl    $0x1,%eax
mov    %eax,-0x74(%rbp)
movsd  0xbec(%rip),%xmm0
mulsd  -0x38(%rbp),%xmm0
call   1030 <exp@plt>
mov    %eax,%ecx
mov    -0x74(%rbp),%eax
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sho_wf>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x30,%rsp
movapd %xmm0,%xmm3
xorpd  %xmm0,%xmm0
ucomisd %xmm3,%xmm0
ja     1410 <sho_wf+0x2e0>
xorpd  %xmm2,%xmm2
ucomisd %xmm1,%xmm2
ja     1410 <sho_wf+0x2e0>
mov    %esi,%r12d
mov    %edi,%r14d
mov    %esi,%eax
or     %edi,%eax
js     1410 <sho_wf+0x2e0>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm3
jne    11a9 <sho_wf+0x79>
jp     11a9 <sho_wf+0x79>
test   %r12d,%r12d
jne    1410 <sho_wf+0x2e0>
movapd %xmm1,%xmm4
test   %r14d,%r14d
jle    13d9 <sho_wf+0x2a9>
cmp    $0x1,%r14d
jne    1364 <sho_wf+0x234>
movsd  0xe90(%rip),%xmm0
movapd %xmm1,%xmm4
jmp    13bf <sho_wf+0x28f>
movsd  %xmm1,0x8(%rsp)
cvttsd2si %xmm1,%r15d
mov    %r15d,%edi
movsd  %xmm3,(%rsp)
call   1050 <sqrt@plt>
cvtsi2sd %eax,%xmm5
mulsd  (%rsp),%xmm5
movsd  %xmm5,0x10(%rsp)
mulsd  %xmm5,%xmm5
lea    (%r12,%r14,1),%ebp
lea    0x1(%rbp,%rbp,1),%eax
cmp    $0x3,%eax
jl     122c <sho_wf+0xfc>
movapd %xmm5,%xmm0
addsd  %xmm5,%xmm0
mov    $0x3,%ecx
movsd  0xe0a(%rip),%xmm1
movsd  0xe0a(%rip),%xmm3
movsd  0x8(%rsp),%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm0,%xmm4
divsd  %xmm1,%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm1
add    $0x2,%ecx
cmp    %eax,%ecx
jle    1210 <sho_wf+0xe0>
cvttsd2si %xmm2,%r15d
movsd  %xmm5,0x20(%rsp)
mov    %r15d,%edi
call   1050 <sqrt@plt>
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
addsd  0xdc0(%rip),%xmm1
movsd  0x10(%rsp),%xmm0
divsd  %xmm0,%xmm1
subsd  %xmm0,%xmm1
mulsd  %xmm3,%xmm1
lea    -0x1(%rbp),%eax
movapd %xmm1,%xmm2
cmp    %r12d,%eax
je     1272 <sho_wf+0x142>
movapd %xmm3,%xmm2
add    $0xfffffffe,%ebp
cmp    %r12d,%ebp
jl     1316 <sho_wf+0x1e6>
lea    (%r12,%r14,1),%ebx
add    $0xffffffff,%ebx
mov    $0x1,%edi
nopw   0x0(%rax,%rax,1)
movsd  %xmm3,0x28(%rsp)
movsd  %xmm1,0x18(%rsp)
add    $0xffffffff,%ebx
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
addsd  0xd6a(%rip),%xmm1
divsd  %xmm0,%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
lea    0x1(%rdi),%ebp
call   1050 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0x28(%rsp),%xmm0
movsd  (%rsp),%xmm1
mulsd  0x18(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
mov    %ebp,%edi
call   1050 <sqrt@plt>
movsd  (%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
divsd  %xmm0,%xmm2
movsd  0x10(%rsp),%xmm0
mov    %ebp,%edi
movsd  0x18(%rsp),%xmm3
movapd %xmm2,%xmm1
cmp    %r12d,%ebx
jg     1290 <sho_wf+0x160>
movsd  %xmm2,(%rsp)
mov    0x2cc6(%rip),%rax
movsd  0x8(%rsp),%xmm0
divsd  (%rax),%xmm0
movsd  0xcec(%rip),%xmm1
call   1040 <pow@plt>
mov    %eax,%ebx
movsd  0x20(%rsp),%xmm0
mulsd  0xcdf(%rip),%xmm0
call   1030 <exp@plt>
imul   %ebx,%eax
add    %eax,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  (%rsp),%xmm0
jmp    1410 <sho_wf+0x2e0>
mov    %r14d,%eax
and    $0xfffffffe,%eax
mov    $0x1,%ecx
movsd  0xc99(%rip),%xmm0
movapd %xmm1,%xmm4
nopl   0x0(%rax,%rax,1)
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
movapd %xmm3,%xmm2
addsd  %xmm0,%xmm2
divsd  %xmm3,%xmm2
mulsd  %xmm4,%xmm2
lea    0x1(%rcx),%edx
xorps  %xmm3,%xmm3
cvtsi2sd %edx,%xmm3
movapd %xmm3,%xmm4
addsd  %xmm0,%xmm4
divsd  %xmm3,%xmm4
mulsd  %xmm2,%xmm4
add    $0x2,%ecx
cmp    %eax,%edx
jne    1380 <sho_wf+0x250>
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
test   $0x1,%r14b
je     13d9 <sho_wf+0x2a9>
movsd  0xc43(%rip),%xmm2
addsd  %xmm0,%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm2,%xmm4
movsd  %xmm4,(%rsp)
mov    0x2c03(%rip),%rax
divsd  (%rax),%xmm1
cvttsd2si %xmm1,%edi
call   1050 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  (%rsp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sqrt@plt>
add    %eax,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
add    $0x30,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define double @sho_wf(double %0, double %1, i32 %2, i32 %3) {
%5 = fcmp olt double %0, 0.000000e+00
%6 = fcmp olt double %1, 0.000000e+00
%7 = or i1 %5, %6
%8 = or i32 %3, %2
%9 = icmp slt i32 %8, 0
%10 = or i1 %7, %9
br i1 %10, label %128, label %11
11: ; preds = %4
%12 = fcmp oeq double %0, 0.000000e+00
br i1 %12, label %13, label %61
13: ; preds = %11
%14 = icmp eq i32 %3, 0
br i1 %14, label %15, label %128
15: ; preds = %13
%16 = icmp slt i32 %2, 1
br i1 %16, label %49, label %17
17: ; preds = %15
%18 = and i32 %2, 1
%19 = icmp eq i32 %2, 1
br i1 %19, label %40, label %20
20: ; preds = %17
%21 = and i32 %2, -2
br label %22
22: ; preds = %22, %20
%23 = phi double [ %1, %20 ], [ %34, %22 ]
%24 = phi i32 [ 1, %20 ], [ %35, %22 ]
%25 = phi i32 [ 0, %20 ], [ %36, %22 ]
%26 = sitofp i32 %24 to double
%27 = fadd double %26, 5.000000e-01
%28 = fdiv double %27, %26
%29 = fmul double %23, %28
%30 = add nuw i32 %24, 1
%31 = sitofp i32 %30 to double
%32 = fadd double %31, 5.000000e-01
%33 = fdiv double %32, %31
%34 = fmul double %29, %33
%35 = add nuw i32 %24, 2
%36 = add i32 %25, 2
%37 = icmp eq i32 %36, %21
br i1 %37, label %38, label %22
38: ; preds = %22
%39 = sitofp i32 %35 to double
br label %40
40: ; preds = %38, %17
%41 = phi double [ undef, %17 ], [ %34, %38 ]
%42 = phi double [ %1, %17 ], [ %34, %38 ]
%43 = phi double [ 1.000000e+00, %17 ], [ %39, %38 ]
%44 = icmp eq i32 %18, 0
br i1 %44, label %49, label %45
45: ; preds = %40
%46 = fadd double %43, 5.000000e-01
%47 = fdiv double %46, %43
%48 = fmul double %42, %47
br label %49
49: ; preds = %45, %40, %15
%50 = phi double [ %1, %15 ], [ %41, %40 ], [ %48, %45 ]
%51 = load double, double* @M_PI
%52 = fdiv double %1, %51
%53 = fptosi double %52 to i32
%54 = tail call i32 @sqrt(i32 %53)
%55 = sitofp i32 %54 to double
%56 = fmul double %50, %55
%57 = fptosi double %56 to i32
%58 = tail call i32 @sqrt(i32 %57)
%59 = shl nsw i32 %58, 1
%60 = sitofp i32 %59 to double
br label %128
61: ; preds = %11
%62 = fptosi double %1 to i32
%63 = tail call i32 @sqrt(i32 %62)
%64 = sitofp i32 %63 to double
%65 = fmul double %64, %0
%66 = fmul double %65, %65
%67 = add nsw i32 %3, %2
%68 = shl nsw i32 %67, 1
%69 = or i32 %68, 1
%70 = icmp slt i32 %69, 3
br i1 %70, label %83, label %71
71: ; preds = %61
%72 = fmul double %66, 2.000000e+00
br label %73
73: ; preds = %71, %73
%74 = phi double [ %1, %71 ], [ %78, %73 ]
%75 = phi i32 [ 3, %71 ], [ %79, %73 ]
%76 = sitofp i32 %75 to double
%77 = fdiv double %72, %76
%78 = fmul double %74, %77
%79 = add nuw nsw i32 %75, 2
%80 = icmp sgt i32 %79, %69
br i1 %80, label %81, label %73
81: ; preds = %73
%82 = fptosi double %78 to i32
br label %83
83: ; preds = %81, %61
%84 = phi i32 [ %82, %81 ], [ %62, %61 ]
%85 = tail call i32 @sqrt(i32 %84)
%86 = sitofp i32 %85 to double
%87 = sitofp i32 %67 to double
%88 = fadd double %87, 5.000000e-01
%89 = fdiv double %88, %65
%90 = fsub double %89, %65
%91 = fmul double %90, %86
%92 = add nsw i32 %67, -1
%93 = icmp eq i32 %92, %3
%94 = select i1 %93, double %91, double %86
%95 = add nsw i32 %67, -2
%96 = icmp slt i32 %95, %3
br i1 %96, label %117, label %97
97: ; preds = %83, %97
%98 = phi double [ %99, %97 ], [ %86, %83 ]
%99 = phi double [ %114, %97 ], [ %91, %83 ]
%100 = phi i32 [ %115, %97 ], [ %95, %83 ]
%101 = sitofp i32 %100 to double
%102 = fadd double %101, 1.500000e+00
%103 = fdiv double %102, %65
%104 = fsub double %103, %65
%105 = sub nsw i32 %67, %100
%106 = add nsw i32 %105, -1
%107 = tail call i32 @sqrt(i32 %106)
%108 = sitofp i32 %107 to double
%109 = fneg double %98
%110 = fmul double %109, %108
%111 = tail call double @llvm.fmuladd.f64(double %104, double %99, double %110)
%112 = tail call i32 @sqrt(i32 %105)
%113 = sitofp i32 %112 to double
%114 = fdiv double %111, %113
%115 = add nsw i32 %100, -1
%116 = icmp sgt i32 %100, %3
br i1 %116, label %97, label %117
117: ; preds = %97, %83
%118 = phi double [ %94, %83 ], [ %114, %97 ]
%119 = load double, double* @M_PI
%120 = fdiv double %1, %119
%121 = tail call i32 @pow(double %120, double 2.500000e-01)
%122 = shl nsw i32 %121, 1
%123 = fmul double %66, -5.000000e-01
%124 = tail call i32 @exp(double %123)
%125 = mul nsw i32 %122, %124
%126 = sitofp i32 %125 to double
%127 = fmul double %118, %126
br label %128
128: ; preds = %13, %4, %117, %49
%129 = phi double [ %60, %49 ], [ %127, %117 ], [ 0.000000e+00, %4 ], [ 0.000000e+00, %13 ]
ret double %129
}