#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 exp(double);
void gaussian(float *x, int n, float sig) {
  int i, ii;
  for (i = 0; i < n; i++) {
    if (i > n / 2)
      ii = i - n;
    else
      ii = i;
    x[i] = exp(-.5 * ii * ii * sig * sig);
  }
  return;
}


x86 assembly (clang, -O0):
<gaussian>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11bd <gaussian+0xad>
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x1c(%rbp),%eax
cmp    %ecx,%eax
jle    1163 <gaussian+0x53>
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1169 <gaussian+0x59>
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xe8a(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movss  -0x10(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
movss  -0x10(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
call   1030 <exp@plt>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    112b <gaussian+0x1b>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gaussian>:
push   %r15
push   %r14
push   %r12
push   %rbx
push   %rax
test   %esi,%esi
jle    1188 <gaussian+0x78>
mov    %esi,%r15d
mov    %rdi,%r14
mov    %esi,%r12d
shr    %r15d
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%rsp)
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    %r15,%rbx
mov    $0x0,%eax
cmova  %r12d,%eax
mov    %ebx,%ecx
sub    %eax,%ecx
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
movapd %xmm1,%xmm0
mulsd  0xe9d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  (%rsp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
call   1030 <exp@plt>
movss  %xmm0,(%r14,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%r12
jne    1140 <gaussian+0x30>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @gaussian(float* %0, i32 %1, float %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca float
%7 = alloca i32
%8 = alloca i32
store float* %0, float** %4
store i32 %1, i32* %5
store float %2, float* %6
store i32 0, i32* %7
br label %9
9: ; preds = %42, %3
%10 = load i32, i32* %7
%11 = load i32, i32* %5
%12 = icmp slt i32 %10, %11
br i1 %12, label %13, label %45
13: ; preds = %9
%14 = load i32, i32* %7
%15 = load i32, i32* %5
%16 = sdiv i32 %15, 2
%17 = icmp sgt i32 %14, %16
br i1 %17, label %18, label %22
18: ; preds = %13
%19 = load i32, i32* %7
%20 = load i32, i32* %5
%21 = sub nsw i32 %19, %20
store i32 %21, i32* %8
br label %24
22: ; preds = %13
%23 = load i32, i32* %7
store i32 %23, i32* %8
br label %24
24: ; preds = %22, %18
%25 = load i32, i32* %8
%26 = sitofp i32 %25 to double
%27 = fmul double -5.000000e-01, %26
%28 = load i32, i32* %8
%29 = sitofp i32 %28 to double
%30 = fmul double %27, %29
%31 = load float, float* %6
%32 = fpext float %31 to double
%33 = fmul double %30, %32
%34 = load float, float* %6
%35 = fpext float %34 to double
%36 = fmul double %33, %35
%37 = call float @exp(double %36)
%38 = load float*, float** %4
%39 = load i32, i32* %7
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
store float %37, float* %41
br label %42
42: ; preds = %24
%43 = load i32, i32* %7
%44 = add nsw i32 %43, 1
store i32 %44, i32* %7
br label %9
45: ; preds = %9
ret void
}

LLVM IR (-O2):
define void @gaussian(float* nocapture writeonly %0, i32 %1, float %2) {
%4 = icmp sgt i32 %1, 0
br i1 %4, label %5, label %25
5: ; preds = %3
%6 = lshr i32 %1, 1
%7 = fpext float %2 to double
%8 = zext i32 %6 to i64
%9 = zext i32 %1 to i64
br label %10
10: ; preds = %5, %10
%11 = phi i64 [ 0, %5 ], [ %23, %10 ]
%12 = icmp ugt i64 %11, %8
%13 = select i1 %12, i32 %1, i32 0
%14 = trunc i64 %11 to i32
%15 = sub nsw i32 %14, %13
%16 = sitofp i32 %15 to double
%17 = fmul double %16, -5.000000e-01
%18 = fmul double %17, %16
%19 = fmul double %18, %7
%20 = fmul double %19, %7
%21 = tail call float @exp(double %20)
%22 = getelementptr inbounds float, float* %0, i64 %11
store float %21, float* %22
%23 = add nuw nsw i64 %11, 1
%24 = icmp eq i64 %23, %9
br i1 %24, label %25, label %10
25: ; preds = %10, %3
ret void
}