#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 */
int abs(int);
int /*<<< orphan*/ assert(int);
double pow(double, double);
double sqrt(double);
int wiginit(int s1, int s2, int ntheta, const double *cos_theta, double *out) {
  double s12sign = ((s1 + s2) % 2) ? -1.0 : 1.0;
  double prefactor = 1.0;
  int abs_s1;
  int i, swap;
  if (abs(s1) > abs(s2)) {
    prefactor *= s12sign;
    swap = s1;
    s1 = s2;
    s2 = swap;
  }
  if (s2 < 0) {
    prefactor *= s12sign;
    s1 = -s1;
    s2 = -s2;
  }
  abs_s1 = abs(s1);
  assert(abs_s1 <= s2);
  for (i = 1; i <= s2 - abs_s1; i++)
    prefactor *= sqrt((double)(s2 + abs_s1 + i) / (double)i);
  for (i = 0; i < ntheta; i++)
    out[i] = prefactor *
             pow((1.0 + cos_theta[i]) / 2.0, 0.5 * (double)(s2 + s1)) *
             pow((1.0 - cos_theta[i]) / 2.0, 0.5 * (double)(s2 - s1));
  return s2;
}


x86 assembly (clang, -O0):
<wiginit>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x4(%rbp),%eax
add    -0x8(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
movsd  0xe99(%rip),%xmm0
movsd  0xe89(%rip),%xmm1
movsd  %xmm1,-0x50(%rbp)
cmp    $0x0,%edx
movsd  %xmm0,-0x48(%rbp)
jne    1194 <wiginit+0x54>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe5a(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x4(%rbp),%edi
call   1030 <abs@plt>
mov    %eax,-0x54(%rbp)
mov    -0x8(%rbp),%edi
call   1030 <abs@plt>
mov    %eax,%ecx
mov    -0x54(%rbp),%eax
cmp    %ecx,%eax
jle    11ec <wiginit+0xac>
movsd  -0x28(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x8(%rbp)
cmpl   $0x0,-0x8(%rbp)
jge    1215 <wiginit+0xd5>
movsd  -0x28(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
xor    %eax,%eax
sub    -0x4(%rbp),%eax
mov    %eax,-0x4(%rbp)
xor    %eax,%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x8(%rbp)
mov    -0x4(%rbp),%edi
call   1030 <abs@plt>
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x8(%rbp),%eax
setle  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jg     127e <wiginit+0x13e>
mov    -0x8(%rbp),%eax
add    -0x34(%rbp),%eax
add    -0x38(%rbp),%eax
cvtsi2sd %eax,%xmm0
cvtsi2sdl -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    123a <wiginit+0xfa>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    134b <wiginit+0x20b>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  0xd55(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
movsd  0xd60(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%eax
add    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm2
movsd  0xd42(%rip),%xmm1
mulsd  %xmm2,%xmm1
call   1050 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  0xd08(%rip),%xmm0
subsd  (%rax,%rcx,8),%xmm0
movsd  0xd13(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%eax
sub    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm2
movsd  0xcf5(%rip),%xmm1
mulsd  %xmm2,%xmm1
call   1050 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1285 <wiginit+0x145>
mov    -0x8(%rbp),%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<wiginit>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8,%r15
mov    %rcx,%r14
mov    %edx,0x10(%rsp)
mov    %edi,%r13d
lea    (%rsi,%r13,1),%eax
xor    %ecx,%ecx
test   $0x1,%al
sete   %cl
lea    0xea3(%rip),%rax
movsd  (%rax,%rcx,8),%xmm0
mov    %edi,%eax
neg    %eax
cmovs  %edi,%eax
mov    %esi,%ecx
neg    %ecx
cmovs  %esi,%ecx
movsd  0xe98(%rip),%xmm1
movapd %xmm0,%xmm2
cmp    %ecx,%eax
ja     1184 <wiginit+0x54>
movapd %xmm1,%xmm2
mov    %esi,%ecx
cmova  %r13d,%ecx
cmova  %esi,%r13d
mov    %ecx,%ebx
neg    %ebx
cmovs  %ecx,%ebx
mov    %r13d,%eax
neg    %eax
mov    %r13d,%r12d
cmovns %eax,%r12d
test   %ecx,%ecx
js     11a9 <wiginit+0x79>
movapd %xmm1,%xmm0
cmovs  %eax,%r13d
mulsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%rsp)
xor    %edi,%edi
cmp    %r12d,%ebx
setge  %dil
call   1030 <assert@plt>
mov    %ebx,%eax
sub    %r12d,%eax
mov    %rbx,0x30(%rsp)
jle    11fb <wiginit+0xcb>
lea    (%rbx,%r12,1),%ebp
add    $0x1,%r12d
cmp    %r12d,%ebx
jne    120a <wiginit+0xda>
mov    $0x1,%r12d
movsd  0x8(%rsp),%xmm2
mov    0x10(%rsp),%ebx
test   $0x1,%al
jne    12e3 <wiginit+0x1b3>
jmp    131e <wiginit+0x1ee>
movsd  0x8(%rsp),%xmm2
mov    0x10(%rsp),%ebx
jmp    131e <wiginit+0x1ee>
mov    %r15,0x18(%rsp)
mov    %eax,0x2c(%rsp)
mov    %eax,%r15d
and    $0xfffffffe,%r15d
mov    $0x1,%eax
xorpd  %xmm3,%xmm3
movsd  0x8(%rsp),%xmm2
jmp    125b <wiginit+0x12b>
nopl   0x0(%rax,%rax,1)
movsd  %xmm4,0x20(%rsp)
call   1050 <sqrt@plt>
movsd  0x20(%rsp),%xmm4
xorpd  %xmm3,%xmm3
movsd  0x8(%rsp),%xmm2
mulsd  %xmm4,%xmm2
mulsd  %xmm0,%xmm2
lea    0x1(%rbx),%eax
cmp    %r15d,%ebx
je     12ce <wiginit+0x19e>
mov    %eax,%r12d
lea    (%r12,%rbp,1),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r12d,%xmm1
divsd  %xmm1,%xmm0
ucomisd %xmm3,%xmm0
movsd  %xmm2,0x8(%rsp)
jb     1290 <wiginit+0x160>
xorps  %xmm4,%xmm4
sqrtsd %xmm0,%xmm4
jmp    129d <wiginit+0x16d>
nopw   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
xorpd  %xmm3,%xmm3
movapd %xmm0,%xmm4
lea    (%r12,%rbp,1),%eax
add    $0x1,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
lea    0x1(%r12),%ebx
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
divsd  %xmm1,%xmm0
ucomisd %xmm3,%xmm0
jb     1230 <wiginit+0x100>
sqrtsd %xmm0,%xmm0
jmp    1245 <wiginit+0x115>
add    $0x2,%r12d
mov    0x18(%rsp),%r15
mov    0x10(%rsp),%ebx
mov    0x2c(%rsp),%eax
test   $0x1,%al
je     131e <wiginit+0x1ee>
add    %r12d,%ebp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r12d,%xmm1
divsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1309 <wiginit+0x1d9>
sqrtsd %xmm0,%xmm0
jmp    131a <wiginit+0x1ea>
movsd  %xmm2,0x8(%rsp)
call   1050 <sqrt@plt>
movsd  0x8(%rsp),%xmm2
mulsd  %xmm0,%xmm2
test   %ebx,%ebx
jle    14ae <wiginit+0x37e>
mov    0x30(%rsp),%rcx
lea    (%rcx,%r13,1),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0xcda(%rip),%xmm0
mov    %ecx,%eax
sub    %r13d,%eax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm3
mov    %ebx,%r12d
cmp    $0x1,%ebx
movsd  %xmm3,0x20(%rsp)
movsd  %xmm2,0x8(%rsp)
jne    136d <wiginit+0x23d>
xor    %ebp,%ebp
jmp    144c <wiginit+0x31c>
mov    %r12d,%r13d
and    $0xfffffffe,%r13d
xor    %ebp,%ebp
movsd  %xmm1,0x10(%rsp)
nopl   0x0(%rax)
movsd  (%r14,%rbp,8),%xmm0
movsd  0xc82(%rip),%xmm2
addsd  %xmm2,%xmm0
movsd  0xc7e(%rip),%xmm2
mulsd  %xmm2,%xmm0
call   1040 <pow@plt>
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0xc59(%rip),%xmm0
subsd  (%r14,%rbp,8),%xmm0
mulsd  0xc53(%rip),%xmm0
movsd  0x20(%rsp),%xmm1
call   1040 <pow@plt>
mulsd  0x18(%rsp),%xmm0
movsd  %xmm0,(%r15,%rbp,8)
movsd  0x8(%r14,%rbp,8),%xmm0
addsd  0xc25(%rip),%xmm0
mulsd  0xc25(%rip),%xmm0
movsd  0x10(%rsp),%xmm1
call   1040 <pow@plt>
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  0xbfe(%rip),%xmm0
subsd  0x8(%r14,%rbp,8),%xmm0
mulsd  0xbf7(%rip),%xmm0
movsd  0x20(%rsp),%xmm1
call   1040 <pow@plt>
movsd  0x10(%rsp),%xmm1
mulsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x8(%r15,%rbp,8)
add    $0x2,%rbp
cmp    %rbp,%r13
jne    1380 <wiginit+0x250>
movsd  0xbc4(%rip),%xmm3
test   $0x1,%r12b
je     14ae <wiginit+0x37e>
movsd  (%r14,%rbp,8),%xmm0
addsd  0xba8(%rip),%xmm0
mulsd  %xmm3,%xmm0
call   1040 <pow@plt>
movsd  0x8(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  0xb87(%rip),%xmm0
subsd  (%r14,%rbp,8),%xmm0
mulsd  0xb81(%rip),%xmm0
movsd  0x20(%rsp),%xmm1
call   1040 <pow@plt>
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%r15,%rbp,8)
mov    0x30(%rsp),%rax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @wiginit(i32 %0, i32 %1, i32 %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca double
%12 = alloca double
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double* %3, double** %9
store double* %4, double** %10
%16 = load i32, i32* %6
%17 = load i32, i32* %7
%18 = add nsw i32 %16, %17
%19 = srem i32 %18, 2
%20 = icmp ne i32 %19, 0
%21 = zext i1 %20 to i64
%22 = select i1 %20, double -1.000000e+00, double 1.000000e+00
store double %22, double* %11
store double 1.000000e+00, double* %12
%23 = load i32, i32* %6
%24 = call i32 @abs(i32 %23)
%25 = load i32, i32* %7
%26 = call i32 @abs(i32 %25)
%27 = icmp sgt i32 %24, %26
br i1 %27, label %28, label %35
28: ; preds = %5
%29 = load double, double* %11
%30 = load double, double* %12
%31 = fmul double %30, %29
store double %31, double* %12
%32 = load i32, i32* %6
store i32 %32, i32* %15
%33 = load i32, i32* %7
store i32 %33, i32* %6
%34 = load i32, i32* %15
store i32 %34, i32* %7
br label %35
35: ; preds = %28, %5
%36 = load i32, i32* %7
%37 = icmp slt i32 %36, 0
br i1 %37, label %38, label %46
38: ; preds = %35
%39 = load double, double* %11
%40 = load double, double* %12
%41 = fmul double %40, %39
store double %41, double* %12
%42 = load i32, i32* %6
%43 = sub nsw i32 0, %42
store i32 %43, i32* %6
%44 = load i32, i32* %7
%45 = sub nsw i32 0, %44
store i32 %45, i32* %7
br label %46
46: ; preds = %38, %35
%47 = load i32, i32* %6
%48 = call i32 @abs(i32 %47)
store i32 %48, i32* %13
%49 = load i32, i32* %13
%50 = load i32, i32* %7
%51 = icmp sle i32 %49, %50
%52 = zext i1 %51 to i32
%53 = call i32 @assert(i32 %52)
store i32 1, i32* %14
br label %54
54: ; preds = %73, %46
%55 = load i32, i32* %14
%56 = load i32, i32* %7
%57 = load i32, i32* %13
%58 = sub nsw i32 %56, %57
%59 = icmp sle i32 %55, %58
br i1 %59, label %60, label %76
60: ; preds = %54
%61 = load i32, i32* %7
%62 = load i32, i32* %13
%63 = add nsw i32 %61, %62
%64 = load i32, i32* %14
%65 = add nsw i32 %63, %64
%66 = sitofp i32 %65 to double
%67 = load i32, i32* %14
%68 = sitofp i32 %67 to double
%69 = fdiv double %66, %68
%70 = call double @sqrt(double %69)
%71 = load double, double* %12
%72 = fmul double %71, %70
store double %72, double* %12
br label %73
73: ; preds = %60
%74 = load i32, i32* %14
%75 = add nsw i32 %74, 1
store i32 %75, i32* %14
br label %54
76: ; preds = %54
store i32 0, i32* %14
br label %77
77: ; preds = %115, %76
%78 = load i32, i32* %14
%79 = load i32, i32* %8
%80 = icmp slt i32 %78, %79
br i1 %80, label %81, label %118
81: ; preds = %77
%82 = load double, double* %12
%83 = load double*, double** %9
%84 = load i32, i32* %14
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = fadd double 1.000000e+00, %87
%89 = fdiv double %88, 2.000000e+00
%90 = load i32, i32* %7
%91 = load i32, i32* %6
%92 = add nsw i32 %90, %91
%93 = sitofp i32 %92 to double
%94 = fmul double 5.000000e-01, %93
%95 = call double @pow(double %89, double %94)
%96 = fmul double %82, %95
%97 = load double*, double** %9
%98 = load i32, i32* %14
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %97, i64 %99
%101 = load double, double* %100
%102 = fsub double 1.000000e+00, %101
%103 = fdiv double %102, 2.000000e+00
%104 = load i32, i32* %7
%105 = load i32, i32* %6
%106 = sub nsw i32 %104, %105
%107 = sitofp i32 %106 to double
%108 = fmul double 5.000000e-01, %107
%109 = call double @pow(double %103, double %108)
%110 = fmul double %96, %109
%111 = load double*, double** %10
%112 = load i32, i32* %14
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %111, i64 %113
store double %110, double* %114
br label %115
115: ; preds = %81
%116 = load i32, i32* %14
%117 = add nsw i32 %116, 1
store i32 %117, i32* %14
br label %77
118: ; preds = %77
%119 = load i32, i32* %7
ret i32 %119
}

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