#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 */
scalar_t__ floor(double);
extern double cubic_bspline_2d_interpol_dx(double *coeffs, long width,
                                           long height, double x, double y) {
  int spline_degree = 3;
  long x_index[3], y_index[4];
  double x_weight[3], y_weight[4];
  double interpolated;
  double w;
  long width2 = 2L * width - 2L, height2 = 2L * height - 2L;
  long i, j, k;
  i = (long)floor(x + 1) - (spline_degree - 1) / 2L;
  j = (long)floor(y) - spline_degree / 2L;
  for (k = 0L; k <= spline_degree; k++) {
    if (k < spline_degree) {
      x_index[k] = i++;
    }
    y_index[k] = j++;
  }
  w = x + 0.5 - (double)x_index[1];
  x_weight[1] = 3.0 / 4.0 - w * w;
  x_weight[2] = (1.0 / 2.0) * (w - x_weight[1] + 1.0);
  x_weight[0] = 1.0 - x_weight[1] - x_weight[2];
  w = y - (double)y_index[1];
  y_weight[3] = (1.0 / 6.0) * w * w * w;
  y_weight[0] = (1.0 / 6.0) + (1.0 / 2.0) * w * (w - 1.0) - y_weight[3];
  y_weight[2] = w + y_weight[0] - 2.0 * y_weight[3];
  y_weight[1] = 1.0 - y_weight[0] - y_weight[2] - y_weight[3];
  for (k = 0L; k <= spline_degree; k++) {
    if (k < spline_degree) {
      x_index[k] =
          (width == 1L)
              ? (0L)
              : ((x_index[k] < 0L)
                     ? (-x_index[k] - width2 * ((-x_index[k]) / width2))
                     : (x_index[k] - width2 * (x_index[k] / width2)));
      if (width <= x_index[k]) {
        x_index[k] = width2 - x_index[k];
      }
    }
    y_index[k] =
        (height == 1L)
            ? (0L)
            : ((y_index[k] < 0L)
                   ? (-y_index[k] - height2 * ((-y_index[k]) / height2))
                   : (y_index[k] - height2 * (y_index[k] / height2)));
    if (height <= y_index[k]) {
      y_index[k] = height2 - y_index[k];
    }
  }
  interpolated = 0.0;
  for (i = 0L; i < spline_degree; i++) {
    for (j = 0L; j <= spline_degree; j++) {
      interpolated += (coeffs[x_index[i] * width + y_index[j]] -
                       coeffs[(x_index[i] - 1) * width + y_index[j]]) *
                      x_weight[i] * y_weight[j];
    }
  }
  return (interpolated);
}


