void disconts(u, v, prev_u, prev_v, du, dv, s1, s2, s3, nx, ny, Ix, Iy, It, err,
              discont) float *Ix,
    *Iy, *It, *u, *v, *du, *dv, *prev_u, *prev_v, *err, *discont;
int nx, ny;
float *s1, *s2, *s3;
{
  float image_abs_max(), variance(), mean();
  int index, nx_max, ny_max, i, j;
  float t1;
  float t3;
  inten_const(err, Ix, Iy, It, du, dv, 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;
      t3 = s3[index] * 1.4142136;
      if ((((err[index]) < 0 ? -(err[index]) : (err[index])) >= t1) ||
          (((u[index] - prev_u[index]) < 0
                ? -(u[index] - prev_u[index])
                : (u[index] - prev_u[index])) >= t3) ||
          (((v[index] - prev_v[index]) < 0 ? -(v[index] - prev_v[index])
                                           : (v[index] - prev_v[index])) >= t3))
        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)] = 1.0;
    discont[(i * nx) + nx_max] = 0.0;
  }
  for (j = 0; j < nx; j++) {
    discont[j] = 1.0;
    discont[(ny_max * nx) + j] = 0.0;
  }
}


x86 assembly (clang, -O0):
<disconts>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x58(%rbp),%rax
mov    0x50(%rbp),%rax
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    0x50(%rbp),%rdi
mov    0x38(%rbp),%rsi
mov    0x40(%rbp),%rdx
mov    0x48(%rbp),%rcx
mov    -0x28(%rbp),%r8
mov    -0x30(%rbp),%r9
mov    0x28(%rbp),%r10d
mov    0x30(%rbp),%eax
mov    %r10d,(%rsp)
mov    %eax,0x8(%rsp)
mov    $0x0,%al
call   1030 <inten_const@plt>
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    0x30(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13c1 <disconts+0x2b1>
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    0x28(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13ae <disconts+0x29e>
mov    -0x40(%rbp),%eax
imul   0x28(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe2e(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe08(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    0x50(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomiss (%rax,%rcx,4),%xmm0
jbe    1246 <disconts+0x136>
mov    0x50(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x50(%rbp)
jmp    1258 <disconts+0x148>
mov    0x50(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
ucomiss -0x48(%rbp),%xmm0
jae    1371 <disconts+0x261>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    12be <disconts+0x1ae>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    12dd <disconts+0x1cd>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x54(%rbp),%xmm0
ucomiss -0x4c(%rbp),%xmm0
jae    1371 <disconts+0x261>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    1343 <disconts+0x233>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x58(%rbp)
jmp    1362 <disconts+0x252>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x58(%rbp),%xmm0
ucomiss -0x4c(%rbp),%xmm0
jb     1386 <disconts+0x276>
mov    0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    139b <disconts+0x28b>
mov    0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  0xc6a(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13a0 <disconts+0x290>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11a3 <disconts+0x93>
jmp    13b3 <disconts+0x2a3>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    118b <disconts+0x7b>
mov    0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    0x28(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    0x30(%rbp),%eax
jge    1428 <disconts+0x318>
mov    0x58(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  0xc04(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    0x58(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   0x28(%rbp),%ecx
add    -0x38(%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    13da <disconts+0x2ca>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    1477 <disconts+0x367>
mov    0x58(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  0xbb5(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    0x58(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   0x28(%rbp),%ecx
add    -0x44(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    142f <disconts+0x31f>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<disconts>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %rcx,%r15
mov    %rdx,%rbp
mov    %rsi,%r12
mov    %rdi,%r14
mov    0x80(%rsp),%rdi
mov    0x60(%rsp),%ebx
mov    0x58(%rsp),%r13d
mov    0x68(%rsp),%rsi
mov    0x70(%rsp),%rdx
mov    0x78(%rsp),%rcx
xor    %eax,%eax
push   %rbx
push   %r13
call   1030 <inten_const@plt>
add    $0x10,%rsp
lea    -0x1(%rbx),%eax
mov    %eax,0x4(%rsp)
lea    -0x1(%r13),%eax
mov    %eax,(%rsp)
movslq %r13d,%r10
cmp    $0x3,%ebx
jl     12f5 <disconts+0x1e5>
mov    0x50(%rsp),%rcx
mov    0x40(%rsp),%rax
mov    0x4(%rsp),%r11d
mov    (%rsp),%edx
mov    0x88(%rsp),%rsi
lea    0x4(%rsi,%r10,4),%rdi
lea    0x0(,%r10,4),%r8
lea    0x4(%rax,%r10,4),%rax
lea    0x4(%rcx,%r10,4),%rbx
mov    %rbp,%rsi
mov    0x80(%rsp),%rcx
lea    0x4(%rcx,%r10,4),%rbp
lea    0x4(%r14,%r10,4),%rcx
lea    0x4(%rsi,%r10,4),%rsi
lea    0x4(%r12,%r10,4),%r14
lea    0x4(%r15,%r10,4),%r12
lea    -0x4(,%rdx,4),%rdx
ff
mov    $0x1,%r9d
movsd  0xe26(%rip),%xmm0
movaps 0xe2f(%rip),%xmm1
movss  0xe47(%rip),%xmm2
jmp    1218 <disconts+0x108>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r9
add    %r8,%rdi
add    %r8,%rax
add    %r8,%rbx
add    %r8,%rbp
add    %r8,%rcx
add    %r8,%rsi
add    %r8,%r14
add    %r8,%r12
cmp    %r11,%r9
mov    %r15,%r13
je     12f5 <disconts+0x1e5>
mov    %r13,%r15
cmp    $0x3,%r13d
jl     11f0 <disconts+0xe0>
xor    %r13d,%r13d
jmp    123f <disconts+0x12f>
cs nopw 0x0(%rax,%rax,1)
movss  %xmm3,(%rdi,%r13,1)
add    $0x4,%r13
cmp    %r13,%rdx
je     11f0 <disconts+0xe0>
movss  (%rax,%r13,1),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm0,%xmm3
xorps  %xmm4,%xmm4
cvtsd2ss %xmm3,%xmm4
movss  0x0(%rbp,%r13,1),%xmm5
movaps %xmm5,%xmm6
xorps  %xmm1,%xmm6
xorps  %xmm3,%xmm3
movaps %xmm5,%xmm7
cmpltss %xmm3,%xmm7
andps  %xmm7,%xmm6
andnps %xmm5,%xmm7
orps   %xmm6,%xmm7
ucomiss %xmm4,%xmm7
jae    1230 <disconts+0x120>
movss  (%rbx,%r13,1),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm0,%xmm3
xorps  %xmm4,%xmm4
cvtsd2ss %xmm3,%xmm4
movss  (%rcx,%r13,1),%xmm5
subss  (%rsi,%r13,1),%xmm5
movaps %xmm5,%xmm6
xorps  %xmm1,%xmm6
xorps  %xmm3,%xmm3
movaps %xmm5,%xmm7
cmpltss %xmm3,%xmm7
andps  %xmm7,%xmm6
andnps %xmm5,%xmm7
orps   %xmm6,%xmm7
ucomiss %xmm4,%xmm7
jae    1230 <disconts+0x120>
movss  (%r14,%r13,1),%xmm5
subss  (%r12,%r13,1),%xmm5
movaps %xmm5,%xmm6
xorps  %xmm1,%xmm6
xorps  %xmm3,%xmm3
movaps %xmm5,%xmm7
cmpltss %xmm3,%xmm7
andps  %xmm7,%xmm6
andnps %xmm5,%xmm7
orps   %xmm6,%xmm7
ucomiss %xmm4,%xmm7
jae    1230 <disconts+0x120>
movaps %xmm2,%xmm3
jmp    1230 <disconts+0x120>
mov    0x60(%rsp),%edx
test   %edx,%edx
mov    0x88(%rsp),%r9
jle    1374 <disconts+0x264>
movslq (%rsp),%rax
mov    %edx,%ecx
cmp    $0x1,%edx
jne    1314 <disconts+0x204>
xor    %edx,%edx
jmp    1358 <disconts+0x248>
mov    %ecx,%esi
and    $0xfffffffe,%esi
lea    (%r10,%rax,1),%rdi
lea    0x0(,%r10,8),%rbp
xor    %edx,%edx
mov    %r9,%rbx
nopw   0x0(%rax,%rax,1)
movl   $0x3f800000,(%rbx)
movl   $0x0,(%rbx,%rax,4)
movl   $0x3f800000,(%rbx,%r10,4)
3f
movl   $0x0,(%rbx,%rdi,4)
add    $0x2,%rdx
add    %rbp,%rbx
cmp    %rdx,%rsi
jne    1330 <disconts+0x220>
test   $0x1,%cl
je     1374 <disconts+0x264>
imul   %r10,%rdx
movl   $0x3f800000,(%r9,%rdx,4)
3f
add    %rax,%rdx
movl   $0x0,(%r9,%rdx,4)
test   %r13d,%r13d
jle    1556 <disconts+0x446>
mov    0x4(%rsp),%eax
imul   %r13d,%eax
movslq %eax,%r8
mov    %r13d,%eax
cmp    $0x8,%r13d
jb     13b3 <disconts+0x2a3>
lea    (%r8,%rax,1),%rcx
lea    (%r9,%rcx,4),%rcx
cmp    %r9,%rcx
jbe    1450 <disconts+0x340>
lea    (%r9,%rax,4),%rcx
lea    (%r9,%r8,4),%rdx
cmp    %rcx,%rdx
jae    1450 <disconts+0x340>
xor    %edx,%edx
mov    %rdx,%rcx
not    %rcx
add    %rax,%rcx
mov    %rax,%rsi
and    $0x3,%rsi
je     13e9 <disconts+0x2d9>
lea    (%r9,%r8,4),%rdi
nopl   0x0(%rax,%rax,1)
movl   $0x3f800000,(%r9,%rdx,4)
3f
movl   $0x0,(%rdi,%rdx,4)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rsi
jne    13d0 <disconts+0x2c0>
cmp    $0x3,%rcx
jb     1556 <disconts+0x446>
lea    (%r9,%r8,4),%rcx
add    $0xc,%rcx
nopl   0x0(%rax,%rax,1)
movl   $0x3f800000,(%r9,%rdx,4)
3f
movl   $0x0,-0xc(%rcx,%rdx,4)
movl   $0x3f800000,0x4(%r9,%rdx,4)
80 3f
movl   $0x0,-0x8(%rcx,%rdx,4)
movl   $0x3f800000,0x8(%r9,%rdx,4)
80 3f
movl   $0x0,-0x4(%rcx,%rdx,4)
movl   $0x3f800000,0xc(%r9,%rdx,4)
80 3f
movl   $0x0,(%rcx,%rdx,4)
add    $0x4,%rdx
cmp    %rdx,%rax
jne    1400 <disconts+0x2f0>
jmp    1556 <disconts+0x446>
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    1476 <disconts+0x366>
xor    %ebp,%ebp
jmp    1500 <disconts+0x3f0>
and    $0xfffffffffffffffc,%rdi
lea    0x70(%r9),%rbx
lea    (%r9,%r8,4),%rcx
add    $0x70,%rcx
xor    %ebp,%ebp
movapd 0xb90(%rip),%xmm0
xorps  %xmm1,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd %xmm0,-0x70(%rbx,%rbp,4)
movupd %xmm0,-0x60(%rbx,%rbp,4)
movups %xmm1,-0x70(%rcx,%rbp,4)
movups %xmm1,-0x60(%rcx,%rbp,4)
movupd %xmm0,-0x50(%rbx,%rbp,4)
movupd %xmm0,-0x40(%rbx,%rbp,4)
movups %xmm1,-0x50(%rcx,%rbp,4)
movups %xmm1,-0x40(%rcx,%rbp,4)
movupd %xmm0,-0x30(%rbx,%rbp,4)
movupd %xmm0,-0x20(%rbx,%rbp,4)
movups %xmm1,-0x30(%rcx,%rbp,4)
movups %xmm1,-0x20(%rcx,%rbp,4)
movupd %xmm0,-0x10(%rbx,%rbp,4)
movupd %xmm0,(%rbx,%rbp,4)
movups %xmm1,-0x10(%rcx,%rbp,4)
movups %xmm1,(%rcx,%rbp,4)
add    $0x20,%rbp
add    $0xfffffffffffffffc,%rdi
jne    14a0 <disconts+0x390>
test   %rsi,%rsi
je     154d <disconts+0x43d>
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
movapd 0xaf9(%rip),%xmm0
xorps  %xmm1,%xmm1
nopw   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%rdi,%rbp,1)
movupd %xmm0,(%rdi,%rbp,1)
movups %xmm1,-0x10(%rcx,%rbp,1)
movups %xmm1,(%rcx,%rbp,1)
add    $0x20,%rbp
cmp    %rbp,%rsi
jne    1530 <disconts+0x420>
cmp    %rax,%rdx
jne    13b5 <disconts+0x2a5>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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