#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 */
double ceil(double);
scalar_t__ floor(double);
void subsample_double(double *in, double *out, int dim_in[3], int dim_out[3]) {
  int i, x, y, z;
  double k111, k112, k121, k122, k211, k212, k221, k222;
  double dx1, dx2, dy1, dy2, dz1, dz2, xi, yi, zi, samp[3];
  int off1, off2, xcoord, ycoord, zcoord;
  for (i = 0; i < 3; i++) {
    if (dim_out[i] > dim_in[i])
      samp[i] = ceil((double)dim_out[i] / (double)dim_in[i]);
    else
      samp[i] = 1.0 / (ceil((double)dim_in[i] / (double)dim_out[i]));
  }
  for (z = 0; z < dim_out[2]; z++) {
    zi = 1.0 + (double)z / samp[2];
    for (y = 0; y < dim_out[1]; y++) {
      yi = 1.0 + (double)y / samp[1];
      for (x = 0; x < dim_out[0]; x++) {
        xi = 1.0 + (double)x / samp[0];
        i = z * dim_out[0] * dim_out[1] + y * dim_out[0] + x;
        if (zi >= 0 && zi < dim_in[2] && yi >= 0 && yi < dim_in[1] && xi >= 0 &&
            xi < dim_in[0]) {
          xcoord = (int)floor(xi);
          dx1 = xi - (double)xcoord;
          dx2 = 1.0 - dx1;
          ycoord = (int)floor(yi);
          dy1 = yi - (double)ycoord;
          dy2 = 1.0 - dy1;
          zcoord = (int)floor(zi);
          dz1 = zi - (double)zcoord;
          dz2 = 1.0 - dz1;
          off1 =
              xcoord - 1 + dim_in[0] * (ycoord - 1 + dim_in[1] * (zcoord - 1));
          k222 = (double)in[off1];
          k122 = (double)in[off1 + 1];
          off2 = off1 + dim_in[0];
          k212 = (double)in[off2];
          k112 = (double)in[off2 + 1];
          off1 += dim_in[0] * dim_in[1];
          k221 = (double)in[off1];
          k121 = (double)in[off1 + 1];
          off2 = off1 + dim_in[0];
          k211 = (double)in[off2];
          k111 = (double)in[off2 + 1];
          out[i] = ((((k222 * dx2 + k122 * dx1) * dy2 +
                      (k212 * dx2 + k112 * dx1) * dy1)) *
                        dz2 +
                    (((k221 * dx2 + k121 * dx1) * dy2 +
                      (k211 * dx2 + k111 * dx1) * dy1)) *
                        dz1);
        } else
          out[i] = 0;
      }
    }
  }
}


