compute_weights(weights, w) float weights[21][21];
int w;
{
  int i, j, offset, weighted_size;
  float kernel[21], term1, term2, sum;
  weighted_size = 2 * w + 1;
  if (w == 1) {
    kernel[0] = 0.25;
    kernel[1] = 0.5;
    kernel[2] = 0.25;
    sum = 0.0;
    for (i = 0; i < (weighted_size); i++)
      for (j = 0; j < (weighted_size); j++) {
        weights[i][j] = kernel[i] * kernel[j];
        sum += weights[i][j];
      }
    if (sum != 1.0) {
      printf("Fatal error: sum of Gaussian coefficients is zero\n");
      exit(1);
    }
  } else if (w == 2) {
    kernel[0] = 0.0625;
    kernel[1] = 0.25;
    kernel[2] = 0.375;
    kernel[3] = 0.25;
    kernel[4] = 0.0625;
    sum = 0.0;
    for (i = 0; i < (weighted_size); i++)
      for (j = 0; j < (weighted_size); j++) {
        weights[i][j] = kernel[i] * kernel[j];
        sum += weights[i][j];
      }
    if (sum != 1.0) {
      printf("Fatal error: sum of Gaussian coefficients is zero\n");
      exit(1);
    }
  } else {
    printf("\nFatal error: window size not 1 or 2: w=%d\n", w);
    exit(1);
  }
  printf("\nWeights:\n");
  for (i = 0; i < (weighted_size); i++) {
    for (j = 0; j < (weighted_size); j++)
      printf("%8.5f ", weights[i][j]);
    printf("\n");
  }
}


