void separable_convolution(int *image, int width, int height, int depth,
                           float *gau, int gau_width, float *smx, float *smy) {
  int i, j, k, d, I1, I2;
  float x, y;
  for (d = 0; d < depth; d++) {
    for (i = 0; i < height; i++) {
      for (j = 0; j < width; j++) {
        x = gau[0] * image[((height) * (j) + (i) + (width) * (height) * (d))];
        y = gau[0] * image[((height) * (j) + (i) + (width) * (height) * (d))];
        for (k = 1; k < gau_width; k++) {
          I1 = (i + k) % height;
          I2 = (i - k + height) % height;
          y = y +
              gau[k] *
                  image[((height) * (j) + (I1) + (width) * (height) * (d))] +
              gau[k] *
                  image[((height) * (j) + (I2) + (width) * (height) * (d))];
          I1 = (j + k) % width;
          I2 = (j - k + width) % width;
          x = x +
              gau[k] *
                  image[((height) * (I1) + (i) + (width) * (height) * (d))] +
              gau[k] *
                  image[((height) * (I2) + (i) + (width) * (height) * (d))];
        }
        smx[((height) * (j) + (i) + (width) * (height) * (d))] = x;
        smy[((height) * (j) + (i) + (width) * (height) * (d))] = y;
      }
    }
  }
}