x86 assembly (clang, -O0):
<subsample_double>:
push   %rbp
mov    %rsp,%rbp
sub    $0xf0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x3,-0x24(%rbp)
jge    11ef <subsample_double+0xcf>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x18(%rbp),%rcx
movslq -0x24(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jle    119d <subsample_double+0x7d>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movslq -0x24(%rbp),%rax
movsd  %xmm0,-0xd0(%rbp,%rax,8)
ff ff
jmp    11dc <subsample_double+0xbc>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movaps %xmm0,%xmm1
movsd  0xe35(%rip),%xmm0
divsd  %xmm1,%xmm0
movslq -0x24(%rbp),%rax
movsd  %xmm0,-0xd0(%rbp,%rax,8)
ff ff
jmp    11e1 <subsample_double+0xc1>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1142 <subsample_double+0x22>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    165a <subsample_double+0x53a>
cvtsi2sdl -0x30(%rbp),%xmm1
divsd  -0xc0(%rbp),%xmm1
ff
movsd  0xde5(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    1647 <subsample_double+0x527>
cvtsi2sdl -0x2c(%rbp),%xmm1
divsd  -0xc8(%rbp),%xmm1
ff
movsd  0xdad(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    1634 <subsample_double+0x514>
cvtsi2sdl -0x28(%rbp),%xmm1
divsd  -0xd0(%rbp),%xmm1
ff
movsd  0xd76(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x30(%rbp),%eax
mov    -0x20(%rbp),%rcx
imul   (%rcx),%eax
mov    -0x20(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x2c(%rbp),%ecx
mov    -0x20(%rbp),%rdx
imul   (%rdx),%ecx
add    %ecx,%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1611 <subsample_double+0x4f1>
movsd  -0xb8(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
cvtsi2sdl 0x8(%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    1611 <subsample_double+0x4f1>
movsd  -0xb0(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1611 <subsample_double+0x4f1>
movsd  -0xb0(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
cvtsi2sdl 0x4(%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    1611 <subsample_double+0x4f1>
movsd  -0xa8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1611 <subsample_double+0x4f1>
movsd  -0xa8(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    1611 <subsample_double+0x4f1>
movsd  -0xa8(%rbp),%xmm0
ff
call   1040 <floor@plt>
mov    %eax,-0xdc(%rbp)
movsd  -0xa8(%rbp),%xmm0
ff
cvtsi2sdl -0xdc(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xc83(%rip),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0xb0(%rbp),%xmm0
ff
call   1040 <floor@plt>
mov    %eax,-0xe0(%rbp)
movsd  -0xb0(%rbp),%xmm0
ff
cvtsi2sdl -0xe0(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xc42(%rip),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
call   1040 <floor@plt>
mov    %eax,-0xe4(%rbp)
movsd  -0xb8(%rbp),%xmm0
ff
cvtsi2sdl -0xe4(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xbfb(%rip),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0xdc(%rbp),%eax
sub    $0x1,%eax
mov    -0x18(%rbp),%rcx
mov    (%rcx),%ecx
mov    -0xe0(%rbp),%edx
sub    $0x1,%edx
mov    -0x18(%rbp),%rsi
mov    0x4(%rsi),%esi
mov    -0xe4(%rbp),%edi
sub    $0x1,%edi
imul   %edi,%esi
add    %esi,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0xd4(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xd4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xd4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0xd4(%rbp),%eax
mov    -0x18(%rbp),%rcx
add    (%rcx),%eax
mov    %eax,-0xd8(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xd8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    -0x18(%rbp),%rcx
imul   0x4(%rcx),%eax
add    -0xd4(%rbp),%eax
mov    %eax,-0xd4(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xd4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xd4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0xd4(%rbp),%eax
mov    -0x18(%rbp),%rcx
add    (%rcx),%eax
mov    %eax,-0xd8(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xd8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xd8(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x60(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm4
movsd  -0x40(%rbp),%xmm3
mulsd  -0x78(%rbp),%xmm3
mulsd  %xmm4,%xmm1
addsd  %xmm3,%xmm1
mulsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0x68(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm4
movsd  -0x48(%rbp),%xmm3
mulsd  -0x78(%rbp),%xmm3
mulsd  %xmm4,%xmm1
addsd  %xmm3,%xmm1
movsd  -0x90(%rbp),%xmm4
ff
movsd  -0x58(%rbp),%xmm3
movsd  -0x80(%rbp),%xmm6
movsd  -0x38(%rbp),%xmm5
mulsd  -0x78(%rbp),%xmm5
mulsd  %xmm6,%xmm3
addsd  %xmm5,%xmm3
mulsd  -0x88(%rbp),%xmm3
ff
mulsd  %xmm4,%xmm1
addsd  %xmm3,%xmm1
mulsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1621 <subsample_double+0x501>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1626 <subsample_double+0x506>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1266 <subsample_double+0x146>
jmp    1639 <subsample_double+0x519>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    122e <subsample_double+0x10e>
jmp    164c <subsample_double+0x52c>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11f6 <subsample_double+0xd6>
add    $0xf0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<subsample_double>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rsi,0x38(%rsp)
mov    %rdi,0x48(%rsp)
mov    %rcx,0x8(%rsp)
mov    (%rcx),%ebx
mov    %rdx,(%rsp)
mov    (%rdx),%eax
cmp    %eax,%ebx
jle    1163 <subsample_double+0x43>
cvtsi2sd %ebx,%xmm0
mov    %rax,%r15
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movapd %xmm0,%xmm1
jmp    1183 <subsample_double+0x63>
mov    %rax,%r15
cvtsi2sd %eax,%xmm0
cvtsi2sd %ebx,%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movsd  0xe81(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x30(%rsp)
mov    0x8(%rsp),%rax
mov    0x4(%rax),%ebp
mov    (%rsp),%rax
mov    0x4(%rax),%eax
cmp    %eax,%ebp
jle    11be <subsample_double+0x9e>
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
mov    %eax,%r12d
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movsd  %xmm0,0x28(%rsp)
jmp    11ea <subsample_double+0xca>
mov    %eax,%r12d
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movsd  0xe20(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x28(%rsp)
mov    0x8(%rsp),%rax
mov    0x8(%rax),%r14d
mov    (%rsp),%rax
mov    0x8(%rax),%eax
cmp    %eax,%r14d
jle    121f <subsample_double+0xff>
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movsd  %xmm0,0x20(%rsp)
jmp    1249 <subsample_double+0x129>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
divsd  %xmm1,%xmm0
call   1030 <ceil@plt>
movsd  0xdc1(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
test   %r14d,%r14d
mov    0x8(%rsp),%rcx
mov    (%rsp),%r8
mov    0x38(%rsp),%rsi
movsd  0x30(%rsp),%xmm1
jle    1524 <subsample_double+0x404>
mov    %r15,%r10
mov    %r12d,%r11d
xor    %r15d,%r15d
movsd  0xd88(%rip),%xmm8
xorpd  %xmm2,%xmm2
jmp    1291 <subsample_double+0x171>
xchg   %ax,%ax
mov    0x8(%rcx),%r14d
add    $0x1,%r15d
cmp    %r14d,%r15d
jge    1524 <subsample_double+0x404>
xorps  %xmm6,%xmm6
cvtsi2sd %r15d,%xmm6
divsd  0x20(%rsp),%xmm6
addsd  %xmm8,%xmm6
test   %ebp,%ebp
jle    1284 <subsample_double+0x164>
xor    %r12d,%r12d
movsd  %xmm6,0x40(%rsp)
jmp    12c9 <subsample_double+0x1a9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r12d
cmp    %ebp,%r12d
jge    1280 <subsample_double+0x160>
xorps  %xmm3,%xmm3
cvtsi2sd %r12d,%xmm3
divsd  0x28(%rsp),%xmm3
addsd  %xmm8,%xmm3
test   %ebx,%ebx
jle    12c0 <subsample_double+0x1a0>
xor    %r14d,%r14d
movsd  %xmm3,0x10(%rsp)
jmp    1304 <subsample_double+0x1e4>
nopl   0x0(%rax,%rax,1)
movq   $0x0,(%rsi,%r13,8)
add    $0x1,%r14d
mov    0x4(%rcx),%ebp
cmp    %ebx,%r14d
jge    12c0 <subsample_double+0x1a0>
imul   %r15d,%ebp
add    %r12d,%ebp
imul   %ebx,%ebp
add    %r14d,%ebp
ucomisd %xmm2,%xmm6
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
divsd  %xmm1,%xmm0
addsd  %xmm8,%xmm0
movslq %ebp,%r13
jb     12f0 <subsample_double+0x1d0>
xorps  %xmm4,%xmm4
cvtsi2sdl 0x8(%r8),%xmm4
ucomisd %xmm6,%xmm4
jbe    12f0 <subsample_double+0x1d0>
ucomisd %xmm2,%xmm3
jb     12f0 <subsample_double+0x1d0>
xorps  %xmm4,%xmm4
cvtsi2sd %r11d,%xmm4
ucomisd %xmm3,%xmm4
jbe    12f0 <subsample_double+0x1d0>
ucomisd %xmm2,%xmm0
jb     12f0 <subsample_double+0x1d0>
xorps  %xmm4,%xmm4
cvtsi2sd %r10d,%xmm4
ucomisd %xmm0,%xmm4
jbe    12f0 <subsample_double+0x1d0>
movsd  %xmm0,0x18(%rsp)
call   1040 <floor@plt>
mov    %rax,%rbx
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  0x18(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsp)
movsd  0xc71(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x50(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <floor@plt>
mov    %rax,%rbp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
movsd  0x10(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x58(%rsp)
movsd  0xc3a(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x60(%rsp)
movsd  0x40(%rsp),%xmm0
call   1040 <floor@plt>
movsd  0x40(%rsp),%xmm6
movsd  0xc16(%rip),%xmm8
mov    (%rsp),%r8
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movapd %xmm6,%xmm9
subsd  %xmm1,%xmm9
movapd %xmm8,%xmm1
subsd  %xmm9,%xmm1
mov    (%r8),%r10d
mov    0x4(%r8),%r11d
add    $0xffffffff,%eax
imul   %r11d,%eax
add    %ebp,%eax
add    $0xffffffff,%eax
imul   %r10d,%eax
lea    (%rax,%rbx,1),%ecx
add    $0xffffffff,%ecx
add    %ebx,%eax
movslq %ecx,%rdx
movslq %eax,%rsi
lea    (%r10,%rsi,1),%eax
add    $0xffffffff,%eax
movslq %eax,%rdi
lea    (%rsi,%r10,1),%eax
movslq %eax,%rbp
mov    %r11d,%eax
imul   %r10d,%eax
lea    (%rax,%rsi,1),%ecx
add    $0xffffffff,%ecx
movslq %ecx,%rcx
add    %esi,%eax
movslq %eax,%r9
lea    (%rcx,%r10,1),%eax
cltq
mov    0x48(%rsp),%rbx
movsd  (%rbx,%rsi,8),%xmm3
mov    0x38(%rsp),%rsi
movsd  0x18(%rsp),%xmm4
mulsd  %xmm4,%xmm3
movsd  (%rbx,%rdx,8),%xmm2
movsd  0x50(%rsp),%xmm5
mulsd  %xmm5,%xmm2
addsd  %xmm3,%xmm2
movsd  (%rbx,%rbp,8),%xmm3
mulsd  %xmm4,%xmm3
movapd %xmm4,%xmm7
movsd  (%rbx,%rdi,8),%xmm4
mulsd  %xmm5,%xmm4
movapd %xmm5,%xmm10
addsd  %xmm3,%xmm4
movsd  0x58(%rsp),%xmm5
mulsd  %xmm5,%xmm4
movsd  0x60(%rsp),%xmm0
mulsd  %xmm0,%xmm2
addsd  %xmm4,%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rbx,%r9,8),%xmm1
mulsd  %xmm7,%xmm1
movsd  (%rbx,%rcx,8),%xmm3
mov    0x8(%rsp),%rcx
movapd %xmm10,%xmm4
mulsd  %xmm10,%xmm3
addsd  %xmm1,%xmm3
movsd  0x30(%rsp),%xmm1
mulsd  0x8(%rbx,%rax,8),%xmm7
mulsd  (%rbx,%rax,8),%xmm4
addsd  %xmm7,%xmm4
mulsd  %xmm5,%xmm4
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm9,%xmm3
addsd  %xmm2,%xmm3
xorpd  %xmm2,%xmm2
movsd  %xmm3,(%rsi,%r13,8)
movsd  0x10(%rsp),%xmm3
mov    (%rcx),%ebx
jmp    12f8 <subsample_double+0x1d8>
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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