int Extend(filt, x_dim, y_dim, x_pos, y_pos, result, f_or_e)
register double *filt, *result;
register int x_dim;
int y_dim, x_pos, y_pos, f_or_e;
{
  int filt_sz = x_dim * y_dim;
  register int x_start = 0, y_start = 0, x_stop = x_dim, y_stop = filt_sz;
  register int y_filt, x_filt, y_edge, x_edge;
  register int x_base = (x_pos > 0) ? (x_dim - 1) : 0;
  register int y_base = (y_pos > 0) ? (x_dim * (y_dim - 1)) : 0;
  int x_edge_dist =
      (x_pos > 0) ? (x_pos - x_dim) : ((x_pos < -1) ? (x_pos + 1) : 0);
  int y_edge_dist = x_dim * ((y_pos > 0) ? (y_pos - y_dim)
                                         : ((y_pos < -1) ? (y_pos + 1) : 0));
  int i;
  int mx_pos = (x_dim / 2) + 1;
  int my_pos = (y_dim / 2) + 1;
  for (i = 0; i < filt_sz; i++)
    result[i] = 0.0;
  if (f_or_e == 1) {
    if (x_pos == mx_pos)
      x_stop = (x_dim + 1) / 2;
    else if (x_pos == -mx_pos) {
      x_start = x_dim / 2;
      x_edge_dist = 0;
    }
    if (y_pos == my_pos)
      y_stop = x_dim * ((y_dim + 1) / 2);
    else if (y_pos == -my_pos) {
      y_start = x_dim * (y_dim / 2);
      y_edge_dist = 0;
    }
  }
  for (y_filt = y_start, y_edge = y_edge_dist; y_filt < y_stop;
       y_filt += x_dim, y_edge += x_dim)
    for (x_filt = y_filt + x_start, x_edge = x_edge_dist;
         x_filt < y_filt + x_stop; x_filt++, x_edge++)
      if (((!y_base &&
            ((((y_edge) > 0) ? 1 : (((y_edge) < 0) ? -1 : 0)) == -1)) ||
           (y_base &&
            ((((y_edge) > 0) ? 1 : (((y_edge) < 0) ? -1 : 0)) == 1))) !=
          ((!x_base &&
            ((((x_edge) > 0) ? 1 : (((x_edge) < 0) ? -1 : 0)) == -1)) ||
           (x_base &&
            ((((x_edge) > 0) ? 1 : (((x_edge) < 0) ? -1 : 0)) == 1)))) {
        result[(((y_base - (((y_edge) >= 0) ? (y_edge) : (-(y_edge)))) >= 0)
                    ? (y_base - (((y_edge) >= 0) ? (y_edge) : (-(y_edge))))
                    : (-(y_base -
                         (((y_edge) >= 0) ? (y_edge) : (-(y_edge)))))) +
               (((x_base - (((x_edge) >= 0) ? (x_edge) : (-(x_edge)))) >= 0)
                    ? (x_base - (((x_edge) >= 0) ? (x_edge) : (-(x_edge))))
                    : (-(x_base -
                         (((x_edge) >= 0) ? (x_edge) : (-(x_edge))))))] -=
            filt[x_filt];
        result[(((y_base + y_edge) < (0))
                    ? (0)
                    : (((y_base + y_edge) >= (y_dim)) ? (y_dim - 1)
                                                      : (y_base + y_edge))) +
               (((x_base + x_edge) < (0))
                    ? (0)
                    : (((x_base + x_edge) >= (x_dim)) ? (x_dim - 1)
                                                      : (x_base + x_edge)))] +=
            filt[x_filt] + filt[x_filt];
      } else
        result[(((y_base - (((y_edge) >= 0) ? (y_edge) : (-(y_edge)))) >= 0)
                    ? (y_base - (((y_edge) >= 0) ? (y_edge) : (-(y_edge))))
                    : (-(y_base -
                         (((y_edge) >= 0) ? (y_edge) : (-(y_edge)))))) +
               (((x_base - (((x_edge) >= 0) ? (x_edge) : (-(x_edge)))) >= 0)
                    ? (x_base - (((x_edge) >= 0) ? (x_edge) : (-(x_edge))))
                    : (-(x_base -
                         (((x_edge) >= 0) ? (x_edge) : (-(x_edge))))))] +=
            filt[x_filt];
  return (0);
}


