void overlapped_adaptive_threshold(const unsigned char *source,
                                   int source_stride, unsigned char *dest,
                                   int width, int height, int tile_size) {
  int x_step = tile_size;
  int y_step = tile_size;
  int x, y;
  for (y = 0; y < height; y += y_step / 2) {
    for (x = 0; x < width; x += x_step / 2) {
      int xxmax, yymax, sum, count, yy, xx;
      int mean;
      double fmean, sd, fsd;
      int ys, ymax, xs, xmax;
      xxmax = x + x_step;
      if (xxmax > width)
        xxmax = width;
      yymax = y + y_step;
      if (yymax > height)
        yymax = height;
      sum = 0;
      count = (yymax - y) * (xxmax - x);
      for (yy = y; yy < yymax; ++yy) {
        const unsigned char *p = source + (yy * width + x) * source_stride;
        for (xx = x; xx < xxmax; ++xx) {
          sum += *p;
          p += source_stride;
        }
      }
      mean = sum / count;
      fmean = (double)sum / (double)count;
      sd = 0;
      for (yy = y; yy < yymax; ++yy) {
        const unsigned char *p = source + (yy * width + x) * source_stride;
        for (xx = x; xx < xxmax; ++xx) {
          double z = *p - fmean;
          sd += z * z;
          p += source_stride;
        }
      }
      fsd = (double)sd / (double)count;
      if (fsd < 10 * 10)
        mean -= 20;
      ys = y + y_step / 4;
      ymax = (y + y_step) - y_step / 4;
      if (ymax > yymax)
        ymax = yymax;
      if (y == 0)
        ys = 0;
      xs = x + x_step / 4;
      xmax = (x + x_step) - x_step / 4;
      if (xmax > xxmax)
        xmax = xxmax;
      if (x == 0)
        xs = 0;
      for (yy = ys; yy < ymax; ++yy) {
        const unsigned char *p = source + (yy * width + xs) * source_stride;
        for (xx = xs; xx < xmax; ++xx) {
          if (*p > mean) {
            dest[(yy * width + xx)] = 255;
          } else {
            dest[(yy * width + xx)] = 0;
          }
          p += source_stride;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<overlapped_adaptive_threshold>:
push   %rbp
mov    %rsp,%rbp
sub    $0x28,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    14ed <overlapped_adaptive_threshold+0x3ed>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    14d2 <overlapped_adaptive_threshold+0x3d2>
mov    -0x30(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jle    116b <overlapped_adaptive_threshold+0x6b>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jle    1186 <overlapped_adaptive_threshold+0x86>
mov    -0x20(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    -0x38(%rbp),%ecx
sub    -0x30(%rbp),%ecx
imul   %ecx,%eax
mov    %eax,-0x44(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    121e <overlapped_adaptive_threshold+0x11e>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x80(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    120b <overlapped_adaptive_threshold+0x10b>
mov    -0x80(%rbp),%rax
movzbl (%rax),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x80(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x80(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11d3 <overlapped_adaptive_threshold+0xd3>
jmp    1210 <overlapped_adaptive_threshold+0x110>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    11a5 <overlapped_adaptive_threshold+0xa5>
mov    -0x40(%rbp),%eax
cltd
idivl  -0x44(%rbp)
mov    %eax,-0x50(%rbp)
cvtsi2sdl -0x40(%rbp),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    12fb <overlapped_adaptive_threshold+0x1fb>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x88(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    12e8 <overlapped_adaptive_threshold+0x1e8>
mov    -0x88(%rbp),%rax
movzbl (%rax),%eax
cvtsi2sd %eax,%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x88(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x88(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    127a <overlapped_adaptive_threshold+0x17a>
jmp    12ed <overlapped_adaptive_threshold+0x1ed>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1249 <overlapped_adaptive_threshold+0x149>
movsd  -0x60(%rbp),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0xcea(%rip),%xmm0
ucomisd -0x68(%rbp),%xmm0
jbe    132a <overlapped_adaptive_threshold+0x22a>
mov    -0x50(%rbp),%eax
sub    $0x14,%eax
mov    %eax,-0x50(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0xa0(%rbp)
mov    -0x2c(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xa0(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x6c(%rbp)
mov    -0x34(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x9c(%rbp)
mov    -0x2c(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x9c(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jle    1381 <overlapped_adaptive_threshold+0x281>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x70(%rbp)
cmpl   $0x0,-0x34(%rbp)
jne    1392 <overlapped_adaptive_threshold+0x292>
movl   $0x0,-0x6c(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0xa8(%rbp)
mov    -0x28(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xa8(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x74(%rbp)
mov    -0x30(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0xa4(%rbp)
mov    -0x28(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xa4(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x78(%rbp)
mov    -0x78(%rbp),%eax
cmp    -0x38(%rbp),%eax
jle    13e9 <overlapped_adaptive_threshold+0x2e9>
mov    -0x38(%rbp),%eax
mov    %eax,-0x78(%rbp)
cmpl   $0x0,-0x30(%rbp)
jne    13fa <overlapped_adaptive_threshold+0x2fa>
movl   $0x0,-0x74(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x70(%rbp),%eax
jge    14b7 <overlapped_adaptive_threshold+0x3b7>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x74(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x98(%rbp)
mov    -0x74(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x78(%rbp),%eax
jge    14a4 <overlapped_adaptive_threshold+0x3a4>
mov    -0x98(%rbp),%rax
movzbl (%rax),%eax
cmp    -0x50(%rbp),%eax
jle    146a <overlapped_adaptive_threshold+0x36a>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movb   $0xff,(%rax,%rcx,1)
jmp    147f <overlapped_adaptive_threshold+0x37f>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movb   $0x0,(%rax,%rcx,1)
mov    -0xc(%rbp),%ecx
mov    -0x98(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x98(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1431 <overlapped_adaptive_threshold+0x331>
jmp    14a9 <overlapped_adaptive_threshold+0x3a9>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1400 <overlapped_adaptive_threshold+0x300>
jmp    14bc <overlapped_adaptive_threshold+0x3bc>
mov    -0x28(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    1144 <overlapped_adaptive_threshold+0x44>
jmp    14d7 <overlapped_adaptive_threshold+0x3d7>
mov    -0x2c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
jmp    1131 <overlapped_adaptive_threshold+0x31>
add    $0x28,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<overlapped_adaptive_threshold>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9,-0x58(%rsp)
mov    %ecx,-0x7c(%rsp)
mov    %rdx,-0x40(%rsp)
mov    %rdi,-0x70(%rsp)
mov    %r8d,-0x78(%rsp)
test   %r8d,%r8d
jle    14a6 <overlapped_adaptive_threshold+0x3a6>
movslq %esi,%r11
mov    -0x58(%rsp),%rcx
lea    0x3(%rcx),%edx
test   %ecx,%ecx
cmovns %ecx,%edx
sar    $0x2,%edx
mov    %ecx,%eax
shr    $0x1f,%eax
add    %ecx,%eax
sar    %eax
cltq
mov    %rax,-0x48(%rsp)
movslq -0x7c(%rsp),%r12
movslq %ecx,%rax
mov    %rax,-0x10(%rsp)
mov    %rdx,-0x50(%rsp)
movslq %edx,%rax
mov    %rax,-0x18(%rsp)
mov    -0x40(%rsp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rsp)
movsd  0xe7e(%rip),%xmm0
xor    %r14d,%r14d
jmp    11a0 <overlapped_adaptive_threshold+0xa0>
nopw   0x0(%rax,%rax,1)
add    -0x48(%rsp),%r14
cmp    -0x78(%rsp),%r14d
jge    14a6 <overlapped_adaptive_threshold+0x3a6>
cmpl   $0x0,-0x7c(%rsp)
jle    1190 <overlapped_adaptive_threshold+0x90>
mov    -0x58(%rsp),%rax
lea    (%r14,%rax,1),%ecx
mov    -0x78(%rsp),%eax
cmp    %eax,%ecx
mov    %ecx,%edx
cmovg  %eax,%edx
mov    %edx,%eax
sub    %r14d,%eax
mov    %eax,-0x60(%rsp)
mov    -0x50(%rsp),%rax
lea    (%rax,%r14,1),%esi
sub    %eax,%ecx
cmp    %edx,%ecx
cmovg  %edx,%ecx
test   %r14d,%r14d
cmove  %r14d,%esi
mov    %edx,-0x74(%rsp)
movslq %edx,%r10
mov    %esi,-0x64(%rsp)
movslq %esi,%rax
mov    %ecx,-0x5c(%rsp)
movslq %ecx,%rcx
mov    %rcx,0x20(%rsp)
mov    %r12,%rcx
mov    %rax,-0x20(%rsp)
imul   %rax,%rcx
add    -0x30(%rsp),%rcx
mov    %rcx,-0x28(%rsp)
xor    %ebx,%ebx
mov    %r14,(%rsp)
mov    %r10,-0x8(%rsp)
jmp    123c <overlapped_adaptive_threshold+0x13c>
nopl   0x0(%rax,%rax,1)
mov    0x8(%rsp),%rbx
add    -0x48(%rsp),%rbx
cmp    %r12,%rbx
mov    (%rsp),%r14
mov    -0x8(%rsp),%r10
jge    1190 <overlapped_adaptive_threshold+0x90>
mov    -0x10(%rsp),%rax
add    %rbx,%rax
cmp    %r12,%rax
mov    %rax,%r13
mov    %eax,%ecx
cmovg  -0x7c(%rsp),%ecx
mov    %ecx,%r9d
sub    %ebx,%r9d
imul   -0x60(%rsp),%r9d
cmp    %r14d,-0x74(%rsp)
jle    12c0 <overlapped_adaptive_threshold+0x1c0>
movslq %ecx,%r8
xor    %esi,%esi
mov    %r14,%rax
jmp    1279 <overlapped_adaptive_threshold+0x179>
xchg   %ax,%ax
add    $0x1,%rax
cmp    %r10,%rax
jge    12d0 <overlapped_adaptive_threshold+0x1d0>
cmp    %r8,%rbx
jge    1270 <overlapped_adaptive_threshold+0x170>
mov    %rax,%rdx
imul   %r12,%rdx
add    %rbx,%rdx
imul   %r11,%rdx
add    -0x70(%rsp),%rdx
mov    %ebx,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movzbl (%rdx),%ebp
add    %ebp,%esi
add    %r11,%rdx
add    $0x1,%edi
cmp    %ecx,%edi
jl     12a0 <overlapped_adaptive_threshold+0x1a0>
jmp    1270 <overlapped_adaptive_threshold+0x170>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sd %r9d,%xmm1
xor    %eax,%eax
jmp    1360 <overlapped_adaptive_threshold+0x260>
nop
mov    %esi,%eax
cltd
idiv   %r9d
xorps  %xmm1,%xmm1
cvtsi2sd %r9d,%xmm1
cmp    %r14d,-0x74(%rsp)
jle    1360 <overlapped_adaptive_threshold+0x260>
xorps  %xmm3,%xmm3
cvtsi2sd %esi,%xmm3
divsd  %xmm1,%xmm3
xorpd  %xmm2,%xmm2
mov    %r14,%rdx
jmp    1309 <overlapped_adaptive_threshold+0x209>
nopl   0x0(%rax)
add    $0x1,%rdx
cmp    %r10,%rdx
jge    1367 <overlapped_adaptive_threshold+0x267>
cmp    %r8,%rbx
jge    1300 <overlapped_adaptive_threshold+0x200>
mov    %rdx,%rsi
imul   %r12,%rsi
add    %rbx,%rsi
imul   %r11,%rsi
add    -0x70(%rsp),%rsi
mov    %ebx,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movzbl (%rsi),%ebp
xorps  %xmm4,%xmm4
cvtsi2sd %ebp,%xmm4
subsd  %xmm3,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm4,%xmm2
add    %r11,%rsi
add    $0x1,%edi
cmp    %ecx,%edi
jl     1330 <overlapped_adaptive_threshold+0x230>
jmp    1300 <overlapped_adaptive_threshold+0x200>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %ecx,%r8
xorpd  %xmm2,%xmm2
divsd  %xmm1,%xmm2
lea    -0x14(%rax),%edx
ucomisd %xmm2,%xmm0
cmovbe %eax,%edx
sub    -0x18(%rsp),%r13
cmp    %r8,%r13
cmovg  %ecx,%r13d
mov    %r13,-0x38(%rsp)
mov    -0x50(%rsp),%rax
lea    (%rax,%rbx,1),%r15d
mov    %rbx,0x8(%rsp)
test   %rbx,%rbx
mov    $0x0,%ecx
cmove  %ecx,%r15d
mov    -0x64(%rsp),%ecx
cmp    -0x5c(%rsp),%ecx
jge    1220 <overlapped_adaptive_threshold+0x120>
movslq %r15d,%rbp
movslq -0x38(%rsp),%rsi
mov    %rsi,%r10
sub    %rbp,%r10
mov    %rbp,%rax
not    %rax
lea    0x1(%rbp),%rcx
mov    %rcx,0x18(%rsp)
mov    %rsi,%r8
neg    %r8
mov    -0x28(%rsp),%r9
mov    -0x20(%rsp),%r14
mov    %r8,0x10(%rsp)
jmp    1402 <overlapped_adaptive_threshold+0x302>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r14
add    %r12,%r9
cmp    0x20(%rsp),%r14
jge    1220 <overlapped_adaptive_threshold+0x120>
cmp    -0x38(%rsp),%r15d
jge    13f0 <overlapped_adaptive_threshold+0x2f0>
mov    %r14,%rcx
imul   %r12,%rcx
add    %rbp,%rcx
mov    %rcx,%r13
imul   %r11,%r13
add    -0x70(%rsp),%r13
test   $0x1,%r10b
jne    1430 <overlapped_adaptive_threshold+0x330>
mov    %rbp,%rcx
cmp    %r8,%rax
jne    1470 <overlapped_adaptive_threshold+0x370>
jmp    13f0 <overlapped_adaptive_threshold+0x2f0>
nop
mov    %rbp,%r8
mov    %r12,%rbp
movzbl 0x0(%r13),%r12d
cmp    %r12d,%edx
setl   %bl
neg    %bl
mov    -0x40(%rsp),%rdi
mov    %bl,(%rdi,%rcx,1)
mov    %rbp,%r12
mov    %r8,%rbp
mov    0x10(%rsp),%r8
add    %r11,%r13
mov    0x18(%rsp),%rcx
cmp    %r8,%rax
je     13f0 <overlapped_adaptive_threshold+0x2f0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movzbl 0x0(%r13),%ebx
cmp    %ebx,%edx
setl   %bl
neg    %bl
mov    %bl,-0x1(%r9,%rcx,1)
movzbl 0x0(%r13,%r11,1),%ebx
add    %r11,%r13
cmp    %ebx,%edx
setl   %bl
neg    %bl
mov    %bl,(%r9,%rcx,1)
add    %r11,%r13
add    $0x2,%rcx
cmp    %rsi,%rcx
jl     1470 <overlapped_adaptive_threshold+0x370>
jmp    13f0 <overlapped_adaptive_threshold+0x2f0>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @overlapped_adaptive_threshold(i8* nocapture %0, i32 %1, i8* nocapture writeonly %2, i32 %3, i32 %4, i32 %5) {
%7 = icmp sgt i32 %4, 0
br i1 %7, label %8, label %183
8: ; preds = %6
%9 = icmp sgt i32 %3, 0
%10 = sext i32 %1 to i64
%11 = sdiv i32 %5, 4
%12 = sdiv i32 %5, 2
%13 = sext i32 %12 to i64
%14 = sext i32 %3 to i64
%15 = sext i32 %5 to i64
%16 = sext i32 %11 to i64
br label %17
17: ; preds = %8, %179
%18 = phi i64 [ 0, %8 ], [ %180, %179 ]
%19 = trunc i64 %18 to i32
br i1 %9, label %20, label %179
20: ; preds = %17
%21 = add nsw i32 %19, %5
%22 = icmp sgt i32 %21, %4
%23 = select i1 %22, i32 %4, i32 %21
%24 = sub nsw i32 %23, %19
%25 = icmp sgt i32 %23, %19
%26 = add nsw i32 %11, %19
%27 = sub nsw i32 %21, %11
%28 = icmp sgt i32 %27, %23
%29 = select i1 %28, i32 %23, i32 %27
%30 = icmp eq i32 %19, 0
%31 = select i1 %30, i32 0, i32 %26
%32 = icmp slt i32 %31, %29
%33 = sext i32 %23 to i64
%34 = sext i32 %31 to i64
%35 = sext i32 %29 to i64
br label %36
36: ; preds = %20, %176
%37 = phi i64 [ 0, %20 ], [ %177, %176 ]
%38 = add nsw i64 %37, %15
%39 = icmp sgt i64 %38, %14
%40 = trunc i64 %38 to i32
%41 = select i1 %39, i32 %3, i32 %40
%42 = trunc i64 %37 to i32
%43 = sub nsw i32 %41, %42
%44 = mul nsw i32 %43, %24
br i1 %25, label %47, label %45
45: ; preds = %36
%46 = sitofp i32 %44 to double
br label %77
47: ; preds = %36
%48 = sext i32 %41 to i64
%49 = icmp slt i64 %37, %48
br label %50
50: ; preds = %47, %68
%51 = phi i64 [ %18, %47 ], [ %70, %68 ]
%52 = phi i32 [ 0, %47 ], [ %69, %68 ]
br i1 %49, label %53, label %68
53: ; preds = %50
%54 = mul nsw i64 %51, %14
%55 = add nsw i64 %54, %37
%56 = mul nsw i64 %55, %10
%57 = getelementptr inbounds i8, i8* %0, i64 %56
br label %58
58: ; preds = %53, %58
%59 = phi i8* [ %65, %58 ], [ %57, %53 ]
%60 = phi i32 [ %66, %58 ], [ %42, %53 ]
%61 = phi i32 [ %64, %58 ], [ %52, %53 ]
%62 = load i8, i8* %59
%63 = zext i8 %62 to i32
%64 = add nsw i32 %61, %63
%65 = getelementptr inbounds i8, i8* %59, i64 %10
%66 = add nsw i32 %60, 1
%67 = icmp slt i32 %66, %41
br i1 %67, label %58, label %68
68: ; preds = %58, %50
%69 = phi i32 [ %52, %50 ], [ %64, %58 ]
%70 = add nsw i64 %51, 1
%71 = icmp slt i64 %70, %33
br i1 %71, label %50, label %72
72: ; preds = %68
%73 = sdiv i32 %69, %44
%74 = sitofp i32 %69 to double
%75 = sitofp i32 %44 to double
%76 = fdiv double %74, %75
br i1 %25, label %81, label %77
77: ; preds = %45, %72
%78 = phi double [ %46, %45 ], [ %75, %72 ]
%79 = phi i32 [ 0, %45 ], [ %73, %72 ]
%80 = sext i32 %41 to i64
br label %107
81: ; preds = %72
%82 = sext i32 %41 to i64
%83 = icmp slt i64 %37, %82
br label %84
84: ; preds = %81, %103
%85 = phi i64 [ %18, %81 ], [ %105, %103 ]
%86 = phi double [ 0.000000e+00, %81 ], [ %104, %103 ]
br i1 %83, label %87, label %103
87: ; preds = %84
%88 = mul nsw i64 %85, %14
%89 = add nsw i64 %88, %37
%90 = mul nsw i64 %89, %10
%91 = getelementptr inbounds i8, i8* %0, i64 %90
br label %92
92: ; preds = %87, %92
%93 = phi i8* [ %100, %92 ], [ %91, %87 ]
%94 = phi double [ %99, %92 ], [ %86, %87 ]
%95 = phi i32 [ %101, %92 ], [ %42, %87 ]
%96 = load i8, i8* %93
%97 = uitofp i8 %96 to double
%98 = fsub double %97, %76
%99 = tail call double @llvm.fmuladd.f64(double %98, double %98, double %94)
%100 = getelementptr inbounds i8, i8* %93, i64 %10
%101 = add nsw i32 %95, 1
%102 = icmp slt i32 %101, %41
br i1 %102, label %92, label %103
103: ; preds = %92, %84
%104 = phi double [ %86, %84 ], [ %99, %92 ]
%105 = add nsw i64 %85, 1
%106 = icmp slt i64 %105, %33
br i1 %106, label %84, label %107
107: ; preds = %103, %77
%108 = phi double [ %78, %77 ], [ %75, %103 ]
%109 = phi i32 [ %79, %77 ], [ %73, %103 ]
%110 = phi i64 [ %80, %77 ], [ %82, %103 ]
%111 = phi double [ 0.000000e+00, %77 ], [ %104, %103 ]
%112 = fdiv double %111, %108
%113 = fcmp olt double %112, 1.000000e+02
%114 = add nsw i32 %109, -20
%115 = select i1 %113, i32 %114, i32 %109
%116 = sub nsw i64 %38, %16
%117 = icmp sgt i64 %116, %110
%118 = trunc i64 %116 to i32
%119 = select i1 %117, i32 %41, i32 %118
%120 = icmp eq i64 %37, 0
%121 = trunc i64 %37 to i32
%122 = add i32 %11, %121
%123 = select i1 %120, i32 0, i32 %122
br i1 %32, label %124, label %176
124: ; preds = %107
%125 = icmp slt i32 %123, %119
%126 = sext i32 %123 to i64
%127 = sext i32 %119 to i64
%128 = sub nsw i64 %127, %126
%129 = xor i64 %126, -1
%130 = and i64 %128, 1
%131 = icmp eq i64 %130, 0
%132 = add nsw i64 %126, 1
%133 = sub nsw i64 0, %127
%134 = icmp eq i64 %129, %133
br label %135
135: ; preds = %124, %173
%136 = phi i64 [ %34, %124 ], [ %174, %173 ]
%137 = mul nsw i64 %136, %14
br i1 %125, label %138, label %173
138: ; preds = %135
%139 = add nsw i64 %137, %126
%140 = mul nsw i64 %139, %10
%141 = getelementptr inbounds i8, i8* %0, i64 %140
br i1 %131, label %150, label %142
142: ; preds = %138
%143 = load i8, i8* %141
%144 = zext i8 %143 to i32
%145 = icmp slt i32 %115, %144
%146 = add nsw i64 %137, %126
%147 = getelementptr inbounds i8, i8* %2, i64 %146
%148 = sext i1 %145 to i8
store i8 %148, i8* %147
%149 = getelementptr inbounds i8, i8* %141, i64 %10
br label %150
150: ; preds = %142, %138
%151 = phi i64 [ %126, %138 ], [ %132, %142 ]
%152 = phi i8* [ %141, %138 ], [ %149, %142 ]
br i1 %134, label %173, label %153
153: ; preds = %150, %153
%154 = phi i64 [ %171, %153 ], [ %151, %150 ]
%155 = phi i8* [ %170, %153 ], [ %152, %150 ]
%156 = load i8, i8* %155
%157 = zext i8 %156 to i32
%158 = icmp slt i32 %115, %157
%159 = add nsw i64 %154, %137
%160 = getelementptr inbounds i8, i8* %2, i64 %159
%161 = sext i1 %158 to i8
store i8 %161, i8* %160
%162 = getelementptr inbounds i8, i8* %155, i64 %10
%163 = add nsw i64 %154, 1
%164 = load i8, i8* %162
%165 = zext i8 %164 to i32
%166 = icmp slt i32 %115, %165
%167 = add nsw i64 %163, %137
%168 = getelementptr inbounds i8, i8* %2, i64 %167
%169 = sext i1 %166 to i8
store i8 %169, i8* %168
%170 = getelementptr inbounds i8, i8* %162, i64 %10
%171 = add nsw i64 %154, 2
%172 = icmp slt i64 %171, %127
br i1 %172, label %153, label %173
173: ; preds = %150, %153, %135
%174 = add nsw i64 %136, 1
%175 = icmp slt i64 %174, %35
br i1 %175, label %135, label %176
176: ; preds = %173, %107
%177 = add nsw i64 %37, %13
%178 = icmp slt i64 %177, %14
br i1 %178, label %36, label %179
179: ; preds = %176, %17
%180 = add nsw i64 %18, %13
%181 = trunc i64 %180 to i32
%182 = icmp slt i32 %181, %4
br i1 %182, label %17, label %183
183: ; preds = %179, %6
ret void
}