void smooth_double(double *vol, int dims[3], double separations[3], double s[3],
                   int use_mask) {
  int i;
  double xsum, ysum, zsum;
  double *x, *y, *z;
  int xyz[3], nvol, sum_mask;
  double *mask;
  unsigned char *mask2;
  nvol = dims[0] * dims[1] * dims[2];
  for (i = 0; i < 3; i++) {
    s[i] /= separations[i];
    if (s[i] < 1.0)
      s[i] = 1.0;
    s[i] /= sqrt(8.0 * log(2.0));
    xyz[i] = (int)floor((6.0 * s[i]) + 0.5);
  }
  x = (double *)malloc(sizeof(double) * ((2 * xyz[0]) + 1));
  y = (double *)malloc(sizeof(double) * ((2 * xyz[1]) + 1));
  z = (double *)malloc(sizeof(double) * ((2 * xyz[2]) + 1));
  if (use_mask) {
    mask = (double *)malloc(sizeof(double) * nvol);
    mask2 = (unsigned char *)malloc(sizeof(unsigned char) * nvol);
    sum_mask = 0;
    for (i = 0; i < nvol; i++) {
      if (vol[i] > 0.0) {
        mask[i] = 1.0;
        mask2[i] = 1;
        sum_mask++;
      } else {
        mask[i] = 0.0;
        mask2[i] = 0;
      }
    }
  }
  for (i = -xyz[0]; i <= xyz[0]; i++)
    x[i + xyz[0]] = (double)i;
  for (i = -xyz[1]; i <= xyz[1]; i++)
    y[i + xyz[1]] = (double)i;
  for (i = -xyz[2]; i <= xyz[2]; i++)
    z[i + xyz[2]] = (double)i;
  xsum = 0.0;
  ysum = 0.0;
  zsum = 0.0;
  for (i = 0; i < ((2 * xyz[0]) + 1); i++) {
    x[i] = exp(-pow(x[i], 2) / (2.0 * pow(s[0], 2)));
    xsum += x[i];
  }
  for (i = 0; i < ((2 * xyz[1]) + 1); i++) {
    y[i] = exp(-pow(y[i], 2) / (2.0 * pow(s[1], 2)));
    ysum += y[i];
  }
  for (i = 0; i < ((2 * xyz[2]) + 1); i++) {
    z[i] = exp(-pow(z[i], 2) / (2.0 * pow(s[2], 2)));
    zsum += z[i];
  }
  for (i = 0; i < ((2 * xyz[0]) + 1); i++)
    x[i] /= xsum;
  for (i = 0; i < ((2 * xyz[1]) + 1); i++)
    y[i] /= ysum;
  for (i = 0; i < ((2 * xyz[2]) + 1); i++)
    z[i] /= zsum;
  convxyz_double(vol, x, y, z, ((2 * xyz[0]) + 1), ((2 * xyz[1]) + 1),
                 ((2 * xyz[2]) + 1), -xyz[0], -xyz[1], -xyz[2], vol, dims);
  if (use_mask) {
    if (sum_mask > 0) {
      convxyz_double(mask, x, y, z, ((2 * xyz[0]) + 1), ((2 * xyz[1]) + 1),
                     ((2 * xyz[2]) + 1), -xyz[0], -xyz[1], -xyz[2], mask, dims);
      for (i = 0; i < nvol; i++) {
        if (mask2[i] > 0)
          vol[i] /= (double)mask[i];
        else
          vol[i] = 0.0;
      }
    }
    free(mask);
    free(mask2);
  }
  free(x);
  free(y);
  free(z);
}