x86 assembly (clang, -O0):
<separable_convolution>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1395 <separable_convolution+0x295>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1382 <separable_convolution+0x282>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    136f <separable_convolution+0x26f>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1311 <separable_convolution+0x211>
mov    -0x28(%rbp),%eax
add    -0x30(%rbp),%eax
cltd
idivl  -0x10(%rbp)
mov    %edx,-0x38(%rbp)
mov    -0x28(%rbp),%eax
sub    -0x30(%rbp),%eax
add    -0x10(%rbp),%eax
cltd
idivl  -0x10(%rbp)
mov    %edx,-0x3c(%rbp)
movss  -0x44(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x38(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x3c(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x30(%rbp),%eax
cltd
idivl  -0xc(%rbp)
mov    %edx,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
sub    -0x30(%rbp),%eax
add    -0xc(%rbp),%eax
cltd
idivl  -0xc(%rbp)
mov    %edx,-0x3c(%rbp)
movss  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x38(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x3c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11c9 <separable_convolution+0xc9>
movss  -0x40(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x44(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    114e <separable_convolution+0x4e>
jmp    1374 <separable_convolution+0x274>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113b <separable_convolution+0x3b>
jmp    1387 <separable_convolution+0x287>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1128 <separable_convolution+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<separable_convolution>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x58(%rsp)
mov    %ecx,-0x40(%rsp)
test   %ecx,%ecx
jle    1354 <separable_convolution+0x254>
mov    %r9d,%r11d
mov    %edx,%ebp
mov    0x40(%rsp),%r9
mov    0x38(%rsp),%r14
mov    %edx,%eax
mov    -0x58(%rsp),%rcx
imul   %ecx,%eax
mov    %eax,-0x44(%rsp)
mov    %edx,%eax
mov    %rax,-0x18(%rsp)
mov    %ecx,%r15d
mov    %r11d,%esi
lea    -0x1(%rbp),%eax
mov    %eax,-0x48(%rsp)
lea    -0x1(%rcx),%eax
mov    %eax,-0x3c(%rsp)
add    $0xffffffffffffffff,%rsi
xor    %eax,%eax
mov    %rax,-0x68(%rsp)
mov    %r11d,-0x38(%rsp)
mov    %r15,-0x10(%rsp)
jmp    1187 <separable_convolution+0x87>
nopw   0x0(%rax,%rax,1)
mov    -0x68(%rsp),%rax
add    $0x1,%eax
mov    %rax,-0x68(%rsp)
cmp    -0x40(%rsp),%eax
je     1354 <separable_convolution+0x254>
test   %ebp,%ebp
jle    1170 <separable_convolution+0x70>
mov    -0x44(%rsp),%eax
imul   -0x68(%rsp),%eax
mov    %rax,-0x30(%rsp)
mov    $0x1,%ecx
mov    -0x48(%rsp),%eax
mov    %eax,-0x6c(%rsp)
xor    %eax,%eax
mov    %rax,-0x60(%rsp)
jmp    11d1 <separable_convolution+0xd1>
nop
mov    -0x60(%rsp),%rdx
add    $0x1,%rdx
addl   $0x1,-0x6c(%rsp)
add    $0x1,%rcx
mov    %rdx,%rax
mov    %rdx,-0x60(%rsp)
cmp    -0x18(%rsp),%rdx
je     1170 <separable_convolution+0x70>
cmpl   $0x0,-0x58(%rsp)
jle    11b0 <separable_convolution+0xb0>
mov    -0x30(%rsp),%rax
mov    -0x60(%rsp),%rdx
lea    (%rax,%rdx,1),%r13d
mov    -0x3c(%rsp),%edx
xor    %ebx,%ebx
mov    %r13,-0x28(%rsp)
jmp    121c <separable_convolution+0x11c>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  %xmm0,(%r14,%rax,4)
shufps $0x55,%xmm0,%xmm0
movss  %xmm0,(%r9,%rax,4)
add    $0x1,%rbx
add    $0x1,%edx
cmp    %r15,%rbx
je     11b0 <separable_convolution+0xb0>
mov    %ebx,%eax
imul   %ebp,%eax
mov    %rax,-0x50(%rsp)
add    %r13d,%eax
cltq
xorps  %xmm0,%xmm0
cvtsi2ssl (%rdi,%rax,4),%xmm0
mulss  (%r8),%xmm0
shufps $0x0,%xmm0,%xmm0
cmp    $0x2,%r11d
jl     1200 <separable_convolution+0x100>
mov    %rax,-0x8(%rsp)
mov    -0x50(%rsp),%rax
add    -0x30(%rsp),%eax
mov    %rax,-0x50(%rsp)
xor    %r12d,%r12d
mov    %edx,-0x34(%rsp)
mov    %edx,%r11d
mov    -0x6c(%rsp),%r14d
mov    %rbx,-0x20(%rsp)
nopl   0x0(%rax)
lea    (%rcx,%r12,1),%eax
cltd
idiv   %ebp
mov    %rsi,%r9
mov    %rbp,%rsi
mov    %edx,%r15d
mov    %r14d,%eax
cltd
idiv   %esi
mov    %edx,%ebx
mov    -0x20(%rsp),%rax
add    %r12d,%eax
add    $0x1,%eax
cltd
mov    -0x58(%rsp),%r10
idiv   %r10d
mov    %edx,%ebp
mov    %r11d,%eax
cltd
idiv   %r10d
imul   %esi,%ebp
mov    -0x28(%rsp),%r10
add    %r10d,%ebp
movslq %ebp,%rax
movss  (%rdi,%rax,4),%xmm1
mov    -0x50(%rsp),%r13
add    %r13d,%r15d
movslq %r15d,%rax
mov    %rsi,%rbp
mov    %r9,%rsi
movss  (%rdi,%rax,4),%xmm2
add    %r13d,%ebx
movslq %ebx,%rax
movss  (%rdi,%rax,4),%xmm3
imul   %ebp,%edx
add    %r10d,%edx
movslq %edx,%rax
movss  (%rdi,%rax,4),%xmm4
unpcklps %xmm2,%xmm1
cvtdq2ps %xmm1,%xmm1
movss  0x4(%r8,%r12,4),%xmm5
shufps $0x0,%xmm5,%xmm5
mulps  %xmm5,%xmm1
addps  %xmm0,%xmm1
unpcklps %xmm3,%xmm4
cvtdq2ps %xmm4,%xmm2
mulps  %xmm5,%xmm2
movaps %xmm2,%xmm0
addps  %xmm1,%xmm0
add    $0xffffffff,%r14d
add    $0xffffffff,%r11d
add    $0x1,%r12
cmp    %r12,%r9
jne    1270 <separable_convolution+0x170>
addps  %xmm1,%xmm2
movaps %xmm2,%xmm0
mov    -0x38(%rsp),%r11d
mov    0x40(%rsp),%r9
mov    0x38(%rsp),%r14
mov    -0x10(%rsp),%r15
mov    -0x28(%rsp),%r13
mov    -0x34(%rsp),%edx
mov    -0x20(%rsp),%rbx
mov    -0x8(%rsp),%rax
jmp    1200 <separable_convolution+0x100>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @separable_convolution(i32* %0, i32 %1, i32 %2, i32 %3, float* %4, i32 %5, float* %6, float* %7) {
%9 = alloca i32*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca float*
%14 = alloca i32
%15 = alloca float*
%16 = alloca float*
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca float
%24 = alloca float
store i32* %0, i32** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store float* %4, float** %13
store i32 %5, i32* %14
store float* %6, float** %15
store float* %7, float** %16
store i32 0, i32* %20
br label %25
25: ; preds = %241, %8
%26 = load i32, i32* %20
%27 = load i32, i32* %12
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %244
29: ; preds = %25
store i32 0, i32* %17
br label %30
30: ; preds = %237, %29
%31 = load i32, i32* %17
%32 = load i32, i32* %11
%33 = icmp slt i32 %31, %32
br i1 %33, label %34, label %240
34: ; preds = %30
store i32 0, i32* %18
br label %35
35: ; preds = %233, %34
%36 = load i32, i32* %18
%37 = load i32, i32* %10
%38 = icmp slt i32 %36, %37
br i1 %38, label %39, label %236
39: ; preds = %35
%40 = load float*, float** %13
%41 = getelementptr inbounds float, float* %40, i64 0
%42 = load float, float* %41
%43 = load i32*, i32** %9
%44 = load i32, i32* %11
%45 = load i32, i32* %18
%46 = mul nsw i32 %44, %45
%47 = load i32, i32* %17
%48 = add nsw i32 %46, %47
%49 = load i32, i32* %10
%50 = load i32, i32* %11
%51 = mul nsw i32 %49, %50
%52 = load i32, i32* %20
%53 = mul nsw i32 %51, %52
%54 = add nsw i32 %48, %53
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds i32, i32* %43, i64 %55
%57 = load i32, i32* %56
%58 = sitofp i32 %57 to float
%59 = fmul float %42, %58
store float %59, float* %23
%60 = load float*, float** %13
%61 = getelementptr inbounds float, float* %60, i64 0
%62 = load float, float* %61
%63 = load i32*, i32** %9
%64 = load i32, i32* %11
%65 = load i32, i32* %18
%66 = mul nsw i32 %64, %65
%67 = load i32, i32* %17
%68 = add nsw i32 %66, %67
%69 = load i32, i32* %10
%70 = load i32, i32* %11
%71 = mul nsw i32 %69, %70
%72 = load i32, i32* %20
%73 = mul nsw i32 %71, %72
%74 = add nsw i32 %68, %73
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds i32, i32* %63, i64 %75
%77 = load i32, i32* %76
%78 = sitofp i32 %77 to float
%79 = fmul float %62, %78
store float %79, float* %24
store i32 1, i32* %19
br label %80
80: ; preds = %199, %39
%81 = load i32, i32* %19
%82 = load i32, i32* %14
%83 = icmp slt i32 %81, %82
br i1 %83, label %84, label %202
84: ; preds = %80
%85 = load i32, i32* %17
%86 = load i32, i32* %19
%87 = add nsw i32 %85, %86
%88 = load i32, i32* %11
%89 = srem i32 %87, %88
store i32 %89, i32* %21
%90 = load i32, i32* %17
%91 = load i32, i32* %19
%92 = sub nsw i32 %90, %91
%93 = load i32, i32* %11
%94 = add nsw i32 %92, %93
%95 = load i32, i32* %11
%96 = srem i32 %94, %95
store i32 %96, i32* %22
%97 = load float, float* %24
%98 = load float*, float** %13
%99 = load i32, i32* %19
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %98, i64 %100
%102 = load float, float* %101
%103 = load i32*, i32** %9
%104 = load i32, i32* %11
%105 = load i32, i32* %18
%106 = mul nsw i32 %104, %105
%107 = load i32, i32* %21
%108 = add nsw i32 %106, %107
%109 = load i32, i32* %10
%110 = load i32, i32* %11
%111 = mul nsw i32 %109, %110
%112 = load i32, i32* %20
%113 = mul nsw i32 %111, %112
%114 = add nsw i32 %108, %113
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds i32, i32* %103, i64 %115
%117 = load i32, i32* %116
%118 = sitofp i32 %117 to float
%119 = call float @llvm.fmuladd.f32(float %102, float %118, float %97)
%120 = load float*, float** %13
%121 = load i32, i32* %19
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds float, float* %120, i64 %122
%124 = load float, float* %123
%125 = load i32*, i32** %9
%126 = load i32, i32* %11
%127 = load i32, i32* %18
%128 = mul nsw i32 %126, %127
%129 = load i32, i32* %22
%130 = add nsw i32 %128, %129
%131 = load i32, i32* %10
%132 = load i32, i32* %11
%133 = mul nsw i32 %131, %132
%134 = load i32, i32* %20
%135 = mul nsw i32 %133, %134
%136 = add nsw i32 %130, %135
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds i32, i32* %125, i64 %137
%139 = load i32, i32* %138
%140 = sitofp i32 %139 to float
%141 = call float @llvm.fmuladd.f32(float %124, float %140, float %119)
store float %141, float* %24
%142 = load i32, i32* %18
%143 = load i32, i32* %19
%144 = add nsw i32 %142, %143
%145 = load i32, i32* %10
%146 = srem i32 %144, %145
store i32 %146, i32* %21
%147 = load i32, i32* %18
%148 = load i32, i32* %19
%149 = sub nsw i32 %147, %148
%150 = load i32, i32* %10
%151 = add nsw i32 %149, %150
%152 = load i32, i32* %10
%153 = srem i32 %151, %152
store i32 %153, i32* %22
%154 = load float, float* %23
%155 = load float*, float** %13
%156 = load i32, i32* %19
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds float, float* %155, i64 %157
%159 = load float, float* %158
%160 = load i32*, i32** %9
%161 = load i32, i32* %11
%162 = load i32, i32* %21
%163 = mul nsw i32 %161, %162
%164 = load i32, i32* %17
%165 = add nsw i32 %163, %164
%166 = load i32, i32* %10
%167 = load i32, i32* %11
%168 = mul nsw i32 %166, %167
%169 = load i32, i32* %20
%170 = mul nsw i32 %168, %169
%171 = add nsw i32 %165, %170
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds i32, i32* %160, i64 %172
%174 = load i32, i32* %173
%175 = sitofp i32 %174 to float
%176 = call float @llvm.fmuladd.f32(float %159, float %175, float %154)
%177 = load float*, float** %13
%178 = load i32, i32* %19
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds float, float* %177, i64 %179
%181 = load float, float* %180
%182 = load i32*, i32** %9
%183 = load i32, i32* %11
%184 = load i32, i32* %22
%185 = mul nsw i32 %183, %184
%186 = load i32, i32* %17
%187 = add nsw i32 %185, %186
%188 = load i32, i32* %10
%189 = load i32, i32* %11
%190 = mul nsw i32 %188, %189
%191 = load i32, i32* %20
%192 = mul nsw i32 %190, %191
%193 = add nsw i32 %187, %192
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds i32, i32* %182, i64 %194
%196 = load i32, i32* %195
%197 = sitofp i32 %196 to float
%198 = call float @llvm.fmuladd.f32(float %181, float %197, float %176)
store float %198, float* %23
br label %199
199: ; preds = %84
%200 = load i32, i32* %19
%201 = add nsw i32 %200, 1
store i32 %201, i32* %19
br label %80
202: ; preds = %80
%203 = load float, float* %23
%204 = load float*, float** %15
%205 = load i32, i32* %11
%206 = load i32, i32* %18
%207 = mul nsw i32 %205, %206
%208 = load i32, i32* %17
%209 = add nsw i32 %207, %208
%210 = load i32, i32* %10
%211 = load i32, i32* %11
%212 = mul nsw i32 %210, %211
%213 = load i32, i32* %20
%214 = mul nsw i32 %212, %213
%215 = add nsw i32 %209, %214
%216 = sext i32 %215 to i64
%217 = getelementptr inbounds float, float* %204, i64 %216
store float %203, float* %217
%218 = load float, float* %24
%219 = load float*, float** %16
%220 = load i32, i32* %11
%221 = load i32, i32* %18
%222 = mul nsw i32 %220, %221
%223 = load i32, i32* %17
%224 = add nsw i32 %222, %223
%225 = load i32, i32* %10
%226 = load i32, i32* %11
%227 = mul nsw i32 %225, %226
%228 = load i32, i32* %20
%229 = mul nsw i32 %227, %228
%230 = add nsw i32 %224, %229
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds float, float* %219, i64 %231
store float %218, float* %232
br label %233
233: ; preds = %202
%234 = load i32, i32* %18
%235 = add nsw i32 %234, 1
store i32 %235, i32* %18
br label %35
236: ; preds = %35
br label %237
237: ; preds = %236
%238 = load i32, i32* %17
%239 = add nsw i32 %238, 1
store i32 %239, i32* %17
br label %30
240: ; preds = %30
br label %241
241: ; preds = %240
%242 = load i32, i32* %20
%243 = add nsw i32 %242, 1
store i32 %243, i32* %20
br label %25
244: ; preds = %25
ret void
}

LLVM IR (-O2):
define void @separable_convolution(i32* nocapture %0, i32 %1, i32 %2, i32 %3, float* nocapture %4, i32 %5, float* nocapture writeonly %6, float* nocapture writeonly %7) {
%9 = icmp sgt i32 %3, 0
br i1 %9, label %10, label %110
10: ; preds = %8
%11 = icmp sgt i32 %2, 0
%12 = icmp sgt i32 %1, 0
%13 = mul nsw i32 %2, %1
%14 = icmp sgt i32 %5, 1
%15 = zext i32 %2 to i64
%16 = zext i32 %1 to i64
%17 = zext i32 %5 to i64
br label %18
18: ; preds = %10, %107
%19 = phi i32 [ 0, %10 ], [ %108, %107 ]
br i1 %11, label %20, label %107
20: ; preds = %18
%21 = mul nsw i32 %13, %19
br label %22
22: ; preds = %20, %104
%23 = phi i64 [ 0, %20 ], [ %105, %104 ]
br i1 %12, label %24, label %104
24: ; preds = %22
%25 = trunc i64 %23 to i32
%26 = add i32 %21, %25
%27 = add i32 %25, %2
br label %28
28: ; preds = %24, %96
%29 = phi i64 [ 0, %24 ], [ %102, %96 ]
%30 = load float, float* %4
%31 = trunc i64 %29 to i32
%32 = mul i32 %31, %2
%33 = add i32 %26, %32
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds i32, i32* %0, i64 %34
%36 = load i32, i32* %35
%37 = sitofp i32 %36 to float
%38 = fmul float %30, %37
%39 = insertelement <2 x float> poison, float %38, i64 0
%40 = shufflevector <2 x float> %39, <2 x float> poison, <2 x i32> zeroinitializer
br i1 %14, label %41, label %96
41: ; preds = %28
%42 = add i32 %32, %21
%43 = trunc i64 %29 to i32
%44 = add i32 %43, %1
br label %45
45: ; preds = %41, %45
%46 = phi i64 [ 1, %41 ], [ %90, %45 ]
%47 = phi i32 [ -1, %41 ], [ %91, %45 ]
%48 = phi <2 x float> [ %40, %41 ], [ %89, %45 ]
%49 = trunc i64 %46 to i32
%50 = add nuw nsw i64 %46, %23
%51 = trunc i64 %50 to i32
%52 = srem i32 %51, %2
%53 = add i32 %27, %47
%54 = srem i32 %53, %2
%55 = getelementptr inbounds float, float* %4, i64 %46
%56 = load float, float* %55
%57 = add i32 %42, %52
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds i32, i32* %0, i64 %58
%60 = load i32, i32* %59
%61 = add i32 %42, %54
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds i32, i32* %0, i64 %62
%64 = load i32, i32* %63
%65 = add nuw nsw i64 %46, %29
%66 = trunc i64 %65 to i32
%67 = srem i32 %66, %1
%68 = add i32 %44, %47
%69 = srem i32 %68, %1
%70 = mul nsw i32 %67, %2
%71 = add i32 %26, %70
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds i32, i32* %0, i64 %72
%74 = load i32, i32* %73
%75 = insertelement <2 x i32> poison, i32 %74, i64 0
%76 = insertelement <2 x i32> %75, i32 %60, i64 1
%77 = sitofp <2 x i32> %76 to <2 x float>
%78 = insertelement <2 x float> poison, float %56, i64 0
%79 = shufflevector <2 x float> %78, <2 x float> poison, <2 x i32> zeroinitializer
%80 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %79, <2 x float> %77, <2 x float> %48)
%81 = mul nsw i32 %69, %2
%82 = add i32 %26, %81
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds i32, i32* %0, i64 %83
%85 = load i32, i32* %84
%86 = insertelement <2 x i32> poison, i32 %85, i64 0
%87 = insertelement <2 x i32> %86, i32 %64, i64 1
%88 = sitofp <2 x i32> %87 to <2 x float>
%89 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %79, <2 x float> %88, <2 x float> %80)
%90 = add nuw nsw i64 %46, 1
%91 = xor i32 %49, -1
%92 = icmp eq i64 %90, %17
br i1 %92, label %93, label %45
93: ; preds = %45
%94 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %79, <2 x float> %77, <2 x float> %48)
%95 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %79, <2 x float> %88, <2 x float> %94)
br label %96
96: ; preds = %93, %28
%97 = phi <2 x float> [ %40, %28 ], [ %95, %93 ]
%98 = getelementptr inbounds float, float* %6, i64 %34
%99 = extractelement <2 x float> %97, i64 0
store float %99, float* %98
%100 = getelementptr inbounds float, float* %7, i64 %34
%101 = extractelement <2 x float> %97, i64 1
store float %101, float* %100
%102 = add nuw nsw i64 %29, 1
%103 = icmp eq i64 %102, %16
br i1 %103, label %104, label %28
104: ; preds = %96, %22
%105 = add nuw nsw i64 %23, 1
%106 = icmp eq i64 %105, %15
br i1 %106, label %107, label %22
107: ; preds = %104, %18
%108 = add nuw nsw i32 %19, 1
%109 = icmp eq i32 %108, %3
br i1 %109, label %110, label %18
110: ; preds = %107, %8
ret void
}