#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 */
float expf(float);
float CalcHistAt(float pt, int bins, float *out_means, float *out_amps,
                 float *out_sigmas, int cur_dists) {
  int i;
  float mark = 0.0;
  for (i = 0; i < cur_dists; i++) {
    float amp = out_amps[i];
    float mean = out_means[i];
    float var = out_sigmas[i];
    float x = mean - pt;
    if (x < -bins / 2)
      x += bins;
    if (x > bins / 2)
      x -= bins;
    float nrm = amp / (var * 2.506628253) * expf(-(x * x) / (2 * var * var));
    mark += nrm;
  }
  return mark;
}


x86 assembly (clang, -O0):
<CalcHistAt>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movss  %xmm0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    127d <CalcHistAt+0x16d>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x34(%rbp),%xmm0
subss  -0x4(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm1
xor    %eax,%eax
sub    -0x8(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cvtsi2ss %eax,%xmm0
ucomiss %xmm1,%xmm0
jbe    11be <CalcHistAt+0xae>
cvtsi2ssl -0x8(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
mov    -0x8(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cvtsi2ss %eax,%xmm1
ucomiss %xmm1,%xmm0
jbe    11ee <CalcHistAt+0xde>
cvtsi2ssl -0x8(%rbp),%xmm1
movss  -0x3c(%rbp),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x38(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xe00(%rip),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movss  -0x3c(%rbp),%xmm0
mulss  -0x3c(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  0xdcc(%rip),%xmm1
mulss  -0x38(%rbp),%xmm1
mulss  -0x38(%rbp),%xmm1
divss  %xmm1,%xmm0
call   1030 <expf@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113f <CalcHistAt+0x2f>
movss  -0x2c(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<CalcHistAt>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x30,%rsp
test   %r8d,%r8d
jle    1250 <CalcHistAt+0x140>
mov    %rcx,%r14
mov    %rdx,%r15
mov    %rsi,%r12
movaps %xmm0,%xmm3
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
sar    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  %xmm0,0x10(%rsp)
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  %xmm0,0xc(%rsp)
xorps  %xmm0,%xmm0
cvtsi2ss %edi,%xmm0
movaps %xmm0,0x20(%rsp)
mov    %r8d,%r13d
xorps  %xmm0,%xmm0
xor    %ebx,%ebx
movss  %xmm3,0x14(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  %xmm0,0x8(%rsp)
movss  (%r15,%rbx,4),%xmm0
movss  (%r12,%rbx,4),%xmm1
movss  (%r14,%rbx,4),%xmm2
subss  %xmm3,%xmm1
movaps %xmm1,%xmm3
cmpltss 0xc(%rsp),%xmm3
movaps %xmm3,%xmm4
movaps 0x20(%rsp),%xmm5
andps  %xmm5,%xmm4
movss  0xe47(%rip),%xmm6
andnps %xmm6,%xmm3
orps   %xmm4,%xmm3
addss  %xmm1,%xmm3
movss  0x10(%rsp),%xmm1
cmpltss %xmm3,%xmm1
andps  %xmm5,%xmm1
subss  %xmm1,%xmm3
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mulsd  0xe1d(%rip),%xmm0
divsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsp)
movaps %xmm3,%xmm0
xorps  0xe11(%rip),%xmm0
mulss  %xmm3,%xmm0
movaps %xmm2,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
call   1030 <expf@plt>
movss  0x14(%rsp),%xmm3
cvtss2sd %xmm0,%xmm0
mulsd  0x18(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  0x8(%rsp),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
movss  0x8(%rsp),%xmm0
add    $0x1,%rbx
cmp    %rbx,%r13
jne    1180 <CalcHistAt+0x70>
jmp    1253 <CalcHistAt+0x143>
xorps  %xmm0,%xmm0
add    $0x30,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define float @CalcHistAt(float %0, i32 %1, float* nocapture %2, float* nocapture %3, float* nocapture %4, i32 %5) {
%7 = icmp sgt i32 %5, 0
br i1 %7, label %8, label %47
8: ; preds = %6
%9 = sdiv i32 %1, -2
%10 = sitofp i32 %9 to float
%11 = sitofp i32 %1 to float
%12 = sdiv i32 %1, 2
%13 = sitofp i32 %12 to float
%14 = zext i32 %5 to i64
br label %15
15: ; preds = %8, %15
%16 = phi i64 [ 0, %8 ], [ %45, %15 ]
%17 = phi float [ 0.000000e+00, %8 ], [ %44, %15 ]
%18 = getelementptr inbounds float, float* %3, i64 %16
%19 = load float, float* %18
%20 = getelementptr inbounds float, float* %2, i64 %16
%21 = load float, float* %20
%22 = getelementptr inbounds float, float* %4, i64 %16
%23 = load float, float* %22
%24 = fsub float %21, %0
%25 = fcmp olt float %24, %10
%26 = select i1 %25, float %11, float -0.000000e+00
%27 = fadd float %24, %26
%28 = fcmp ogt float %27, %13
%29 = select i1 %28, float %11, float 0.000000e+00
%30 = fsub float %27, %29
%31 = fpext float %19 to double
%32 = fpext float %23 to double
%33 = fmul double %32, 0x40040D931D0EEABB
%34 = fdiv double %31, %33
%35 = fneg float %30
%36 = fmul float %30, %35
%37 = fmul float %23, 2.000000e+00
%38 = fmul float %23, %37
%39 = fdiv float %36, %38
%40 = tail call float @expf(float %39)
%41 = fpext float %40 to double
%42 = fmul double %34, %41
%43 = fptrunc double %42 to float
%44 = fadd float %17, %43
%45 = add nuw nsw i64 %16, 1
%46 = icmp eq i64 %45, %14
br i1 %46, label %47, label %15
47: ; preds = %15, %6
%48 = phi float [ 0.000000e+00, %6 ], [ %44, %15 ]
ret float %48
}