void spatial_disconts(u, v, s2, nx, ny, discont) float *u, *v, *discont;
int nx, ny;
float *s2;
{
  float image_abs_max(), variance(), mean();
  int index, nx_max, ny_max, i, j;
  float t2;
  for (i = 1; i < ny - 1; i++) {
    for (j = 1; j < nx - 1; j++) {
      index = i * nx + j;
      t2 = s2[index] * 1.4142136;
      if ((((u[index] - u[index - nx]) < 0
                ? -(u[index] - u[index - nx])
                : (u[index] - u[index - nx])) >= t2) ||
          (((u[index] - u[index - 1]) < 0 ? -(u[index] - u[index - 1])
                                          : (u[index] - u[index - 1])) >= t2) ||
          (((v[index] - v[index - nx]) < 0
                ? -(v[index] - v[index - nx])
                : (v[index] - v[index - nx])) >= t2) ||
          (((v[index] - v[index - 1]) < 0 ? -(v[index] - v[index - 1])
                                          : (v[index] - v[index - 1])) >= t2))
        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):
<spatial_disconts>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    141e <spatial_disconts+0x31e>
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    140b <spatial_disconts+0x30b>
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe8f(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    11df <spatial_disconts+0xdf>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x44(%rbp)
jmp    1203 <spatial_disconts+0x103>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
ucomiss -0x40(%rbp),%xmm0
jae    13ce <spatial_disconts+0x2ce>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    1273 <spatial_disconts+0x173>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x48(%rbp)
jmp    1297 <spatial_disconts+0x197>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
ucomiss -0x40(%rbp),%xmm0
jae    13ce <spatial_disconts+0x2ce>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    1307 <spatial_disconts+0x207>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x4c(%rbp)
jmp    132b <spatial_disconts+0x22b>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x4c(%rbp),%xmm0
ucomiss -0x40(%rbp),%xmm0
jae    13ce <spatial_disconts+0x2ce>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    139b <spatial_disconts+0x29b>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x50(%rbp)
jmp    13bf <spatial_disconts+0x2bf>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
ucomiss -0x40(%rbp),%xmm0
jb     13e3 <spatial_disconts+0x2e3>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13f8 <spatial_disconts+0x2f8>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  0xc15(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13fd <spatial_disconts+0x2fd>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    113a <spatial_disconts+0x3a>
jmp    1410 <spatial_disconts+0x310>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1122 <spatial_disconts+0x22>
mov    -0x20(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1480 <spatial_disconts+0x380>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1437 <spatial_disconts+0x337>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    14ca <spatial_disconts+0x3ca>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x3c(%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    1487 <spatial_disconts+0x387>
pop    %rbp
ret

x86 assembly (clang, -O2):
<spatial_disconts>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    -0x1(%r8),%eax
mov    %eax,-0x4(%rsp)
lea    -0x1(%rcx),%eax
mov    %eax,-0x8(%rsp)
movslq %ecx,%r14
cmp    $0x3,%r8d
jl     129b <spatial_disconts+0x19b>
mov    -0x4(%rsp),%r15d
mov    -0x8(%rsp),%ebp
lea    (%r9,%r14,4),%rax
add    $0x4,%rax
lea    0x0(,%r14,4),%r12
lea    (%rdi,%r14,4),%rbx
add    $0x4,%rbx
lea    (%rsi,%r14,4),%r13
add    $0x4,%r13
lea    (%rdx,%r14,4),%rdx
add    $0x4,%rdx
add    $0x4,%rdi
add    $0x4,%rsi
lea    -0x4(,%rbp,4),%rbp
ff
mov    $0x1,%r11d
movsd  0xe8b(%rip),%xmm0
movaps 0xe94(%rip),%xmm1
movss  0xe9b(%rip),%xmm8
jmp    11af <spatial_disconts+0xaf>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r11
add    %r12,%rax
add    %r12,%rbx
add    %r12,%r13
add    %r12,%rdx
add    %r12,%rdi
add    %r12,%rsi
cmp    %r15,%r11
je     129b <spatial_disconts+0x19b>
cmp    $0x3,%ecx
jl     1190 <spatial_disconts+0x90>
xor    %r10d,%r10d
jmp    11cf <spatial_disconts+0xcf>
nopl   0x0(%rax)
movss  %xmm4,(%rax,%r10,1)
add    $0x4,%r10
cmp    %r10,%rbp
je     1190 <spatial_disconts+0x90>
movss  (%rdx,%r10,1),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm0,%xmm3
cvtsd2ss %xmm3,%xmm3
movss  (%rbx,%r10,1),%xmm2
movaps %xmm2,%xmm5
subss  (%rdi,%r10,1),%xmm5
movaps %xmm5,%xmm7
xorps  %xmm1,%xmm7
xorps  %xmm4,%xmm4
movaps %xmm5,%xmm6
cmpltss %xmm4,%xmm6
andps  %xmm6,%xmm7
andnps %xmm5,%xmm6
orps   %xmm7,%xmm6
ucomiss %xmm3,%xmm6
jae    11c0 <spatial_disconts+0xc0>
subss  -0x4(%rbx,%r10,1),%xmm2
movaps %xmm2,%xmm5
xorps  %xmm1,%xmm5
xorps  %xmm4,%xmm4
movaps %xmm2,%xmm6
cmpltss %xmm4,%xmm6
andps  %xmm6,%xmm5
andnps %xmm2,%xmm6
orps   %xmm5,%xmm6
ucomiss %xmm3,%xmm6
jae    11c0 <spatial_disconts+0xc0>
movss  0x0(%r13,%r10,1),%xmm5
movaps %xmm5,%xmm2
subss  (%rsi,%r10,1),%xmm2
movaps %xmm2,%xmm6
xorps  %xmm1,%xmm6
xorps  %xmm4,%xmm4
movaps %xmm2,%xmm7
cmpltss %xmm4,%xmm7
andps  %xmm7,%xmm6
andnps %xmm2,%xmm7
orps   %xmm6,%xmm7
ucomiss %xmm3,%xmm7
jae    11c0 <spatial_disconts+0xc0>
subss  -0x4(%r13,%r10,1),%xmm5
movaps %xmm5,%xmm2
xorps  %xmm1,%xmm2
xorps  %xmm4,%xmm4
movaps %xmm5,%xmm6
cmpltss %xmm4,%xmm6
andps  %xmm6,%xmm2
andnps %xmm5,%xmm6
orps   %xmm2,%xmm6
ucomiss %xmm3,%xmm6
jae    11c0 <spatial_disconts+0xc0>
movaps %xmm8,%xmm4
jmp    11c0 <spatial_disconts+0xc0>
test   %r8d,%r8d
jle    1315 <spatial_disconts+0x215>
movslq -0x8(%rsp),%rax
mov    %r8d,%r10d
cmp    $0x1,%r8d
jne    12b2 <spatial_disconts+0x1b2>
xor    %esi,%esi
jmp    12f8 <spatial_disconts+0x1f8>
mov    %r10d,%edi
and    $0xfffffffe,%edi
lea    (%r14,%rax,1),%rbp
lea    0x0(,%r14,8),%rbx
xor    %esi,%esi
mov    %r9,%rdx
nopl   0x0(%rax)
movl   $0x0,(%rdx)
movl   $0x0,(%rdx,%rax,4)
movl   $0x0,(%rdx,%r14,4)
movl   $0x0,(%rdx,%rbp,4)
add    $0x2,%rsi
add    %rbx,%rdx
cmp    %rsi,%rdi
jne    12d0 <spatial_disconts+0x1d0>
test   $0x1,%r10b
je     1315 <spatial_disconts+0x215>
imul   %r14,%rsi
movl   $0x0,(%r9,%rsi,4)
add    %rax,%rsi
movl   $0x0,(%r9,%rsi,4)
test   %ecx,%ecx
jle    14e8 <spatial_disconts+0x3e8>
mov    -0x4(%rsp),%eax
imul   %ecx,%eax
movslq %eax,%r8
mov    %ecx,%eax
cmp    $0x8,%ecx
jb     1350 <spatial_disconts+0x250>
lea    (%r8,%rax,1),%rcx
lea    (%r9,%rcx,4),%rcx
cmp    %r9,%rcx
jbe    13f0 <spatial_disconts+0x2f0>
lea    (%r9,%rax,4),%rcx
lea    (%r9,%r8,4),%rdx
cmp    %rcx,%rdx
jae    13f0 <spatial_disconts+0x2f0>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rsi
and    $0x3,%rsi
je     1389 <spatial_disconts+0x289>
lea    (%r9,%r8,4),%rdi
nopl   0x0(%rax,%rax,1)
movl   $0x0,(%r9,%rcx,4)
movl   $0x0,(%rdi,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1370 <spatial_disconts+0x270>
cmp    $0x3,%rdx
jb     14e8 <spatial_disconts+0x3e8>
lea    (%r9,%r8,4),%rdx
add    $0xc,%rdx
nopl   0x0(%rax,%rax,1)
movl   $0x0,(%r9,%rcx,4)
movl   $0x0,-0xc(%rdx,%rcx,4)
movl   $0x0,0x4(%r9,%rcx,4)
movl   $0x0,-0x8(%rdx,%rcx,4)
movl   $0x0,0x8(%r9,%rcx,4)
movl   $0x0,-0x4(%rdx,%rcx,4)
movl   $0x0,0xc(%r9,%rcx,4)
movl   $0x0,(%rdx,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    13a0 <spatial_disconts+0x2a0>
jmp    14e8 <spatial_disconts+0x3e8>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rdx
mov    %rdx,%rdi
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %edi,%esi
and    $0x3,%esi
cmp    $0x18,%rdx
jae    1416 <spatial_disconts+0x316>
xor    %ebp,%ebp
jmp    1498 <spatial_disconts+0x398>
and    $0xfffffffffffffffc,%rdi
lea    0x70(%r9),%rbx
lea    (%r9,%r8,4),%rdx
add    $0x70,%rdx
xor    %ebp,%ebp
xorpd  %xmm0,%xmm0
nopl   0x0(%rax)
movupd %xmm0,-0x70(%rbx,%rbp,4)
movupd %xmm0,-0x60(%rbx,%rbp,4)
movupd %xmm0,-0x70(%rdx,%rbp,4)
movupd %xmm0,-0x60(%rdx,%rbp,4)
movupd %xmm0,-0x50(%rbx,%rbp,4)
movupd %xmm0,-0x40(%rbx,%rbp,4)
movupd %xmm0,-0x50(%rdx,%rbp,4)
movupd %xmm0,-0x40(%rdx,%rbp,4)
movupd %xmm0,-0x30(%rbx,%rbp,4)
movupd %xmm0,-0x20(%rbx,%rbp,4)
movupd %xmm0,-0x30(%rdx,%rbp,4)
movupd %xmm0,-0x20(%rdx,%rbp,4)
movupd %xmm0,-0x10(%rbx,%rbp,4)
movupd %xmm0,(%rbx,%rbp,4)
movupd %xmm0,-0x10(%rdx,%rbp,4)
movupd %xmm0,(%rdx,%rbp,4)
add    $0x20,%rbp
add    $0xfffffffffffffffc,%rdi
jne    1430 <spatial_disconts+0x330>
test   %rsi,%rsi
je     14df <spatial_disconts+0x3df>
lea    (%r8,%rbp,1),%rdx
lea    (%r9,%rdx,4),%rdx
add    $0x10,%rdx
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(%rdx,%rbp,1)
movupd %xmm0,(%rdx,%rbp,1)
add    $0x20,%rbp
cmp    %rbp,%rsi
jne    14c0 <spatial_disconts+0x3c0>
cmp    %rax,%rcx
jne    1352 <spatial_disconts+0x252>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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