#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 /*<<< orphan*/ self_adj_index_search(int *, double, double *, int);
void self_adj_theta_update(double *fx, int N_population, double *theta,
                           double *grid_points, int grid_size, double *freq_ref,
                           double gain, double *freq_loc) {
  int n;
  int k;
  for (k = 1; k <= grid_size + 1; k++)
    freq_loc[k] = 0.0;
  for (n = 1; n <= N_population; n++) {
    self_adj_index_search(&k, fx[n], grid_points, grid_size);
    freq_loc[k] += 1.0;
  }
  for (k = 1; k <= grid_size + 1; k++) {
    freq_loc[k] /= N_population;
    theta[k] += gain * (freq_loc[k] - freq_ref[k]);
  }
}


x86 assembly (clang, -O0):
<self_adj_theta_update>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     116e <self_adj_theta_update+0x5e>
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    113f <self_adj_theta_update+0x2f>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11c6 <self_adj_theta_update+0xb6>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rsi
mov    -0x24(%rbp),%edx
lea    -0x40(%rbp),%rdi
call   1030 <self_adj_index_search@plt>
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  0xe52(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1175 <self_adj_theta_update+0x65>
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1240 <self_adj_theta_update+0x130>
cvtsi2sdl -0xc(%rbp),%xmm1
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11cd <self_adj_theta_update+0xbd>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<self_adj_theta_update>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
movapd %xmm0,0x30(%rsp)
mov    %r8d,%r12d
mov    %rcx,%r13
mov    %esi,%ebx
mov    %rdi,%rbp
mov    0x80(%rsp),%rax
mov    %r8d,%ecx
test   %r8d,%r8d
mov    %r9,0x28(%rsp)
mov    %rdx,0x20(%rsp)
mov    %rcx,0x18(%rsp)
js     117b <self_adj_theta_update+0x5b>
mov    %rax,%r14
lea    0x8(%rax),%rdi
lea    0x8(,%rcx,8),%rdx
xor    %esi,%esi
call   1040 <memset@plt>
lea    0x2(%r12),%eax
jmp    1183 <self_adj_theta_update+0x63>
mov    %rax,%r14
mov    $0x1,%eax
mov    %eax,0xc(%rsp)
mov    %rbx,0x10(%rsp)
test   %ebx,%ebx
mov    %r14,%rbx
jle    11e7 <self_adj_theta_update+0xc7>
mov    0x10(%rsp),%rax
lea    0x1(%rax),%r15d
mov    $0x1,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%r14,8),%xmm0
lea    0xc(%rsp),%rdi
mov    %r13,%rsi
mov    %r12d,%edx
call   1030 <self_adj_index_search@plt>
movslq 0xc(%rsp),%rax
movsd  (%rbx,%rax,8),%xmm0
addsd  0xe27(%rip),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
add    $0x1,%r14
cmp    %r14,%r15
jne    11b0 <self_adj_theta_update+0x90>
test   %r12d,%r12d
movapd 0x30(%rsp),%xmm5
mov    0x28(%rsp),%rdi
mov    0x20(%rsp),%rbp
js     130c <self_adj_theta_update+0x1ec>
xorps  %xmm0,%xmm0
cvtsi2sdl 0x10(%rsp),%xmm0
mov    0x18(%rsp),%r12
lea    0x2(%r12),%rax
add    $0x1,%r12
mov    $0x1,%ecx
cmp    $0x2,%r12
jb     12e0 <self_adj_theta_update+0x1c0>
lea    0x8(%rbx),%rdx
lea    (%rbx,%rax,8),%r15
lea    0x8(%rbp),%r8
lea    0x0(,%rax,8),%r9
add    %rbp,%r9
lea    0x8(%rdi),%r14
lea    (%rdi,%rax,8),%rsi
cmp    %r9,%rdx
setb   %r10b
cmp    %r15,%r8
setb   %r11b
cmp    %rsi,%rdx
setb   %dl
cmp    %r15,%r14
setb   %r15b
cmp    %rsi,%r8
setb   %r8b
cmp    %r9,%r14
setb   %r9b
test   %r11b,%r10b
jne    12e0 <self_adj_theta_update+0x1c0>
and    %r15b,%dl
jne    12e0 <self_adj_theta_update+0x1c0>
and    %r9b,%r8b
jne    12e0 <self_adj_theta_update+0x1c0>
mov    %r12,%r8
mov    %r12,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %r12,%rcx
or     $0x1,%rcx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
movapd %xmm5,%xmm2
unpcklpd %xmm5,%xmm2
xor    %esi,%esi
movupd 0x8(%rbx,%rsi,8),%xmm3
divpd  %xmm1,%xmm3
movupd %xmm3,0x8(%rbx,%rsi,8)
movupd 0x8(%rdi,%rsi,8),%xmm4
subpd  %xmm4,%xmm3
mulpd  %xmm2,%xmm3
movupd 0x8(%rbp,%rsi,8),%xmm4
addpd  %xmm3,%xmm4
movupd %xmm4,0x8(%rbp,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%rdx
jne    12a0 <self_adj_theta_update+0x180>
cmp    %rdx,%r8
je     130c <self_adj_theta_update+0x1ec>
nopl   0x0(%rax)
movsd  (%rbx,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
subsd  (%rdi,%rcx,8),%xmm1
mulsd  %xmm5,%xmm1
addsd  0x0(%rbp,%rcx,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12e0 <self_adj_theta_update+0x1c0>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @self_adj_theta_update(double* %0, i32 %1, double* %2, double* %3, i32 %4, double* %5, double %6, double* %7) {
%9 = alloca double*
%10 = alloca i32
%11 = alloca double*
%12 = alloca double*
%13 = alloca i32
%14 = alloca double*
%15 = alloca double
%16 = alloca double*
%17 = alloca i32
%18 = alloca i32
store double* %0, double** %9
store i32 %1, i32* %10
store double* %2, double** %11
store double* %3, double** %12
store i32 %4, i32* %13
store double* %5, double** %14
store double %6, double* %15
store double* %7, double** %16
store i32 1, i32* %18
br label %19
19: ; preds = %29, %8
%20 = load i32, i32* %18
%21 = load i32, i32* %13
%22 = add nsw i32 %21, 1
%23 = icmp sle i32 %20, %22
br i1 %23, label %24, label %32
24: ; preds = %19
%25 = load double*, double** %16
%26 = load i32, i32* %18
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
store double 0.000000e+00, double* %28
br label %29
29: ; preds = %24
%30 = load i32, i32* %18
%31 = add nsw i32 %30, 1
store i32 %31, i32* %18
br label %19
32: ; preds = %19
store i32 1, i32* %17
br label %33
33: ; preds = %52, %32
%34 = load i32, i32* %17
%35 = load i32, i32* %10
%36 = icmp sle i32 %34, %35
br i1 %36, label %37, label %55
37: ; preds = %33
%38 = load double*, double** %9
%39 = load i32, i32* %17
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
%42 = load double, double* %41
%43 = load double*, double** %12
%44 = load i32, i32* %13
%45 = call i32 @self_adj_index_search(i32* %18, double %42, double* %43, i32 %44)
%46 = load double*, double** %16
%47 = load i32, i32* %18
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
%50 = load double, double* %49
%51 = fadd double %50, 1.000000e+00
store double %51, double* %49
br label %52
52: ; preds = %37
%53 = load i32, i32* %17
%54 = add nsw i32 %53, 1
store i32 %54, i32* %17
br label %33
55: ; preds = %33
store i32 1, i32* %18
br label %56
56: ; preds = %88, %55
%57 = load i32, i32* %18
%58 = load i32, i32* %13
%59 = add nsw i32 %58, 1
%60 = icmp sle i32 %57, %59
br i1 %60, label %61, label %91
61: ; preds = %56
%62 = load i32, i32* %10
%63 = sitofp i32 %62 to double
%64 = load double*, double** %16
%65 = load i32, i32* %18
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %64, i64 %66
%68 = load double, double* %67
%69 = fdiv double %68, %63
store double %69, double* %67
%70 = load double, double* %15
%71 = load double*, double** %16
%72 = load i32, i32* %18
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %71, i64 %73
%75 = load double, double* %74
%76 = load double*, double** %14
%77 = load i32, i32* %18
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %76, i64 %78
%80 = load double, double* %79
%81 = fsub double %75, %80
%82 = load double*, double** %11
%83 = load i32, i32* %18
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %82, i64 %84
%86 = load double, double* %85
%87 = call double @llvm.fmuladd.f64(double %70, double %81, double %86)
store double %87, double* %85
br label %88
88: ; preds = %61
%89 = load i32, i32* %18
%90 = add nsw i32 %89, 1
store i32 %90, i32* %18
br label %56
91: ; preds = %56
ret void
}

LLVM IR (-O2):
define void @self_adj_theta_update(double* nocapture %0, i32 %1, double* nocapture %2, double* %3, i32 %4, double* nocapture %5, double %6, double* nocapture %7) {
%9 = alloca i32
%10 = bitcast i32* %9 to i8*
%11 = icmp slt i32 %4, 0
br i1 %11, label %19, label %12
12: ; preds = %8
%13 = getelementptr double, double* %7, i64 1
%14 = bitcast double* %13 to i8*
%15 = add nuw nsw i32 %4, 1
%16 = zext i32 %15 to i64
%17 = shl nuw nsw i64 %16, 3
%18 = add nuw i32 %4, 2
br label %19
19: ; preds = %12, %8
%20 = phi i32 [ 1, %8 ], [ %18, %12 ]
store i32 %20, i32* %9
%21 = icmp slt i32 %1, 1
br i1 %21, label %25, label %22
22: ; preds = %19
%23 = add nuw i32 %1, 1
%24 = zext i32 %23 to i64
br label %80
25: ; preds = %80, %19
br i1 %11, label %105, label %26
26: ; preds = %25
%27 = sitofp i32 %1 to double
%28 = add nuw i32 %4, 2
%29 = zext i32 %28 to i64
%30 = add nsw i64 %29, -1
%31 = icmp ult i64 %30, 2
br i1 %31, label %78, label %32
32: ; preds = %26
%33 = getelementptr double, double* %7, i64 1
%34 = getelementptr double, double* %7, i64 %29
%35 = getelementptr double, double* %2, i64 1
%36 = getelementptr double, double* %2, i64 %29
%37 = getelementptr double, double* %5, i64 1
%38 = getelementptr double, double* %5, i64 %29
%39 = icmp ult double* %33, %36
%40 = icmp ult double* %35, %34
%41 = and i1 %39, %40
%42 = icmp ult double* %33, %38
%43 = icmp ult double* %37, %34
%44 = and i1 %42, %43
%45 = or i1 %41, %44
%46 = icmp ult double* %35, %38
%47 = icmp ult double* %37, %36
%48 = and i1 %46, %47
%49 = or i1 %45, %48
br i1 %49, label %78, label %50
50: ; preds = %32
%51 = and i64 %30, -2
%52 = or i64 %30, 1
%53 = insertelement <2 x double> poison, double %27, i64 0
%54 = shufflevector <2 x double> %53, <2 x double> poison, <2 x i32> zeroinitializer
%55 = insertelement <2 x double> poison, double %6, i64 0
%56 = shufflevector <2 x double> %55, <2 x double> poison, <2 x i32> zeroinitializer
br label %57
57: ; preds = %57, %50
%58 = phi i64 [ 0, %50 ], [ %74, %57 ]
%59 = or i64 %58, 1
%60 = getelementptr inbounds double, double* %7, i64 %59
%61 = bitcast double* %60 to <2 x double>*
%62 = load <2 x double>, <2 x double>* %61
%63 = fdiv <2 x double> %62, %54
%64 = bitcast double* %60 to <2 x double>*
store <2 x double> %63, <2 x double>* %64
%65 = getelementptr inbounds double, double* %5, i64 %59
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = fsub <2 x double> %63, %67
%69 = getelementptr inbounds double, double* %2, i64 %59
%70 = bitcast double* %69 to <2 x double>*
%71 = load <2 x double>, <2 x double>* %70
%72 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %56, <2 x double> %68, <2 x double> %71)
%73 = bitcast double* %69 to <2 x double>*
store <2 x double> %72, <2 x double>* %73
%74 = add nuw i64 %58, 2
%75 = icmp eq i64 %74, %51
br i1 %75, label %76, label %57
76: ; preds = %57
%77 = icmp eq i64 %30, %51
br i1 %77, label %105, label %78
78: ; preds = %32, %26, %76
%79 = phi i64 [ 1, %32 ], [ 1, %26 ], [ %52, %76 ]
br label %92
80: ; preds = %22, %80
%81 = phi i64 [ 1, %22 ], [ %90, %80 ]
%82 = getelementptr inbounds double, double* %0, i64 %81
%83 = load double, double* %82
%84 = call i32 @self_adj_index_search(i32* %9, double %83, double* %3, i32 %4)
%85 = load i32, i32* %9
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %7, i64 %86
%88 = load double, double* %87
%89 = fadd double %88, 1.000000e+00
store double %89, double* %87
%90 = add nuw nsw i64 %81, 1
%91 = icmp eq i64 %90, %24
br i1 %91, label %25, label %80
92: ; preds = %78, %92
%93 = phi i64 [ %103, %92 ], [ %79, %78 ]
%94 = getelementptr inbounds double, double* %7, i64 %93
%95 = load double, double* %94
%96 = fdiv double %95, %27
store double %96, double* %94
%97 = getelementptr inbounds double, double* %5, i64 %93
%98 = load double, double* %97
%99 = fsub double %96, %98
%100 = getelementptr inbounds double, double* %2, i64 %93
%101 = load double, double* %100
%102 = call double @llvm.fmuladd.f64(double %6, double %99, double %101)
store double %102, double* %100
%103 = add nuw nsw i64 %93, 1
%104 = icmp eq i64 %103, %29
br i1 %104, label %105, label %92
105: ; preds = %92, %76, %25
ret void
}