#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 */
scalar_t__ exp(double);
double log(double);
void self_adj_theta_norm(double *theta, double *freq_ref, double *freq_est,
                         int grid_size, double norm_const) {
  int i;
  int N_empty;
  double sumv;
  double maxv;
  double pr_e;
  pr_e = 0.0;
  N_empty = 0;
  for (i = 1; i <= grid_size + 1; i++)
    if (freq_est[i] <= 0.0) {
      pr_e += freq_ref[i];
      N_empty++;
    }
  pr_e = pr_e / (grid_size + 1 - N_empty);
  for (i = 1; i <= grid_size + 1; i++)
    theta[i] += log(freq_ref[i] + pr_e);
  if (norm_const > 0.0) {
    maxv = theta[1];
    for (i = 1; i <= grid_size + 1; i++)
      if (theta[i] > maxv)
        maxv = theta[i];
    sumv = 0.0;
    for (i = 1; i <= grid_size + 1; i++) {
      theta[i] -= maxv;
      sumv += exp(theta[i]);
    }
    for (i = 1; i <= grid_size + 1; i++)
      theta[i] += (-log(sumv) + log(norm_const));
  }
}


x86 assembly (clang, -O0):
<self_adj_theta_norm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     11ac <self_adj_theta_norm+0x8c>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jb     1199 <self_adj_theta_norm+0x79>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    119e <self_adj_theta_norm+0x7e>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1152 <self_adj_theta_norm+0x32>
movsd  -0x48(%rbp),%xmm0
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
sub    -0x30(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1216 <self_adj_theta_norm+0xf6>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
call   1040 <log@plt>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ce <self_adj_theta_norm+0xae>
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    136e <self_adj_theta_norm+0x24e>
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     128b <self_adj_theta_norm+0x16b>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1278 <self_adj_theta_norm+0x158>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    127d <self_adj_theta_norm+0x15d>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    123d <self_adj_theta_norm+0x11d>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     12f5 <self_adj_theta_norm+0x1d5>
movsd  -0x40(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <exp@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    129a <self_adj_theta_norm+0x17a>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1369 <self_adj_theta_norm+0x249>
movsd  -0x38(%rbp),%xmm0
call   1040 <log@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1040 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12fc <self_adj_theta_norm+0x1dc>
jmp    136e <self_adj_theta_norm+0x24e>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<self_adj_theta_norm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
movsd  %xmm0,0x10(%rsp)
mov    %ecx,%r14d
mov    %rdi,%rbx
lea    0x1(%r14),%r12d
test   %ecx,%ecx
js     1227 <self_adj_theta_norm+0x107>
mov    %rsi,%r15
lea    0x2(%r14),%r13d
lea    -0x1(%r13),%rcx
test   %r14d,%r14d
je     11b6 <self_adj_theta_norm+0x96>
mov    %rcx,%rdi
and    $0xfffffffffffffffe,%rdi
neg    %rdi
xor    %eax,%eax
xorpd  %xmm0,%xmm0
mov    $0x1,%esi
xorpd  %xmm1,%xmm1
jmp    1192 <self_adj_theta_norm+0x72>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%rdi,%rsi,1),%rbp
add    $0x2,%rbp
add    $0x2,%rsi
cmp    $0x1,%rbp
je     11c1 <self_adj_theta_norm+0xa1>
ucomisd (%rdx,%rsi,8),%xmm0
jb     11a2 <self_adj_theta_norm+0x82>
addsd  (%r15,%rsi,8),%xmm1
add    $0x1,%eax
ucomisd 0x8(%rdx,%rsi,8),%xmm0
jb     1180 <self_adj_theta_norm+0x60>
addsd  0x8(%r15,%rsi,8),%xmm1
add    $0x1,%eax
jmp    1180 <self_adj_theta_norm+0x60>
xor    %eax,%eax
xorpd  %xmm1,%xmm1
mov    $0x1,%esi
test   $0x1,%cl
je     11da <self_adj_theta_norm+0xba>
xorpd  %xmm0,%xmm0
ucomisd (%rdx,%rsi,8),%xmm0
jb     11da <self_adj_theta_norm+0xba>
addsd  (%r15,%rsi,8),%xmm1
add    $0x1,%eax
test   %r14d,%r14d
js     1227 <self_adj_theta_norm+0x107>
mov    %r12d,%ecx
sub    %eax,%ecx
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
divsd  %xmm0,%xmm1
mov    $0x1,%ebp
movsd  %xmm1,(%rsp)
nopl   0x0(%rax)
movsd  (%r15,%rbp,8),%xmm0
addsd  %xmm1,%xmm0
call   1040 <log@plt>
movsd  (%rsp),%xmm1
addsd  (%rbx,%rbp,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r13
jne    1200 <self_adj_theta_norm+0xe0>
xorpd  %xmm0,%xmm0
movsd  0x10(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jbe    13b4 <self_adj_theta_norm+0x294>
test   %r14d,%r14d
js     13b4 <self_adj_theta_norm+0x294>
movsd  0x8(%rbx),%xmm0
movapd %xmm0,%xmm2
je     12f6 <self_adj_theta_norm+0x1d6>
mov    %r12d,%edx
lea    -0x1(%rdx),%rcx
add    $0xfffffffffffffffe,%rdx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x3,%rdx
jae    1274 <self_adj_theta_norm+0x154>
mov    $0x2,%edx
movapd %xmm0,%xmm1
jmp    12c9 <self_adj_theta_norm+0x1a9>
and    $0xfffffffffffffffc,%rcx
neg    %rcx
mov    $0x2,%edx
movapd %xmm0,%xmm1
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx,%rdx,8),%xmm2
maxsd  %xmm1,%xmm2
movsd  0x8(%rbx,%rdx,8),%xmm1
maxsd  %xmm2,%xmm1
movsd  0x10(%rbx,%rdx,8),%xmm2
maxsd  %xmm1,%xmm2
movsd  0x18(%rbx,%rdx,8),%xmm1
maxsd  %xmm2,%xmm1
lea    (%rcx,%rdx,1),%rsi
add    $0x4,%rsi
add    $0x4,%rdx
cmp    $0x2,%rsi
jne    1290 <self_adj_theta_norm+0x170>
movapd %xmm1,%xmm2
test   %rax,%rax
je     12f6 <self_adj_theta_norm+0x1d6>
lea    (%rbx,%rdx,8),%rcx
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rdx,8),%xmm2
maxsd  %xmm1,%xmm2
add    $0x1,%rdx
movapd %xmm2,%xmm1
cmp    %rdx,%rax
jne    12e0 <self_adj_theta_norm+0x1c0>
test   %r14d,%r14d
js     13b4 <self_adj_theta_norm+0x294>
lea    0x2(%r14),%r15d
movsd  %xmm2,0x8(%rsp)
subsd  %xmm2,%xmm0
movsd  %xmm0,0x8(%rbx)
call   1030 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
test   %r14d,%r14d
je     136e <self_adj_theta_norm+0x24e>
mov    $0x2,%ebp
nopl   0x0(%rax)
movsd  %xmm0,(%rsp)
movsd  (%rbx,%rbp,8),%xmm0
subsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
call   1030 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  (%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  (%rsp),%xmm0
add    $0x1,%rbp
cmp    %rbp,%r15
jne    1330 <self_adj_theta_norm+0x210>
test   %r14d,%r14d
js     13b4 <self_adj_theta_norm+0x294>
mov    $0x1,%ebp
movsd  %xmm0,(%rsp)
nopl   (%rax)
movsd  (%rsp),%xmm0
call   1040 <log@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <log@plt>
subsd  0x8(%rsp),%xmm0
addsd  (%rbx,%rbp,8),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r15
jne    1380 <self_adj_theta_norm+0x260>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @self_adj_theta_norm(double* %0, double* %1, double* %2, i32 %3, double %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca double
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double
%15 = alloca double
store double* %0, double** %6
store double* %1, double** %7
store double* %2, double** %8
store i32 %3, i32* %9
store double %4, double* %10
store double 0.000000e+00, double* %15
store i32 0, i32* %12
store i32 1, i32* %11
br label %16
16: ; preds = %39, %5
%17 = load i32, i32* %11
%18 = load i32, i32* %9
%19 = add nsw i32 %18, 1
%20 = icmp sle i32 %17, %19
br i1 %20, label %21, label %42
21: ; preds = %16
%22 = load double*, double** %8
%23 = load i32, i32* %11
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
%26 = load double, double* %25
%27 = fcmp ole double %26, 0.000000e+00
br i1 %27, label %28, label %38
28: ; preds = %21
%29 = load double*, double** %7
%30 = load i32, i32* %11
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %29, i64 %31
%33 = load double, double* %32
%34 = load double, double* %15
%35 = fadd double %34, %33
store double %35, double* %15
%36 = load i32, i32* %12
%37 = add nsw i32 %36, 1
store i32 %37, i32* %12
br label %38
38: ; preds = %28, %21
br label %39
39: ; preds = %38
%40 = load i32, i32* %11
%41 = add nsw i32 %40, 1
store i32 %41, i32* %11
br label %16
42: ; preds = %16
%43 = load double, double* %15
%44 = load i32, i32* %9
%45 = add nsw i32 %44, 1
%46 = load i32, i32* %12
%47 = sub nsw i32 %45, %46
%48 = sitofp i32 %47 to double
%49 = fdiv double %43, %48
store double %49, double* %15
store i32 1, i32* %11
br label %50
50: ; preds = %70, %42
%51 = load i32, i32* %11
%52 = load i32, i32* %9
%53 = add nsw i32 %52, 1
%54 = icmp sle i32 %51, %53
br i1 %54, label %55, label %73
55: ; preds = %50
%56 = load double*, double** %7
%57 = load i32, i32* %11
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %56, i64 %58
%60 = load double, double* %59
%61 = load double, double* %15
%62 = fadd double %60, %61
%63 = call double @log(double %62)
%64 = load double*, double** %6
%65 = load i32, i32* %11
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %64, i64 %66
%68 = load double, double* %67
%69 = fadd double %68, %63
store double %69, double* %67
br label %70
70: ; preds = %55
%71 = load i32, i32* %11
%72 = add nsw i32 %71, 1
store i32 %72, i32* %11
br label %50
73: ; preds = %50
%74 = load double, double* %10
%75 = fcmp ogt double %74, 0.000000e+00
br i1 %75, label %76, label %152
76: ; preds = %73
%77 = load double*, double** %6
%78 = getelementptr inbounds double, double* %77, i64 1
%79 = load double, double* %78
store double %79, double* %14
store i32 1, i32* %11
br label %80
80: ; preds = %100, %76
%81 = load i32, i32* %11
%82 = load i32, i32* %9
%83 = add nsw i32 %82, 1
%84 = icmp sle i32 %81, %83
br i1 %84, label %85, label %103
85: ; preds = %80
%86 = load double*, double** %6
%87 = load i32, i32* %11
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %86, i64 %88
%90 = load double, double* %89
%91 = load double, double* %14
%92 = fcmp ogt double %90, %91
br i1 %92, label %93, label %99
93: ; preds = %85
%94 = load double*, double** %6
%95 = load i32, i32* %11
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %94, i64 %96
%98 = load double, double* %97
store double %98, double* %14
br label %99
99: ; preds = %93, %85
br label %100
100: ; preds = %99
%101 = load i32, i32* %11
%102 = add nsw i32 %101, 1
store i32 %102, i32* %11
br label %80
103: ; preds = %80
store double 0.000000e+00, double* %13
store i32 1, i32* %11
br label %104
104: ; preds = %126, %103
%105 = load i32, i32* %11
%106 = load i32, i32* %9
%107 = add nsw i32 %106, 1
%108 = icmp sle i32 %105, %107
br i1 %108, label %109, label %129
109: ; preds = %104
%110 = load double, double* %14
%111 = load double*, double** %6
%112 = load i32, i32* %11
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %111, i64 %113
%115 = load double, double* %114
%116 = fsub double %115, %110
store double %116, double* %114
%117 = load double*, double** %6
%118 = load i32, i32* %11
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %117, i64 %119
%121 = load double, double* %120
%122 = call i64 @exp(double %121)
%123 = sitofp i64 %122 to double
%124 = load double, double* %13
%125 = fadd double %124, %123
store double %125, double* %13
br label %126
126: ; preds = %109
%127 = load i32, i32* %11
%128 = add nsw i32 %127, 1
store i32 %128, i32* %11
br label %104
129: ; preds = %104
store i32 1, i32* %11
br label %130
130: ; preds = %148, %129
%131 = load i32, i32* %11
%132 = load i32, i32* %9
%133 = add nsw i32 %132, 1
%134 = icmp sle i32 %131, %133
br i1 %134, label %135, label %151
135: ; preds = %130
%136 = load double, double* %13
%137 = call double @log(double %136)
%138 = fneg double %137
%139 = load double, double* %10
%140 = call double @log(double %139)
%141 = fadd double %138, %140
%142 = load double*, double** %6
%143 = load i32, i32* %11
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double, double* %142, i64 %144
%146 = load double, double* %145
%147 = fadd double %146, %141
store double %147, double* %145
br label %148
148: ; preds = %135
%149 = load i32, i32* %11
%150 = add nsw i32 %149, 1
store i32 %150, i32* %11
br label %130
151: ; preds = %130
br label %152
152: ; preds = %151, %73
ret void
}

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