x86 assembly (clang, -O0):
<Extend>:
push   %rbp
mov    %rsp,%rbp
sub    $0xe0,%rsp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    -0xc(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x14(%rbp)
jle    115f <Extend+0x5f>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1169 <Extend+0x69>
xor    %eax,%eax
mov    %eax,-0x64(%rbp)
jmp    1169 <Extend+0x69>
mov    -0x64(%rbp),%eax
mov    %eax,-0x48(%rbp)
cmpl   $0x0,-0x18(%rbp)
jle    118d <Extend+0x8d>
mov    -0xc(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x68(%rbp)
jmp    1197 <Extend+0x97>
xor    %eax,%eax
mov    %eax,-0x68(%rbp)
jmp    1197 <Extend+0x97>
mov    -0x68(%rbp),%eax
mov    %eax,-0x4c(%rbp)
cmpl   $0x0,-0x14(%rbp)
jle    11b5 <Extend+0xb5>
mov    -0x14(%rbp),%eax
sub    -0xc(%rbp),%eax
mov    %eax,-0x6c(%rbp)
jmp    11dd <Extend+0xdd>
cmpl   $0xffffffff,-0x14(%rbp)
jge    11cd <Extend+0xcd>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    11d7 <Extend+0xd7>
xor    %eax,%eax
mov    %eax,-0x70(%rbp)
jmp    11d7 <Extend+0xd7>
mov    -0x70(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x74(%rbp)
cmpl   $0x0,-0x18(%rbp)
jle    1201 <Extend+0x101>
mov    -0x18(%rbp),%eax
sub    -0x10(%rbp),%eax
mov    %eax,-0x78(%rbp)
jmp    1229 <Extend+0x129>
cmpl   $0xffffffff,-0x18(%rbp)
jge    1219 <Extend+0x119>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    1223 <Extend+0x123>
xor    %eax,%eax
mov    %eax,-0x7c(%rbp)
jmp    1223 <Extend+0x123>
mov    -0x7c(%rbp),%eax
mov    %eax,-0x78(%rbp)
mov    -0x74(%rbp),%eax
mov    -0x78(%rbp),%ecx
imul   %ecx,%eax
mov    %eax,-0x54(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x10(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1288 <Extend+0x188>
mov    -0x20(%rbp),%rax
movslq -0x58(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    125e <Extend+0x15e>
cmpl   $0x1,0x10(%rbp)
jne    1351 <Extend+0x251>
mov    -0x14(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jne    12b4 <Extend+0x1b4>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x30(%rbp)
jmp    12de <Extend+0x1de>
mov    -0x14(%rbp),%eax
xor    %ecx,%ecx
sub    -0x5c(%rbp),%ecx
cmp    %ecx,%eax
jne    12d9 <Extend+0x1d9>
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x50(%rbp)
jmp    12de <Extend+0x1de>
mov    -0x18(%rbp),%eax
cmp    -0x60(%rbp),%eax
jne    130e <Extend+0x20e>
mov    -0xc(%rbp),%eax
mov    %eax,-0x80(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x80(%rbp),%eax
imul   %ecx,%eax
mov    %eax,-0x34(%rbp)
jmp    134c <Extend+0x24c>
mov    -0x18(%rbp),%eax
xor    %ecx,%ecx
sub    -0x60(%rbp),%ecx
cmp    %ecx,%eax
jne    1347 <Extend+0x247>
mov    -0xc(%rbp),%eax
mov    %eax,-0x84(%rbp)
mov    -0x10(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x84(%rbp),%eax
imul   %ecx,%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x54(%rbp)
jmp    134c <Extend+0x24c>
jmp    1351 <Extend+0x251>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1a3e <Extend+0x93e>
mov    -0x38(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
cmp    %ecx,%eax
jge    1a22 <Extend+0x922>
cmpl   $0x0,-0x4c(%rbp)
jne    13da <Extend+0x2da>
cmpl   $0x0,-0x40(%rbp)
jle    13ad <Extend+0x2ad>
mov    $0x1,%eax
mov    %eax,-0x88(%rbp)
jmp    13c3 <Extend+0x2c3>
mov    -0x40(%rbp),%edx
xor    %eax,%eax
mov    $0xffffffff,%ecx
cmp    $0x0,%edx
cmovl  %ecx,%eax
mov    %eax,-0x88(%rbp)
mov    -0x88(%rbp),%ecx
mov    $0x1,%al
cmp    $0xffffffff,%ecx
mov    %al,-0x89(%rbp)
je     143a <Extend+0x33a>
xor    %eax,%eax
cmpl   $0x0,-0x4c(%rbp)
mov    %al,-0x8a(%rbp)
je     142e <Extend+0x32e>
cmpl   $0x0,-0x40(%rbp)
jle    1406 <Extend+0x306>
mov    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    141c <Extend+0x31c>
mov    -0x40(%rbp),%edx
xor    %eax,%eax
mov    $0xffffffff,%ecx
cmp    $0x0,%edx
cmovl  %ecx,%eax
mov    %eax,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    $0x1,%eax
sete   %al
mov    %al,-0x8a(%rbp)
mov    -0x8a(%rbp),%al
mov    %al,-0x89(%rbp)
mov    -0x89(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x94(%rbp)
cmpl   $0x0,-0x48(%rbp)
jne    149c <Extend+0x39c>
cmpl   $0x0,-0x44(%rbp)
jle    146f <Extend+0x36f>
mov    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    1485 <Extend+0x385>
mov    -0x44(%rbp),%edx
xor    %eax,%eax
mov    $0xffffffff,%ecx
cmp    $0x0,%edx
cmovl  %ecx,%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%ecx
mov    $0x1,%al
cmp    $0xffffffff,%ecx
mov    %al,-0x99(%rbp)
je     14fc <Extend+0x3fc>
xor    %eax,%eax
cmpl   $0x0,-0x48(%rbp)
mov    %al,-0x9a(%rbp)
je     14f0 <Extend+0x3f0>
cmpl   $0x0,-0x44(%rbp)
jle    14c8 <Extend+0x3c8>
mov    $0x1,%eax
mov    %eax,-0xa0(%rbp)
jmp    14de <Extend+0x3de>
mov    -0x44(%rbp),%edx
xor    %eax,%eax
mov    $0xffffffff,%ecx
cmp    $0x0,%edx
cmovl  %ecx,%eax
mov    %eax,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
cmp    $0x1,%eax
sete   %al
mov    %al,-0x9a(%rbp)
mov    -0x9a(%rbp),%al
mov    %al,-0x99(%rbp)
mov    -0x94(%rbp),%eax
mov    -0x99(%rbp),%cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
je     1818 <Extend+0x718>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x20(%rbp),%rax
mov    %rax,-0xb0(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0xa4(%rbp)
cmpl   $0x0,-0x40(%rbp)
jl     1557 <Extend+0x457>
mov    -0x40(%rbp),%eax
mov    %eax,-0xbc(%rbp)
jmp    1562 <Extend+0x462>
xor    %eax,%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0xbc(%rbp)
mov    -0xa4(%rbp),%eax
mov    -0xbc(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jl     15be <Extend+0x4be>
mov    -0x4c(%rbp),%eax
mov    %eax,-0xc0(%rbp)
cmpl   $0x0,-0x40(%rbp)
jl     159a <Extend+0x49a>
mov    -0x40(%rbp),%eax
mov    %eax,-0xc4(%rbp)
jmp    15a5 <Extend+0x4a5>
xor    %eax,%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0xc4(%rbp)
mov    -0xc0(%rbp),%eax
mov    -0xc4(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0xc8(%rbp)
jmp    1602 <Extend+0x502>
mov    -0x4c(%rbp),%eax
mov    %eax,-0xcc(%rbp)
cmpl   $0x0,-0x40(%rbp)
jl     15df <Extend+0x4df>
mov    -0x40(%rbp),%eax
mov    %eax,-0xd0(%rbp)
jmp    15ea <Extend+0x4ea>
xor    %eax,%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0xd0(%rbp)
mov    -0xcc(%rbp),%ecx
mov    -0xd0(%rbp),%eax
sub    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0xc8(%rbp)
mov    -0xc8(%rbp),%eax
mov    %eax,-0xd8(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0xd4(%rbp)
cmpl   $0x0,-0x44(%rbp)
jl     162f <Extend+0x52f>
mov    -0x44(%rbp),%eax
mov    %eax,-0xdc(%rbp)
jmp    163a <Extend+0x53a>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0xdc(%rbp)
mov    -0xd4(%rbp),%eax
mov    -0xdc(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jl     1696 <Extend+0x596>
mov    -0x48(%rbp),%eax
mov    %eax,-0xe0(%rbp)
cmpl   $0x0,-0x44(%rbp)
jl     1672 <Extend+0x572>
mov    -0x44(%rbp),%eax
mov    %eax,-0xe4(%rbp)
jmp    167d <Extend+0x57d>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0xe4(%rbp)
mov    -0xe0(%rbp),%eax
mov    -0xe4(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0xe8(%rbp)
jmp    16da <Extend+0x5da>
mov    -0x48(%rbp),%eax
mov    %eax,-0xec(%rbp)
cmpl   $0x0,-0x44(%rbp)
jl     16b7 <Extend+0x5b7>
mov    -0x44(%rbp),%eax
mov    %eax,-0xf0(%rbp)
jmp    16c2 <Extend+0x5c2>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0xf0(%rbp)
mov    -0xec(%rbp),%ecx
mov    -0xf0(%rbp),%eax
sub    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0xe8(%rbp)
mov    -0xb0(%rbp),%rax
movsd  -0xb8(%rbp),%xmm1
ff
mov    -0xd8(%rbp),%ecx
mov    -0xe8(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
mov    -0x20(%rbp),%rax
mov    %rax,-0xf8(%rbp)
mov    -0x4c(%rbp),%eax
add    -0x40(%rbp),%eax
cmp    $0x0,%eax
jge    1751 <Extend+0x651>
xor    %eax,%eax
mov    %eax,-0x104(%rbp)
jmp    1789 <Extend+0x689>
mov    -0x4c(%rbp),%eax
add    -0x40(%rbp),%eax
cmp    -0x10(%rbp),%eax
jl     1771 <Extend+0x671>
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x108(%rbp)
jmp    177d <Extend+0x67d>
mov    -0x4c(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x108(%rbp)
mov    -0x108(%rbp),%eax
mov    %eax,-0x104(%rbp)
mov    -0x104(%rbp),%eax
mov    %eax,-0x10c(%rbp)
mov    -0x48(%rbp),%eax
add    -0x44(%rbp),%eax
cmp    $0x0,%eax
jge    17b1 <Extend+0x6b1>
xor    %eax,%eax
mov    %eax,-0x110(%rbp)
jmp    17e9 <Extend+0x6e9>
mov    -0x48(%rbp),%eax
add    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     17d1 <Extend+0x6d1>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x114(%rbp)
jmp    17dd <Extend+0x6dd>
mov    -0x48(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x114(%rbp)
mov    -0x114(%rbp),%eax
mov    %eax,-0x110(%rbp)
mov    -0xf8(%rbp),%rax
movsd  -0x100(%rbp),%xmm0
ff
mov    -0x10c(%rbp),%ecx
mov    -0x110(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1a06 <Extend+0x906>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
mov    -0x20(%rbp),%rax
mov    %rax,-0x120(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x118(%rbp)
cmpl   $0x0,-0x40(%rbp)
jl     1859 <Extend+0x759>
mov    -0x40(%rbp),%eax
mov    %eax,-0x12c(%rbp)
jmp    1864 <Extend+0x764>
xor    %eax,%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0x12c(%rbp)
mov    -0x118(%rbp),%eax
mov    -0x12c(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jl     18c0 <Extend+0x7c0>
mov    -0x4c(%rbp),%eax
mov    %eax,-0x130(%rbp)
cmpl   $0x0,-0x40(%rbp)
jl     189c <Extend+0x79c>
mov    -0x40(%rbp),%eax
mov    %eax,-0x134(%rbp)
jmp    18a7 <Extend+0x7a7>
xor    %eax,%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0x134(%rbp)
mov    -0x130(%rbp),%eax
mov    -0x134(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x138(%rbp)
jmp    1904 <Extend+0x804>
mov    -0x4c(%rbp),%eax
mov    %eax,-0x13c(%rbp)
cmpl   $0x0,-0x40(%rbp)
jl     18e1 <Extend+0x7e1>
mov    -0x40(%rbp),%eax
mov    %eax,-0x140(%rbp)
jmp    18ec <Extend+0x7ec>
xor    %eax,%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0x140(%rbp)
mov    -0x13c(%rbp),%ecx
mov    -0x140(%rbp),%eax
sub    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x138(%rbp)
mov    -0x138(%rbp),%eax
mov    %eax,-0x148(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x144(%rbp)
cmpl   $0x0,-0x44(%rbp)
jl     1931 <Extend+0x831>
mov    -0x44(%rbp),%eax
mov    %eax,-0x14c(%rbp)
jmp    193c <Extend+0x83c>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0x14c(%rbp)
mov    -0x144(%rbp),%eax
mov    -0x14c(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jl     1998 <Extend+0x898>
mov    -0x48(%rbp),%eax
mov    %eax,-0x150(%rbp)
cmpl   $0x0,-0x44(%rbp)
jl     1974 <Extend+0x874>
mov    -0x44(%rbp),%eax
mov    %eax,-0x154(%rbp)
jmp    197f <Extend+0x87f>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0x154(%rbp)
mov    -0x150(%rbp),%eax
mov    -0x154(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x158(%rbp)
jmp    19dc <Extend+0x8dc>
mov    -0x48(%rbp),%eax
mov    %eax,-0x15c(%rbp)
cmpl   $0x0,-0x44(%rbp)
jl     19b9 <Extend+0x8b9>
mov    -0x44(%rbp),%eax
mov    %eax,-0x160(%rbp)
jmp    19c4 <Extend+0x8c4>
xor    %eax,%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0x160(%rbp)
mov    -0x15c(%rbp),%ecx
mov    -0x160(%rbp),%eax
sub    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x158(%rbp)
mov    -0x120(%rbp),%rax
movsd  -0x128(%rbp),%xmm0
ff
mov    -0x148(%rbp),%ecx
mov    -0x158(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1a0b <Extend+0x90b>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1378 <Extend+0x278>
jmp    1a27 <Extend+0x927>
mov    -0xc(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
jmp    135d <Extend+0x25d>
xor    %eax,%eax
add    $0xe0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Extend>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %ecx,%r13d
mov    %edx,%r14d
imul   %esi,%r14d
lea    -0x1(%rsi),%r10d
lea    -0x1(%rdx),%eax
mov    %eax,0x30(%rsp)
imul   %esi,%eax
mov    %ecx,%ebx
mov    %rsi,0x28(%rsp)
sub    %esi,%ebx
lea    0x1(%r13),%esi
xor    %r11d,%r11d
cmp    $0xffffffff,%ecx
cmovge %r11d,%esi
test   %ecx,%ecx
mov    $0x0,%ebp
mov    %r10d,0x34(%rsp)
cmovg  %r10d,%ebp
mov    %r8d,%r15d
cmovg  %ebx,%esi
mov    %esi,0x14(%rsp)
mov    %r8d,%esi
mov    %rdx,0x20(%rsp)
sub    %edx,%esi
lea    0x1(%r15),%edx
cmp    $0xffffffff,%r8d
cmovge %r11d,%edx
mov    %r9,%rbx
mov    %rdi,%r12
test   %r8d,%r8d
cmovg  %esi,%edx
mov    %rdx,0x8(%rsp)
mov    0x90(%rsp),%esi
cmovle %r11d,%eax
mov    %rax,0x40(%rsp)
test   %r14d,%r14d
jle    11be <Extend+0xae>
mov    %r14d,%edx
shl    $0x3,%rdx
mov    %rbx,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x90(%rsp),%esi
xor    %r11d,%r11d
mov    0x28(%rsp),%rdx
mov    0x8(%rsp),%rax
imul   %edx,%eax
mov    %rax,0x8(%rsp)
cmp    $0x1,%esi
jne    1216 <Extend+0x106>
mov    %rdx,%r9
shr    $0x1f,%edx
mov    0x20(%rsp),%rsi
mov    %esi,%eax
shr    $0x1f,%eax
add    %r9d,%edx
add    %esi,%eax
sar    %edx
sar    %eax
lea    0x1(%rdx),%esi
lea    0x1(%rax),%r8d
cmp    %r13d,%esi
jne    1231 <Extend+0x121>
lea    0x1(%r9),%edx
shr    $0x1f,%edx
add    %r9d,%edx
add    $0x1,%edx
sar    %edx
mov    %edx,0x4(%rsp)
xor    %edx,%edx
mov    %rdx,0x18(%rsp)
jmp    1255 <Extend+0x145>
mov    %edx,%eax
mov    %edx,0x4(%rsp)
xor    %eax,%eax
mov    %rax,0x18(%rsp)
cmp    %r14d,%r11d
jl     12ae <Extend+0x19e>
jmp    1409 <Extend+0x2f9>
mov    %edx,%esi
not    %esi
xor    %edi,%edi
cmp    %r13d,%esi
mov    0x14(%rsp),%esi
cmove  %edi,%esi
mov    %esi,0x14(%rsp)
cmove  %edx,%edi
mov    %rdi,0x18(%rsp)
mov    %r9d,%edx
mov    %r9d,0x4(%rsp)
mov    %r9,%rdx
cmp    %r15d,%r8d
jne    1281 <Extend+0x171>
mov    0x20(%rsp),%rcx
lea    0x1(%rcx),%eax
shr    $0x1f,%eax
lea    (%rcx,%rax,1),%r14d
add    $0x1,%r14d
sar    %r14d
imul   %edx,%r14d
cmp    %r14d,%r11d
jl     12ae <Extend+0x19e>
jmp    1409 <Extend+0x2f9>
mov    %eax,%ecx
not    %ecx
imul   %edx,%eax
xor    %edx,%edx
cmp    %r15d,%ecx
mov    0x8(%rsp),%rcx
cmove  %edx,%ecx
mov    %rcx,0x8(%rsp)
cmovne %edx,%eax
mov    0x28(%rsp),%rdx
mov    %eax,%r11d
cmp    %r14d,%r11d
jge    1409 <Extend+0x2f9>
mov    0x18(%rsp),%rax
lea    (%rax,%r11,1),%r15d
movslq %r11d,%rsi
movslq %edx,%rax
mov    %rax,0x50(%rsp)
movslq 0x4(%rsp),%rax
mov    %rax,0x38(%rsp)
movslq %r14d,%rax
mov    %rax,0x48(%rsp)
jmp    12ff <Extend+0x1ef>
cs nopw 0x0(%rax,%rax,1)
add    0x50(%rsp),%rsi
mov    0x8(%rsp),%rax
add    %edx,%eax
mov    %rax,0x8(%rsp)
add    %edx,%r15d
cmp    0x48(%rsp),%rsi
jge    1409 <Extend+0x2f9>
mov    0x18(%rsp),%rax
cmp    0x4(%rsp),%eax
jge    12e0 <Extend+0x1d0>
xor    %eax,%eax
xor    %edi,%edi
mov    0x8(%rsp),%rcx
test   %ecx,%ecx
sets   %al
setg   %dil
mov    0x40(%rsp),%r10
test   %r10d,%r10d
cmove  %eax,%edi
mov    0x38(%rsp),%rax
lea    (%rsi,%rax,1),%r8
mov    %ecx,%eax
neg    %eax
cmovns %ecx,%eax
add    %r10d,%eax
mov    %eax,%r9d
neg    %r9d
cmovs  %eax,%r9d
lea    (%rcx,%r10,1),%eax
cmp    0x20(%rsp),%eax
mov    0x30(%rsp),%r10d
cmovl  %eax,%r10d
movslq %r15d,%r11
test   %eax,%eax
mov    $0x0,%eax
cmovs  %eax,%r10d
mov    0x14(%rsp),%eax
jmp    138a <Extend+0x27a>
nopl   0x0(%rax)
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx,%r14,8)
add    $0x1,%r11
add    $0x1,%eax
cmp    %r8,%r11
jge    12e0 <Extend+0x1d0>
xor    %ecx,%ecx
xor    %r13d,%r13d
test   %eax,%eax
sets   %cl
setg   %r13b
test   %ebp,%ebp
cmove  %ecx,%r13d
mov    %eax,%ecx
neg    %ecx
cmovns %eax,%ecx
add    %ebp,%ecx
mov    %ecx,%r14d
neg    %r14d
cmovs  %ecx,%r14d
add    %r9d,%r14d
movsd  (%r12,%r11,8),%xmm0
movsd  (%rbx,%r14,8),%xmm1
cmp    %r13b,%dil
je     1370 <Extend+0x260>
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%r14,8)
movsd  (%r12,%r11,8),%xmm0
addsd  %xmm0,%xmm0
lea    (%rax,%rbp,1),%ecx
cmp    %edx,%ecx
mov    0x34(%rsp),%edx
cmovl  %ecx,%edx
test   %ecx,%ecx
mov    $0x0,%ecx
cmovs  %ecx,%edx
add    %r10d,%edx
movslq %edx,%rcx
mov    0x28(%rsp),%rdx
addsd  (%rbx,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
jmp    137a <Extend+0x26a>
xor    %eax,%eax
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @Extend(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, double* %5, i32 %6) {
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
store double* %0, double** %8
store i32 %1, i32* %9
store i32 %2, i32* %10
store i32 %3, i32* %11
store i32 %4, i32* %12
store double* %5, double** %13
store i32 %6, i32* %14
%31 = load i32, i32* %9
%32 = load i32, i32* %10
%33 = mul nsw i32 %31, %32
store i32 %33, i32* %15
store i32 0, i32* %16
store i32 0, i32* %17
%34 = load i32, i32* %9
store i32 %34, i32* %18
%35 = load i32, i32* %15
store i32 %35, i32* %19
%36 = load i32, i32* %11
%37 = icmp sgt i32 %36, 0
br i1 %37, label %38, label %41
38: ; preds = %7
%39 = load i32, i32* %9
%40 = sub nsw i32 %39, 1
br label %42
41: ; preds = %7
br label %42
42: ; preds = %41, %38
%43 = phi i32 [ %40, %38 ], [ 0, %41 ]
store i32 %43, i32* %24
%44 = load i32, i32* %12
%45 = icmp sgt i32 %44, 0
br i1 %45, label %46, label %51
46: ; preds = %42
%47 = load i32, i32* %9
%48 = load i32, i32* %10
%49 = sub nsw i32 %48, 1
%50 = mul nsw i32 %47, %49
br label %52
51: ; preds = %42
br label %52
52: ; preds = %51, %46
%53 = phi i32 [ %50, %46 ], [ 0, %51 ]
store i32 %53, i32* %25
%54 = load i32, i32* %11
%55 = icmp sgt i32 %54, 0
br i1 %55, label %56, label %60
56: ; preds = %52
%57 = load i32, i32* %11
%58 = load i32, i32* %9
%59 = sub nsw i32 %57, %58
br label %69
60: ; preds = %52
%61 = load i32, i32* %11
%62 = icmp slt i32 %61, -1
br i1 %62, label %63, label %66
63: ; preds = %60
%64 = load i32, i32* %11
%65 = add nsw i32 %64, 1
br label %67
66: ; preds = %60
br label %67
67: ; preds = %66, %63
%68 = phi i32 [ %65, %63 ], [ 0, %66 ]
br label %69
69: ; preds = %67, %56
%70 = phi i32 [ %59, %56 ], [ %68, %67 ]
store i32 %70, i32* %26
%71 = load i32, i32* %9
%72 = load i32, i32* %12
%73 = icmp sgt i32 %72, 0
br i1 %73, label %74, label %78
74: ; preds = %69
%75 = load i32, i32* %12
%76 = load i32, i32* %10
%77 = sub nsw i32 %75, %76
br label %87
78: ; preds = %69
%79 = load i32, i32* %12
%80 = icmp slt i32 %79, -1
br i1 %80, label %81, label %84
81: ; preds = %78
%82 = load i32, i32* %12
%83 = add nsw i32 %82, 1
br label %85
84: ; preds = %78
br label %85
85: ; preds = %84, %81
%86 = phi i32 [ %83, %81 ], [ 0, %84 ]
br label %87
87: ; preds = %85, %74
%88 = phi i32 [ %77, %74 ], [ %86, %85 ]
%89 = mul nsw i32 %71, %88
store i32 %89, i32* %27
%90 = load i32, i32* %9
%91 = sdiv i32 %90, 2
%92 = add nsw i32 %91, 1
store i32 %92, i32* %29
%93 = load i32, i32* %10
%94 = sdiv i32 %93, 2
%95 = add nsw i32 %94, 1
store i32 %95, i32* %30
store i32 0, i32* %28
br label %96
96: ; preds = %105, %87
%97 = load i32, i32* %28
%98 = load i32, i32* %15
%99 = icmp slt i32 %97, %98
br i1 %99, label %100, label %108
100: ; preds = %96
%101 = load double*, double** %13
%102 = load i32, i32* %28
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %101, i64 %103
store double 0.000000e+00, double* %104
br label %105
105: ; preds = %100
%106 = load i32, i32* %28
%107 = add nsw i32 %106, 1
store i32 %107, i32* %28
br label %96
108: ; preds = %96
%109 = load i32, i32* %14
%110 = icmp eq i32 %109, 1
br i1 %110, label %111, label %150
111: ; preds = %108
%112 = load i32, i32* %11
%113 = load i32, i32* %29
%114 = icmp eq i32 %112, %113
br i1 %114, label %115, label %119
115: ; preds = %111
%116 = load i32, i32* %9
%117 = add nsw i32 %116, 1
%118 = sdiv i32 %117, 2
store i32 %118, i32* %18
br label %128
119: ; preds = %111
%120 = load i32, i32* %11
%121 = load i32, i32* %29
%122 = sub nsw i32 0, %121
%123 = icmp eq i32 %120, %122
br i1 %123, label %124, label %127
124: ; preds = %119
%125 = load i32, i32* %9
%126 = sdiv i32 %125, 2
store i32 %126, i32* %16
store i32 0, i32* %26
br label %127
127: ; preds = %124, %119
br label %128
128: ; preds = %127, %115
%129 = load i32, i32* %12
%130 = load i32, i32* %30
%131 = icmp eq i32 %129, %130
br i1 %131, label %132, label %138
132: ; preds = %128
%133 = load i32, i32* %9
%134 = load i32, i32* %10
%135 = add nsw i32 %134, 1
%136 = sdiv i32 %135, 2
%137 = mul nsw i32 %133, %136
store i32 %137, i32* %19
br label %149
138: ; preds = %128
%139 = load i32, i32* %12
%140 = load i32, i32* %30
%141 = sub nsw i32 0, %140
%142 = icmp eq i32 %139, %141
br i1 %142, label %143, label %148
143: ; preds = %138
%144 = load i32, i32* %9
%145 = load i32, i32* %10
%146 = sdiv i32 %145, 2
%147 = mul nsw i32 %144, %146
store i32 %147, i32* %17
store i32 0, i32* %27
br label %148
148: ; preds = %143, %138
br label %149
149: ; preds = %148, %132
br label %150
150: ; preds = %149, %108
%151 = load i32, i32* %17
store i32 %151, i32* %20
%152 = load i32, i32* %27
store i32 %152, i32* %22
br label %153
153: ; preds = %486, %150
%154 = load i32, i32* %20
%155 = load i32, i32* %19
%156 = icmp slt i32 %154, %155
br i1 %156, label %157, label %493
157: ; preds = %153
%158 = load i32, i32* %20
%159 = load i32, i32* %16
%160 = add nsw i32 %158, %159
store i32 %160, i32* %21
%161 = load i32, i32* %26
store i32 %161, i32* %23
br label %162
162: ; preds = %480, %157
%163 = load i32, i32* %21
%164 = load i32, i32* %20
%165 = load i32, i32* %18
%166 = add nsw i32 %164, %165
%167 = icmp slt i32 %163, %166
br i1 %167, label %168, label %485
168: ; preds = %162
%169 = load i32, i32* %25
%170 = icmp ne i32 %169, 0
br i1 %170, label %183, label %171
171: ; preds = %168
%172 = load i32, i32* %22
%173 = icmp sgt i32 %172, 0
br i1 %173, label %174, label %175
174: ; preds = %171
br label %180
175: ; preds = %171
%176 = load i32, i32* %22
%177 = icmp slt i32 %176, 0
%178 = zext i1 %177 to i64
%179 = select i1 %177, i32 -1, i32 0
br label %180
180: ; preds = %175, %174
%181 = phi i32 [ 1, %174 ], [ %179, %175 ]
%182 = icmp eq i32 %181, -1
br i1 %182, label %200, label %183
183: ; preds = %180, %168
%184 = load i32, i32* %25
%185 = icmp ne i32 %184, 0
br i1 %185, label %186, label %198
186: ; preds = %183
%187 = load i32, i32* %22
%188 = icmp sgt i32 %187, 0
br i1 %188, label %189, label %190
189: ; preds = %186
br label %195
190: ; preds = %186
%191 = load i32, i32* %22
%192 = icmp slt i32 %191, 0
%193 = zext i1 %192 to i64
%194 = select i1 %192, i32 -1, i32 0
br label %195
195: ; preds = %190, %189
%196 = phi i32 [ 1, %189 ], [ %194, %190 ]
%197 = icmp eq i32 %196, 1
br label %198
198: ; preds = %195, %183
%199 = phi i1 [ false, %183 ], [ %197, %195 ]
br label %200
200: ; preds = %198, %180
%201 = phi i1 [ true, %180 ], [ %199, %198 ]
%202 = zext i1 %201 to i32
%203 = load i32, i32* %24
%204 = icmp ne i32 %203, 0
br i1 %204, label %217, label %205
205: ; preds = %200
%206 = load i32, i32* %23
%207 = icmp sgt i32 %206, 0
br i1 %207, label %208, label %209
208: ; preds = %205
br label %214
209: ; preds = %205
%210 = load i32, i32* %23
%211 = icmp slt i32 %210, 0
%212 = zext i1 %211 to i64
%213 = select i1 %211, i32 -1, i32 0
br label %214
214: ; preds = %209, %208
%215 = phi i32 [ 1, %208 ], [ %213, %209 ]
%216 = icmp eq i32 %215, -1
br i1 %216, label %234, label %217
217: ; preds = %214, %200
%218 = load i32, i32* %24
%219 = icmp ne i32 %218, 0
br i1 %219, label %220, label %232
220: ; preds = %217
%221 = load i32, i32* %23
%222 = icmp sgt i32 %221, 0
br i1 %222, label %223, label %224
223: ; preds = %220
br label %229
224: ; preds = %220
%225 = load i32, i32* %23
%226 = icmp slt i32 %225, 0
%227 = zext i1 %226 to i64
%228 = select i1 %226, i32 -1, i32 0
br label %229
229: ; preds = %224, %223
%230 = phi i32 [ 1, %223 ], [ %228, %224 ]
%231 = icmp eq i32 %230, 1
br label %232
232: ; preds = %229, %217
%233 = phi i1 [ false, %217 ], [ %231, %229 ]
br label %234
234: ; preds = %232, %214
%235 = phi i1 [ true, %214 ], [ %233, %232 ]
%236 = zext i1 %235 to i32
%237 = icmp ne i32 %202, %236
br i1 %237, label %238, label %389
238: ; preds = %234
%239 = load double*, double** %8
%240 = load i32, i32* %21
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds double, double* %239, i64 %241
%243 = load double, double* %242
%244 = load double*, double** %13
%245 = load i32, i32* %25
%246 = load i32, i32* %22
%247 = icmp sge i32 %246, 0
br i1 %247, label %248, label %250
248: ; preds = %238
%249 = load i32, i32* %22
br label %253
250: ; preds = %238
%251 = load i32, i32* %22
%252 = sub nsw i32 0, %251
br label %253
253: ; preds = %250, %248
%254 = phi i32 [ %249, %248 ], [ %252, %250 ]
%255 = sub nsw i32 %245, %254
%256 = icmp sge i32 %255, 0
br i1 %256, label %257, label %269
257: ; preds = %253
%258 = load i32, i32* %25
%259 = load i32, i32* %22
%260 = icmp sge i32 %259, 0
br i1 %260, label %261, label %263
261: ; preds = %257
%262 = load i32, i32* %22
br label %266
263: ; preds = %257
%264 = load i32, i32* %22
%265 = sub nsw i32 0, %264
br label %266
266: ; preds = %263, %261
%267 = phi i32 [ %262, %261 ], [ %265, %263 ]
%268 = sub nsw i32 %258, %267
br label %282
269: ; preds = %253
%270 = load i32, i32* %25
%271 = load i32, i32* %22
%272 = icmp sge i32 %271, 0
br i1 %272, label %273, label %275
273: ; preds = %269
%274 = load i32, i32* %22
br label %278
275: ; preds = %269
%276 = load i32, i32* %22
%277 = sub nsw i32 0, %276
br label %278
278: ; preds = %275, %273
%279 = phi i32 [ %274, %273 ], [ %277, %275 ]
%280 = sub nsw i32 %270, %279
%281 = sub nsw i32 0, %280
br label %282
282: ; preds = %278, %266
%283 = phi i32 [ %268, %266 ], [ %281, %278 ]
%284 = load i32, i32* %24
%285 = load i32, i32* %23
%286 = icmp sge i32 %285, 0
br i1 %286, label %287, label %289
287: ; preds = %282
%288 = load i32, i32* %23
br label %292
289: ; preds = %282
%290 = load i32, i32* %23
%291 = sub nsw i32 0, %290
br label %292
292: ; preds = %289, %287
%293 = phi i32 [ %288, %287 ], [ %291, %289 ]
%294 = sub nsw i32 %284, %293
%295 = icmp sge i32 %294, 0
br i1 %295, label %296, label %308
296: ; preds = %292
%297 = load i32, i32* %24
%298 = load i32, i32* %23
%299 = icmp sge i32 %298, 0
br i1 %299, label %300, label %302
300: ; preds = %296
%301 = load i32, i32* %23
br label %305
302: ; preds = %296
%303 = load i32, i32* %23
%304 = sub nsw i32 0, %303
br label %305
305: ; preds = %302, %300
%306 = phi i32 [ %301, %300 ], [ %304, %302 ]
%307 = sub nsw i32 %297, %306
br label %321
308: ; preds = %292
%309 = load i32, i32* %24
%310 = load i32, i32* %23
%311 = icmp sge i32 %310, 0
br i1 %311, label %312, label %314
312: ; preds = %308
%313 = load i32, i32* %23
br label %317
314: ; preds = %308
%315 = load i32, i32* %23
%316 = sub nsw i32 0, %315
br label %317
317: ; preds = %314, %312
%318 = phi i32 [ %313, %312 ], [ %316, %314 ]
%319 = sub nsw i32 %309, %318
%320 = sub nsw i32 0, %319
br label %321
321: ; preds = %317, %305
%322 = phi i32 [ %307, %305 ], [ %320, %317 ]
%323 = add nsw i32 %283, %322
%324 = sext i32 %323 to i64
%325 = getelementptr inbounds double, double* %244, i64 %324
%326 = load double, double* %325
%327 = fsub double %326, %243
store double %327, double* %325
%328 = load double*, double** %8
%329 = load i32, i32* %21
%330 = sext i32 %329 to i64
%331 = getelementptr inbounds double, double* %328, i64 %330
%332 = load double, double* %331
%333 = load double*, double** %8
%334 = load i32, i32* %21
%335 = sext i32 %334 to i64
%336 = getelementptr inbounds double, double* %333, i64 %335
%337 = load double, double* %336
%338 = fadd double %332, %337
%339 = load double*, double** %13
%340 = load i32, i32* %25
%341 = load i32, i32* %22
%342 = add nsw i32 %340, %341
%343 = icmp slt i32 %342, 0
br i1 %343, label %344, label %345
344: ; preds = %321
br label %360
345: ; preds = %321
%346 = load i32, i32* %25
%347 = load i32, i32* %22
%348 = add nsw i32 %346, %347
%349 = load i32, i32* %10
%350 = icmp sge i32 %348, %349
br i1 %350, label %351, label %354
351: ; preds = %345
%352 = load i32, i32* %10
%353 = sub nsw i32 %352, 1
br label %358
354: ; preds = %345
%355 = load i32, i32* %25
%356 = load i32, i32* %22
%357 = add nsw i32 %355, %356
br label %358
358: ; preds = %354, %351
%359 = phi i32 [ %353, %351 ], [ %357, %354 ]
br label %360
360: ; preds = %358, %344
%361 = phi i32 [ 0, %344 ], [ %359, %358 ]
%362 = load i32, i32* %24
%363 = load i32, i32* %23
%364 = add nsw i32 %362, %363
%365 = icmp slt i32 %364, 0
br i1 %365, label %366, label %367
366: ; preds = %360
br label %382
367: ; preds = %360
%368 = load i32, i32* %24
%369 = load i32, i32* %23
%370 = add nsw i32 %368, %369
%371 = load i32, i32* %9
%372 = icmp sge i32 %370, %371
br i1 %372, label %373, label %376
373: ; preds = %367
%374 = load i32, i32* %9
%375 = sub nsw i32 %374, 1
br label %380
376: ; preds = %367
%377 = load i32, i32* %24
%378 = load i32, i32* %23
%379 = add nsw i32 %377, %378
br label %380
380: ; preds = %376, %373
%381 = phi i32 [ %375, %373 ], [ %379, %376 ]
br label %382
382: ; preds = %380, %366
%383 = phi i32 [ 0, %366 ], [ %381, %380 ]
%384 = add nsw i32 %361, %383
%385 = sext i32 %384 to i64
%386 = getelementptr inbounds double, double* %339, i64 %385
%387 = load double, double* %386
%388 = fadd double %387, %338
store double %388, double* %386
br label %479
389: ; preds = %234
%390 = load double*, double** %8
%391 = load i32, i32* %21
%392 = sext i32 %391 to i64
%393 = getelementptr inbounds double, double* %390, i64 %392
%394 = load double, double* %393
%395 = load double*, double** %13
%396 = load i32, i32* %25
%397 = load i32, i32* %22
%398 = icmp sge i32 %397, 0
br i1 %398, label %399, label %401
399: ; preds = %389
%400 = load i32, i32* %22
br label %404
401: ; preds = %389
%402 = load i32, i32* %22
%403 = sub nsw i32 0, %402
br label %404
404: ; preds = %401, %399
%405 = phi i32 [ %400, %399 ], [ %403, %401 ]
%406 = sub nsw i32 %396, %405
%407 = icmp sge i32 %406, 0
br i1 %407, label %408, label %420
408: ; preds = %404
%409 = load i32, i32* %25
%410 = load i32, i32* %22
%411 = icmp sge i32 %410, 0
br i1 %411, label %412, label %414
412: ; preds = %408
%413 = load i32, i32* %22
br label %417
414: ; preds = %408
%415 = load i32, i32* %22
%416 = sub nsw i32 0, %415
br label %417
417: ; preds = %414, %412
%418 = phi i32 [ %413, %412 ], [ %416, %414 ]
%419 = sub nsw i32 %409, %418
br label %433
420: ; preds = %404
%421 = load i32, i32* %25
%422 = load i32, i32* %22
%423 = icmp sge i32 %422, 0
br i1 %423, label %424, label %426
424: ; preds = %420
%425 = load i32, i32* %22
br label %429
426: ; preds = %420
%427 = load i32, i32* %22
%428 = sub nsw i32 0, %427
br label %429
429: ; preds = %426, %424
%430 = phi i32 [ %425, %424 ], [ %428, %426 ]
%431 = sub nsw i32 %421, %430
%432 = sub nsw i32 0, %431
br label %433
433: ; preds = %429, %417
%434 = phi i32 [ %419, %417 ], [ %432, %429 ]
%435 = load i32, i32* %24
%436 = load i32, i32* %23
%437 = icmp sge i32 %436, 0
br i1 %437, label %438, label %440
438: ; preds = %433
%439 = load i32, i32* %23
br label %443
440: ; preds = %433
%441 = load i32, i32* %23
%442 = sub nsw i32 0, %441
br label %443
443: ; preds = %440, %438
%444 = phi i32 [ %439, %438 ], [ %442, %440 ]
%445 = sub nsw i32 %435, %444
%446 = icmp sge i32 %445, 0
br i1 %446, label %447, label %459
447: ; preds = %443
%448 = load i32, i32* %24
%449 = load i32, i32* %23
%450 = icmp sge i32 %449, 0
br i1 %450, label %451, label %453
451: ; preds = %447
%452 = load i32, i32* %23
br label %456
453: ; preds = %447
%454 = load i32, i32* %23
%455 = sub nsw i32 0, %454
br label %456
456: ; preds = %453, %451
%457 = phi i32 [ %452, %451 ], [ %455, %453 ]
%458 = sub nsw i32 %448, %457
br label %472
459: ; preds = %443
%460 = load i32, i32* %24
%461 = load i32, i32* %23
%462 = icmp sge i32 %461, 0
br i1 %462, label %463, label %465
463: ; preds = %459
%464 = load i32, i32* %23
br label %468
465: ; preds = %459
%466 = load i32, i32* %23
%467 = sub nsw i32 0, %466
br label %468
468: ; preds = %465, %463
%469 = phi i32 [ %464, %463 ], [ %467, %465 ]
%470 = sub nsw i32 %460, %469
%471 = sub nsw i32 0, %470
br label %472
472: ; preds = %468, %456
%473 = phi i32 [ %458, %456 ], [ %471, %468 ]
%474 = add nsw i32 %434, %473
%475 = sext i32 %474 to i64
%476 = getelementptr inbounds double, double* %395, i64 %475
%477 = load double, double* %476
%478 = fadd double %477, %394
store double %478, double* %476
br label %479
479: ; preds = %472, %382
br label %480
480: ; preds = %479
%481 = load i32, i32* %21
%482 = add nsw i32 %481, 1
store i32 %482, i32* %21
%483 = load i32, i32* %23
%484 = add nsw i32 %483, 1
store i32 %484, i32* %23
br label %162
485: ; preds = %162
br label %486
486: ; preds = %485
%487 = load i32, i32* %9
%488 = load i32, i32* %20
%489 = add nsw i32 %488, %487
store i32 %489, i32* %20
%490 = load i32, i32* %9
%491 = load i32, i32* %22
%492 = add nsw i32 %491, %490
store i32 %492, i32* %22
br label %153
493: ; preds = %153
ret i32 0
}

LLVM IR (-O2):
define i32 @Extend(double* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, double* nocapture %5, i32 %6) {
%8 = mul i32 %2, %1
%9 = icmp sgt i32 %3, 0
%10 = add nsw i32 %1, -1
%11 = select i1 %9, i32 %10, i32 0
%12 = icmp sgt i32 %4, 0
%13 = add nsw i32 %2, -1
%14 = mul nsw i32 %13, %1
%15 = select i1 %12, i32 %14, i32 0
%16 = sub nsw i32 %3, %1
%17 = icmp slt i32 %3, -1
%18 = add nsw i32 %3, 1
%19 = select i1 %17, i32 %18, i32 0
%20 = select i1 %9, i32 %16, i32 %19
%21 = sub nsw i32 %4, %2
%22 = icmp slt i32 %4, -1
%23 = add nsw i32 %4, 1
%24 = select i1 %22, i32 %23, i32 0
%25 = select i1 %12, i32 %21, i32 %24
%26 = mul nsw i32 %25, %1
%27 = sdiv i32 %1, 2
%28 = add nsw i32 %27, 1
%29 = sdiv i32 %2, 2
%30 = add nsw i32 %29, 1
%31 = icmp sgt i32 %8, 0
br i1 %31, label %32, label %36
32: ; preds = %7
%33 = bitcast double* %5 to i8*
%34 = zext i32 %8 to i64
%35 = shl nuw nsw i64 %34, 3
br label %36
36: ; preds = %32, %7
%37 = icmp eq i32 %6, 1
br i1 %37, label %38, label %63
38: ; preds = %36
%39 = icmp eq i32 %28, %3
br i1 %39, label %40, label %43
40: ; preds = %38
%41 = add nsw i32 %1, 1
%42 = sdiv i32 %41, 2
br label %48
43: ; preds = %38
%44 = xor i32 %27, -1
%45 = icmp eq i32 %44, %3
%46 = select i1 %45, i32 0, i32 %20
%47 = select i1 %45, i32 %27, i32 0
br label %48
48: ; preds = %43, %40
%49 = phi i32 [ %42, %40 ], [ %1, %43 ]
%50 = phi i32 [ %20, %40 ], [ %46, %43 ]
%51 = phi i32 [ 0, %40 ], [ %47, %43 ]
%52 = icmp eq i32 %30, %4
br i1 %52, label %53, label %57
53: ; preds = %48
%54 = add nsw i32 %2, 1
%55 = sdiv i32 %54, 2
%56 = mul nsw i32 %55, %1
br label %63
57: ; preds = %48
%58 = xor i32 %29, -1
%59 = icmp eq i32 %58, %4
%60 = mul nsw i32 %29, %1
%61 = select i1 %59, i32 0, i32 %26
%62 = select i1 %59, i32 %60, i32 0
br label %63
63: ; preds = %57, %53, %36
%64 = phi i32 [ %56, %53 ], [ %8, %36 ], [ %8, %57 ]
%65 = phi i32 [ %49, %53 ], [ %1, %36 ], [ %49, %57 ]
%66 = phi i32 [ %50, %53 ], [ %20, %36 ], [ %50, %57 ]
%67 = phi i32 [ %26, %53 ], [ %26, %36 ], [ %61, %57 ]
%68 = phi i32 [ 0, %53 ], [ 0, %36 ], [ %62, %57 ]
%69 = phi i32 [ %51, %53 ], [ 0, %36 ], [ %51, %57 ]
%70 = icmp slt i32 %68, %64
br i1 %70, label %71, label %139
71: ; preds = %63
%72 = icmp eq i32 %15, 0
%73 = icmp eq i32 %11, 0
%74 = add i32 %69, %68
%75 = sext i32 %68 to i64
%76 = sext i32 %1 to i64
%77 = sext i32 %65 to i64
%78 = sext i32 %64 to i64
%79 = icmp slt i32 %69, %65
br label %80
80: ; preds = %71, %134
%81 = phi i64 [ %75, %71 ], [ %135, %134 ]
%82 = phi i32 [ %74, %71 ], [ %138, %134 ]
%83 = phi i32 [ %67, %71 ], [ %136, %134 ]
%84 = add nsw i64 %81, %77
br i1 %79, label %85, label %134
85: ; preds = %80
%86 = sext i32 %82 to i64
%87 = icmp slt i32 %83, 0
%88 = icmp sgt i32 %83, 0
%89 = select i1 %72, i1 %87, i1 %88
%90 = tail call i32 @llvm.abs.i32(i32 %83, i1 true)
%91 = sub nsw i32 %15, %90
%92 = tail call i32 @llvm.abs.i32(i32 %91, i1 true)
%93 = add nsw i32 %83, %15
%94 = icmp slt i32 %93, 0
%95 = icmp slt i32 %93, %2
%96 = select i1 %95, i32 %93, i32 %13
%97 = select i1 %94, i32 0, i32 %96
br label %98
98: ; preds = %85, %130
%99 = phi i64 [ %86, %85 ], [ %131, %130 ]
%100 = phi i32 [ %66, %85 ], [ %132, %130 ]
%101 = icmp slt i32 %100, 0
%102 = icmp sgt i32 %100, 0
%103 = select i1 %73, i1 %101, i1 %102
%104 = xor i1 %89, %103
%105 = getelementptr inbounds double, double* %0, i64 %99
%106 = load double, double* %105
%107 = tail call i32 @llvm.abs.i32(i32 %100, i1 true)
%108 = sub nsw i32 %11, %107
%109 = tail call i32 @llvm.abs.i32(i32 %108, i1 true)
%110 = add nuw nsw i32 %109, %92
%111 = zext i32 %110 to i64
%112 = getelementptr inbounds double, double* %5, i64 %111
%113 = load double, double* %112
br i1 %104, label %114, label %128
114: ; preds = %98
%115 = fsub double %113, %106
store double %115, double* %112
%116 = load double, double* %105
%117 = fadd double %116, %116
%118 = add nsw i32 %100, %11
%119 = icmp slt i32 %118, 0
%120 = icmp slt i32 %118, %1
%121 = select i1 %120, i32 %118, i32 %10
%122 = select i1 %119, i32 0, i32 %121
%123 = add nsw i32 %122, %97
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %5, i64 %124
%126 = load double, double* %125
%127 = fadd double %126, %117
store double %127, double* %125
br label %130
128: ; preds = %98
%129 = fadd double %106, %113
store double %129, double* %112
br label %130
130: ; preds = %114, %128
%131 = add nsw i64 %99, 1
%132 = add nsw i32 %100, 1
%133 = icmp slt i64 %131, %84
br i1 %133, label %98, label %134
134: ; preds = %130, %80
%135 = add i64 %81, %76
%136 = add nsw i32 %83, %1
%137 = icmp slt i64 %135, %78
%138 = add i32 %82, %1
br i1 %137, label %80, label %139
139: ; preds = %134, %63
ret i32 0
}