x86 assembly (clang, -O0):
<smooth_double>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %rbx
sub    $0xc8,%rsp
mov    %rdi,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
mov    %rcx,-0x38(%rbp)
mov    %r8d,-0x3c(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    -0x28(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x28(%rbp),%rcx
imul   0x8(%rcx),%eax
mov    %eax,-0x80(%rbp)
movl   $0x0,-0x40(%rbp)
cmpl   $0x3,-0x40(%rbp)
jge    129f <smooth_double+0x11f>
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  0xe07(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1221 <smooth_double+0xa1>
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  0xdec(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  0xdd7(%rip),%xmm0
call   1070 <log@plt>
movaps %xmm0,%xmm1
movsd  0xde7(%rip),%xmm0
mulsd  %xmm1,%xmm0
call   1090 <sqrt@plt>
movaps %xmm0,%xmm1
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xda0(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd9c(%rip),%xmm1
addsd  %xmm1,%xmm0
call   10a0 <floor@plt>
cvttsd2si %xmm0,%ecx
movslq -0x40(%rbp),%rax
mov    %ecx,-0x7c(%rbp,%rax,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11c4 <smooth_double+0x44>
mov    -0x7c(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x60(%rbp)
mov    -0x78(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x68(%rbp)
mov    -0x74(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x70(%rbp)
cmpl   $0x0,-0x3c(%rbp)
je     13c8 <smooth_double+0x248>
movslq -0x80(%rbp),%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x90(%rbp)
movslq -0x80(%rbp),%rdi
shl    $0x0,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x98(%rbp)
movl   $0x0,-0x84(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x80(%rbp),%eax
jge    13c3 <smooth_double+0x243>
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    138e <smooth_double+0x20e>
mov    -0x90(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  0xca2(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x98(%rbp),%rax
movslq -0x40(%rbp),%rcx
movb   $0x1,(%rax,%rcx,1)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    13b0 <smooth_double+0x230>
mov    -0x90(%rbp),%rax
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x98(%rbp),%rax
movslq -0x40(%rbp),%rcx
movb   $0x0,(%rax,%rcx,1)
jmp    13b5 <smooth_double+0x235>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    132d <smooth_double+0x1ad>
jmp    13c8 <smooth_double+0x248>
xor    %eax,%eax
sub    -0x7c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x7c(%rbp),%eax
jg     1401 <smooth_double+0x281>
cvtsi2sdl -0x40(%rbp),%xmm0
mov    -0x60(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x7c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    13d0 <smooth_double+0x250>
xor    %eax,%eax
sub    -0x78(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x78(%rbp),%eax
jg     143a <smooth_double+0x2ba>
cvtsi2sdl -0x40(%rbp),%xmm0
mov    -0x68(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x78(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1409 <smooth_double+0x289>
xor    %eax,%eax
sub    -0x74(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x74(%rbp),%eax
jg     1473 <smooth_double+0x2f3>
cvtsi2sdl -0x40(%rbp),%xmm0
mov    -0x70(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x74(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1442 <smooth_double+0x2c2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x7c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1546 <smooth_double+0x3c6>
mov    -0x60(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xb45(%rip),%xmm1
call   1080 <pow@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xb11(%rip),%xmm1
call   1080 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0xa0(%rbp),%xmm0
ff
movsd  0xaf9(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1060 <exp@plt>
mov    -0x60(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1492 <smooth_double+0x312>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x78(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1602 <smooth_double+0x482>
mov    -0x68(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xa8a(%rip),%xmm1
call   1080 <pow@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x38(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  0xa55(%rip),%xmm1
call   1080 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0xa8(%rbp),%xmm0
ff
movsd  0xa3d(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1060 <exp@plt>
mov    -0x68(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    154d <smooth_double+0x3cd>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x74(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    16be <smooth_double+0x53e>
mov    -0x70(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0x9ce(%rip),%xmm1
call   1080 <pow@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x38(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  0x999(%rip),%xmm1
call   1080 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0x981(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1060 <exp@plt>
mov    -0x70(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x70(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1609 <smooth_double+0x489>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x7c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1702 <smooth_double+0x582>
movsd  -0x48(%rbp),%xmm1
mov    -0x60(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    16c5 <smooth_double+0x545>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x78(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1746 <smooth_double+0x5c6>
movsd  -0x50(%rbp),%xmm1
mov    -0x68(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1709 <smooth_double+0x589>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x74(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    178a <smooth_double+0x60a>
movsd  -0x58(%rbp),%xmm1
mov    -0x70(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    174d <smooth_double+0x5cd>
mov    -0x20(%rbp),%rdi
mov    -0x60(%rbp),%rsi
mov    -0x68(%rbp),%rdx
mov    -0x70(%rbp),%rcx
mov    -0x7c(%rbp),%r8d
shl    $0x1,%r8d
add    $0x1,%r8d
mov    -0x78(%rbp),%r9d
shl    $0x1,%r9d
add    $0x1,%r9d
mov    -0x74(%rbp),%r15d
shl    $0x1,%r15d
add    $0x1,%r15d
xor    %r14d,%r14d
sub    -0x7c(%rbp),%r14d
xor    %ebx,%ebx
sub    -0x78(%rbp),%ebx
xor    %r11d,%r11d
sub    -0x74(%rbp),%r11d
mov    -0x20(%rbp),%r10
mov    -0x28(%rbp),%rax
mov    %r15d,(%rsp)
mov    %r14d,0x8(%rsp)
mov    %ebx,0x10(%rsp)
mov    %r11d,0x18(%rsp)
mov    %r10,0x20(%rsp)
mov    %rax,0x28(%rsp)
mov    $0x0,%al
call   1040 <convxyz_double@plt>
cmpl   $0x0,-0x3c(%rbp)
je     1921 <smooth_double+0x7a1>
cmpl   $0x0,-0x84(%rbp)
jle    1909 <smooth_double+0x789>
mov    -0x90(%rbp),%rdi
mov    -0x60(%rbp),%rsi
mov    -0x68(%rbp),%rdx
mov    -0x70(%rbp),%rcx
mov    -0x7c(%rbp),%r8d
shl    $0x1,%r8d
add    $0x1,%r8d
mov    -0x78(%rbp),%r9d
shl    $0x1,%r9d
add    $0x1,%r9d
mov    -0x74(%rbp),%r15d
shl    $0x1,%r15d
add    $0x1,%r15d
xor    %r14d,%r14d
sub    -0x7c(%rbp),%r14d
xor    %ebx,%ebx
sub    -0x78(%rbp),%ebx
xor    %r11d,%r11d
sub    -0x74(%rbp),%r11d
mov    -0x90(%rbp),%r10
mov    -0x28(%rbp),%rax
mov    %r15d,(%rsp)
mov    %r14d,0x8(%rsp)
mov    %ebx,0x10(%rsp)
mov    %r11d,0x18(%rsp)
mov    %r10,0x20(%rsp)
mov    %rax,0x28(%rsp)
mov    $0x0,%al
call   1040 <convxyz_double@plt>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x80(%rbp),%eax
jge    1904 <smooth_double+0x784>
mov    -0x98(%rbp),%rax
movslq -0x40(%rbp),%rcx
movzbl (%rax,%rcx,1),%eax
cmp    $0x0,%eax
jle    18e1 <smooth_double+0x761>
mov    -0x90(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    18f1 <smooth_double+0x771>
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    18f6 <smooth_double+0x776>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1892 <smooth_double+0x712>
jmp    1909 <smooth_double+0x789>
mov    -0x90(%rbp),%rdi
call   1030 <free@plt>
mov    -0x98(%rbp),%rdi
call   1030 <free@plt>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
mov    -0x68(%rbp),%rdi
call   1030 <free@plt>
mov    -0x70(%rbp),%rdi
call   1030 <free@plt>
add    $0xc8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

x86 assembly (clang, -O2):
<smooth_double>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
mov    %r8d,%r14d
mov    %rcx,%rbx
mov    %rdx,%r15
mov    %rdi,0x40(%rsp)
mov    0x4(%rsi),%r12d
movsd  (%rcx),%xmm0
movsd  0x8(%rcx),%xmm1
movsd  %xmm1,0x50(%rsp)
divsd  (%rdx),%xmm0
movsd  0xe72(%rip),%xmm1
maxsd  %xmm0,%xmm1
imul   (%rsi),%r12d
divsd  0xe6a(%rip),%xmm1
mov    %rsi,0xb8(%rsp)
imul   0x8(%rsi),%r12d
movsd  %xmm1,(%rcx)
movsd  %xmm1,0x20(%rsp)
movapd %xmm1,%xmm0
mulsd  0xe4f(%rip),%xmm0
addsd  0xe4f(%rip),%xmm0
call   1070 <floor@plt>
cvttsd2si %xmm0,%r13d
movsd  0x50(%rsp),%xmm1
divsd  0x8(%r15),%xmm1
movsd  0xe19(%rip),%xmm0
maxsd  %xmm1,%xmm0
divsd  0xe15(%rip),%xmm0
movsd  %xmm0,0x8(%rbx)
mulsd  0xe10(%rip),%xmm0
addsd  0xe10(%rip),%xmm0
call   1070 <floor@plt>
cvttsd2si %xmm0,%ebp
movsd  0x10(%rbx),%xmm1
divsd  0x10(%r15),%xmm1
movsd  0xddc(%rip),%xmm0
maxsd  %xmm1,%xmm0
divsd  0xdd8(%rip),%xmm0
mov    %rbx,0xa0(%rsp)
movsd  %xmm0,0x10(%rbx)
mulsd  0xdcb(%rip),%xmm0
addsd  0xdcb(%rip),%xmm0
call   1070 <floor@plt>
cvttsd2si %xmm0,%ebx
mov    %r13,0x68(%rsp)
lea    0x1(%r13,%r13,1),%eax
mov    %eax,0xc(%rsp)
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,0x10(%rsp)
mov    %rbp,0x70(%rsp)
lea    0x1(%rbp,%rbp,1),%eax
mov    %eax,0x8(%rsp)
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,%r15
lea    (%rbx,%rbx,1),%eax
add    $0x1,%eax
mov    %rax,0x48(%rsp)
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,%r13
mov    %r14d,0x7c(%rsp)
test   %r14d,%r14d
mov    %r12d,0x78(%rsp)
je     130e <smooth_double+0x1be>
movslq %r12d,%r14
lea    0x0(,%r14,8),%rdi
call   1050 <malloc@plt>
mov    %rax,0x38(%rsp)
mov    %r14,%rdi
call   1050 <malloc@plt>
test   %r14d,%r14d
mov    %rax,0x80(%rsp)
jle    1310 <smooth_double+0x1c0>
mov    %r12d,%r8d
mov    %rbx,%r9
cmp    $0x4,%r12d
jae    1328 <smooth_double+0x1d8>
xor    %ecx,%ecx
xor    %edi,%edi
mov    0x10(%rsp),%rsi
mov    0x38(%rsp),%rbp
mov    0x40(%rsp),%r14
movsd  0x20(%rsp),%xmm7
jmp    152c <smooth_double+0x3dc>
jmp    1318 <smooth_double+0x1c8>
movl   $0x0,0x1c(%rsp)
mov    0x10(%rsp),%rsi
movsd  0x20(%rsp),%xmm7
jmp    143f <smooth_double+0x2ef>
mov    %r8d,%ecx
and    $0xfffffffc,%ecx
xorpd  %xmm0,%xmm0
xor    %edx,%edx
movapd 0xce3(%rip),%xmm8
movdqa 0xcea(%rip),%xmm9
xorpd  %xmm1,%xmm1
mov    0x38(%rsp),%rbp
mov    0x40(%rsp),%r14
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%r14,%rdx,8),%xmm2
movupd 0x10(%r14,%rdx,8),%xmm3
xorpd  %xmm4,%xmm4
cmpltpd %xmm2,%xmm4
movapd %xmm4,%xmm2
andpd  %xmm8,%xmm2
xorpd  %xmm5,%xmm5
cmpltpd %xmm3,%xmm5
movapd %xmm5,%xmm3
andpd  %xmm8,%xmm3
movq   %xmm4,%rsi
movzbl %sil,%esi
movdqa %xmm4,%xmm6
punpckhqdq %xmm4,%xmm6
movq   %xmm6,%rdi
shl    $0x8,%edi
or     %esi,%edi
movd   %edi,%xmm6
pand   %xmm9,%xmm6
movq   %xmm5,%rsi
movzbl %sil,%esi
movdqa %xmm5,%xmm7
punpckhqdq %xmm5,%xmm7
movq   %xmm7,%rdi
shl    $0x8,%edi
or     %esi,%edi
movd   %edi,%xmm7
pand   %xmm9,%xmm7
shufps $0xe8,%xmm4,%xmm4
psubd  %xmm4,%xmm0
shufps $0xe8,%xmm5,%xmm5
psubd  %xmm5,%xmm1
movupd %xmm2,0x0(%rbp,%rdx,8)
movupd %xmm3,0x10(%rbp,%rdx,8)
movd   %xmm6,%esi
mov    %si,(%rax,%rdx,1)
movd   %xmm7,%esi
mov    %si,0x2(%rax,%rdx,1)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1360 <smooth_double+0x210>
paddd  %xmm0,%xmm1
pshufd $0x55,%xmm1,%xmm0
paddd  %xmm1,%xmm0
movd   %xmm0,%edi
cmp    %r8,%rcx
mov    0x10(%rsp),%rsi
movsd  0x20(%rsp),%xmm7
jne    152c <smooth_double+0x3dc>
mov    %edi,0x1c(%rsp)
mov    %r9,%rbx
mov    0x68(%rsp),%rax
mov    %eax,%r11d
neg    %r11d
test   %eax,%eax
js     15c4 <smooth_double+0x474>
mov    0x68(%rsp),%rcx
lea    (%rcx,%rcx,1),%eax
movslq %r11d,%r10
mov    %ecx,%r9d
lea    0x1(%rcx),%r8d
cmp    $0x4,%eax
jae    1471 <smooth_double+0x321>
mov    %r10,%rax
jmp    15a4 <smooth_double+0x454>
and    $0xfffffffc,%eax
movd   %r11d,%xmm0
pshufd $0x50,%xmm0,%xmm0
paddd  0xbba(%rip),%xmm0
lea    -0x4(%rax),%rdx
mov    %rdx,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rbx,%rbp
test   %rdx,%rdx
je     1576 <smooth_double+0x426>
mov    %rcx,%rdi
and    $0xfffffffffffffffe,%rdi
lea    (%r10,%r9,1),%rdx
lea    (%rsi,%rdx,8),%rbx
add    $0x30,%rbx
xor    %edx,%edx
movdqa 0xb92(%rip),%xmm1
movdqa 0xb9a(%rip),%xmm2
movdqa 0xba2(%rip),%xmm3
movdqa 0xbaa(%rip),%xmm4
cs nopw 0x0(%rax,%rax,1)
movdqa %xmm0,%xmm5
cvtdq2pd %xmm0,%xmm6
paddd  %xmm1,%xmm5
cvtdq2pd %xmm5,%xmm5
movups %xmm6,-0x30(%rbx,%rdx,8)
movups %xmm5,-0x20(%rbx,%rdx,8)
movdqa %xmm0,%xmm5
paddd  %xmm2,%xmm5
movdqa %xmm0,%xmm6
cvtdq2pd %xmm5,%xmm5
paddd  %xmm3,%xmm6
cvtdq2pd %xmm6,%xmm6
movups %xmm5,-0x10(%rbx,%rdx,8)
movdqu %xmm6,(%rbx,%rdx,8)
add    $0x8,%rdx
paddd  %xmm4,%xmm0
add    $0xfffffffffffffffe,%rdi
jne    14e0 <smooth_double+0x390>
jmp    1578 <smooth_double+0x428>
movq   0xacc(%rip),%xmm2
pxor   %xmm1,%xmm1
jmp    155b <smooth_double+0x40b>
nopw   0x0(%rax,%rax,1)
mov    %bl,%dl
seta   (%rax,%rcx,1)
movq   %xmm0,0x0(%rbp,%rcx,8)
add    %edx,%edi
add    $0x1,%rcx
cmp    %rcx,%r8
je     1438 <smooth_double+0x2e8>
movsd  (%r14,%rcx,8),%xmm0
xor    %edx,%edx
ucomisd %xmm1,%xmm0
seta   %bl
movdqa %xmm2,%xmm0
ja     1540 <smooth_double+0x3f0>
pxor   %xmm0,%xmm0
jmp    1540 <smooth_double+0x3f0>
xor    %edx,%edx
add    %r10,%rax
test   $0x1,%cl
mov    %rbp,%rbx
je     15a4 <smooth_double+0x454>
cvtdq2pd %xmm0,%xmm1
paddd  0xac1(%rip),%xmm0
cvtdq2pd %xmm0,%xmm0
add    %r10,%rdx
add    %r9,%rdx
movupd %xmm1,(%rsi,%rdx,8)
movdqu %xmm0,0x10(%rsi,%rdx,8)
lea    (%rsi,%r9,8),%rcx
mov    %r8d,%edx
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rax
cmp    %eax,%edx
jne    15b0 <smooth_double+0x460>
mov    0x70(%rsp),%rax
mov    %eax,%ebp
neg    %ebp
test   %eax,%eax
js     1704 <smooth_double+0x5b4>
mov    0x70(%rsp),%rcx
lea    (%rcx,%rcx,1),%eax
movslq %ebp,%r10
mov    %ecx,%r9d
lea    0x1(%rcx),%r8d
cmp    $0x4,%eax
jae    15f4 <smooth_double+0x4a4>
mov    %r10,%rax
jmp    16dc <smooth_double+0x58c>
and    $0xfffffffc,%eax
movd   %ebp,%xmm0
pshufd $0x50,%xmm0,%xmm0
paddd  0xa38(%rip),%xmm0
lea    -0x4(%rax),%rdx
mov    %rdx,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rbx,%r12
test   %rdx,%rdx
je     16ac <smooth_double+0x55c>
mov    %rcx,%rdi
and    $0xfffffffffffffffe,%rdi
lea    (%r10,%r9,1),%rdx
lea    (%r15,%rdx,8),%rbx
add    $0x30,%rbx
xor    %edx,%edx
movdqa 0xa10(%rip),%xmm1
movdqa 0xa18(%rip),%xmm2
movdqa 0xa20(%rip),%xmm3
movdqa 0xa28(%rip),%xmm4
nopl   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm5
cvtdq2pd %xmm0,%xmm6
paddd  %xmm1,%xmm5
cvtdq2pd %xmm5,%xmm5
movups %xmm6,-0x30(%rbx,%rdx,8)
movups %xmm5,-0x20(%rbx,%rdx,8)
movdqa %xmm0,%xmm5
paddd  %xmm2,%xmm5
movdqa %xmm0,%xmm6
cvtdq2pd %xmm5,%xmm5
paddd  %xmm3,%xmm6
cvtdq2pd %xmm6,%xmm6
movups %xmm5,-0x10(%rbx,%rdx,8)
movdqu %xmm6,(%rbx,%rdx,8)
add    $0x8,%rdx
paddd  %xmm4,%xmm0
add    $0xfffffffffffffffe,%rdi
jne    1660 <smooth_double+0x510>
jmp    16ae <smooth_double+0x55e>
xor    %edx,%edx
add    %r10,%rax
test   $0x1,%cl
mov    %r12,%rbx
je     16dc <smooth_double+0x58c>
cvtdq2pd %xmm0,%xmm1
paddd  0x98b(%rip),%xmm0
cvtdq2pd %xmm0,%xmm0
add    %r10,%rdx
add    %r9,%rdx
movupd %xmm1,(%r15,%rdx,8)
movdqu %xmm0,0x10(%r15,%rdx,8)
lea    (%r15,%r9,8),%rcx
mov    %r8d,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rax
cmp    %eax,%edx
jne    16f0 <smooth_double+0x5a0>
mov    %ebx,%r12d
neg    %r12d
test   %ebx,%ebx
mov    %rbx,0xa8(%rsp)
js     1844 <smooth_double+0x6f4>
lea    (%rbx,%rbx,1),%eax
movslq %r12d,%r10
mov    %ebx,%r9d
lea    0x1(%rbx),%r8d
cmp    $0x4,%eax
jae    1734 <smooth_double+0x5e4>
mov    %r10,%rax
jmp    181a <smooth_double+0x6ca>
and    $0xfffffffc,%eax
movd   %r12d,%xmm0
pshufd $0x50,%xmm0,%xmm0
paddd  0x8f7(%rip),%xmm0
lea    -0x4(%rax),%rdx
mov    %rdx,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     17ec <smooth_double+0x69c>
mov    %rcx,%rdi
and    $0xfffffffffffffffe,%rdi
lea    (%r10,%r9,1),%rdx
lea    0x30(,%rdx,8),%rbx
add    %r13,%rbx
xor    %edx,%edx
movdqa 0x8cf(%rip),%xmm1
movdqa 0x8d7(%rip),%xmm2
movdqa 0x8df(%rip),%xmm3
movdqa 0x8e7(%rip),%xmm4
nopl   0x0(%rax)
movdqa %xmm0,%xmm5
cvtdq2pd %xmm0,%xmm6
paddd  %xmm1,%xmm5
cvtdq2pd %xmm5,%xmm5
movups %xmm6,-0x30(%rbx,%rdx,8)
movups %xmm5,-0x20(%rbx,%rdx,8)
movdqa %xmm0,%xmm5
paddd  %xmm2,%xmm5
movdqa %xmm0,%xmm6
cvtdq2pd %xmm5,%xmm5
paddd  %xmm3,%xmm6
cvtdq2pd %xmm6,%xmm6
movups %xmm5,-0x10(%rbx,%rdx,8)
movdqu %xmm6,(%rbx,%rdx,8)
add    $0x8,%rdx
paddd  %xmm4,%xmm0
add    $0xfffffffffffffffe,%rdi
jne    17a0 <smooth_double+0x650>
jmp    17ee <smooth_double+0x69e>
xor    %edx,%edx
add    %r10,%rax
test   $0x1,%cl
je     181a <smooth_double+0x6ca>
cvtdq2pd %xmm0,%xmm1
paddd  0x84e(%rip),%xmm0
cvtdq2pd %xmm0,%xmm0
add    %r10,%rdx
add    %r9,%rdx
movupd %xmm1,0x0(%r13,%rdx,8)
movdqu %xmm0,0x10(%r13,%rdx,8)
lea    0x0(,%r9,8),%rcx
add    %r13,%rcx
mov    %r8d,%edx
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rax
cmp    %eax,%edx
jne    1830 <smooth_double+0x6e0>
mov    %rbp,0x88(%rsp)
mov    %r11,0x90(%rsp)
mov    0xc(%rsp),%eax
mov    %eax,%ecx
mov    %rcx,0x98(%rsp)
xorpd  %xmm1,%xmm1
test   %eax,%eax
jle    18e7 <smooth_double+0x797>
mov    0x98(%rsp),%rax
lea    -0x1(%rax),%rbx
pxor   %xmm0,%xmm0
movdqa %xmm0,0x20(%rsp)
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rsi,%rbp,8),%xmm1
movapd %xmm1,%xmm0
xorpd  0x7ef(%rip),%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm7,%xmm7
addsd  %xmm7,%xmm7
divsd  %xmm7,%xmm0
call   1060 <exp@plt>
mov    0x10(%rsp),%rsi
movsd  %xmm0,(%rsi,%rbp,8)
movapd 0x20(%rsp),%xmm1
addsd  %xmm0,%xmm1
cmp    %rbp,%rbx
je     18e7 <smooth_double+0x797>
movapd %xmm1,0x20(%rsp)
mov    0xa0(%rsp),%rax
movsd  (%rax),%xmm7
add    $0x1,%rbp
jmp    1890 <smooth_double+0x740>
movapd %xmm1,0x20(%rsp)
mov    %r12,0xb0(%rsp)
mov    0x8(%rsp),%eax
mov    %eax,%r14d
test   %eax,%eax
mov    0xa0(%rsp),%rbp
pxor   %xmm0,%xmm0
jle    1977 <smooth_double+0x827>
pxor   %xmm0,%xmm0
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movdqa %xmm0,0x50(%rsp)
movsd  (%r15,%rbx,8),%xmm1
movapd %xmm1,%xmm0
xorpd  0x758(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
call   1060 <exp@plt>
movsd  %xmm0,(%r15,%rbx,8)
movapd 0x50(%rsp),%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,0x50(%rsp)
movdqa 0x50(%rsp),%xmm0
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1920 <smooth_double+0x7d0>
movdqa %xmm0,0x50(%rsp)
mov    0x48(%rsp),%rax
mov    %eax,%r12d
xorps  %xmm4,%xmm4
test   %eax,%eax
jle    19e2 <smooth_double+0x892>
xor    %ebx,%ebx
xchg   %ax,%ax
movaps %xmm4,0xc0(%rsp)
movsd  0x0(%r13,%rbx,8),%xmm1
movapd %xmm1,%xmm0
xorpd  0x6e5(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
call   1060 <exp@plt>
movapd 0xc0(%rsp),%xmm4
movsd  %xmm0,0x0(%r13,%rbx,8)
addsd  %xmm0,%xmm4
add    $0x1,%rbx
cmp    %rbx,%r12
jne    1990 <smooth_double+0x840>
cmpl   $0x0,0xc(%rsp)
mov    0x10(%rsp),%rsi
mov    0x40(%rsp),%rdi
mov    0x90(%rsp),%r10
mov    0x88(%rsp),%r11
movapd 0x50(%rsp),%xmm3
mov    0x98(%rsp),%rbp
movapd 0x20(%rsp),%xmm5
jle    1ab7 <smooth_double+0x967>
cmpl   $0x0,0x68(%rsp)
je     1a7a <smooth_double+0x92a>
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
lea    -0x3(%rbp),%rcx
mov    %rcx,%rax
shr    %rax
add    $0x1,%rax
test   %rcx,%rcx
je     1a7e <smooth_double+0x92e>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movupd (%rsi,%rcx,8),%xmm1
movupd 0x10(%rsi,%rcx,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,(%rsi,%rcx,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%rsi,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rdx
jne    1a50 <smooth_double+0x900>
jmp    1a80 <smooth_double+0x930>
xor    %edx,%edx
jmp    1aa0 <smooth_double+0x950>
xor    %ecx,%ecx
lea    -0x1(%rbp),%rdx
test   $0x1,%al
je     1aa0 <smooth_double+0x950>
movupd (%rsi,%rcx,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%rsi,%rcx,8)
cs nopw 0x0(%rax,%rax,1)
movsd  (%rsi,%rdx,8),%xmm0
divsd  %xmm5,%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rbp
jne    1aa0 <smooth_double+0x950>
cmpl   $0x0,0x8(%rsp)
jle    1b59 <smooth_double+0xa09>
cmpl   $0x0,0x70(%rsp)
je     1b1e <smooth_double+0x9ce>
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x3(%r14),%rcx
mov    %rcx,%rax
shr    %rax
add    $0x1,%rax
test   %rcx,%rcx
je     1b22 <smooth_double+0x9d2>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %ecx,%ecx
nopl   (%rax)
movupd (%r15,%rcx,8),%xmm1
movupd 0x10(%r15,%rcx,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rcx,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%r15,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rdx
jne    1af0 <smooth_double+0x9a0>
jmp    1b24 <smooth_double+0x9d4>
xor    %edx,%edx
jmp    1b40 <smooth_double+0x9f0>
xor    %ecx,%ecx
lea    -0x1(%r14),%rdx
test   $0x1,%al
je     1b40 <smooth_double+0x9f0>
movupd (%r15,%rcx,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rcx,8)
nopl   0x0(%rax)
movsd  (%r15,%rdx,8),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%r15,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r14
jne    1b40 <smooth_double+0x9f0>
cmpl   $0x0,0x48(%rsp)
jle    1c1b <smooth_double+0xacb>
cmpl   $0x0,0xa8(%rsp)
je     1bd0 <smooth_double+0xa80>
movapd %xmm4,%xmm0
unpcklpd %xmm4,%xmm0
lea    -0x3(%r12),%rcx
mov    %rcx,%rax
shr    %rax
add    $0x1,%rax
test   %rcx,%rcx
je     1bd4 <smooth_double+0xa84>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd 0x0(%r13,%rcx,8),%xmm1
movupd 0x10(%r13,%rcx,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,0x0(%r13,%rcx,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%r13,%rcx,8)
add    $0x4,%rcx
add    $0xfffffffffffffffe,%rdx
jne    1ba0 <smooth_double+0xa50>
jmp    1bd6 <smooth_double+0xa86>
xor    %edx,%edx
jmp    1c00 <smooth_double+0xab0>
xor    %ecx,%ecx
lea    -0x1(%r12),%rdx
test   $0x1,%al
je     1c00 <smooth_double+0xab0>
movupd 0x0(%r13,%rcx,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,0x0(%r13,%rcx,8)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rdx,8),%xmm0
divsd  %xmm4,%xmm0
movsd  %xmm0,0x0(%r13,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r12
jne    1c00 <smooth_double+0xab0>
mov    %r15,%rdx
mov    %r13,%rcx
mov    0xc(%rsp),%r8d
mov    0x8(%rsp),%r9d
xor    %eax,%eax
mov    0xb8(%rsp),%rbx
push   %rbx
push   %rdi
mov    0xc0(%rsp),%rbp
push   %rbp
push   %r11
push   %r10
push   0x70(%rsp)
call   1040 <convxyz_double@plt>
add    $0x30,%rsp
cmpl   $0x0,0x7c(%rsp)
je     1d0b <smooth_double+0xbbb>
cmpl   $0x0,0x1c(%rsp)
mov    0x78(%rsp),%r14d
mov    0x38(%rsp),%rdi
mov    0x80(%rsp),%rdx
jle    1cfb <smooth_double+0xbab>
mov    0x10(%rsp),%rsi
mov    %r15,%rdx
mov    %r13,%rcx
mov    0xc(%rsp),%r8d
mov    0x8(%rsp),%r9d
xor    %eax,%eax
push   %rbx
push   %rdi
push   %rbp
push   0xa0(%rsp)
push   0xb0(%rsp)
push   0x70(%rsp)
call   1040 <convxyz_double@plt>
mov    0xb0(%rsp),%rdx
mov    0x70(%rsp),%rbp
mov    0x68(%rsp),%rdi
add    $0x30,%rsp
test   %r14d,%r14d
jle    1cfb <smooth_double+0xbab>
mov    %r14d,%eax
xor    %ecx,%ecx
jmp    1ce2 <smooth_double+0xb92>
nopl   0x0(%rax)
movq   $0x0,0x0(%rbp,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
je     1cfb <smooth_double+0xbab>
cmpb   $0x0,(%rdx,%rcx,1)
je     1cd0 <smooth_double+0xb80>
movsd  0x0(%rbp,%rcx,8),%xmm0
divsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
jmp    1cd9 <smooth_double+0xb89>
mov    %rdx,%rbx
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x10(%rsp),%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

LLVM IR (-O2):
define void @smooth_double(double* %0, i32* %1, double* nocapture %2, double* nocapture %3, i32 %4) {
%6 = load i32, i32* %1
%7 = getelementptr inbounds i32, i32* %1, i64 1
%8 = load i32, i32* %7
%9 = getelementptr inbounds i32, i32* %1, i64 2
%10 = load i32, i32* %9
%11 = load double, double* %2
%12 = load double, double* %3
%13 = fdiv double %12, %11
%14 = fcmp olt double %13, 1.000000e+00
%15 = select i1 %14, double 1.000000e+00, double %13
%16 = fdiv double %15, 0x4002D6ABE44AFC43
store double %16, double* %3
%17 = tail call double @llvm.fmuladd.f64(double %16, double 6.000000e+00, double 5.000000e-01)
%18 = tail call double @llvm.floor.f64(double %17)
%19 = fptosi double %18 to i32
%20 = getelementptr inbounds double, double* %2, i64 1
%21 = load double, double* %20
%22 = getelementptr inbounds double, double* %3, i64 1
%23 = load double, double* %22
%24 = fdiv double %23, %21
%25 = fcmp olt double %24, 1.000000e+00
%26 = select i1 %25, double 1.000000e+00, double %24
%27 = fdiv double %26, 0x4002D6ABE44AFC43
store double %27, double* %22
%28 = tail call double @llvm.fmuladd.f64(double %27, double 6.000000e+00, double 5.000000e-01)
%29 = tail call double @llvm.floor.f64(double %28)
%30 = fptosi double %29 to i32
%31 = getelementptr inbounds double, double* %2, i64 2
%32 = load double, double* %31
%33 = getelementptr inbounds double, double* %3, i64 2
%34 = load double, double* %33
%35 = fdiv double %34, %32
%36 = fcmp olt double %35, 1.000000e+00
%37 = select i1 %36, double 1.000000e+00, double %35
%38 = fdiv double %37, 0x4002D6ABE44AFC43
store double %38, double* %33
%39 = tail call double @llvm.fmuladd.f64(double %38, double 6.000000e+00, double 5.000000e-01)
%40 = tail call double @llvm.floor.f64(double %39)
%41 = fptosi double %40 to i32
%42 = mul nsw i32 %8, %6
%43 = mul i32 %42, %10
%44 = shl nsw i32 %19, 1
%45 = or i32 %44, 1
%46 = sext i32 %45 to i64
%47 = shl nsw i64 %46, 3
%48 = tail call i8* @malloc(i64 %47)
%49 = bitcast i8* %48 to double*
%50 = shl nsw i32 %30, 1
%51 = or i32 %50, 1
%52 = sext i32 %51 to i64
%53 = shl nsw i64 %52, 3
%54 = tail call i8* @malloc(i64 %53)
%55 = bitcast i8* %54 to double*
%56 = shl nsw i32 %41, 1
%57 = or i32 %56, 1
%58 = sext i32 %57 to i64
%59 = shl nsw i64 %58, 3
%60 = tail call i8* @malloc(i64 %59)
%61 = bitcast i8* %60 to double*
%62 = icmp eq i32 %4, 0
br i1 %62, label %126, label %63
63: ; preds = %5
%64 = sext i32 %43 to i64
%65 = shl nsw i64 %64, 3
%66 = tail call i8* @malloc(i64 %65)
%67 = bitcast i8* %66 to double*
%68 = tail call i8* @malloc(i64 %64)
%69 = icmp sgt i32 %43, 0
br i1 %69, label %70, label %126
70: ; preds = %63
%71 = zext i32 %43 to i64
%72 = icmp ult i32 %43, 4
br i1 %72, label %109, label %73
73: ; preds = %70
%74 = and i64 %71, 4294967292
br label %75
75: ; preds = %75, %73
%76 = phi i64 [ 0, %73 ], [ %103, %75 ]
%77 = phi <2 x i32> [ zeroinitializer, %73 ], [ %93, %75 ]
%78 = phi <2 x i32> [ zeroinitializer, %73 ], [ %94, %75 ]
%79 = getelementptr inbounds double, double* %0, i64 %76
%80 = bitcast double* %79 to <2 x double>*
%81 = load <2 x double>, <2 x double>* %80
%82 = getelementptr inbounds double, double* %79, i64 2
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = fcmp ogt <2 x double> %81, zeroinitializer
%86 = fcmp ogt <2 x double> %84, zeroinitializer
%87 = select <2 x i1> %85, <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double> zeroinitializer
%88 = select <2 x i1> %86, <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double> zeroinitializer
%89 = zext <2 x i1> %85 to <2 x i8>
%90 = zext <2 x i1> %86 to <2 x i8>
%91 = zext <2 x i1> %85 to <2 x i32>
%92 = zext <2 x i1> %86 to <2 x i32>
%93 = add <2 x i32> %77, %91
%94 = add <2 x i32> %78, %92
%95 = getelementptr inbounds double, double* %67, i64 %76
%96 = bitcast double* %95 to <2 x double>*
store <2 x double> %87, <2 x double>* %96
%97 = getelementptr inbounds double, double* %95, i64 2
%98 = bitcast double* %97 to <2 x double>*
store <2 x double> %88, <2 x double>* %98
%99 = getelementptr inbounds i8, i8* %68, i64 %76
%100 = bitcast i8* %99 to <2 x i8>*
store <2 x i8> %89, <2 x i8>* %100
%101 = getelementptr inbounds i8, i8* %99, i64 2
%102 = bitcast i8* %101 to <2 x i8>*
store <2 x i8> %90, <2 x i8>* %102
%103 = add nuw i64 %76, 4
%104 = icmp eq i64 %103, %74
br i1 %104, label %105, label %75
105: ; preds = %75
%106 = add <2 x i32> %94, %93
%107 = call i32 @llvm.vector.reduce.add.v2i32(<2 x i32> %106)
%108 = icmp eq i64 %74, %71
br i1 %108, label %126, label %109
109: ; preds = %70, %105
%110 = phi i64 [ 0, %70 ], [ %74, %105 ]
%111 = phi i32 [ 0, %70 ], [ %107, %105 ]
br label %112
112: ; preds = %109, %112
%113 = phi i64 [ %124, %112 ], [ %110, %109 ]
%114 = phi i32 [ %121, %112 ], [ %111, %109 ]
%115 = getelementptr inbounds double, double* %0, i64 %113
%116 = load double, double* %115
%117 = fcmp ogt double %116, 0.000000e+00
%118 = select i1 %117, double 1.000000e+00, double 0.000000e+00
%119 = zext i1 %117 to i8
%120 = zext i1 %117 to i32
%121 = add nuw nsw i32 %114, %120
%122 = getelementptr inbounds double, double* %67, i64 %113
store double %118, double* %122
%123 = getelementptr inbounds i8, i8* %68, i64 %113
store i8 %119, i8* %123
%124 = add nuw nsw i64 %113, 1
%125 = icmp eq i64 %124, %71
br i1 %125, label %126, label %112
126: ; preds = %112, %105, %63, %5
%127 = phi i32 [ undef, %5 ], [ 0, %63 ], [ %107, %105 ], [ %121, %112 ]
%128 = phi double* [ undef, %5 ], [ %67, %63 ], [ %67, %105 ], [ %67, %112 ]
%129 = phi i8* [ undef, %5 ], [ %68, %63 ], [ %68, %105 ], [ %68, %112 ]
%130 = sub i32 0, %19
%131 = icmp slt i32 %19, 0
br i1 %131, label %204, label %132
132: ; preds = %126
%133 = sext i32 %130 to i64
%134 = zext i32 %19 to i64
%135 = add nuw i32 %19, 1
%136 = icmp ult i32 %44, 4
br i1 %136, label %193, label %137
137: ; preds = %132
%138 = and i32 %44, -4
%139 = zext i32 %138 to i64
%140 = add nsw i64 %139, %133
%141 = insertelement <2 x i32> poison, i32 %130, i64 0
%142 = shufflevector <2 x i32> %141, <2 x i32> poison, <2 x i32> zeroinitializer
%143 = add <2 x i32> %142, <i32 0, i32 1>
%144 = add nsw i64 %139, -4
%145 = lshr exact i64 %144, 2
%146 = add nuw nsw i64 %145, 1
%147 = and i64 %146, 1
%148 = icmp eq i64 %144, 0
br i1 %148, label %179, label %149
149: ; preds = %137
%150 = and i64 %146, 9223372036854775806
br label %151
151: ; preds = %151, %149
%152 = phi i64 [ 0, %149 ], [ %175, %151 ]
%153 = phi <2 x i32> [ %143, %149 ], [ %176, %151 ]
%154 = phi i64 [ 0, %149 ], [ %177, %151 ]
%155 = add i64 %152, %133
%156 = add <2 x i32> %153, <i32 2, i32 2>
%157 = sitofp <2 x i32> %153 to <2 x double>
%158 = sitofp <2 x i32> %156 to <2 x double>
%159 = add nsw i64 %155, %134
%160 = getelementptr inbounds double, double* %49, i64 %159
%161 = bitcast double* %160 to <2 x double>*
store <2 x double> %157, <2 x double>* %161
%162 = getelementptr inbounds double, double* %160, i64 2
%163 = bitcast double* %162 to <2 x double>*
store <2 x double> %158, <2 x double>* %163
%164 = or i64 %152, 4
%165 = add <2 x i32> %153, <i32 4, i32 4>
%166 = add i64 %164, %133
%167 = add <2 x i32> %153, <i32 6, i32 6>
%168 = sitofp <2 x i32> %165 to <2 x double>
%169 = sitofp <2 x i32> %167 to <2 x double>
%170 = add nsw i64 %166, %134
%171 = getelementptr inbounds double, double* %49, i64 %170
%172 = bitcast double* %171 to <2 x double>*
store <2 x double> %168, <2 x double>* %172
%173 = getelementptr inbounds double, double* %171, i64 2
%174 = bitcast double* %173 to <2 x double>*
store <2 x double> %169, <2 x double>* %174
%175 = add nuw i64 %152, 8
%176 = add <2 x i32> %153, <i32 8, i32 8>
%177 = add i64 %154, 2
%178 = icmp eq i64 %177, %150
br i1 %178, label %179, label %151
179: ; preds = %151, %137
%180 = phi i64 [ 0, %137 ], [ %175, %151 ]
%181 = phi <2 x i32> [ %143, %137 ], [ %176, %151 ]
%182 = icmp eq i64 %147, 0
br i1 %182, label %193, label %183
183: ; preds = %179
%184 = add i64 %180, %133
%185 = add <2 x i32> %181, <i32 2, i32 2>
%186 = sitofp <2 x i32> %181 to <2 x double>
%187 = sitofp <2 x i32> %185 to <2 x double>
%188 = add nsw i64 %184, %134
%189 = getelementptr inbounds double, double* %49, i64 %188
%190 = bitcast double* %189 to <2 x double>*
store <2 x double> %186, <2 x double>* %190
%191 = getelementptr inbounds double, double* %189, i64 2
%192 = bitcast double* %191 to <2 x double>*
store <2 x double> %187, <2 x double>* %192
br label %193
193: ; preds = %183, %179, %132
%194 = phi i64 [ %133, %132 ], [ %140, %179 ], [ %140, %183 ]
br label %195
195: ; preds = %193, %195
%196 = phi i64 [ %201, %195 ], [ %194, %193 ]
%197 = trunc i64 %196 to i32
%198 = sitofp i32 %197 to double
%199 = add nsw i64 %196, %134
%200 = getelementptr inbounds double, double* %49, i64 %199
store double %198, double* %200
%201 = add nsw i64 %196, 1
%202 = trunc i64 %201 to i32
%203 = icmp eq i32 %135, %202
br i1 %203, label %204, label %195
204: ; preds = %195, %126
%205 = sub i32 0, %30
%206 = icmp slt i32 %30, 0
br i1 %206, label %279, label %207
207: ; preds = %204
%208 = sext i32 %205 to i64
%209 = zext i32 %30 to i64
%210 = add nuw i32 %30, 1
%211 = icmp ult i32 %50, 4
br i1 %211, label %268, label %212
212: ; preds = %207
%213 = and i32 %50, -4
%214 = zext i32 %213 to i64
%215 = add nsw i64 %214, %208
%216 = insertelement <2 x i32> poison, i32 %205, i64 0
%217 = shufflevector <2 x i32> %216, <2 x i32> poison, <2 x i32> zeroinitializer
%218 = add <2 x i32> %217, <i32 0, i32 1>
%219 = add nsw i64 %214, -4
%220 = lshr exact i64 %219, 2
%221 = add nuw nsw i64 %220, 1
%222 = and i64 %221, 1
%223 = icmp eq i64 %219, 0
br i1 %223, label %254, label %224
224: ; preds = %212
%225 = and i64 %221, 9223372036854775806
br label %226
226: ; preds = %226, %224
%227 = phi i64 [ 0, %224 ], [ %250, %226 ]
%228 = phi <2 x i32> [ %218, %224 ], [ %251, %226 ]
%229 = phi i64 [ 0, %224 ], [ %252, %226 ]
%230 = add i64 %227, %208
%231 = add <2 x i32> %228, <i32 2, i32 2>
%232 = sitofp <2 x i32> %228 to <2 x double>
%233 = sitofp <2 x i32> %231 to <2 x double>
%234 = add nsw i64 %230, %209
%235 = getelementptr inbounds double, double* %55, i64 %234
%236 = bitcast double* %235 to <2 x double>*
store <2 x double> %232, <2 x double>* %236
%237 = getelementptr inbounds double, double* %235, i64 2
%238 = bitcast double* %237 to <2 x double>*
store <2 x double> %233, <2 x double>* %238
%239 = or i64 %227, 4
%240 = add <2 x i32> %228, <i32 4, i32 4>
%241 = add i64 %239, %208
%242 = add <2 x i32> %228, <i32 6, i32 6>
%243 = sitofp <2 x i32> %240 to <2 x double>
%244 = sitofp <2 x i32> %242 to <2 x double>
%245 = add nsw i64 %241, %209
%246 = getelementptr inbounds double, double* %55, i64 %245
%247 = bitcast double* %246 to <2 x double>*
store <2 x double> %243, <2 x double>* %247
%248 = getelementptr inbounds double, double* %246, i64 2
%249 = bitcast double* %248 to <2 x double>*
store <2 x double> %244, <2 x double>* %249
%250 = add nuw i64 %227, 8
%251 = add <2 x i32> %228, <i32 8, i32 8>
%252 = add i64 %229, 2
%253 = icmp eq i64 %252, %225
br i1 %253, label %254, label %226
254: ; preds = %226, %212
%255 = phi i64 [ 0, %212 ], [ %250, %226 ]
%256 = phi <2 x i32> [ %218, %212 ], [ %251, %226 ]
%257 = icmp eq i64 %222, 0
br i1 %257, label %268, label %258
258: ; preds = %254
%259 = add i64 %255, %208
%260 = add <2 x i32> %256, <i32 2, i32 2>
%261 = sitofp <2 x i32> %256 to <2 x double>
%262 = sitofp <2 x i32> %260 to <2 x double>
%263 = add nsw i64 %259, %209
%264 = getelementptr inbounds double, double* %55, i64 %263
%265 = bitcast double* %264 to <2 x double>*
store <2 x double> %261, <2 x double>* %265
%266 = getelementptr inbounds double, double* %264, i64 2
%267 = bitcast double* %266 to <2 x double>*
store <2 x double> %262, <2 x double>* %267
br label %268
268: ; preds = %258, %254, %207
%269 = phi i64 [ %208, %207 ], [ %215, %254 ], [ %215, %258 ]
br label %270
270: ; preds = %268, %270
%271 = phi i64 [ %276, %270 ], [ %269, %268 ]
%272 = trunc i64 %271 to i32
%273 = sitofp i32 %272 to double
%274 = add nsw i64 %271, %209
%275 = getelementptr inbounds double, double* %55, i64 %274
store double %273, double* %275
%276 = add nsw i64 %271, 1
%277 = trunc i64 %276 to i32
%278 = icmp eq i32 %210, %277
br i1 %278, label %279, label %270
279: ; preds = %270, %204
%280 = sub i32 0, %41
%281 = icmp slt i32 %41, 0
br i1 %281, label %345, label %282
282: ; preds = %279
%283 = sext i32 %280 to i64
%284 = zext i32 %41 to i64
%285 = add nuw i32 %41, 1
%286 = icmp ult i32 %56, 4
br i1 %286, label %301, label %303
287: ; preds = %317, %303
%288 = phi i64 [ 0, %303 ], [ %341, %317 ]
%289 = phi <2 x i32> [ %309, %303 ], [ %342, %317 ]
%290 = icmp eq i64 %313, 0
br i1 %290, label %301, label %291
291: ; preds = %287
%292 = add i64 %288, %283
%293 = add <2 x i32> %289, <i32 2, i32 2>
%294 = sitofp <2 x i32> %289 to <2 x double>
%295 = sitofp <2 x i32> %293 to <2 x double>
%296 = add nsw i64 %292, %284
%297 = getelementptr inbounds double, double* %61, i64 %296
%298 = bitcast double* %297 to <2 x double>*
store <2 x double> %294, <2 x double>* %298
%299 = getelementptr inbounds double, double* %297, i64 2
%300 = bitcast double* %299 to <2 x double>*
store <2 x double> %295, <2 x double>* %300
br label %301
301: ; preds = %291, %287, %282
%302 = phi i64 [ %283, %282 ], [ %306, %287 ], [ %306, %291 ]
br label %349
303: ; preds = %282
%304 = and i32 %56, -4
%305 = zext i32 %304 to i64
%306 = add nsw i64 %305, %283
%307 = insertelement <2 x i32> poison, i32 %280, i64 0
%308 = shufflevector <2 x i32> %307, <2 x i32> poison, <2 x i32> zeroinitializer
%309 = add <2 x i32> %308, <i32 0, i32 1>
%310 = add nsw i64 %305, -4
%311 = lshr exact i64 %310, 2
%312 = add nuw nsw i64 %311, 1
%313 = and i64 %312, 1
%314 = icmp eq i64 %310, 0
br i1 %314, label %287, label %315
315: ; preds = %303
%316 = and i64 %312, 9223372036854775806
br label %317
317: ; preds = %317, %315
%318 = phi i64 [ 0, %315 ], [ %341, %317 ]
%319 = phi <2 x i32> [ %309, %315 ], [ %342, %317 ]
%320 = phi i64 [ 0, %315 ], [ %343, %317 ]
%321 = add i64 %318, %283
%322 = add <2 x i32> %319, <i32 2, i32 2>
%323 = sitofp <2 x i32> %319 to <2 x double>
%324 = sitofp <2 x i32> %322 to <2 x double>
%325 = add nsw i64 %321, %284
%326 = getelementptr inbounds double, double* %61, i64 %325
%327 = bitcast double* %326 to <2 x double>*
store <2 x double> %323, <2 x double>* %327
%328 = getelementptr inbounds double, double* %326, i64 2
%329 = bitcast double* %328 to <2 x double>*
store <2 x double> %324, <2 x double>* %329
%330 = or i64 %318, 4
%331 = add <2 x i32> %319, <i32 4, i32 4>
%332 = add i64 %330, %283
%333 = add <2 x i32> %319, <i32 6, i32 6>
%334 = sitofp <2 x i32> %331 to <2 x double>
%335 = sitofp <2 x i32> %333 to <2 x double>
%336 = add nsw i64 %332, %284
%337 = getelementptr inbounds double, double* %61, i64 %336
%338 = bitcast double* %337 to <2 x double>*
store <2 x double> %334, <2 x double>* %338
%339 = getelementptr inbounds double, double* %337, i64 2
%340 = bitcast double* %339 to <2 x double>*
store <2 x double> %335, <2 x double>* %340
%341 = add nuw i64 %318, 8
%342 = add <2 x i32> %319, <i32 8, i32 8>
%343 = add i64 %320, 2
%344 = icmp eq i64 %343, %316
br i1 %344, label %287, label %317
345: ; preds = %349, %279
%346 = icmp sgt i32 %45, 0
br i1 %346, label %347, label %358
347: ; preds = %345
%348 = zext i32 %45 to i64
br label %363
349: ; preds = %301, %349
%350 = phi i64 [ %355, %349 ], [ %302, %301 ]
%351 = trunc i64 %350 to i32
%352 = sitofp i32 %351 to double
%353 = add nsw i64 %350, %284
%354 = getelementptr inbounds double, double* %61, i64 %353
store double %352, double* %354
%355 = add nsw i64 %350, 1
%356 = trunc i64 %355 to i32
%357 = icmp eq i32 %285, %356
br i1 %357, label %345, label %349
358: ; preds = %363, %345
%359 = phi double [ 0.000000e+00, %345 ], [ %375, %363 ]
%360 = icmp sgt i32 %51, 0
br i1 %360, label %361, label %380
361: ; preds = %358
%362 = zext i32 %51 to i64
br label %385
363: ; preds = %378, %347
%364 = phi double [ %16, %347 ], [ %379, %378 ]
%365 = phi i64 [ 0, %347 ], [ %376, %378 ]
%366 = phi double [ 0.000000e+00, %347 ], [ %375, %378 ]
%367 = getelementptr inbounds double, double* %49, i64 %365
%368 = load double, double* %367
%369 = fneg double %368
%370 = fmul double %368, %369
%371 = fmul double %364, %364
%372 = fmul double %371, 2.000000e+00
%373 = fdiv double %370, %372
%374 = tail call double @exp(double %373)
store double %374, double* %367
%375 = fadd double %366, %374
%376 = add nuw nsw i64 %365, 1
%377 = icmp eq i64 %376, %348
br i1 %377, label %358, label %378
378: ; preds = %363
%379 = load double, double* %3
br label %363
380: ; preds = %385, %358
%381 = phi double [ 0.000000e+00, %358 ], [ %397, %385 ]
%382 = icmp sgt i32 %57, 0
br i1 %382, label %383, label %400
383: ; preds = %380
%384 = zext i32 %57 to i64
br label %444
385: ; preds = %361, %385
%386 = phi i64 [ 0, %361 ], [ %398, %385 ]
%387 = phi double [ 0.000000e+00, %361 ], [ %397, %385 ]
%388 = getelementptr inbounds double, double* %55, i64 %386
%389 = load double, double* %388
%390 = fneg double %389
%391 = fmul double %389, %390
%392 = load double, double* %22
%393 = fmul double %392, %392
%394 = fmul double %393, 2.000000e+00
%395 = fdiv double %391, %394
%396 = tail call double @exp(double %395)
store double %396, double* %388
%397 = fadd double %387, %396
%398 = add nuw nsw i64 %386, 1
%399 = icmp eq i64 %398, %362
br i1 %399, label %380, label %385
400: ; preds = %444, %380
%401 = phi double [ 0.000000e+00, %380 ], [ %456, %444 ]
br i1 %346, label %402, label %459
402: ; preds = %400
%403 = zext i32 %45 to i64
%404 = icmp eq i32 %19, 0
br i1 %404, label %414, label %416
405: ; preds = %427, %416
%406 = phi i64 [ 0, %416 ], [ %441, %427 ]
%407 = icmp eq i64 %423, 0
br i1 %407, label %414, label %408
408: ; preds = %405
%409 = getelementptr inbounds double, double* %49, i64 %406
%410 = bitcast double* %409 to <2 x double>*
%411 = load <2 x double>, <2 x double>* %410
%412 = fdiv <2 x double> %411, %419
%413 = bitcast double* %409 to <2 x double>*
store <2 x double> %412, <2 x double>* %413
br label %414
414: ; preds = %408, %405, %402
%415 = phi i64 [ 0, %402 ], [ %417, %405 ], [ %417, %408 ]
br label %502
416: ; preds = %402
%417 = add nsw i64 %403, -1
%418 = insertelement <2 x double> poison, double %359, i64 0
%419 = shufflevector <2 x double> %418, <2 x double> poison, <2 x i32> zeroinitializer
%420 = add nsw i64 %403, -3
%421 = lshr exact i64 %420, 1
%422 = add nuw i64 %421, 1
%423 = and i64 %422, 1
%424 = icmp eq i64 %420, 0
br i1 %424, label %405, label %425
425: ; preds = %416
%426 = and i64 %422, -2
br label %427
427: ; preds = %427, %425
%428 = phi i64 [ 0, %425 ], [ %441, %427 ]
%429 = phi i64 [ 0, %425 ], [ %442, %427 ]
%430 = getelementptr inbounds double, double* %49, i64 %428
%431 = bitcast double* %430 to <2 x double>*
%432 = load <2 x double>, <2 x double>* %431
%433 = fdiv <2 x double> %432, %419
%434 = bitcast double* %430 to <2 x double>*
store <2 x double> %433, <2 x double>* %434
%435 = or i64 %428, 2
%436 = getelementptr inbounds double, double* %49, i64 %435
%437 = bitcast double* %436 to <2 x double>*
%438 = load <2 x double>, <2 x double>* %437
%439 = fdiv <2 x double> %438, %419
%440 = bitcast double* %436 to <2 x double>*
store <2 x double> %439, <2 x double>* %440
%441 = add nuw i64 %428, 4
%442 = add nuw i64 %429, 2
%443 = icmp eq i64 %442, %426
br i1 %443, label %405, label %427
444: ; preds = %383, %444
%445 = phi i64 [ 0, %383 ], [ %457, %444 ]
%446 = phi double [ 0.000000e+00, %383 ], [ %456, %444 ]
%447 = getelementptr inbounds double, double* %61, i64 %445
%448 = load double, double* %447
%449 = fneg double %448
%450 = fmul double %448, %449
%451 = load double, double* %33
%452 = fmul double %451, %451
%453 = fmul double %452, 2.000000e+00
%454 = fdiv double %450, %453
%455 = tail call double @exp(double %454)
store double %455, double* %447
%456 = fadd double %446, %455
%457 = add nuw nsw i64 %445, 1
%458 = icmp eq i64 %457, %384
br i1 %458, label %400, label %444
459: ; preds = %502, %400
br i1 %360, label %460, label %509
460: ; preds = %459
%461 = zext i32 %51 to i64
%462 = icmp eq i32 %30, 0
br i1 %462, label %472, label %474
463: ; preds = %485, %474
%464 = phi i64 [ 0, %474 ], [ %499, %485 ]
%465 = icmp eq i64 %481, 0
br i1 %465, label %472, label %466
466: ; preds = %463
%467 = getelementptr inbounds double, double* %55, i64 %464
%468 = bitcast double* %467 to <2 x double>*
%469 = load <2 x double>, <2 x double>* %468
%470 = fdiv <2 x double> %469, %477
%471 = bitcast double* %467 to <2 x double>*
store <2 x double> %470, <2 x double>* %471
br label %472
472: ; preds = %466, %463, %460
%473 = phi i64 [ 0, %460 ], [ %475, %463 ], [ %475, %466 ]
br label %552
474: ; preds = %460
%475 = add nsw i64 %461, -1
%476 = insertelement <2 x double> poison, double %381, i64 0
%477 = shufflevector <2 x double> %476, <2 x double> poison, <2 x i32> zeroinitializer
%478 = add nsw i64 %461, -3
%479 = lshr exact i64 %478, 1
%480 = add nuw i64 %479, 1
%481 = and i64 %480, 1
%482 = icmp eq i64 %478, 0
br i1 %482, label %463, label %483
483: ; preds = %474
%484 = and i64 %480, -2
br label %485
485: ; preds = %485, %483
%486 = phi i64 [ 0, %483 ], [ %499, %485 ]
%487 = phi i64 [ 0, %483 ], [ %500, %485 ]
%488 = getelementptr inbounds double, double* %55, i64 %486
%489 = bitcast double* %488 to <2 x double>*
%490 = load <2 x double>, <2 x double>* %489
%491 = fdiv <2 x double> %490, %477
%492 = bitcast double* %488 to <2 x double>*
store <2 x double> %491, <2 x double>* %492
%493 = or i64 %486, 2
%494 = getelementptr inbounds double, double* %55, i64 %493
%495 = bitcast double* %494 to <2 x double>*
%496 = load <2 x double>, <2 x double>* %495
%497 = fdiv <2 x double> %496, %477
%498 = bitcast double* %494 to <2 x double>*
store <2 x double> %497, <2 x double>* %498
%499 = add nuw i64 %486, 4
%500 = add nuw i64 %487, 2
%501 = icmp eq i64 %500, %484
br i1 %501, label %463, label %485
502: ; preds = %414, %502
%503 = phi i64 [ %507, %502 ], [ %415, %414 ]
%504 = getelementptr inbounds double, double* %49, i64 %503
%505 = load double, double* %504
%506 = fdiv double %505, %359
store double %506, double* %504
%507 = add nuw nsw i64 %503, 1
%508 = icmp eq i64 %507, %403
br i1 %508, label %459, label %502
509: ; preds = %552, %459
br i1 %382, label %510, label %566
510: ; preds = %509
%511 = zext i32 %57 to i64
%512 = icmp eq i32 %41, 0
br i1 %512, label %522, label %524
513: ; preds = %535, %524
%514 = phi i64 [ 0, %524 ], [ %549, %535 ]
%515 = icmp eq i64 %531, 0
br i1 %515, label %522, label %516
516: ; preds = %513
%517 = getelementptr inbounds double, double* %61, i64 %514
%518 = bitcast double* %517 to <2 x double>*
%519 = load <2 x double>, <2 x double>* %518
%520 = fdiv <2 x double> %519, %527
%521 = bitcast double* %517 to <2 x double>*
store <2 x double> %520, <2 x double>* %521
br label %522
522: ; preds = %516, %513, %510
%523 = phi i64 [ 0, %510 ], [ %525, %513 ], [ %525, %516 ]
br label %559
524: ; preds = %510
%525 = add nsw i64 %511, -1
%526 = insertelement <2 x double> poison, double %401, i64 0
%527 = shufflevector <2 x double> %526, <2 x double> poison, <2 x i32> zeroinitializer
%528 = add nsw i64 %511, -3
%529 = lshr exact i64 %528, 1
%530 = add nuw i64 %529, 1
%531 = and i64 %530, 1
%532 = icmp eq i64 %528, 0
br i1 %532, label %513, label %533
533: ; preds = %524
%534 = and i64 %530, -2
br label %535
535: ; preds = %535, %533
%536 = phi i64 [ 0, %533 ], [ %549, %535 ]
%537 = phi i64 [ 0, %533 ], [ %550, %535 ]
%538 = getelementptr inbounds double, double* %61, i64 %536
%539 = bitcast double* %538 to <2 x double>*
%540 = load <2 x double>, <2 x double>* %539
%541 = fdiv <2 x double> %540, %527
%542 = bitcast double* %538 to <2 x double>*
store <2 x double> %541, <2 x double>* %542
%543 = or i64 %536, 2
%544 = getelementptr inbounds double, double* %61, i64 %543
%545 = bitcast double* %544 to <2 x double>*
%546 = load <2 x double>, <2 x double>* %545
%547 = fdiv <2 x double> %546, %527
%548 = bitcast double* %544 to <2 x double>*
store <2 x double> %547, <2 x double>* %548
%549 = add nuw i64 %536, 4
%550 = add nuw i64 %537, 2
%551 = icmp eq i64 %550, %534
br i1 %551, label %513, label %535
552: ; preds = %472, %552
%553 = phi i64 [ %557, %552 ], [ %473, %472 ]
%554 = getelementptr inbounds double, double* %55, i64 %553
%555 = load double, double* %554
%556 = fdiv double %555, %381
store double %556, double* %554
%557 = add nuw nsw i64 %553, 1
%558 = icmp eq i64 %557, %461
br i1 %558, label %509, label %552
559: ; preds = %522, %559
%560 = phi i64 [ %564, %559 ], [ %523, %522 ]
%561 = getelementptr inbounds double, double* %61, i64 %560
%562 = load double, double* %561
%563 = fdiv double %562, %401
store double %563, double* %561
%564 = add nuw nsw i64 %560, 1
%565 = icmp eq i64 %564, %511
br i1 %565, label %566, label %559
566: ; preds = %559, %509
%567 = tail call i32 (double*, double*, double*, double*, i32, i32, i32, i32, i32, i32, double*, i32*, ...) bitcast (i32 (...)* @convxyz_double to i32 (double*, double*, double*, double*, i32, i32, i32, i32, i32, i32, double*, i32*, ...)*)(double* %0, double* %49, double* %55, double* %61, i32 %45, i32 %51, i32 %57, i32 %130, i32 %205, i32 %280, double* %0, i32* %1)
br i1 %62, label %593, label %568
568: ; preds = %566
%569 = icmp sgt i32 %127, 0
br i1 %569, label %570, label %591
570: ; preds = %568
%571 = tail call i32 (double*, double*, double*, double*, i32, i32, i32, i32, i32, i32, double*, i32*, ...) bitcast (i32 (...)* @convxyz_double to i32 (double*, double*, double*, double*, i32, i32, i32, i32, i32, i32, double*, i32*, ...)*)(double* %128, double* %49, double* %55, double* %61, i32 %45, i32 %51, i32 %57, i32 %130, i32 %205, i32 %280, double* %128, i32* %1)
%572 = icmp sgt i32 %43, 0
br i1 %572, label %573, label %591
573: ; preds = %570
%574 = zext i32 %43 to i64
br label %575
575: ; preds = %573, %588
%576 = phi i64 [ 0, %573 ], [ %589, %588 ]
%577 = getelementptr inbounds i8, i8* %129, i64 %576
%578 = load i8, i8* %577
%579 = icmp eq i8 %578, 0
br i1 %579, label %586, label %580
580: ; preds = %575
%581 = getelementptr inbounds double, double* %128, i64 %576
%582 = load double, double* %581
%583 = getelementptr inbounds double, double* %0, i64 %576
%584 = load double, double* %583
%585 = fdiv double %584, %582
store double %585, double* %583
br label %588
586: ; preds = %575
%587 = getelementptr inbounds double, double* %0, i64 %576
store double 0.000000e+00, double* %587
br label %588
588: ; preds = %580, %586
%589 = add nuw nsw i64 %576, 1
%590 = icmp eq i64 %589, %574
br i1 %590, label %591, label %575
591: ; preds = %588, %570, %568
%592 = bitcast double* %128 to i8*
tail call void @free(i8* %592)
tail call void @free(i8* %129)
br label %593
593: ; preds = %591, %566
tail call void @free(i8* %48)
tail call void @free(i8* %54)
tail call void @free(i8* %60)
ret void
}