#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(int);
void set_kernel(float sigma, int *radius, float *kernel) {
  if (sigma == 0.0f) {
    return;
  }
  int rad = (int)(sigma * 3.0f + 0.5f);
  if (rad < 1) {
    rad = 1;
  }
  float sum = 0.0f;
  for (int i = -rad; i <= rad; i++) {
    float weight = expf(-(i * i) / (2.0f * sigma * sigma));
    kernel[i + rad] = weight;
    sum += weight;
  }
  for (int i = 0, len = rad * 2 + 1; i < len; kernel[i++] /= sum)
    ;
  *radius = rad;
}


x86 assembly (clang, -O0):
<set_kernel>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movss  -0x4(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1141 <set_kernel+0x31>
jp     1141 <set_kernel+0x31>
jmp    1250 <set_kernel+0x140>
movss  -0x4(%rbp),%xmm0
movss  0xeb2(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  0xeaa(%rip),%xmm1
addss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x1,-0x1c(%rbp)
jge    1176 <set_kernel+0x66>
movl   $0x1,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
xor    %eax,%eax
sub    -0x1c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     11f9 <set_kernel+0xe9>
mov    -0x24(%rbp),%ecx
imul   -0x24(%rbp),%ecx
xor    %eax,%eax
sub    %ecx,%eax
cvtsi2ss %eax,%xmm0
movss  0xe5f(%rip),%xmm1
mulss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   1030 <expf@plt>
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x28(%rbp),%xmm0
addss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1186 <set_kernel+0x76>
movl   $0x0,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    1247 <set_kernel+0x137>
jmp    121d <set_kernel+0x10d>
movss  -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x2c(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    120c <set_kernel+0xfc>
mov    -0x1c(%rbp),%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<set_kernel>:
xorps  %xmm3,%xmm3
ucomiss %xmm3,%xmm0
jne    111b <set_kernel+0xb>
jp     111b <set_kernel+0xb>
ret
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rsi,%rbx
mov    %rdi,%r14
movss  0xec9(%rip),%xmm1
mulss  %xmm0,%xmm1
addss  0xec1(%rip),%xmm1
cvttss2si %xmm1,%eax
test   %eax,%eax
mov    $0x1,%r15d
cmovg  %eax,%r15d
movaps %xmm0,%xmm1
addss  %xmm0,%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0xc(%rsp)
mov    %r15d,%r12d
neg    %r12d
lea    (%r15,%r15,1),%r13
add    $0x1,%r13
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movaps %xmm3,0x10(%rsp)
lea    (%r12,%rbp,1),%eax
imul   %eax,%eax
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
divss  0xc(%rsp),%xmm0
cvttss2si %xmm0,%edi
call   1030 <expf@plt>
movaps 0x10(%rsp),%xmm3
movss  %xmm0,(%rbx,%rbp,4)
addss  %xmm0,%xmm3
add    $0x1,%rbp
cmp    %ebp,%r13d
jne    1180 <set_kernel+0x70>
lea    (%r15,%r15,1),%eax
add    $0x1,%eax
cmp    $0x4,%eax
jae    11cb <set_kernel+0xbb>
xor    %ecx,%ecx
jmp    1240 <set_kernel+0x130>
mov    %eax,%ecx
and    $0x7ffffffc,%ecx
movaps %xmm3,%xmm0
shufps $0x0,%xmm3,%xmm0
lea    -0x4(%rcx),%rsi
mov    %rsi,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
test   %rsi,%rsi
je     1229 <set_kernel+0x119>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movups (%rbx,%rsi,4),%xmm1
movups 0x10(%rbx,%rsi,4),%xmm2
divps  %xmm0,%xmm1
movups %xmm1,(%rbx,%rsi,4)
divps  %xmm0,%xmm2
movups %xmm2,0x10(%rbx,%rsi,4)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1200 <set_kernel+0xf0>
test   $0x1,%dl
je     1240 <set_kernel+0x130>
jmp    1230 <set_kernel+0x120>
xor    %esi,%esi
test   $0x1,%dl
je     1240 <set_kernel+0x130>
movups (%rbx,%rsi,4),%xmm1
divps  %xmm0,%xmm1
movups %xmm1,(%rbx,%rsi,4)
nopl   0x0(%rax,%rax,1)
movss  (%rbx,%rcx,4),%xmm0
divss  %xmm3,%xmm0
movss  %xmm0,(%rbx,%rcx,4)
lea    0x1(%rcx),%rdx
mov    %rdx,%rcx
cmp    %rdx,%rax
jne    1240 <set_kernel+0x130>
mov    %r15d,(%r14)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @set_kernel(float %0, i32* %1, float* %2) {
%4 = alloca float
%5 = alloca i32*
%6 = alloca float*
%7 = alloca i32
%8 = alloca float
%9 = alloca i32
%10 = alloca float
%11 = alloca i32
%12 = alloca i32
store float %0, float* %4
store i32* %1, i32** %5
store float* %2, float** %6
%13 = load float, float* %4
%14 = fcmp oeq float %13, 0.000000e+00
br i1 %14, label %15, label %16
15: ; preds = %3
br label %77
16: ; preds = %3
%17 = load float, float* %4
%18 = call float @llvm.fmuladd.f32(float %17, float 3.000000e+00, float 5.000000e-01)
%19 = fptosi float %18 to i32
store i32 %19, i32* %7
%20 = load i32, i32* %7
%21 = icmp slt i32 %20, 1
br i1 %21, label %22, label %23
22: ; preds = %16
store i32 1, i32* %7
br label %23
23: ; preds = %22, %16
store float 0.000000e+00, float* %8
%24 = load i32, i32* %7
%25 = sub nsw i32 0, %24
store i32 %25, i32* %9
br label %26
26: ; preds = %53, %23
%27 = load i32, i32* %9
%28 = load i32, i32* %7
%29 = icmp sle i32 %27, %28
br i1 %29, label %30, label %56
30: ; preds = %26
%31 = load i32, i32* %9
%32 = load i32, i32* %9
%33 = mul nsw i32 %31, %32
%34 = sub nsw i32 0, %33
%35 = sitofp i32 %34 to float
%36 = load float, float* %4
%37 = fmul float 2.000000e+00, %36
%38 = load float, float* %4
%39 = fmul float %37, %38
%40 = fdiv float %35, %39
%41 = fptosi float %40 to i32
%42 = call float @expf(i32 %41)
store float %42, float* %10
%43 = load float, float* %10
%44 = load float*, float** %6
%45 = load i32, i32* %9
%46 = load i32, i32* %7
%47 = add nsw i32 %45, %46
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %44, i64 %48
store float %43, float* %49
%50 = load float, float* %10
%51 = load float, float* %8
%52 = fadd float %51, %50
store float %52, float* %8
br label %53
53: ; preds = %30
%54 = load i32, i32* %9
%55 = add nsw i32 %54, 1
store i32 %55, i32* %9
br label %26
56: ; preds = %26
store i32 0, i32* %11
%57 = load i32, i32* %7
%58 = mul nsw i32 %57, 2
%59 = add nsw i32 %58, 1
store i32 %59, i32* %12
br label %60
60: ; preds = %65, %56
%61 = load i32, i32* %11
%62 = load i32, i32* %12
%63 = icmp slt i32 %61, %62
br i1 %63, label %64, label %74
64: ; preds = %60
br label %65
65: ; preds = %64
%66 = load float, float* %8
%67 = load float*, float** %6
%68 = load i32, i32* %11
%69 = add nsw i32 %68, 1
store i32 %69, i32* %11
%70 = sext i32 %68 to i64
%71 = getelementptr inbounds float, float* %67, i64 %70
%72 = load float, float* %71
%73 = fdiv float %72, %66
store float %73, float* %71
br label %60
74: ; preds = %60
%75 = load i32, i32* %7
%76 = load i32*, i32** %5
store i32 %75, i32* %76
br label %77
77: ; preds = %74, %15
ret void
}

LLVM IR (-O2):
define void @set_kernel(float %0, i32* nocapture writeonly %1, float* nocapture %2) {
%4 = fcmp oeq float %0, 0.000000e+00
br i1 %4, label %84, label %5
5: ; preds = %3
%6 = tail call float @llvm.fmuladd.f32(float %0, float 3.000000e+00, float 5.000000e-01)
%7 = fptosi float %6 to i32
%8 = icmp sgt i32 %7, 1
%9 = select i1 %8, i32 %7, i32 1
%10 = sub nsw i32 0, %9
%11 = fmul float %0, 2.000000e+00
%12 = fmul float %11, %0
%13 = sext i32 %10 to i64
%14 = zext i32 %9 to i64
%15 = add nuw i32 %9, 1
br label %60
16: ; preds = %60
%17 = shl nuw nsw i32 %9, 1
%18 = or i32 %17, 1
%19 = zext i32 %18 to i64
%20 = icmp ult i32 %18, 4
br i1 %20, label %30, label %32
21: ; preds = %43, %32
%22 = phi i64 [ 0, %32 ], [ %57, %43 ]
%23 = icmp eq i64 %39, 0
br i1 %23, label %30, label %24
24: ; preds = %21
%25 = getelementptr inbounds float, float* %2, i64 %22
%26 = bitcast float* %25 to <4 x float>*
%27 = load <4 x float>, <4 x float>* %26
%28 = fdiv <4 x float> %27, %35
%29 = bitcast float* %25 to <4 x float>*
store <4 x float> %28, <4 x float>* %29
br label %30
30: ; preds = %24, %21, %16
%31 = phi i64 [ 0, %16 ], [ %33, %21 ], [ %33, %24 ]
br label %77
32: ; preds = %16
%33 = and i64 %19, 2147483644
%34 = insertelement <4 x float> poison, float %72, i64 0
%35 = shufflevector <4 x float> %34, <4 x float> poison, <4 x i32> zeroinitializer
%36 = add nsw i64 %33, -4
%37 = lshr exact i64 %36, 2
%38 = add nuw nsw i64 %37, 1
%39 = and i64 %38, 1
%40 = icmp eq i64 %36, 0
br i1 %40, label %21, label %41
41: ; preds = %32
%42 = and i64 %38, 9223372036854775806
br label %43
43: ; preds = %43, %41
%44 = phi i64 [ 0, %41 ], [ %57, %43 ]
%45 = phi i64 [ 0, %41 ], [ %58, %43 ]
%46 = getelementptr inbounds float, float* %2, i64 %44
%47 = bitcast float* %46 to <4 x float>*
%48 = load <4 x float>, <4 x float>* %47
%49 = fdiv <4 x float> %48, %35
%50 = bitcast float* %46 to <4 x float>*
store <4 x float> %49, <4 x float>* %50
%51 = or i64 %44, 4
%52 = getelementptr inbounds float, float* %2, i64 %51
%53 = bitcast float* %52 to <4 x float>*
%54 = load <4 x float>, <4 x float>* %53
%55 = fdiv <4 x float> %54, %35
%56 = bitcast float* %52 to <4 x float>*
store <4 x float> %55, <4 x float>* %56
%57 = add nuw i64 %44, 8
%58 = add i64 %45, 2
%59 = icmp eq i64 %58, %42
br i1 %59, label %21, label %43
60: ; preds = %5, %60
%61 = phi i64 [ %13, %5 ], [ %73, %60 ]
%62 = phi float [ 0.000000e+00, %5 ], [ %72, %60 ]
%63 = trunc i64 %61 to i32
%64 = mul i32 %63, %63
%65 = sub i32 0, %64
%66 = sitofp i32 %65 to float
%67 = fdiv float %66, %12
%68 = fptosi float %67 to i32
%69 = tail call float @expf(i32 %68)
%70 = add nsw i64 %61, %14
%71 = getelementptr inbounds float, float* %2, i64 %70
store float %69, float* %71
%72 = fadd float %62, %69
%73 = add nsw i64 %61, 1
%74 = trunc i64 %73 to i32
%75 = icmp eq i32 %15, %74
br i1 %75, label %16, label %60
76: ; preds = %77
store i32 %9, i32* %1
br label %84
77: ; preds = %30, %77
%78 = phi i64 [ %79, %77 ], [ %31, %30 ]
%79 = add nuw nsw i64 %78, 1
%80 = getelementptr inbounds float, float* %2, i64 %78
%81 = load float, float* %80
%82 = fdiv float %81, %72
store float %82, float* %80
%83 = icmp eq i64 %79, %19
br i1 %83, label %76, label %77
84: ; preds = %3, %76
ret void
}