x86 assembly (clang, -O0):
<compute_weights>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x1,-0x14(%rbp)
jne    125b <compute_weights+0x13b>
movss  0xeb4(%rip),%xmm0
movss  %xmm0,-0x80(%rbp)
movss  0xeaf(%rip),%xmm0
movss  %xmm0,-0x7c(%rbp)
movss  0xe9a(%rip),%xmm0
movss  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1215 <compute_weights+0xf5>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1202 <compute_weights+0xe2>
movslq -0x18(%rbp),%rax
movss  -0x80(%rbp,%rax,4),%xmm0
movslq -0x1c(%rbp),%rax
mulss  -0x80(%rbp,%rax,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x8c(%rbp),%xmm0
ff
movss  %xmm0,-0x8c(%rbp)
ff
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1194 <compute_weights+0x74>
jmp    1207 <compute_weights+0xe7>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1181 <compute_weights+0x61>
movss  -0x8c(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
movsd  0xde7(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    123e <compute_weights+0x11e>
jp     123e <compute_weights+0x11e>
jmp    1256 <compute_weights+0x136>
lea    0xdd3(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>
jmp    13b2 <compute_weights+0x292>
cmpl   $0x2,-0x14(%rbp)
jne    1392 <compute_weights+0x272>
movss  0xd93(%rip),%xmm0
movss  %xmm0,-0x80(%rbp)
movss  0xd8a(%rip),%xmm0
movss  %xmm0,-0x7c(%rbp)
movss  0xd81(%rip),%xmm0
movss  %xmm0,-0x78(%rbp)
movss  0xd70(%rip),%xmm0
movss  %xmm0,-0x74(%rbp)
movss  0xd5f(%rip),%xmm0
movss  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    134c <compute_weights+0x22c>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1339 <compute_weights+0x219>
movslq -0x18(%rbp),%rax
movss  -0x80(%rbp,%rax,4),%xmm0
movslq -0x1c(%rbp),%rax
mulss  -0x80(%rbp,%rax,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x8c(%rbp),%xmm0
ff
movss  %xmm0,-0x8c(%rbp)
ff
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    12cb <compute_weights+0x1ab>
jmp    133e <compute_weights+0x21e>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    12b8 <compute_weights+0x198>
movss  -0x8c(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
movsd  0xcb0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1375 <compute_weights+0x255>
jp     1375 <compute_weights+0x255>
jmp    138d <compute_weights+0x26d>
lea    0xc9c(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>
jmp    13ad <compute_weights+0x28d>
mov    -0x14(%rbp),%esi
lea    0xcaf(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>
jmp    13b2 <compute_weights+0x292>
lea    0xcbe(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    143a <compute_weights+0x31a>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    141e <compute_weights+0x2fe>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc79(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13da <compute_weights+0x2ba>
lea    0xc5b(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    13c7 <compute_weights+0x2a7>
mov    -0x4(%rbp),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_weights>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rdi,%r14
lea    (%rsi,%rsi,1),%ebp
add    $0x1,%ebp
mov    %ebp,%r15d
cmp    $0x2,%esi
je     1211 <compute_weights+0xd1>
cmp    $0x1,%esi
jne    1349 <compute_weights+0x209>
movabs $0x3f0000003e800000,%rax
mov    %rax,(%rsp)
movl   $0x3e800000,0x8(%rsp)
3e
test   %ebp,%ebp
jle    1333 <compute_weights+0x1f3>
lea    0x4(%r14),%rax
lea    -0x1(%r15),%rcx
xorps  %xmm0,%xmm0
xor    %edx,%edx
jmp    11cc <compute_weights+0x8c>
nopl   0x0(%rax)
mulss  (%rsp,%rdi,4),%xmm1
movaps %xmm0,%xmm2
imul   $0x54,%rdx,%rbx
add    %r14,%rbx
movss  %xmm1,(%rbx,%rdi,4)
movaps %xmm1,%xmm0
addss  %xmm2,%xmm0
add    $0x1,%rdx
add    $0x54,%rax
cmp    %r15,%rdx
je     12b1 <compute_weights+0x171>
movss  (%rsp,%rdx,4),%xmm1
xor    %edi,%edi
test   %esi,%esi
je     11a0 <compute_weights+0x60>
nopw   0x0(%rax,%rax,1)
movss  (%rsp,%rdi,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,-0x4(%rax,%rdi,4)
addss  %xmm0,%xmm2
movss  0x4(%rsp,%rdi,4),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rdi,4)
addss  %xmm2,%xmm0
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    11e0 <compute_weights+0xa0>
jmp    11a0 <compute_weights+0x60>
movaps 0xde8(%rip),%xmm0
movaps %xmm0,(%rsp)
movl   $0x3d800000,0x10(%rsp)
3d
test   %ebp,%ebp
jle    1333 <compute_weights+0x1f3>
lea    0x4(%r14),%rax
lea    -0x1(%r15),%rcx
xorps  %xmm0,%xmm0
xor    %edx,%edx
jmp    1268 <compute_weights+0x128>
nopl   0x0(%rax,%rax,1)
mulss  (%rsp,%rdi,4),%xmm1
movaps %xmm0,%xmm2
imul   $0x54,%rdx,%rbx
add    %r14,%rbx
movss  %xmm1,(%rbx,%rdi,4)
movaps %xmm1,%xmm0
addss  %xmm2,%xmm0
add    $0x1,%rdx
add    $0x54,%rax
cmp    %r15,%rdx
je     12b1 <compute_weights+0x171>
movss  (%rsp,%rdx,4),%xmm1
xor    %edi,%edi
test   %esi,%esi
je     1240 <compute_weights+0x100>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rsp,%rdi,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,-0x4(%rax,%rdi,4)
addss  %xmm0,%xmm2
movss  0x4(%rsp,%rdi,4),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rdi,4)
addss  %xmm2,%xmm0
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    1280 <compute_weights+0x140>
jmp    1240 <compute_weights+0x100>
ucomiss 0xd58(%rip),%xmm0
jne    1333 <compute_weights+0x1f3>
jp     1333 <compute_weights+0x1f3>
lea    0xd84(%rip),%rdi
call   1040 <puts@plt>
test   %ebp,%ebp
jle    1324 <compute_weights+0x1e4>
xor    %r12d,%r12d
lea    0xd6a(%rip),%r13
cs nopw 0x0(%rax,%rax,1)
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%r14,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1050 <printf@plt>
add    $0x1,%rbx
cmp    %rbx,%r15
jne    12f0 <compute_weights+0x1b0>
mov    $0xa,%edi
call   1030 <putchar@plt>
add    $0x1,%r12
add    $0x54,%r14
cmp    %r15,%r12
jne    12e0 <compute_weights+0x1a0>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0xd17(%rip),%rdi
call   1040 <puts@plt>
mov    $0x1,%edi
call   1060 <exit@plt>
lea    0xcc4(%rip),%rdi
xor    %eax,%eax
call   1050 <printf@plt>
mov    $0x1,%edi
call   1060 <exit@plt>

LLVM IR (-O0):
@.str = [51 x i8] c"Fatal error: sum of Gaussian coefficients is zero\0A\00"
@.str.1 = [44 x i8] c"\0AFatal error: window size not 1 or 2: w=%d\0A\00"
@.str.2 = [11 x i8] c"\0AWeights:\0A\00"
@.str.3 = [7 x i8] c"%8.5f \00"
@.str.4 = [2 x i8] c"\0A\00"
define i32 @compute_weights([21 x float]* %0, i32 %1) {
%3 = alloca i32
%4 = alloca [21 x float]*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca [21 x float]
%11 = alloca float
%12 = alloca float
%13 = alloca float
store [21 x float]* %0, [21 x float]** %4
store i32 %1, i32* %5
%14 = load i32, i32* %5
%15 = mul nsw i32 2, %14
%16 = add nsw i32 %15, 1
store i32 %16, i32* %9
%17 = load i32, i32* %5
%18 = icmp eq i32 %17, 1
br i1 %18, label %19, label %73
19: ; preds = %2
%20 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 0
store float 2.500000e-01, float* %20
%21 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 1
store float 5.000000e-01, float* %21
%22 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 2
store float 2.500000e-01, float* %22
store float 0.000000e+00, float* %13
store i32 0, i32* %6
br label %23
23: ; preds = %63, %19
%24 = load i32, i32* %6
%25 = load i32, i32* %9
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %66
27: ; preds = %23
store i32 0, i32* %7
br label %28
28: ; preds = %59, %27
%29 = load i32, i32* %7
%30 = load i32, i32* %9
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %62
32: ; preds = %28
%33 = load i32, i32* %6
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 %34
%36 = load float, float* %35
%37 = load i32, i32* %7
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 %38
%40 = load float, float* %39
%41 = fmul float %36, %40
%42 = load [21 x float]*, [21 x float]** %4
%43 = load i32, i32* %6
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds [21 x float], [21 x float]* %42, i64 %44
%46 = load i32, i32* %7
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds [21 x float], [21 x float]* %45, i64 0, i64 %47
store float %41, float* %48
%49 = load [21 x float]*, [21 x float]** %4
%50 = load i32, i32* %6
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [21 x float], [21 x float]* %49, i64 %51
%53 = load i32, i32* %7
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds [21 x float], [21 x float]* %52, i64 0, i64 %54
%56 = load float, float* %55
%57 = load float, float* %13
%58 = fadd float %57, %56
store float %58, float* %13
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
br label %63
63: ; preds = %62
%64 = load i32, i32* %6
%65 = add nsw i32 %64, 1
store i32 %65, i32* %6
br label %23
66: ; preds = %23
%67 = load float, float* %13
%68 = fpext float %67 to double
%69 = fcmp une double %68, 1.000000e+00
br i1 %69, label %70, label %72
70: ; preds = %66
%71 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str, i64 0, i64 0))
unreachable
72: ; preds = %66
br label %136
73: ; preds = %2
%74 = load i32, i32* %5
%75 = icmp eq i32 %74, 2
br i1 %75, label %76, label %132
76: ; preds = %73
%77 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 0
store float 6.250000e-02, float* %77
%78 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 1
store float 2.500000e-01, float* %78
%79 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 2
store float 3.750000e-01, float* %79
%80 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 3
store float 2.500000e-01, float* %80
%81 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 4
store float 6.250000e-02, float* %81
store float 0.000000e+00, float* %13
store i32 0, i32* %6
br label %82
82: ; preds = %122, %76
%83 = load i32, i32* %6
%84 = load i32, i32* %9
%85 = icmp slt i32 %83, %84
br i1 %85, label %86, label %125
86: ; preds = %82
store i32 0, i32* %7
br label %87
87: ; preds = %118, %86
%88 = load i32, i32* %7
%89 = load i32, i32* %9
%90 = icmp slt i32 %88, %89
br i1 %90, label %91, label %121
91: ; preds = %87
%92 = load i32, i32* %6
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 %93
%95 = load float, float* %94
%96 = load i32, i32* %7
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds [21 x float], [21 x float]* %10, i64 0, i64 %97
%99 = load float, float* %98
%100 = fmul float %95, %99
%101 = load [21 x float]*, [21 x float]** %4
%102 = load i32, i32* %6
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds [21 x float], [21 x float]* %101, i64 %103
%105 = load i32, i32* %7
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds [21 x float], [21 x float]* %104, i64 0, i64 %106
store float %100, float* %107
%108 = load [21 x float]*, [21 x float]** %4
%109 = load i32, i32* %6
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds [21 x float], [21 x float]* %108, i64 %110
%112 = load i32, i32* %7
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds [21 x float], [21 x float]* %111, i64 0, i64 %113
%115 = load float, float* %114
%116 = load float, float* %13
%117 = fadd float %116, %115
store float %117, float* %13
br label %118
118: ; preds = %91
%119 = load i32, i32* %7
%120 = add nsw i32 %119, 1
store i32 %120, i32* %7
br label %87
121: ; preds = %87
br label %122
122: ; preds = %121
%123 = load i32, i32* %6
%124 = add nsw i32 %123, 1
store i32 %124, i32* %6
br label %82
125: ; preds = %82
%126 = load float, float* %13
%127 = fpext float %126 to double
%128 = fcmp une double %127, 1.000000e+00
br i1 %128, label %129, label %131
129: ; preds = %125
%130 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str, i64 0, i64 0))
unreachable
131: ; preds = %125
br label %135
132: ; preds = %73
%133 = load i32, i32* %5
%134 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([44 x i8], [44 x i8]* @.str.1, i64 0, i64 0), i32 %133)
unreachable
135: ; preds = %131
br label %136
136: ; preds = %135, %72
%137 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str.2, i64 0, i64 0))
store i32 0, i32* %6
br label %138
138: ; preds = %163, %136
%139 = load i32, i32* %6
%140 = load i32, i32* %9
%141 = icmp slt i32 %139, %140
br i1 %141, label %142, label %166
142: ; preds = %138
store i32 0, i32* %7
br label %143
143: ; preds = %158, %142
%144 = load i32, i32* %7
%145 = load i32, i32* %9
%146 = icmp slt i32 %144, %145
br i1 %146, label %147, label %161
147: ; preds = %143
%148 = load [21 x float]*, [21 x float]** %4
%149 = load i32, i32* %6
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds [21 x float], [21 x float]* %148, i64 %150
%152 = load i32, i32* %7
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds [21 x float], [21 x float]* %151, i64 0, i64 %153
%155 = load float, float* %154
%156 = fpext float %155 to double
%157 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([7 x i8], [7 x i8]* @.str.3, i64 0, i64 0), double %156)
br label %158
158: ; preds = %147
%159 = load i32, i32* %7
%160 = add nsw i32 %159, 1
store i32 %160, i32* %7
br label %143
161: ; preds = %143
%162 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.4, i64 0, i64 0))
br label %163
163: ; preds = %161
%164 = load i32, i32* %6
%165 = add nsw i32 %164, 1
store i32 %165, i32* %6
br label %138
166: ; preds = %138
%167 = load i32, i32* %3
ret i32 %167
}

