#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ assert(int);
int /*<<< orphan*/ central2d_correct_sd(float *, float *, float *, float *,
                                        float const *, float const *,
                                        float const *, float, float, int, int);
int /*<<< orphan*/ limited_deriv1(float *, float const *, int);
int /*<<< orphan*/ limited_derivk(float *, float const *, int, int);
void central2d_correct(float *restrict v, float *restrict scratch,
                       const float *restrict u, const float *restrict f,
                       const float *restrict g, float dtcdx2, float dtcdy2,
                       int xlo, int xhi, int ylo, int yhi, int nx, int ny,
                       int nfield) {
  assert(0 <= xlo && xlo < xhi && xhi <= nx);
  assert(0 <= ylo && ylo < yhi && yhi <= ny);
  float *restrict ux = scratch;
  float *restrict uy = scratch + nx;
  float *restrict s0 = scratch + 2 * nx;
  float *restrict d0 = scratch + 3 * nx;
  float *restrict s1 = scratch + 4 * nx;
  float *restrict d1 = scratch + 5 * nx;
  for (int k = 0; k < nfield; ++k) {
    float *restrict vk = v + k * ny * nx;
    const float *restrict uk = u + k * ny * nx;
    const float *restrict fk = f + k * ny * nx;
    const float *restrict gk = g + k * ny * nx;
    limited_deriv1(ux + 1, uk + ylo * nx + 1, nx - 2);
    limited_derivk(uy + 1, uk + ylo * nx + 1, nx - 2, nx);
    central2d_correct_sd(s1, d1, ux, uy, uk + ylo * nx, fk + ylo * nx,
                         gk + ylo * nx, dtcdx2, dtcdy2, xlo, xhi);
    for (int iy = ylo; iy < yhi; ++iy) {
      float *tmp;
      tmp = s0;
      s0 = s1;
      s1 = tmp;
      tmp = d0;
      d0 = d1;
      d1 = tmp;
      limited_deriv1(ux + 1, uk + (iy + 1) * nx + 1, nx - 2);
      limited_derivk(uy + 1, uk + (iy + 1) * nx + 1, nx - 2, nx);
      central2d_correct_sd(s1, d1, ux, uy, uk + (iy + 1) * nx,
                           fk + (iy + 1) * nx, gk + (iy + 1) * nx, dtcdx2,
                           dtcdy2, xlo, xhi);
      for (int ix = xlo; ix < xhi; ++ix)
        vk[iy * nx + ix] = (s1[ix] + s0[ix]) - (d1[ix] - d0[ix]);
    }
  }
}


