#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*/ assert(int);
float exp(double);
float *xmalloc(int);
float *malloc_gaussian_symm_kernel(float sigma, int rad) {
  assert(sigma >= 0);
  float *gker = xmalloc((rad + 1) * sizeof(float));
  gker[0] = 1.;
  if (sigma > 0) {
    float sum = gker[0];
    for (int i = 1; i <= rad; i++) {
      gker[i] = exp(-0.5 * (float)i * (float)i / sigma / sigma);
      sum += 2 * gker[i];
    }
    for (int i = 0; i <= rad; i++)
      gker[i] /= sum;
  } else {
    for (int i = 1; i <= rad; i++) {
      gker[i] = 0.0;
    }
  }
  return gker;
}


x86 assembly (clang, -O0):
<malloc_gaussian_symm_kernel>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movss  %xmm0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
movss  -0x4(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cltq
shl    $0x2,%rax
mov    %eax,%edi
call   1030 <xmalloc@plt>
mov    %rax,-0x10(%rbp)
mov    -0x10(%rbp),%rax
movss  0xe85(%rip),%xmm0
movss  %xmm0,(%rax)
movss  -0x4(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    126c <malloc_gaussian_symm_kernel+0x13c>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x14(%rbp)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     122b <malloc_gaussian_symm_kernel+0xfb>
cvtsi2ssl -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe47(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsi2ssl -0x18(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
movss  -0x4(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
movss  -0x4(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
call   1050 <exp@plt>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x14(%rbp),%xmm1
addss  %xmm0,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x14(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11a4 <malloc_gaussian_symm_kernel+0x74>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     1267 <malloc_gaussian_symm_kernel+0x137>
movss  -0x14(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1232 <malloc_gaussian_symm_kernel+0x102>
jmp    12a2 <malloc_gaussian_symm_kernel+0x172>
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     129d <malloc_gaussian_symm_kernel+0x16d>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1273 <malloc_gaussian_symm_kernel+0x143>
jmp    12a2 <malloc_gaussian_symm_kernel+0x172>
mov    -0x10(%rbp),%rax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<malloc_gaussian_symm_kernel>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %edi,%r14d
movss  %xmm0,(%rsp)
xorps  %xmm1,%xmm1
xor    %edi,%edi
ucomiss %xmm1,%xmm0
setae  %dil
call   1050 <assert@plt>
lea    0x4(,%r14,4),%edi
call   1040 <xmalloc@plt>
movss  (%rsp),%xmm0
mov    %rax,%r12
ucomiss 0xe84(%rip),%xmm0
movl   $0x3f800000,(%rax)
jbe    1209 <malloc_gaussian_symm_kernel+0xc9>
test   %r14d,%r14d
jle    122c <malloc_gaussian_symm_kernel+0xec>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
lea    0x1(%r14),%r15d
movss  0xe55(%rip),%xmm3
mov    $0x1,%ebx
movaps %xmm3,(%rsp)
xorps  %xmm0,%xmm0
cvtsi2ss %ebx,%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movaps %xmm1,%xmm0
mulsd  0xe3b(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0x10(%rsp),%xmm1
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
call   1060 <exp@plt>
movaps (%rsp),%xmm3
movss  %xmm0,(%r12,%rbx,4)
addss  %xmm0,%xmm0
addss  %xmm0,%xmm3
add    $0x1,%rbx
cmp    %rbx,%r15
jne    11b0 <malloc_gaussian_symm_kernel+0x70>
test   %r14d,%r14d
jns    123d <malloc_gaussian_symm_kernel+0xfd>
jmp    12d9 <malloc_gaussian_symm_kernel+0x199>
test   %r14d,%r14d
jle    12d9 <malloc_gaussian_symm_kernel+0x199>
mov    %r12,%rdi
add    $0x4,%rdi
mov    %r14d,%edx
shl    $0x2,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
jmp    12d9 <malloc_gaussian_symm_kernel+0x199>
movss  0xdcc(%rip),%xmm3
test   %r14d,%r14d
js     12d9 <malloc_gaussian_symm_kernel+0x199>
lea    0x1(%r14),%eax
cmp    $0x3,%r14d
jae    124b <malloc_gaussian_symm_kernel+0x10b>
xor    %ecx,%ecx
jmp    12c0 <malloc_gaussian_symm_kernel+0x180>
mov    %eax,%ecx
and    $0xfffffffc,%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     12e8 <malloc_gaussian_symm_kernel+0x1a8>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movups (%r12,%rsi,4),%xmm1
movups 0x10(%r12,%rsi,4),%xmm2
divps  %xmm0,%xmm1
movups %xmm1,(%r12,%rsi,4)
divps  %xmm0,%xmm2
movups %xmm2,0x10(%r12,%rsi,4)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1280 <malloc_gaussian_symm_kernel+0x140>
test   $0x1,%dl
je     12b8 <malloc_gaussian_symm_kernel+0x178>
movups (%r12,%rsi,4),%xmm1
divps  %xmm0,%xmm1
movups %xmm1,(%r12,%rsi,4)
cmp    %rax,%rcx
je     12d9 <malloc_gaussian_symm_kernel+0x199>
nopl   (%rax)
movss  (%r12,%rcx,4),%xmm0
divss  %xmm3,%xmm0
movss  %xmm0,(%r12,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12c0 <malloc_gaussian_symm_kernel+0x180>
mov    %r12,%rax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
xor    %esi,%esi
test   $0x1,%dl
jne    12ab <malloc_gaussian_symm_kernel+0x16b>
jmp    12b8 <malloc_gaussian_symm_kernel+0x178>

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

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