LLVM IR (-O2):
@.str.1 = [44 x i8] c"\0AFatal error: window size not 1 or 2: w=%d\0A\00"
@.str.3 = [7 x i8] c"%8.5f \00"
@str = [10 x i8] c"\0AWeights:\00"
@str.6 = [50 x i8] c"Fatal error: sum of Gaussian coefficients is zero\00"
define i32 @compute_weights([21 x float]* nocapture %0, i32 %1) {
%3 = alloca [21 x float]
%4 = bitcast [21 x float]* %3 to i8*
%5 = shl nsw i32 %1, 1
%6 = or i32 %5, 1
switch i32 %1, label %100 [
i32 1, label %7
i32 2, label %54
]
7: ; preds = %2
%8 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 0
store float 2.500000e-01, float* %8
%9 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 1
store float 5.000000e-01, float* %9
%10 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 2
store float 2.500000e-01, float* %10
%11 = icmp sgt i32 %6, 0
br i1 %11, label %12, label %52
12: ; preds = %7
%13 = zext i32 %6 to i64
%14 = icmp eq i32 %1, 0
%15 = add nsw i64 %13, -1
br label %16
16: ; preds = %39, %12
%17 = phi i64 [ 0, %12 ], [ %48, %39 ]
%18 = phi float [ 0.000000e+00, %12 ], [ %47, %39 ]
%19 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %17
%20 = load float, float* %19
br i1 %14, label %39, label %21
21: ; preds = %16, %21
%22 = phi i64 [ %36, %21 ], [ 0, %16 ]
%23 = phi float [ %35, %21 ], [ %18, %16 ]
%24 = phi i64 [ %37, %21 ], [ 0, %16 ]
%25 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %22
%26 = load float, float* %25
%27 = fmul float %20, %26
%28 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %17, i64 %22
store float %27, float* %28
%29 = fadd float %23, %27
%30 = or i64 %22, 1
%31 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %30
%32 = load float, float* %31
%33 = fmul float %20, %32
%34 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %17, i64 %30
store float %33, float* %34
%35 = fadd float %29, %33
%36 = add nuw nsw i64 %22, 2
%37 = add i64 %24, 2
%38 = icmp eq i64 %37, %15
br i1 %38, label %39, label %21
39: ; preds = %21, %16
%40 = phi float [ undef, %16 ], [ %35, %21 ]
%41 = phi i64 [ 0, %16 ], [ %36, %21 ]
%42 = phi float [ %18, %16 ], [ %35, %21 ]
%43 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %41
%44 = load float, float* %43
%45 = fmul float %20, %44
%46 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %17, i64 %41
store float %45, float* %46
%47 = fadd float %42, %45
%48 = add nuw nsw i64 %17, 1
%49 = icmp eq i64 %48, %13
br i1 %49, label %50, label %16
50: ; preds = %39
%51 = fcmp une float %47, 1.000000e+00
br i1 %51, label %52, label %102
52: ; preds = %7, %50
%53 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([50 x i8], [50 x i8]* @str.6, i64 0, i64 0))
tail call void @exit(i32 1)
unreachable
54: ; preds = %2
%55 = bitcast [21 x float]* %3 to <4 x float>*
store <4 x float> <float 6.250000e-02, float 2.500000e-01, float 3.750000e-01, float 2.500000e-01>, <4 x float>* %55
%56 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 4
store float 6.250000e-02, float* %56
%57 = icmp sgt i32 %6, 0
br i1 %57, label %58, label %98
58: ; preds = %54
%59 = zext i32 %6 to i64
%60 = icmp eq i32 %1, 0
%61 = add nsw i64 %59, -1
br label %62
62: ; preds = %85, %58
%63 = phi i64 [ 0, %58 ], [ %94, %85 ]
%64 = phi float [ 0.000000e+00, %58 ], [ %93, %85 ]
%65 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %63
%66 = load float, float* %65
br i1 %60, label %85, label %67
67: ; preds = %62, %67
%68 = phi i64 [ %82, %67 ], [ 0, %62 ]
%69 = phi float [ %81, %67 ], [ %64, %62 ]
%70 = phi i64 [ %83, %67 ], [ 0, %62 ]
%71 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %68
%72 = load float, float* %71
%73 = fmul float %66, %72
%74 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %63, i64 %68
store float %73, float* %74
%75 = fadd float %69, %73
%76 = or i64 %68, 1
%77 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %76
%78 = load float, float* %77
%79 = fmul float %66, %78
%80 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %63, i64 %76
store float %79, float* %80
%81 = fadd float %75, %79
%82 = add nuw nsw i64 %68, 2
%83 = add i64 %70, 2
%84 = icmp eq i64 %83, %61
br i1 %84, label %85, label %67
85: ; preds = %67, %62
%86 = phi float [ undef, %62 ], [ %81, %67 ]
%87 = phi i64 [ 0, %62 ], [ %82, %67 ]
%88 = phi float [ %64, %62 ], [ %81, %67 ]
%89 = getelementptr inbounds [21 x float], [21 x float]* %3, i64 0, i64 %87
%90 = load float, float* %89
%91 = fmul float %66, %90
%92 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %63, i64 %87
store float %91, float* %92
%93 = fadd float %88, %91
%94 = add nuw nsw i64 %63, 1
%95 = icmp eq i64 %94, %59
br i1 %95, label %96, label %62
96: ; preds = %85
%97 = fcmp une float %93, 1.000000e+00
br i1 %97, label %98, label %102
98: ; preds = %54, %96
%99 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([50 x i8], [50 x i8]* @str.6, i64 0, i64 0))
tail call void @exit(i32 1)
unreachable
100: ; preds = %2
%101 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([44 x i8], [44 x i8]* @.str.1, i64 0, i64 0), i32 %1)
tail call void @exit(i32 1)
unreachable
102: ; preds = %96, %50
%103 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([10 x i8], [10 x i8]* @str, i64 0, i64 0))
%104 = icmp sgt i32 %6, 0
br i1 %104, label %105, label %121
105: ; preds = %102
%106 = zext i32 %6 to i64
br label %107
107: ; preds = %117, %105
%108 = phi i64 [ 0, %105 ], [ %119, %117 ]
br label %109
109: ; preds = %107, %109
%110 = phi i64 [ 0, %107 ], [ %115, %109 ]
%111 = getelementptr inbounds [21 x float], [21 x float]* %0, i64 %108, i64 %110
%112 = load float, float* %111
%113 = fpext float %112 to double
%114 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.3, i64 0, i64 0), double %113)
%115 = add nuw nsw i64 %110, 1
%116 = icmp eq i64 %115, %106
br i1 %116, label %117, label %109
117: ; preds = %109
%118 = tail call i32 @putchar(i32 10)
%119 = add nuw nsw i64 %108, 1
%120 = icmp eq i64 %119, %106
br i1 %120, label %121, label %107
121: ; preds = %117, %102
ret i32 undef
}