void data_disconts(image1, image2, u, v, s1, nx, ny, warp, It,
                   discont) float *image1,
    *image2, *It, *u, *v, *warp, *discont;
int nx, ny;
float *s1;
{
  int index, nx_max, ny_max, i, j;
  float t1;
  warp_image(image1, warp, u, v, nx, ny);
  horn_dt(warp, image2, It, nx, ny);
  for (i = 1; i < ny - 1; i++) {
    for (j = 1; j < nx - 1; j++) {
      index = i * nx + j;
      t1 = s1[index] * 1.4142136;
      if (((It[index]) < 0 ? -(It[index]) : (It[index])) >= t1)
        discont[index] = 0.0;
      else
        discont[index] = 1.0;
    }
  }
  ny_max = ny - 1;
  nx_max = nx - 1;
  for (i = 0; i < ny; i++) {
    discont[i * nx] = 0.0;
    discont[i * nx + nx_max] = 0.0;
  }
  for (j = 0; j < nx; j++) {
    discont[j] = 0.0;
    discont[ny_max * nx + j] = 0.0;
  }
}


x86 assembly (clang, -O0):
<data_disconts>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x8(%rbp),%rdi
mov    0x18(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    -0x2c(%rbp),%r8d
mov    0x10(%rbp),%r9d
mov    $0x0,%al
call   1030 <warp_image@plt>
mov    0x18(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    0x20(%rbp),%rdx
mov    -0x2c(%rbp),%ecx
mov    0x10(%rbp),%r8d
mov    $0x0,%al
call   1040 <horn_dt@plt>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1295 <data_disconts+0x175>
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1282 <data_disconts+0x162>
mov    -0x3c(%rbp),%eax
imul   -0x2c(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe22(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomiss (%rax,%rcx,4),%xmm0
jbe    1224 <data_disconts+0x104>
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x48(%rbp)
jmp    1236 <data_disconts+0x116>
mov    0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
ucomiss -0x44(%rbp),%xmm0
jb     125a <data_disconts+0x13a>
mov    0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    126f <data_disconts+0x14f>
mov    0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  0xd9e(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1274 <data_disconts+0x154>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11a7 <data_disconts+0x87>
jmp    1287 <data_disconts+0x167>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    118f <data_disconts+0x6f>
mov    0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    12f7 <data_disconts+0x1d7>
mov    0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12ae <data_disconts+0x18e>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1341 <data_disconts+0x221>
mov    0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12fe <data_disconts+0x1de>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<data_disconts>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9d,%ebx
mov    %r8,0x30(%rsp)
mov    %rsi,%r13
mov    0xc0(%rsp),%r15
mov    0xb0(%rsp),%ebp
mov    0xb8(%rsp),%r14
mov    %r14,%rsi
mov    %r9d,%r8d
mov    %ebp,%r9d
xor    %eax,%eax
call   1030 <warp_image@plt>
mov    %r14,%rdi
mov    %r13,%rsi
mov    %r15,%rdx
mov    %ebx,%ecx
mov    %ebp,%r8d
xor    %eax,%eax
call   1040 <horn_dt@plt>
mov    %rbx,%rcx
lea    -0x1(%rbp),%eax
mov    %eax,0x4(%rsp)
lea    -0x1(%rbx),%eax
mov    %eax,(%rsp)
movslq %ecx,%rax
mov    %rax,0x10(%rsp)
cmp    $0x3,%ebp
jl     140c <data_disconts+0x2ec>
mov    0x4(%rsp),%ebp
mov    (%rsp),%edi
mov    0x10(%rsp),%rsi
lea    0x1(%rsi),%rax
mov    %rax,0x58(%rsp)
lea    (%rsi,%rdi,1),%rax
mov    %rax,0x50(%rsp)
lea    -0x1(%rdi),%rax
mov    %rax,0x38(%rsp)
and    $0xfffffffffffffffc,%rax
mov    %rax,0x48(%rsp)
add    $0x1,%rax
mov    %rax,0x40(%rsp)
mov    %rcx,0x8(%rsp)
mov    0xc8(%rsp),%rcx
lea    0x4(%rcx,%rsi,4),%rdx
lea    0x0(,%rsi,4),%r8
mov    0xc0(%rsp),%rax
lea    0x4(%rax,%rsi,4),%r9
mov    0x30(%rsp),%rbx
lea    0x4(%rbx,%rsi,4),%r11
lea    (%rcx,%rsi,4),%r10
mov    0x8(%rsp),%rcx
lea    (%rbx,%rsi,4),%r13
lea    (%rax,%rsi,4),%r12
mov    %rbp,%rax
mov    $0x1,%r15d
movsd  0xe0e(%rip),%xmm0
movaps 0xde7(%rip),%xmm1
xorps  %xmm2,%xmm2
movss  0xe04(%rip),%xmm3
movapd 0xdc3(%rip),%xmm10
xorps  %xmm8,%xmm8
movaps 0xdd7(%rip),%xmm9
xor    %ebp,%ebp
mov    %rax,0x28(%rsp)
jmp    1283 <data_disconts+0x163>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r15
add    $0x1,%rbp
add    %r8,%rdx
add    %r8,%r9
add    %r8,%r11
add    %r8,%r10
add    %r8,%r13
add    %r8,%r12
cmp    %rax,%r15
je     140c <data_disconts+0x2ec>
cmp    $0x3,%ecx
jl     1260 <data_disconts+0x140>
mov    $0x1,%r14d
cmpq   $0x4,0x38(%rsp)
jb     13c0 <data_disconts+0x2a0>
mov    %rbp,%rcx
imul   0x10(%rsp),%rcx
mov    0x58(%rsp),%rax
lea    (%rax,%rcx,1),%rbx
add    0x50(%rsp),%rcx
mov    0xc8(%rsp),%rsi
lea    (%rsi,%rbx,4),%rax
mov    %rax,0x60(%rsp)
lea    (%rsi,%rcx,4),%rax
mov    %rax,0x68(%rsp)
mov    0x30(%rsp),%rax
lea    (%rax,%rbx,4),%rsi
mov    %rsi,0x18(%rsp)
lea    (%rax,%rcx,4),%rax
mov    %rax,0x20(%rsp)
mov    0xc0(%rsp),%rax
lea    (%rax,%rbx,4),%rbx
mov    %rbx,0x70(%rsp)
lea    (%rax,%rcx,4),%rcx
mov    0x60(%rsp),%rsi
cmp    0x20(%rsp),%rsi
setb   0x20(%rsp)
mov    0x68(%rsp),%rax
cmp    %rax,0x18(%rsp)
setb   0x18(%rsp)
cmp    %rcx,%rsi
setb   %bl
cmp    %rax,0x70(%rsp)
setb   %cl
mov    0x18(%rsp),%al
test   %al,0x20(%rsp)
jne    13b3 <data_disconts+0x293>
and    %cl,%bl
mov    0x8(%rsp),%rcx
mov    0x28(%rsp),%rax
jne    13c0 <data_disconts+0x2a0>
xor    %ebx,%ebx
mov    0x48(%rsp),%rsi
nopw   0x0(%rax,%rax,1)
cvtps2pd 0x8(%r11,%rbx,4),%xmm7
cvtps2pd (%r11,%rbx,4),%xmm5
mulpd  %xmm10,%xmm5
mulpd  %xmm10,%xmm7
cvtpd2ps %xmm7,%xmm7
cvtpd2ps %xmm5,%xmm5
unpcklpd %xmm7,%xmm5
movups (%r9,%rbx,4),%xmm7
movaps %xmm7,%xmm6
cmpltps %xmm8,%xmm6
movaps %xmm6,%xmm4
andnps %xmm7,%xmm4
xorps  %xmm1,%xmm7
andps  %xmm6,%xmm7
orps   %xmm4,%xmm7
cmpnleps %xmm7,%xmm5
andps  %xmm9,%xmm5
movups %xmm5,(%rdx,%rbx,4)
add    $0x4,%rbx
cmp    %rbx,%rsi
jne    1350 <data_disconts+0x230>
mov    0x40(%rsp),%r14
cmp    %rsi,0x38(%rsp)
jne    13c0 <data_disconts+0x2a0>
jmp    1260 <data_disconts+0x140>
mov    0x8(%rsp),%rcx
mov    0x28(%rsp),%rax
nopl   (%rax)
movss  0x0(%r13,%r14,4),%xmm4
cvtss2sd %xmm4,%xmm4
mulsd  %xmm0,%xmm4
cvtsd2ss %xmm4,%xmm4
movss  (%r12,%r14,4),%xmm5
movaps %xmm5,%xmm6
xorps  %xmm1,%xmm6
movaps %xmm5,%xmm7
cmpltss %xmm2,%xmm7
andps  %xmm7,%xmm6
andnps %xmm5,%xmm7
orps   %xmm6,%xmm7
cmpnless %xmm7,%xmm4
andps  %xmm3,%xmm4
movss  %xmm4,(%r10,%r14,4)
add    $0x1,%r14
cmp    %r14,%rdi
jne    13c0 <data_disconts+0x2a0>
jmp    1260 <data_disconts+0x140>
mov    %rcx,%r10
mov    0xb0(%rsp),%edx
test   %edx,%edx
mov    0xc8(%rsp),%r9
mov    0x10(%rsp),%rcx
jle    1494 <data_disconts+0x374>
movslq (%rsp),%rax
mov    %edx,%r8d
cmp    $0x1,%edx
jne    1437 <data_disconts+0x317>
xor    %esi,%esi
jmp    1477 <data_disconts+0x357>
mov    %r8d,%edi
and    $0xfffffffe,%edi
lea    (%rcx,%rax,1),%rbp
lea    0x0(,%rcx,8),%rbx
xor    %esi,%esi
mov    %r9,%rdx
xchg   %ax,%ax
movl   $0x0,(%rdx)
movl   $0x0,(%rdx,%rax,4)
movl   $0x0,(%rdx,%rcx,4)
movl   $0x0,(%rdx,%rbp,4)
add    $0x2,%rsi
add    %rbx,%rdx
cmp    %rsi,%rdi
jne    1450 <data_disconts+0x330>
test   $0x1,%r8b
je     1494 <data_disconts+0x374>
imul   %rcx,%rsi
movl   $0x0,(%r9,%rsi,4)
add    %rax,%rsi
movl   $0x0,(%r9,%rsi,4)
test   %r10d,%r10d
jle    1668 <data_disconts+0x548>
mov    0x4(%rsp),%eax
imul   %r10d,%eax
movslq %eax,%r8
mov    %r10d,%eax
cmp    $0x8,%r10d
jb     14d3 <data_disconts+0x3b3>
lea    (%r8,%rax,1),%rcx
lea    (%r9,%rcx,4),%rcx
cmp    %r9,%rcx
jbe    1570 <data_disconts+0x450>
lea    (%r9,%rax,4),%rcx
lea    (%r9,%r8,4),%rdx
cmp    %rcx,%rdx
jae    1570 <data_disconts+0x450>
xor    %edx,%edx
mov    %rdx,%rcx
not    %rcx
add    %rax,%rcx
mov    %rax,%rsi
and    $0x3,%rsi
je     1509 <data_disconts+0x3e9>
lea    (%r9,%r8,4),%rdi
nopl   0x0(%rax,%rax,1)
movl   $0x0,(%r9,%rdx,4)
movl   $0x0,(%rdi,%rdx,4)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rsi
jne    14f0 <data_disconts+0x3d0>
cmp    $0x3,%rcx
jb     1668 <data_disconts+0x548>
lea    (%r9,%r8,4),%rcx
add    $0xc,%rcx
nopl   0x0(%rax,%rax,1)
movl   $0x0,(%r9,%rdx,4)
movl   $0x0,-0xc(%rcx,%rdx,4)
movl   $0x0,0x4(%r9,%rdx,4)
movl   $0x0,-0x8(%rcx,%rdx,4)
movl   $0x0,0x8(%r9,%rdx,4)
movl   $0x0,-0x4(%rcx,%rdx,4)
movl   $0x0,0xc(%r9,%rdx,4)
movl   $0x0,(%rcx,%rdx,4)
add    $0x4,%rdx
cmp    %rdx,%rax
jne    1520 <data_disconts+0x400>
jmp    1668 <data_disconts+0x548>
mov    %eax,%edx
and    $0xfffffff8,%edx
lea    -0x8(%rdx),%rcx
mov    %rcx,%rdi
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %edi,%esi
and    $0x3,%esi
cmp    $0x18,%rcx
jae    1596 <data_disconts+0x476>
xor    %ebp,%ebp
jmp    1618 <data_disconts+0x4f8>
and    $0xfffffffffffffffc,%rdi
lea    0x70(%r9),%rbx
lea    (%r9,%r8,4),%rcx
add    $0x70,%rcx
xor    %ebp,%ebp
xorpd  %xmm0,%xmm0
nopl   0x0(%rax)
movupd %xmm0,-0x70(%rbx,%rbp,4)
movupd %xmm0,-0x60(%rbx,%rbp,4)
movupd %xmm0,-0x70(%rcx,%rbp,4)
movupd %xmm0,-0x60(%rcx,%rbp,4)
movupd %xmm0,-0x50(%rbx,%rbp,4)
movupd %xmm0,-0x40(%rbx,%rbp,4)
movupd %xmm0,-0x50(%rcx,%rbp,4)
movupd %xmm0,-0x40(%rcx,%rbp,4)
movupd %xmm0,-0x30(%rbx,%rbp,4)
movupd %xmm0,-0x20(%rbx,%rbp,4)
movupd %xmm0,-0x30(%rcx,%rbp,4)
movupd %xmm0,-0x20(%rcx,%rbp,4)
movupd %xmm0,-0x10(%rbx,%rbp,4)
movupd %xmm0,(%rbx,%rbp,4)
movupd %xmm0,-0x10(%rcx,%rbp,4)
movupd %xmm0,(%rcx,%rbp,4)
add    $0x20,%rbp
add    $0xfffffffffffffffc,%rdi
jne    15b0 <data_disconts+0x490>
test   %rsi,%rsi
je     165f <data_disconts+0x53f>
lea    (%r8,%rbp,1),%rcx
lea    (%r9,%rcx,4),%rcx
add    $0x10,%rcx
lea    (%r9,%rbp,4),%rdi
add    $0x10,%rdi
shl    $0x5,%rsi
xor    %ebp,%ebp
xorpd  %xmm0,%xmm0
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%rdi,%rbp,1)
movupd %xmm0,(%rdi,%rbp,1)
movupd %xmm0,-0x10(%rcx,%rbp,1)
movupd %xmm0,(%rcx,%rbp,1)
add    $0x20,%rbp
cmp    %rbp,%rsi
jne    1640 <data_disconts+0x520>
cmp    %rax,%rdx
jne    14d5 <data_disconts+0x3b5>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @data_disconts(float* %0, float* %1, float* %2, float* %3, float* %4, i32 %5, i32 %6, float* %7, float* %8, float* %9) {
%11 = alloca float*
%12 = alloca float*
%13 = alloca float*
%14 = alloca float*
%15 = alloca float*
%16 = alloca i32
%17 = alloca i32
%18 = alloca float*
%19 = alloca float*
%20 = alloca float*
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca float
store float* %0, float** %11
store float* %1, float** %12
store float* %2, float** %13
store float* %3, float** %14
store float* %4, float** %15
store i32 %5, i32* %16
store i32 %6, i32* %17
store float* %7, float** %18
store float* %8, float** %19
store float* %9, float** %20
%27 = load float*, float** %11
%28 = load float*, float** %18
%29 = load float*, float** %13
%30 = load float*, float** %14
%31 = load i32, i32* %16
%32 = load i32, i32* %17
%33 = call i32 (float*, float*, float*, float*, i32, i32, ...) bitcast (i32 (...)* @warp_image to i32 (float*, float*, float*, float*, i32, i32, ...)*)(float* %27, float* %28, float* %29, float* %30, i32 %31, i32 %32)
%34 = load float*, float** %18
%35 = load float*, float** %12
%36 = load float*, float** %19
%37 = load i32, i32* %16
%38 = load i32, i32* %17
%39 = call i32 (float*, float*, float*, i32, i32, ...) bitcast (i32 (...)* @horn_dt to i32 (float*, float*, float*, i32, i32, ...)*)(float* %34, float* %35, float* %36, i32 %37, i32 %38)
store i32 1, i32* %24
br label %40
40: ; preds = %103, %10
%41 = load i32, i32* %24
%42 = load i32, i32* %17
%43 = sub nsw i32 %42, 1
%44 = icmp slt i32 %41, %43
br i1 %44, label %45, label %106
45: ; preds = %40
store i32 1, i32* %25
br label %46
46: ; preds = %99, %45
%47 = load i32, i32* %25
%48 = load i32, i32* %16
%49 = sub nsw i32 %48, 1
%50 = icmp slt i32 %47, %49
br i1 %50, label %51, label %102
51: ; preds = %46
%52 = load i32, i32* %24
%53 = load i32, i32* %16
%54 = mul nsw i32 %52, %53
%55 = load i32, i32* %25
%56 = add nsw i32 %54, %55
store i32 %56, i32* %21
%57 = load float*, float** %15
%58 = load i32, i32* %21
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
%61 = load float, float* %60
%62 = fpext float %61 to double
%63 = fmul double %62, 0x3FF6A09E7098EF50
%64 = fptrunc double %63 to float
store float %64, float* %26
%65 = load float*, float** %19
%66 = load i32, i32* %21
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds float, float* %65, i64 %67
%69 = load float, float* %68
%70 = fcmp olt float %69, 0.000000e+00
br i1 %70, label %71, label %78
71: ; preds = %51
%72 = load float*, float** %19
%73 = load i32, i32* %21
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds float, float* %72, i64 %74
%76 = load float, float* %75
%77 = fneg float %76
br label %84
78: ; preds = %51
%79 = load float*, float** %19
%80 = load i32, i32* %21
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %79, i64 %81
%83 = load float, float* %82
br label %84
84: ; preds = %78, %71
%85 = phi float [ %77, %71 ], [ %83, %78 ]
%86 = load float, float* %26
%87 = fcmp oge float %85, %86
br i1 %87, label %88, label %93
88: ; preds = %84
%89 = load float*, float** %20
%90 = load i32, i32* %21
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds float, float* %89, i64 %91
store float 0.000000e+00, float* %92
br label %98
93: ; preds = %84
%94 = load float*, float** %20
%95 = load i32, i32* %21
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds float, float* %94, i64 %96
store float 1.000000e+00, float* %97
br label %98
98: ; preds = %93, %88
br label %99
99: ; preds = %98
%100 = load i32, i32* %25
%101 = add nsw i32 %100, 1
store i32 %101, i32* %25
br label %46
102: ; preds = %46
br label %103
103: ; preds = %102
%104 = load i32, i32* %24
%105 = add nsw i32 %104, 1
store i32 %105, i32* %24
br label %40
106: ; preds = %40
%107 = load i32, i32* %17
%108 = sub nsw i32 %107, 1
store i32 %108, i32* %23
%109 = load i32, i32* %16
%110 = sub nsw i32 %109, 1
store i32 %110, i32* %22
store i32 0, i32* %24
br label %111
111: ; preds = %130, %106
%112 = load i32, i32* %24
%113 = load i32, i32* %17
%114 = icmp slt i32 %112, %113
br i1 %114, label %115, label %133
115: ; preds = %111
%116 = load float*, float** %20
%117 = load i32, i32* %24
%118 = load i32, i32* %16
%119 = mul nsw i32 %117, %118
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %116, i64 %120
store float 0.000000e+00, float* %121
%122 = load float*, float** %20
%123 = load i32, i32* %24
%124 = load i32, i32* %16
%125 = mul nsw i32 %123, %124
%126 = load i32, i32* %22
%127 = add nsw i32 %125, %126
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds float, float* %122, i64 %128
store float 0.000000e+00, float* %129
br label %130
130: ; preds = %115
%131 = load i32, i32* %24
%132 = add nsw i32 %131, 1
store i32 %132, i32* %24
br label %111
133: ; preds = %111
store i32 0, i32* %25
br label %134
134: ; preds = %151, %133
%135 = load i32, i32* %25
%136 = load i32, i32* %16
%137 = icmp slt i32 %135, %136
br i1 %137, label %138, label %154
138: ; preds = %134
%139 = load float*, float** %20
%140 = load i32, i32* %25
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds float, float* %139, i64 %141
store float 0.000000e+00, float* %142
%143 = load float*, float** %20
%144 = load i32, i32* %23
%145 = load i32, i32* %16
%146 = mul nsw i32 %144, %145
%147 = load i32, i32* %25
%148 = add nsw i32 %146, %147
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds float, float* %143, i64 %149
store float 0.000000e+00, float* %150
br label %151
151: ; preds = %138
%152 = load i32, i32* %25
%153 = add nsw i32 %152, 1
store i32 %153, i32* %25
br label %134
154: ; preds = %134
ret void
}

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