x86 assembly (clang, -O0):
<central2d_correct>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movss  %xmm0,-0x2c(%rbp)
movss  %xmm1,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
xor    %eax,%eax
xor    %ecx,%ecx
cmp    -0x34(%rbp),%ecx
mov    %al,-0xa5(%rbp)
jg     11b5 <central2d_correct+0x75>
mov    -0x34(%rbp),%ecx
xor    %eax,%eax
cmp    0x10(%rbp),%ecx
mov    %al,-0xa5(%rbp)
jge    11b5 <central2d_correct+0x75>
mov    0x10(%rbp),%eax
cmp    0x28(%rbp),%eax
setle  %al
mov    %al,-0xa5(%rbp)
mov    -0xa5(%rbp),%al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
xor    %eax,%eax
xor    %ecx,%ecx
cmp    0x18(%rbp),%ecx
mov    %al,-0xa6(%rbp)
jg     11fb <central2d_correct+0xbb>
mov    0x18(%rbp),%ecx
xor    %eax,%eax
cmp    0x20(%rbp),%ecx
mov    %al,-0xa6(%rbp)
jge    11fb <central2d_correct+0xbb>
mov    0x20(%rbp),%eax
cmp    0x30(%rbp),%eax
setle  %al
mov    %al,-0xa6(%rbp)
mov    -0xa6(%rbp),%al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x10(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq 0x28(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    0x28(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
imul   $0x3,0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    0x28(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x10(%rbp),%rax
imul   $0x5,0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x68(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    0x38(%rbp),%eax
jge    15d5 <central2d_correct+0x495>
mov    -0x8(%rbp),%rax
mov    -0x6c(%rbp),%ecx
imul   0x30(%rbp),%ecx
imul   0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x78(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
imul   0x30(%rbp),%ecx
imul   0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x80(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x6c(%rbp),%ecx
imul   0x30(%rbp),%ecx
imul   0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x88(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x6c(%rbp),%ecx
imul   0x30(%rbp),%ecx
imul   0x28(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x90(%rbp)
mov    -0x40(%rbp),%rdi
add    $0x4,%rdi
mov    -0x80(%rbp),%rsi
mov    0x18(%rbp),%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rsi
add    $0x4,%rsi
mov    0x28(%rbp),%edx
sub    $0x2,%edx
call   1030 <limited_deriv1@plt>
mov    -0x48(%rbp),%rdi
add    $0x4,%rdi
mov    -0x80(%rbp),%rsi
mov    0x18(%rbp),%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rsi
add    $0x4,%rsi
mov    0x28(%rbp),%edx
sub    $0x2,%edx
mov    0x28(%rbp),%ecx
call   1060 <limited_derivk@plt>
mov    -0x60(%rbp),%rdi
mov    -0x68(%rbp),%rsi
mov    -0x40(%rbp),%rdx
mov    -0x48(%rbp),%rcx
mov    -0x80(%rbp),%r8
mov    0x18(%rbp),%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%r8
mov    -0x88(%rbp),%r9
mov    0x18(%rbp),%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%r9
mov    -0x90(%rbp),%r11
mov    0x18(%rbp),%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%r11
movss  -0x2c(%rbp),%xmm0
movss  -0x30(%rbp),%xmm1
mov    -0x34(%rbp),%r10d
mov    0x10(%rbp),%eax
mov    %r11,(%rsp)
mov    %r10d,0x8(%rsp)
mov    %eax,0x10(%rsp)
call   1050 <central2d_correct_sd@plt>
mov    0x18(%rbp),%eax
mov    %eax,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    15c2 <central2d_correct+0x482>
mov    -0x50(%rbp),%rax
mov    %rax,-0xa0(%rbp)
mov    -0x60(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0xa0(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x58(%rbp),%rax
mov    %rax,-0xa0(%rbp)
mov    -0x68(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0xa0(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x40(%rbp),%rdi
add    $0x4,%rdi
mov    -0x80(%rbp),%rsi
mov    -0x94(%rbp),%eax
add    $0x1,%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rsi
add    $0x4,%rsi
mov    0x28(%rbp),%edx
sub    $0x2,%edx
call   1030 <limited_deriv1@plt>
mov    -0x48(%rbp),%rdi
add    $0x4,%rdi
mov    -0x80(%rbp),%rsi
mov    -0x94(%rbp),%eax
add    $0x1,%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rsi
add    $0x4,%rsi
mov    0x28(%rbp),%edx
sub    $0x2,%edx
mov    0x28(%rbp),%ecx
call   1060 <limited_derivk@plt>
mov    -0x60(%rbp),%rdi
mov    -0x68(%rbp),%rsi
mov    -0x40(%rbp),%rdx
mov    -0x48(%rbp),%rcx
mov    -0x80(%rbp),%r8
mov    -0x94(%rbp),%eax
add    $0x1,%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%r8
mov    -0x88(%rbp),%r9
mov    -0x94(%rbp),%eax
add    $0x1,%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%r9
mov    -0x90(%rbp),%r11
mov    -0x94(%rbp),%eax
add    $0x1,%eax
imul   0x28(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%r11
movss  -0x2c(%rbp),%xmm0
movss  -0x30(%rbp),%xmm1
mov    -0x34(%rbp),%r10d
mov    0x10(%rbp),%eax
mov    %r11,(%rsp)
mov    %r10d,0x8(%rsp)
mov    %eax,0x10(%rsp)
call   1050 <central2d_correct_sd@plt>
mov    -0x34(%rbp),%eax
mov    %eax,-0xa4(%rbp)
mov    -0xa4(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    15a9 <central2d_correct+0x469>
mov    -0x60(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0xa4(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x68(%rbp),%rax
movslq -0xa4(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x58(%rbp),%rax
movslq -0xa4(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
mov    -0x78(%rbp),%rax
mov    -0x94(%rbp),%ecx
imul   0x28(%rbp),%ecx
add    -0xa4(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1526 <central2d_correct+0x3e6>
jmp    15ae <central2d_correct+0x46e>
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    13e6 <central2d_correct+0x2a6>
jmp    15c7 <central2d_correct+0x487>
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1289 <central2d_correct+0x149>
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<central2d_correct>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x168,%rsp
movss  %xmm1,0x2c(%rsp)
movss  %xmm0,0x28(%rsp)
mov    %r8,0x100(%rsp)
mov    %rcx,0xf8(%rsp)
mov    %rdx,0xf0(%rsp)
mov    %rsi,0x50(%rsp)
mov    %rdi,0x48(%rsp)
mov    0x1c8(%rsp),%r15d
mov    0x1c0(%rsp),%ebp
mov    0x1b8(%rsp),%edx
mov    0x1b0(%rsp),%r14d
mov    0x1a8(%rsp),%ebx
mov    0x1a0(%rsp),%esi
test   %r9d,%r9d
setns  %al
mov    %r9,0x58(%rsp)
cmp    %esi,%r9d
setl   %cl
and    %al,%cl
cmp    %edx,%esi
setle  %al
and    %cl,%al
movzbl %al,%edi
call   1040 <assert@plt>
test   %ebx,%ebx
setns  %al
cmp    %r14d,%ebx
setl   %cl
and    %al,%cl
cmp    %ebp,%r14d
setle  %al
and    %cl,%al
movzbl %al,%edi
call   1040 <assert@plt>
test   %r15d,%r15d
jle    18a0 <central2d_correct+0x760>
mov    0x1b8(%rsp),%edx
mov    %rdx,%r9
movslq %r9d,%r12
lea    (%rdx,%rdx,4),%eax
cltq
mov    0x50(%rsp),%rcx
lea    (%rcx,%rax,4),%r11
lea    0x0(,%rdx,4),%eax
cltq
lea    (%rcx,%rax,4),%rax
mov    %rax,0x10(%rsp)
lea    (%rdx,%rdx,2),%eax
cltq
lea    (%rcx,%rax,4),%rax
mov    %rax,0x8(%rsp)
lea    (%rdx,%rdx,1),%eax
cltq
lea    (%rcx,%rax,4),%r14
imul   %r9d,%ebp
mov    %r9d,%eax
mov    0x1a8(%rsp),%r8d
imul   %r8d,%eax
cltq
mov    %rax,0xe8(%rsp)
movslq 0x58(%rsp),%rsi
movslq %r8d,%rdx
mov    %r15d,%eax
mov    %rax,0xe0(%rsp)
mov    0x1a0(%rsp),%eax
movslq %eax,%r10
mov    %r12,%rax
mov    %rdx,0xd8(%rsp)
imul   %rdx,%rax
lea    (%rax,%rsi,1),%r8
mov    %r10,%rdx
sub    %rsi,%rdx
mov    %rdx,0x78(%rsp)
and    $0xfffffffffffffff8,%rdx
mov    %r10,%rdi
neg    %rdi
mov    %rdi,0x140(%rsp)
mov    0x48(%rsp),%rdi
lea    (%rdi,%rax,4),%rbx
add    $0x4,%rbx
mov    %rbx,0xc0(%rsp)
add    %r10,%rax
mov    %rax,0xd0(%rsp)
lea    (%rcx,%r12,4),%rax
add    $0x4,%rcx
mov    %rcx,0x88(%rsp)
lea    -0x2(%r9),%ecx
mov    %ecx,0x24(%rsp)
mov    %rax,0x90(%rsp)
lea    0x4(%rax),%rax
mov    %rax,0x80(%rsp)
movslq 0x1b0(%rsp),%rax
mov    %rax,0x150(%rsp)
mov    %rdx,0x120(%rsp)
lea    (%rdx,%rsi,1),%rax
mov    %rax,0x118(%rsp)
mov    %rsi,0x18(%rsp)
lea    0x4(%rsi),%rax
mov    %rax,0x110(%rsp)
mov    %r8,0xc8(%rsp)
lea    0x10(%rdi,%r8,4),%rax
mov    %rax,0xb8(%rsp)
mov    %r12,0x60(%rsp)
lea    0x0(,%r12,4),%rax
mov    %rax,0x148(%rsp)
xor    %eax,%eax
mov    %rax,0x40(%rsp)
xor    %ecx,%ecx
mov    %ebp,0x74(%rsp)
mov    %r10,0x30(%rsp)
jmp    13a9 <central2d_correct+0x269>
mov    0x8(%rsp),%r15
mov    %r14,%r12
mov    0x38(%rsp),%r11
mov    0x10(%rsp),%rbp
mov    %rbp,0x10(%rsp)
mov    0x108(%rsp),%rcx
add    $0x1,%rcx
mov    0x74(%rsp),%ebp
mov    0x40(%rsp),%rax
add    %ebp,%eax
mov    %rax,0x40(%rsp)
mov    %r12,%r14
mov    %r15,0x8(%rsp)
cmp    0xe0(%rsp),%rcx
je     18a0 <central2d_correct+0x760>
mov    %ebp,%eax
mov    %rcx,0x108(%rsp)
imul   %ecx,%eax
movslq %eax,%r15
mov    0xf0(%rsp),%rax
lea    (%rax,%r15,4),%rcx
mov    0xf8(%rsp),%rax
lea    (%rax,%r15,4),%rax
mov    %rax,0xa8(%rsp)
mov    0x100(%rsp),%rax
lea    (%rax,%r15,4),%rax
mov    %rax,0xa0(%rsp)
mov    0xe8(%rsp),%rbp
lea    (%rcx,%rbp,4),%rax
mov    %rax,0x68(%rsp)
mov    %rcx,0x98(%rsp)
lea    0x4(%rcx,%rbp,4),%r12
mov    0x88(%rsp),%rdi
mov    %r12,%rsi
mov    0x24(%rsp),%ebx
mov    %ebx,%edx
mov    %r11,%r13
call   1030 <limited_deriv1@plt>
mov    0x80(%rsp),%rdi
mov    %r12,%rsi
mov    0xa8(%rsp),%r12
mov    %ebx,%edx
mov    0xa0(%rsp),%rbx
mov    0x1b8(%rsp),%ecx
call   1060 <limited_derivk@plt>
lea    (%r12,%rbp,4),%r9
lea    (%rbx,%rbp,4),%rax
sub    $0x8,%rsp
mov    0x18(%rsp),%rdi
mov    %r13,0x40(%rsp)
mov    %r13,%rsi
mov    0x58(%rsp),%rdx
mov    0x98(%rsp),%rcx
mov    0x70(%rsp),%r8
movss  0x30(%rsp),%xmm0
movss  0x34(%rsp),%xmm1
mov    0x1a8(%rsp),%ebp
push   %rbp
push   0x68(%rsp)
push   %rax
call   1050 <central2d_correct_sd@plt>
add    $0x20,%rsp
mov    0x1a8(%rsp),%eax
cmp    0x1b0(%rsp),%eax
jge    1360 <central2d_correct+0x220>
mov    0xc8(%rsp),%rax
add    %r15,%rax
mov    %rax,0x138(%rsp)
mov    0xd0(%rsp),%rax
lea    (%rax,%r15,1),%rax
mov    %rax,0x130(%rsp)
movslq 0x40(%rsp),%rax
mov    0xc0(%rsp),%rcx
lea    (%rcx,%rax,4),%r13
mov    0xb8(%rsp),%rcx
lea    (%rcx,%rax,4),%rcx
mov    0x48(%rsp),%rax
lea    (%rax,%r15,4),%rax
mov    %rax,0x128(%rsp)
xor    %eax,%eax
mov    %rax,0xb0(%rsp)
mov    0xd8(%rsp),%rdx
mov    0x38(%rsp),%r11
mov    0x10(%rsp),%rbp
mov    0x98(%rsp),%rax
mov    0x8(%rsp),%rsi
jmp    156d <central2d_correct+0x42d>
xchg   %ax,%ax
addq   $0x1,0xb0(%rsp)
mov    0x148(%rsp),%rax
mov    0x8(%rsp),%rcx
add    %rax,%rcx
add    %rax,%r13
mov    %r12,%r14
mov    %r15,%rsi
mov    0x68(%rsp),%rdx
cmp    0x150(%rsp),%rdx
mov    0x98(%rsp),%rax
je     1372 <central2d_correct+0x232>
mov    %rcx,0x8(%rsp)
mov    %r11,%r15
mov    %rbp,%r12
mov    %rsi,0x38(%rsp)
mov    %r14,0x10(%rsp)
mov    %rdx,0x158(%rsp)
lea    0x1(%rdx),%rbp
mov    %rbp,0x68(%rsp)
imul   0x60(%rsp),%rbp
lea    (%rax,%rbp,4),%rcx
mov    %rcx,0x160(%rsp)
lea    (%rax,%rbp,4),%r14
add    $0x4,%r14
mov    0x88(%rsp),%rdi
mov    %r14,%rsi
mov    0x24(%rsp),%ebx
mov    %ebx,%edx
call   1030 <limited_deriv1@plt>
mov    0x80(%rsp),%rdi
mov    %r14,%rsi
mov    %ebx,%edx
mov    0x1b8(%rsp),%ecx
call   1060 <limited_derivk@plt>
mov    0xa8(%rsp),%rax
lea    (%rax,%rbp,4),%r9
mov    0xa0(%rsp),%rax
lea    (%rax,%rbp,4),%rax
sub    $0x8,%rsp
mov    0x18(%rsp),%rbx
mov    %rbx,%rdi
mov    0x40(%rsp),%r14
mov    %r14,%rsi
mov    0x58(%rsp),%rdx
mov    0x98(%rsp),%rcx
mov    0x168(%rsp),%r8
movss  0x30(%rsp),%xmm0
movss  0x34(%rsp),%xmm1
mov    0x1a8(%rsp),%ebp
push   %rbp
mov    0x68(%rsp),%rbp
push   %rbp
push   %rax
call   1050 <central2d_correct_sd@plt>
mov    0x50(%rsp),%r10
mov    %r14,%r11
add    $0x20,%rsp
mov    0x1a0(%rsp),%eax
cmp    %eax,%ebp
mov    %rbx,%rbp
jge    1530 <central2d_correct+0x3f0>
mov    0x18(%rsp),%rax
cmpq   $0x8,0x78(%rsp)
jae    1730 <central2d_correct+0x5f0>
mov    0x1a0(%rsp),%edx
sub    %eax,%edx
mov    %rax,%rcx
test   $0x1,%dl
je     16bb <central2d_correct+0x57b>
movss  0x0(%rbp,%rax,4),%xmm0
addss  (%r12,%rax,4),%xmm0
mov    0x158(%rsp),%rdx
imul   0x60(%rsp),%rdx
movss  (%r11,%rax,4),%xmm1
subss  (%r15,%rax,4),%xmm1
subss  %xmm1,%xmm0
add    %rax,%rdx
mov    0x128(%rsp),%rcx
movss  %xmm0,(%rcx,%rdx,4)
lea    0x1(%rax),%rcx
not    %rax
cmp    0x140(%rsp),%rax
je     1530 <central2d_correct+0x3f0>
nopl   0x0(%rax)
movss  0x0(%rbp,%rcx,4),%xmm0
addss  (%r12,%rcx,4),%xmm0
movss  (%r11,%rcx,4),%xmm1
subss  (%r15,%rcx,4),%xmm1
subss  %xmm1,%xmm0
movss  0x4(%rbp,%rcx,4),%xmm1
addss  0x4(%r12,%rcx,4),%xmm1
movss  %xmm0,-0x4(%r13,%rcx,4)
movss  0x4(%r11,%rcx,4),%xmm0
subss  0x4(%r15,%rcx,4),%xmm0
subss  %xmm0,%xmm1
movss  %xmm1,0x0(%r13,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%r10
jne    16d0 <central2d_correct+0x590>
jmp    1530 <central2d_correct+0x3f0>
nopw   0x0(%rax,%rax,1)
mov    0xb0(%rsp),%rax
imul   0x60(%rsp),%rax
mov    0x138(%rsp),%rcx
add    %rax,%rcx
add    0x130(%rsp),%rax
mov    0x48(%rsp),%rdx
lea    (%rdx,%rcx,4),%rsi
lea    (%rdx,%rax,4),%rcx
mov    0x18(%rsp),%rax
lea    (%r15,%rax,4),%rdi
lea    (%r15,%r10,4),%rdx
lea    (%r11,%rax,4),%r8
mov    %rbp,%r9
lea    (%r11,%r10,4),%rbp
cmp    %rdx,%rsi
mov    %r11,%rbx
setb   %r10b
cmp    %rcx,%rdi
lea    (%r12,%rax,4),%rdi
setb   %r11b
cmp    %rbp,%rsi
mov    %r9,%rbp
setb   %r14b
cmp    %rcx,%r8
mov    0x30(%rsp),%rdx
lea    (%r12,%rdx,4),%rdx
setb   %r9b
cmp    %rdx,%rsi
setb   %dl
cmp    %rcx,%rdi
mov    0x30(%rsp),%rdi
lea    0x0(%rbp,%rdi,4),%rdi
setb   %r8b
cmp    %rdi,%rsi
lea    0x0(%rbp,%rax,4),%rsi
setb   %dil
cmp    %rcx,%rsi
setb   %sil
test   %r11b,%r10b
mov    0x30(%rsp),%r10
mov    %rbx,%r11
jne    166c <central2d_correct+0x52c>
mov    0x18(%rsp),%rax
and    %r9b,%r14b
jne    166c <central2d_correct+0x52c>
mov    0x18(%rsp),%rax
and    %r8b,%dl
jne    166c <central2d_correct+0x52c>
mov    0x18(%rsp),%rax
and    %sil,%dil
jne    166c <central2d_correct+0x52c>
mov    0x110(%rsp),%rsi
lea    0x0(,%rsi,4),%rax
add    %rbp,%rax
lea    (%r12,%rsi,4),%rcx
lea    (%r11,%rsi,4),%rdx
lea    (%r15,%rsi,4),%rsi
xor    %edi,%edi
mov    0x120(%rsp),%r8
mov    0x8(%rsp),%rbx
nopl   0x0(%rax,%rax,1)
movups -0x10(%rax,%rdi,4),%xmm0
movups (%rax,%rdi,4),%xmm1
movups -0x10(%rcx,%rdi,4),%xmm2
addps  %xmm0,%xmm2
movups (%rcx,%rdi,4),%xmm0
addps  %xmm1,%xmm0
movups -0x10(%rdx,%rdi,4),%xmm1
movups (%rdx,%rdi,4),%xmm3
movups -0x10(%rsi,%rdi,4),%xmm4
subps  %xmm4,%xmm1
subps  %xmm1,%xmm2
movups (%rsi,%rdi,4),%xmm1
subps  %xmm1,%xmm3
subps  %xmm3,%xmm0
movups %xmm2,-0x10(%rbx,%rdi,4)
movups %xmm0,(%rbx,%rdi,4)
add    $0x8,%rdi
cmp    %rdi,%r8
jne    1840 <central2d_correct+0x700>
mov    0x118(%rsp),%rax
cmp    %r8,0x78(%rsp)
je     1530 <central2d_correct+0x3f0>
jmp    166c <central2d_correct+0x52c>
add    $0x168,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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