x86 assembly (clang, -O0):
<cubic_bspline_2d_interpol_dx>:
push   %rbp
mov    %rsp,%rbp
sub    $0x140,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movl   $0x3,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
shl    $0x1,%rax
sub    $0x2,%rax
mov    %rax,-0xc8(%rbp)
mov    -0x18(%rbp),%rax
shl    $0x1,%rax
sub    $0x2,%rax
mov    %rax,-0xd0(%rbp)
movsd  0xe9a(%rip),%xmm0
addsd  -0x20(%rbp),%xmm0
call   1030 <floor@plt>
mov    %rax,-0xf8(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
cltq
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0xf8(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0xd8(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1030 <floor@plt>
mov    %rax,-0xf0(%rbp)
movslq -0x2c(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0xf0(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0xe0(%rbp)
movq   $0x0,-0xe8(%rbp)
mov    -0xe8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jg     125c <cubic_bspline_2d_interpol_dx+0x14c>
mov    -0xe8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jge    1224 <cubic_bspline_2d_interpol_dx+0x114>
mov    -0xd8(%rbp),%rcx
mov    %rcx,%rax
add    $0x1,%rax
mov    %rax,-0xd8(%rbp)
mov    -0xe8(%rbp),%rax
mov    %rcx,-0x50(%rbp,%rax,8)
mov    -0xe0(%rbp),%rcx
mov    %rcx,%rax
add    $0x1,%rax
mov    %rax,-0xe0(%rbp)
mov    -0xe8(%rbp),%rax
mov    %rcx,-0x70(%rbp,%rax,8)
mov    -0xe8(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0xe8(%rbp)
jmp    11db <cubic_bspline_2d_interpol_dx+0xcb>
movsd  0xdac(%rip),%xmm0
addsd  -0x20(%rbp),%xmm0
cvtsi2sdq -0x48(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xc0(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xd6a(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0xc0(%rbp),%xmm1
ff
subsd  -0x88(%rbp),%xmm1
ff
movsd  0xd2e(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xd32(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  0xd11(%rip),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
cvtsi2sdq -0x68(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  0xce5(%rip),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
mulsd  -0xc0(%rbp),%xmm0
ff
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xcc5(%rip),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xc0(%rbp),%xmm1
ff
movsd  0xc9d(%rip),%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xc95(%rip),%xmm1
addsd  %xmm1,%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
addsd  -0xb0(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm1
ff
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  0xc49(%rip),%xmm0
subsd  -0xb0(%rbp),%xmm0
ff
subsd  -0xa0(%rbp),%xmm0
ff
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movq   $0x0,-0xe8(%rbp)
mov    -0xe8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jg     165d <cubic_bspline_2d_interpol_dx+0x54d>
mov    -0xe8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jge    1528 <cubic_bspline_2d_interpol_dx+0x418>
cmpq   $0x1,-0x10(%rbp)
jne    1423 <cubic_bspline_2d_interpol_dx+0x313>
xor    %eax,%eax
mov    %rax,-0x100(%rbp)
jmp    14db <cubic_bspline_2d_interpol_dx+0x3cb>
mov    -0xe8(%rbp),%rax
cmpq   $0x0,-0x50(%rbp,%rax,8)
jge    1486 <cubic_bspline_2d_interpol_dx+0x376>
mov    -0xe8(%rbp),%rcx
xor    %eax,%eax
sub    -0x50(%rbp,%rcx,8),%rax
mov    %rax,-0x110(%rbp)
mov    -0xc8(%rbp),%rcx
mov    -0xe8(%rbp),%rdx
xor    %eax,%eax
sub    -0x50(%rbp,%rdx,8),%rax
cqto
idivq  -0xc8(%rbp)
mov    %rax,%rdx
mov    -0x110(%rbp),%rax
imul   %rdx,%rcx
sub    %rcx,%rax
mov    %rax,-0x108(%rbp)
jmp    14cd <cubic_bspline_2d_interpol_dx+0x3bd>
mov    -0xe8(%rbp),%rax
mov    -0x50(%rbp,%rax,8),%rax
mov    %rax,-0x118(%rbp)
mov    -0xc8(%rbp),%rcx
mov    -0xe8(%rbp),%rax
mov    -0x50(%rbp,%rax,8),%rax
cqto
idivq  -0xc8(%rbp)
mov    %rax,%rdx
mov    -0x118(%rbp),%rax
imul   %rdx,%rcx
sub    %rcx,%rax
mov    %rax,-0x108(%rbp)
mov    -0x108(%rbp),%rax
mov    %rax,-0x100(%rbp)
mov    -0x100(%rbp),%rcx
mov    -0xe8(%rbp),%rax
mov    %rcx,-0x50(%rbp,%rax,8)
mov    -0x10(%rbp),%rax
mov    -0xe8(%rbp),%rcx
cmp    -0x50(%rbp,%rcx,8),%rax
jg     1523 <cubic_bspline_2d_interpol_dx+0x413>
mov    -0xc8(%rbp),%rcx
mov    -0xe8(%rbp),%rax
sub    -0x50(%rbp,%rax,8),%rcx
mov    -0xe8(%rbp),%rax
mov    %rcx,-0x50(%rbp,%rax,8)
jmp    1528 <cubic_bspline_2d_interpol_dx+0x418>
cmpq   $0x1,-0x18(%rbp)
jne    1541 <cubic_bspline_2d_interpol_dx+0x431>
xor    %eax,%eax
mov    %rax,-0x120(%rbp)
jmp    15f9 <cubic_bspline_2d_interpol_dx+0x4e9>
mov    -0xe8(%rbp),%rax
cmpq   $0x0,-0x70(%rbp,%rax,8)
jge    15a4 <cubic_bspline_2d_interpol_dx+0x494>
mov    -0xe8(%rbp),%rcx
xor    %eax,%eax
sub    -0x70(%rbp,%rcx,8),%rax
mov    %rax,-0x130(%rbp)
mov    -0xd0(%rbp),%rcx
mov    -0xe8(%rbp),%rdx
xor    %eax,%eax
sub    -0x70(%rbp,%rdx,8),%rax
cqto
idivq  -0xd0(%rbp)
mov    %rax,%rdx
mov    -0x130(%rbp),%rax
imul   %rdx,%rcx
sub    %rcx,%rax
mov    %rax,-0x128(%rbp)
jmp    15eb <cubic_bspline_2d_interpol_dx+0x4db>
mov    -0xe8(%rbp),%rax
mov    -0x70(%rbp,%rax,8),%rax
mov    %rax,-0x138(%rbp)
mov    -0xd0(%rbp),%rcx
mov    -0xe8(%rbp),%rax
mov    -0x70(%rbp,%rax,8),%rax
cqto
idivq  -0xd0(%rbp)
mov    %rax,%rdx
mov    -0x138(%rbp),%rax
imul   %rdx,%rcx
sub    %rcx,%rax
mov    %rax,-0x128(%rbp)
mov    -0x128(%rbp),%rax
mov    %rax,-0x120(%rbp)
mov    -0x120(%rbp),%rcx
mov    -0xe8(%rbp),%rax
mov    %rcx,-0x70(%rbp,%rax,8)
mov    -0x18(%rbp),%rax
mov    -0xe8(%rbp),%rcx
cmp    -0x70(%rbp,%rcx,8),%rax
jg     1641 <cubic_bspline_2d_interpol_dx+0x531>
mov    -0xd0(%rbp),%rcx
mov    -0xe8(%rbp),%rax
sub    -0x70(%rbp,%rax,8),%rcx
mov    -0xe8(%rbp),%rax
mov    %rcx,-0x70(%rbp,%rax,8)
jmp    1646 <cubic_bspline_2d_interpol_dx+0x536>
mov    -0xe8(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0xe8(%rbp)
jmp    13e2 <cubic_bspline_2d_interpol_dx+0x2d2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movq   $0x0,-0xd8(%rbp)
mov    -0xd8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jge    1761 <cubic_bspline_2d_interpol_dx+0x651>
movq   $0x0,-0xe0(%rbp)
mov    -0xe0(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jg     1745 <cubic_bspline_2d_interpol_dx+0x635>
mov    -0x8(%rbp),%rax
mov    -0xd8(%rbp),%rcx
mov    -0x50(%rbp,%rcx,8),%rcx
imul   -0x10(%rbp),%rcx
mov    -0xe0(%rbp),%rdx
add    -0x70(%rbp,%rdx,8),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xd8(%rbp),%rcx
mov    -0x50(%rbp,%rcx,8),%rcx
sub    $0x1,%rcx
imul   -0x10(%rbp),%rcx
mov    -0xe0(%rbp),%rdx
add    -0x70(%rbp,%rdx,8),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0xd8(%rbp),%rax
mulsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
mov    -0xe0(%rbp),%rax
movsd  -0xb0(%rbp,%rax,8),%xmm2
ff ff
movsd  -0xb8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0xe0(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0xe0(%rbp)
jmp    1692 <cubic_bspline_2d_interpol_dx+0x582>
jmp    174a <cubic_bspline_2d_interpol_dx+0x63a>
mov    -0xd8(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0xd8(%rbp)
jmp    1673 <cubic_bspline_2d_interpol_dx+0x563>
movsd  -0xb8(%rbp),%xmm0
ff
add    $0x140,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cubic_bspline_2d_interpol_dx>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
movsd  %xmm1,(%rsp)
movsd  %xmm0,0x8(%rsp)
mov    %rdx,%r15
mov    %rsi,%r14
mov    %rdi,%rbx
addsd  0xec7(%rip),%xmm0
call   1030 <floor@plt>
mov    %rax,%r12
movsd  (%rsp),%xmm0
call   1030 <floor@plt>
movsd  (%rsp),%xmm1
movsd  0x8(%rsp),%xmm7
mov    %rax,%r8
lea    (%r14,%r14,1),%r11
add    $0xfffffffffffffffe,%r11
movsd  0xe9f(%rip),%xmm4
addsd  %xmm4,%xmm7
xorps  %xmm0,%xmm0
cvtsi2sd %r12,%xmm0
subsd  %xmm0,%xmm7
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0xe7f(%rip),%xmm2
subsd  %xmm0,%xmm1
movsd  0xe7b(%rip),%xmm5
movapd %xmm1,%xmm8
mulsd  %xmm5,%xmm8
mulsd  %xmm1,%xmm8
mulsd  %xmm1,%xmm8
movapd %xmm1,%xmm6
mulsd  %xmm4,%xmm6
movsd  0xe5f(%rip),%xmm0
addsd  %xmm1,%xmm0
mulsd  %xmm6,%xmm0
addsd  %xmm5,%xmm0
cmp    $0x1,%r14
jne    11df <cubic_bspline_2d_interpol_dx+0xcf>
xor    %eax,%eax
movsd  0xe23(%rip),%xmm6
jmp    1211 <cubic_bspline_2d_interpol_dx+0x101>
lea    -0x1(%r12),%rax
mov    %rax,%rcx
or     %r11,%rcx
shr    $0x20,%rcx
je     11f7 <cubic_bspline_2d_interpol_dx+0xe7>
cqto
idiv   %r11
jmp    11fc <cubic_bspline_2d_interpol_dx+0xec>
xor    %edx,%edx
div    %r11d
movsd  0xdfc(%rip),%xmm6
mov    %rdx,%rax
neg    %rax
test   %r12,%r12
cmovg  %rdx,%rax
subsd  %xmm3,%xmm2
subsd  %xmm8,%xmm0
lea    (%r15,%r15,1),%r10
add    $0xfffffffffffffffe,%r10
mov    %r11,%r13
sub    %rax,%r13
cmp    %r14,%rax
cmovl  %rax,%r13
cmp    $0x1,%r15
jne    1239 <cubic_bspline_2d_interpol_dx+0x129>
xor    %eax,%eax
jmp    1262 <cubic_bspline_2d_interpol_dx+0x152>
lea    -0x1(%r8),%rax
mov    %rax,%rcx
or     %r10,%rcx
shr    $0x20,%rcx
je     1250 <cubic_bspline_2d_interpol_dx+0x140>
cqto
idiv   %r10
jmp    1255 <cubic_bspline_2d_interpol_dx+0x145>
xor    %edx,%edx
div    %r10d
mov    %rdx,%rax
neg    %rax
test   %r8,%r8
cmovg  %rdx,%rax
subsd  %xmm2,%xmm7
addsd  %xmm0,%xmm1
movapd %xmm8,%xmm3
addsd  %xmm8,%xmm3
mov    %r10,%r9
sub    %rax,%r9
cmp    %r15,%rax
cmovl  %rax,%r9
cmp    $0x1,%r14
jne    128b <cubic_bspline_2d_interpol_dx+0x17b>
xor    %eax,%eax
jmp    12b6 <cubic_bspline_2d_interpol_dx+0x1a6>
mov    %r12,%rax
or     %r11,%rax
shr    $0x20,%rax
je     12a1 <cubic_bspline_2d_interpol_dx+0x191>
mov    %r12,%rax
cqto
idiv   %r11
jmp    12a9 <cubic_bspline_2d_interpol_dx+0x199>
mov    %r12d,%eax
xor    %edx,%edx
div    %r11d
mov    %rdx,%rax
neg    %rax
test   %r12,%r12
cmovns %rdx,%rax
addsd  %xmm6,%xmm7
subsd  %xmm3,%xmm1
movapd %xmm6,%xmm3
subsd  %xmm0,%xmm3
mov    %r11,%rcx
sub    %rax,%rcx
cmp    %r14,%rax
cmovl  %rax,%rcx
cmp    $0x1,%r15
jne    12dd <cubic_bspline_2d_interpol_dx+0x1cd>
xor    %eax,%eax
jmp    1308 <cubic_bspline_2d_interpol_dx+0x1f8>
mov    %r8,%rax
or     %r10,%rax
shr    $0x20,%rax
je     12f3 <cubic_bspline_2d_interpol_dx+0x1e3>
mov    %r8,%rax
cqto
idiv   %r10
jmp    12fb <cubic_bspline_2d_interpol_dx+0x1eb>
mov    %r8d,%eax
xor    %edx,%edx
div    %r10d
mov    %rdx,%rax
neg    %rax
test   %r8,%r8
cmovns %rdx,%rax
mulsd  %xmm4,%xmm7
subsd  %xmm2,%xmm6
subsd  %xmm1,%xmm3
mov    %r10,%rdi
sub    %rax,%rdi
cmp    %r15,%rax
cmovl  %rax,%rdi
cmp    $0x1,%r14
jne    132b <cubic_bspline_2d_interpol_dx+0x21b>
xor    %eax,%eax
jmp    1356 <cubic_bspline_2d_interpol_dx+0x246>
lea    0x1(%r12),%rax
mov    %rax,%rdx
or     %r11,%rdx
shr    $0x20,%rdx
je     1343 <cubic_bspline_2d_interpol_dx+0x233>
cqto
idiv   %r11
jmp    1348 <cubic_bspline_2d_interpol_dx+0x238>
xor    %edx,%edx
div    %r11d
mov    %rdx,%rax
neg    %rax
cmp    $0xffffffffffffffff,%r12
cmovge %rdx,%rax
subsd  %xmm7,%xmm6
subsd  %xmm8,%xmm3
sub    %rax,%r11
cmp    %r14,%rax
cmovl  %rax,%r11
cmp    $0x1,%r15
jne    1378 <cubic_bspline_2d_interpol_dx+0x268>
xor    %esi,%esi
xor    %edx,%edx
jmp    13f9 <cubic_bspline_2d_interpol_dx+0x2e9>
lea    0x1(%r8),%rax
mov    %rax,%rdx
or     %r10,%rdx
shr    $0x20,%rdx
je     1392 <cubic_bspline_2d_interpol_dx+0x282>
cqto
idiv   %r10
mov    %rdx,%rsi
jmp    1399 <cubic_bspline_2d_interpol_dx+0x289>
xor    %edx,%edx
div    %r10d
mov    %edx,%esi
lea    0x2(%r8),%rax
cmp    $0xffffffffffffffff,%r8
jge    13ca <cubic_bspline_2d_interpol_dx+0x2ba>
lea    (%r10,%rsi,1),%rdx
neg    %rsi
cmp    %r15,%rsi
cmovge %rdx,%rsi
cmp    $0xfffffffffffffffe,%r8
je     13d7 <cubic_bspline_2d_interpol_dx+0x2c7>
mov    %rax,%rdx
or     %r10,%rdx
shr    $0x20,%rdx
je     13f1 <cubic_bspline_2d_interpol_dx+0x2e1>
cqto
idiv   %r10
jmp    13f6 <cubic_bspline_2d_interpol_dx+0x2e6>
mov    %r10,%rdx
sub    %rsi,%rdx
cmp    %r15,%rsi
cmovge %rdx,%rsi
mov    %rax,%rdx
or     %r10,%rdx
shr    $0x20,%rdx
je     13ea <cubic_bspline_2d_interpol_dx+0x2da>
cqto
idiv   %r10
jmp    13f9 <cubic_bspline_2d_interpol_dx+0x2e9>
xor    %edx,%edx
div    %r10d
jmp    13f9 <cubic_bspline_2d_interpol_dx+0x2e9>
xor    %edx,%edx
div    %r10d
neg    %rdx
sub    %rdx,%r10
cmp    %r15,%rdx
cmovl  %rdx,%r10
mov    %r13,%rax
imul   %r14,%rax
add    $0xffffffffffffffff,%r13
imul   %r14,%r13
lea    (%rax,%r9,1),%rdx
movsd  (%rbx,%rdx,8),%xmm4
lea    (%r9,%r13,1),%rdx
subsd  (%rbx,%rdx,8),%xmm4
mulsd  %xmm6,%xmm4
mulsd  %xmm0,%xmm4
xorpd  %xmm5,%xmm5
addsd  %xmm4,%xmm5
lea    (%rax,%rdi,1),%rdx
movsd  (%rbx,%rdx,8),%xmm4
lea    0x0(%r13,%rdi,1),%rdx
subsd  (%rbx,%rdx,8),%xmm4
mulsd  %xmm6,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm5,%xmm4
lea    (%rax,%rsi,1),%rdx
movsd  (%rbx,%rdx,8),%xmm5
lea    0x0(%r13,%rsi,1),%rdx
subsd  (%rbx,%rdx,8),%xmm5
mulsd  %xmm6,%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm4,%xmm5
add    %r10,%rax
movsd  (%rbx,%rax,8),%xmm4
add    %r10,%r13
subsd  (%rbx,%r13,8),%xmm4
mulsd  %xmm6,%xmm4
mulsd  %xmm8,%xmm4
addsd  %xmm5,%xmm4
mov    %rcx,%rax
imul   %r14,%rax
add    $0xffffffffffffffff,%rcx
imul   %r14,%rcx
lea    (%rax,%r9,1),%rdx
movsd  (%rbx,%rdx,8),%xmm5
lea    (%rcx,%r9,1),%rdx
subsd  (%rbx,%rdx,8),%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
lea    (%rax,%rdi,1),%rdx
movsd  (%rbx,%rdx,8),%xmm6
lea    (%rcx,%rdi,1),%rdx
subsd  (%rbx,%rdx,8),%xmm6
mulsd  %xmm2,%xmm6
mulsd  %xmm3,%xmm6
addsd  %xmm5,%xmm6
lea    (%rax,%rsi,1),%rdx
movsd  (%rbx,%rdx,8),%xmm4
lea    (%rcx,%rsi,1),%rdx
subsd  (%rbx,%rdx,8),%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm6,%xmm4
add    %r10,%rax
movsd  (%rbx,%rax,8),%xmm5
add    %r10,%rcx
subsd  (%rbx,%rcx,8),%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm8,%xmm5
mov    %r11,%rax
imul   %r14,%rax
add    $0xffffffffffffffff,%r11
imul   %r14,%r11
lea    (%rax,%r9,1),%rcx
movsd  (%rbx,%rcx,8),%xmm2
add    %r11,%r9
subsd  (%rbx,%r9,8),%xmm2
addsd  %xmm4,%xmm5
mulsd  %xmm7,%xmm2
mulsd  %xmm0,%xmm2
lea    (%rax,%rdi,1),%rcx
movsd  (%rbx,%rcx,8),%xmm4
add    %r11,%rdi
subsd  (%rbx,%rdi,8),%xmm4
addsd  %xmm5,%xmm2
mulsd  %xmm7,%xmm4
mulsd  %xmm3,%xmm4
lea    (%rax,%rsi,1),%rcx
movsd  (%rbx,%rcx,8),%xmm3
add    %r11,%rsi
subsd  (%rbx,%rsi,8),%xmm3
addsd  %xmm2,%xmm4
mulsd  %xmm7,%xmm3
mulsd  %xmm1,%xmm3
add    %r10,%rax
movsd  (%rbx,%rax,8),%xmm0
add    %r10,%r11
subsd  (%rbx,%r11,8),%xmm0
addsd  %xmm4,%xmm3
mulsd  %xmm7,%xmm0
mulsd  %xmm8,%xmm0
addsd  %xmm3,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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