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

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double pow(double, int);
void cu_compute_vector_stamps(int dp, int ds, int db, int nx, int ny,
                              int nstamps, int stamp_half_width,
                              double *stamp_xpos, double *stamp_ypos,
                              int *kxindex, int *kyindex, int *ext_basis,
                              int nkernel, int kernelRadius, double *V,
                              int BlockDimx, int gridDimx, double *tex0,
                              double *tex1, double *tex2, double *tex3,
                              double *tex4) {
  int idx;
  int np, ns, ki, a, b, d1, i, j, i1, i2, j1, j2;
  int l, m, l1, m1;
  double py, x, y, Bi;
  double temp;
  int blockIdx;
  np = (dp + 1) * (dp + 2) / 2;
  ns = (ds + 1) * (ds + 2) / 2;
  for (blockIdx = 0; blockIdx < gridDimx; blockIdx++) {
    ki = blockIdx < np ? 0 : (blockIdx - np) / ns + 1;
    a = b = 0;
    if (ki < nkernel) {
      a = kxindex[ki];
      b = kyindex[ki];
    }
    if (ki == 0) {
      d1 = dp;
      idx = blockIdx;
    } else if (ki < nkernel) {
      d1 = ds;
      idx = blockIdx - np - (ki - 1) * ns;
    } else {
      d1 = db;
      idx = blockIdx - np - (ki - 1) * ns;
    }
    i = 0;
    for (l = 0; l <= d1; l++) {
      for (m = 0; m <= d1 - l; m++) {
        if (i == idx) {
          l1 = l;
          m1 = m;
        }
        i++;
      }
    }
    Bi = 1.0;
    V[blockIdx] = 0.0;
    for (idx = 0; idx < nstamps; idx++) {
      j1 = ((0) > ((int)stamp_ypos[idx] - stamp_half_width)
                ? (0)
                : ((int)stamp_ypos[idx] - stamp_half_width));
      j2 = ((ny) < ((int)stamp_ypos[idx] + stamp_half_width)
                ? (ny)
                : ((int)stamp_ypos[idx] + stamp_half_width));
      for (j = j1; j < j2; j++) {
        y = (j - 0.5 * (ny - 1)) / (ny - 1);
        py = pow(y, m1);
        i1 = ((0) > ((int)stamp_xpos[idx] - stamp_half_width)
                  ? (0)
                  : ((int)stamp_xpos[idx] - stamp_half_width));
        i2 = ((nx) < ((int)stamp_xpos[idx] + stamp_half_width)
                  ? (nx)
                  : ((int)stamp_xpos[idx] + stamp_half_width));
        for (i = i1; i < i2; i++) {
          x = (i - 0.5 * (nx - 1)) / (nx - 1);
          if (ki == 0) {
            Bi = tex0[i + nx * j];
          } else if (ki < nkernel) {
            if (ext_basis[ki]) {
              Bi = tex1[i + a + nx * (j + b)] - tex0[i + nx * j];
            } else {
              Bi = tex0[i + a + nx * (j + b)] - tex0[i + nx * j];
            }
          } else {
            Bi = 1.0;
          }
          V[blockIdx] += pow(x, l1) * py * Bi * tex2[i + nx * j] *
                         tex3[i + nx * j] * tex4[i + nx * j];
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<cu_compute_vector_stamps>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    0x88(%rbp),%rax
mov    0x80(%rbp),%rax
mov    0x78(%rbp),%rax
mov    0x70(%rbp),%rax
mov    0x68(%rbp),%rax
mov    0x60(%rbp),%eax
mov    0x58(%rbp),%eax
mov    0x50(%rbp),%rax
mov    0x48(%rbp),%eax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9d,-0x18(%rbp)
mov    -0x4(%rbp),%eax
add    $0x1,%eax
mov    -0x4(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    -0x8(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    0x60(%rbp),%eax
jge    16aa <cu_compute_vector_stamps+0x59a>
mov    -0x8c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11d9 <cu_compute_vector_stamps+0xc9>
xor    %eax,%eax
mov    %eax,-0x90(%rbp)
jmp    11ef <cu_compute_vector_stamps+0xdf>
mov    -0x8c(%rbp),%eax
sub    -0x20(%rbp),%eax
cltd
idivl  -0x24(%rbp)
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
mov    -0x90(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
cmp    0x40(%rbp),%eax
jge    122e <cu_compute_vector_stamps+0x11e>
mov    0x28(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    0x30(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x28(%rbp)
jne    124c <cu_compute_vector_stamps+0x13c>
mov    -0x4(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x8c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    129e <cu_compute_vector_stamps+0x18e>
mov    -0x28(%rbp),%eax
cmp    0x40(%rbp),%eax
jge    127b <cu_compute_vector_stamps+0x16b>
mov    -0x8(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x8c(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x24(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x1c(%rbp)
jmp    1299 <cu_compute_vector_stamps+0x189>
mov    -0xc(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x8c(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x24(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x1c(%rbp)
jmp    129e <cu_compute_vector_stamps+0x18e>
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     1312 <cu_compute_vector_stamps+0x202>
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    -0x50(%rbp),%ecx
cmp    %ecx,%eax
jg     12ff <cu_compute_vector_stamps+0x1ef>
mov    -0x38(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jne    12e8 <cu_compute_vector_stamps+0x1d8>
mov    -0x50(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    12bf <cu_compute_vector_stamps+0x1af>
jmp    1304 <cu_compute_vector_stamps+0x1f4>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    12ac <cu_compute_vector_stamps+0x19c>
movsd  0xce6(%rip),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1691 <cu_compute_vector_stamps+0x581>
mov    0x20(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%ecx
sub    0x10(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    136c <cu_compute_vector_stamps+0x25c>
xor    %eax,%eax
mov    %eax,-0x94(%rbp)
jmp    1382 <cu_compute_vector_stamps+0x272>
mov    0x20(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
sub    0x10(%rbp),%eax
mov    %eax,-0x94(%rbp)
mov    -0x94(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x14(%rbp),%eax
mov    0x20(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
cvttsd2si (%rcx,%rdx,8),%ecx
add    0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    13b4 <cu_compute_vector_stamps+0x2a4>
mov    -0x14(%rbp),%eax
mov    %eax,-0x98(%rbp)
jmp    13ca <cu_compute_vector_stamps+0x2ba>
mov    0x20(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
add    0x10(%rbp),%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jge    167e <cu_compute_vector_stamps+0x56e>
cvtsi2sdl -0x3c(%rbp),%xmm1
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xc0c(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
mov    -0x5c(%rbp),%edi
call   1030 <pow@plt>
movsd  %xmm0,-0x68(%rbp)
mov    0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%ecx
sub    0x10(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    1450 <cu_compute_vector_stamps+0x340>
xor    %eax,%eax
mov    %eax,-0x9c(%rbp)
jmp    1466 <cu_compute_vector_stamps+0x356>
mov    0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
sub    0x10(%rbp),%eax
mov    %eax,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x10(%rbp),%eax
mov    0x18(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
cvttsd2si (%rcx,%rdx,8),%ecx
add    0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    1498 <cu_compute_vector_stamps+0x388>
mov    -0x10(%rbp),%eax
mov    %eax,-0xa0(%rbp)
jmp    14ae <cu_compute_vector_stamps+0x39e>
mov    0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
add    0x10(%rbp),%eax
mov    %eax,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    166b <cu_compute_vector_stamps+0x55b>
cvtsi2sdl -0x38(%rbp),%xmm1
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xb28(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
cmpl   $0x0,-0x28(%rbp)
jne    1527 <cu_compute_vector_stamps+0x417>
mov    0x68(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    15db <cu_compute_vector_stamps+0x4cb>
mov    -0x28(%rbp),%eax
cmp    0x40(%rbp),%eax
jge    15c9 <cu_compute_vector_stamps+0x4b9>
mov    0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1587 <cu_compute_vector_stamps+0x477>
mov    0x70(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x2c(%rbp),%ecx
mov    -0x10(%rbp),%edx
mov    -0x3c(%rbp),%esi
add    -0x30(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    15c4 <cu_compute_vector_stamps+0x4b4>
mov    0x68(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x2c(%rbp),%ecx
mov    -0x10(%rbp),%edx
mov    -0x3c(%rbp),%esi
add    -0x30(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    15d6 <cu_compute_vector_stamps+0x4c6>
movsd  0xa2f(%rip),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    15db <cu_compute_vector_stamps+0x4cb>
movsd  -0x70(%rbp),%xmm0
mov    -0x58(%rbp),%edi
call   1030 <pow@plt>
mulsd  -0x68(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
mov    0x78(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x80(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x88(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x3c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x50(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    14bd <cu_compute_vector_stamps+0x3ad>
jmp    1670 <cu_compute_vector_stamps+0x560>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    13d9 <cu_compute_vector_stamps+0x2c9>
jmp    1683 <cu_compute_vector_stamps+0x573>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1339 <cu_compute_vector_stamps+0x229>
jmp    1696 <cu_compute_vector_stamps+0x586>
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    11ae <cu_compute_vector_stamps+0x9e>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cu_compute_vector_stamps>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %r9d,0x40(%rsp)
mov    %r8,0x68(%rsp)
mov    %rcx,0x70(%rsp)
mov    %edx,0x58(%rsp)
mov    %rsi,0x60(%rsp)
mov    %rdi,0x38(%rsp)
mov    0x190(%rsp),%eax
test   %eax,%eax
jle    17ad <cu_compute_vector_stamps+0x69d>
mov    0x38(%rsp),%rdx
lea    0x1(%rdx),%ecx
add    $0x2,%edx
imul   %ecx,%edx
mov    %edx,%edi
shr    $0x1f,%edi
add    %edx,%edi
mov    0x60(%rsp),%rdx
lea    0x1(%rdx),%ecx
add    $0x2,%edx
imul   %ecx,%edx
mov    0x68(%rsp),%rcx
add    $0xffffffff,%ecx
cvtsi2sd %ecx,%xmm0
mov    0x70(%rsp),%rsi
lea    -0x1(%rsi),%ecx
cvtsi2sd %ecx,%xmm1
mov    %edx,%ecx
shr    $0x1f,%ecx
add    %edx,%ecx
sar    %edi
sar    %ecx
mov    %ecx,0x28(%rsp)
mov    %eax,%eax
mov    %rax,0xa0(%rsp)
movsd  0xe57(%rip),%xmm2
movsd  %xmm0,0xc8(%rsp)
mulsd  %xmm0,%xmm2
movsd  %xmm2,0xc0(%rsp)
movsd  0xe39(%rip),%xmm0
movsd  %xmm1,0x100(%rsp)
mulsd  %xmm1,%xmm0
movsd  %xmm0,0xf8(%rsp)
movslq %esi,%rax
mov    %rax,0x48(%rsp)
mov    %edi,0x2c(%rsp)
movslq %edi,%rax
mov    %rax,0x98(%rsp)
mov    0x40(%rsp),%eax
mov    %rax,0xb8(%rsp)
xor    %ebp,%ebp
jmp    1222 <cu_compute_vector_stamps+0x112>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rbp
cmp    0xa0(%rsp),%rbp
je     17ad <cu_compute_vector_stamps+0x69d>
cmp    0x98(%rsp),%rbp
mov    0x44(%rsp),%edi
jge    1270 <cu_compute_vector_stamps+0x160>
cmpl   $0x0,0x170(%rsp)
jle    131f <cu_compute_vector_stamps+0x20f>
movl   $0x0,0x20(%rsp)
mov    0x158(%rsp),%rax
mov    (%rax),%eax
mov    %eax,0x18(%rsp)
mov    0x160(%rsp),%rax
mov    (%rax),%eax
mov    %eax,0x10(%rsp)
mov    $0x1,%al
mov    %eax,0xc(%rsp)
jmp    133f <cu_compute_vector_stamps+0x22f>
nopl   (%rax)
mov    %ebp,%eax
sub    0x2c(%rsp),%eax
cltd
idivl  0x28(%rsp)
lea    0x1(%rax),%ebx
cmp    0x170(%rsp),%ebx
jge    12e0 <cu_compute_vector_stamps+0x1d0>
movslq %ebx,%rsi
mov    0x158(%rsp),%rdx
mov    (%rdx,%rsi,4),%r8d
mov    0x160(%rsp),%rdx
mov    (%rdx,%rsi,4),%esi
mov    $0x1,%dl
mov    %edx,0xc(%rsp)
test   %ebx,%ebx
mov    %r8d,0x18(%rsp)
mov    %esi,0x10(%rsp)
je     1337 <cu_compute_vector_stamps+0x227>
mov    %ebx,0x20(%rsp)
imul   0x28(%rsp),%eax
add    0x2c(%rsp),%eax
mov    %ebp,%r13d
sub    %eax,%r13d
movl   $0x0,0x8(%rsp)
mov    0x60(%rsp),%r14
mov    %r14d,%edx
jmp    1378 <cu_compute_vector_stamps+0x268>
test   %ebx,%ebx
je     1347 <cu_compute_vector_stamps+0x237>
mov    %ebx,0x20(%rsp)
imul   0x28(%rsp),%eax
add    0x2c(%rsp),%eax
mov    %ebp,%r13d
sub    %eax,%r13d
movl   $0x0,0xc(%rsp)
movl   $0x0,0x8(%rsp)
movl   $0x0,0x18(%rsp)
movl   $0x0,0x10(%rsp)
mov    0x58(%rsp),%eax
mov    %eax,%edx
jmp    1378 <cu_compute_vector_stamps+0x268>
movl   $0x0,0x20(%rsp)
mov    $0x1,%al
mov    %eax,0x8(%rsp)
movl   $0x0,0xc(%rsp)
jmp    135d <cu_compute_vector_stamps+0x24d>
movl   $0x0,0x20(%rsp)
mov    $0x1,%al
mov    %eax,0x8(%rsp)
jmp    136d <cu_compute_vector_stamps+0x25d>
mov    $0x1,%al
mov    %eax,0x8(%rsp)
movl   $0x0,0xc(%rsp)
movl   $0x0,0x20(%rsp)
movl   $0x0,0x18(%rsp)
movl   $0x0,0x10(%rsp)
mov    0x38(%rsp),%r14
mov    %r14d,%edx
mov    %ebp,%r13d
test   %edx,%edx
mov    %rbp,0x90(%rsp)
js     1490 <cu_compute_vector_stamps+0x380>
lea    0x1(%rdx),%r11d
mov    %r13d,%r8d
neg    %r8d
xor    %r14d,%r14d
xor    %r15d,%r15d
mov    %rdx,0x30(%rsp)
mov    %r8,0x50(%rsp)
jmp    13c4 <cu_compute_vector_stamps+0x2b4>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
lea    0x1(%r14),%eax
add    $0xffffffff,%r11d
cmp    %r14d,%edx
mov    %eax,%r14d
je     1490 <cu_compute_vector_stamps+0x380>
mov    %edx,%eax
sub    %r14d,%eax
jl     13b0 <cu_compute_vector_stamps+0x2a0>
xor    %ebp,%ebp
cmp    $0x3,%eax
jb     1449 <cu_compute_vector_stamps+0x339>
mov    %r11d,%r12d
and    $0xfffffffc,%r12d
lea    (%r8,%r15,1),%r9d
xor    %ebp,%ebp
mov    %r15d,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    0x1(%rsi),%edx
lea    0x1(%rbp),%r10d
lea    0x2(%rsi),%r8d
mov    %r9d,%eax
add    %ebp,%eax
cmove  %ebp,%ecx
lea    0x2(%rbp),%eax
cmove  %r14d,%edi
cmp    %r13d,%edx
cmove  %r10d,%ecx
cmove  %r14d,%edi
lea    0x3(%rsi),%edx
cmp    %r13d,%r8d
mov    %ecx,%ebx
cmove  %eax,%ebx
cmove  %r14d,%edi
lea    0x3(%rbp),%ecx
cmp    %r13d,%edx
cmove  %r14d,%edi
cmovne %ebx,%ecx
add    $0x4,%esi
add    $0x4,%ebp
cmp    %r12d,%ebp
jne    13f0 <cu_compute_vector_stamps+0x2e0>
add    %ebp,%r15d
mov    0x30(%rsp),%rdx
mov    0x50(%rsp),%r8
test   $0x3,%r11b
je     13b0 <cu_compute_vector_stamps+0x2a0>
mov    %r11d,%eax
and    $0x3,%eax
neg    %eax
lea    (%r8,%r15,1),%ebx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
cmp    %esi,%ebx
cmove  %r14d,%edi
cmove  %ebp,%ecx
add    $0x1,%ebp
add    $0xffffffff,%esi
cmp    %esi,%eax
jne    1470 <cu_compute_vector_stamps+0x360>
sub    %esi,%r15d
jmp    13b0 <cu_compute_vector_stamps+0x2a0>
nopl   0x0(%rax,%rax,1)
mov    %edi,0x44(%rsp)
mov    0x180(%rsp),%rax
mov    0x90(%rsp),%rbp
movq   $0x0,(%rax,%rbp,8)
cmpl   $0x0,0x40(%rsp)
jle    1210 <cu_compute_vector_stamps+0x100>
movslq 0x20(%rsp),%rax
mov    %rax,0xf0(%rsp)
movslq 0x18(%rsp),%rax
mov    %rax,0xb0(%rsp)
movslq 0x10(%rsp),%rax
mov    %rax,0xa8(%rsp)
xor    %ebx,%ebx
mov    %ecx,0x5c(%rsp)
jmp    1502 <cu_compute_vector_stamps+0x3f2>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rbx
cmp    0xb8(%rsp),%rbx
je     1210 <cu_compute_vector_stamps+0x100>
mov    0x150(%rsp),%rax
cvttsd2si (%rax,%rbx,8),%eax
mov    %eax,%edi
mov    0x140(%rsp),%edx
sub    %edx,%edi
mov    $0x0,%esi
cmovs  %esi,%edi
add    %edx,%eax
mov    0x68(%rsp),%rdx
cmp    %edx,%eax
cmovg  %edx,%eax
cmp    %eax,%edi
jge    14f0 <cu_compute_vector_stamps+0x3e0>
mov    %edi,%esi
mov    %eax,%eax
mov    %rax,0xd0(%rsp)
mov    0x48(%rsp),%rax
mov    %rax,%rdx
imul   %rsi,%rdx
mov    %rdx,0x80(%rsp)
mov    0xa8(%rsp),%rdx
add    %rsi,%rdx
imul   %rax,%rdx
add    0xb0(%rsp),%rdx
mov    %rdx,0x78(%rsp)
mov    %rbx,0xd8(%rsp)
jmp    15b8 <cu_compute_vector_stamps+0x4a8>
nopl   0x0(%rax,%rax,1)
mov    0x88(%rsp),%rsi
add    $0x1,%rsi
mov    0x48(%rsp),%rax
add    %rax,0x80(%rsp)
add    %rax,0x78(%rsp)
cmp    0xd0(%rsp),%rsi
mov    0x5c(%rsp),%ecx
mov    0xd8(%rsp),%rbx
je     14f0 <cu_compute_vector_stamps+0x3e0>
mov    %rsi,0x88(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
addsd  0xc0(%rsp),%xmm0
divsd  0xc8(%rsp),%xmm0
mov    %ecx,%edi
call   1030 <pow@plt>
movsd  %xmm0,0x50(%rsp)
mov    0x148(%rsp),%rax
cvttsd2si (%rax,%rbx,8),%eax
mov    %eax,%ecx
mov    0x140(%rsp),%edx
sub    %edx,%ecx
mov    $0x0,%esi
cmovs  %esi,%ecx
add    %edx,%eax
mov    0x70(%rsp),%rdx
cmp    %edx,%eax
cmovg  %edx,%eax
cmp    %eax,%ecx
jge    1580 <cu_compute_vector_stamps+0x470>
mov    0x88(%rsp),%rdx
imul   0x48(%rsp),%rdx
mov    %rdx,0x18(%rsp)
mov    %ecx,%r14d
mov    %eax,%eax
mov    0x80(%rsp),%rcx
lea    (%rcx,%r14,1),%rdx
mov    0x198(%rsp),%rsi
mov    %rdx,0x20(%rsp)
lea    (%rsi,%rdx,8),%rcx
mov    %rcx,0x10(%rsp)
mov    0x78(%rsp),%rcx
add    %r14,%rcx
mov    0x1a0(%rsp),%rdx
lea    (%rdx,%rcx,8),%rdx
mov    %rdx,0xe8(%rsp)
lea    (%rsi,%rcx,8),%rcx
mov    %rcx,0xe0(%rsp)
mov    %r14,%r15
sub    %rax,%r15
xor    %r13d,%r13d
mov    %r14,%rbx
jmp    172c <cu_compute_vector_stamps+0x61c>
nopl   (%rax)
mov    0x18(%rsp),%rcx
lea    (%rbx,%rcx,1),%r12
mov    0x10(%rsp),%rcx
movsd  (%rcx,%r13,8),%xmm0
movsd  %xmm0,0x30(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
addsd  0xf8(%rsp),%xmm0
divsd  0x100(%rsp),%xmm0
call   1030 <pow@plt>
mulsd  0x50(%rsp),%xmm0
mulsd  0x30(%rsp),%xmm0
mov    0x1a8(%rsp),%rax
mulsd  (%rax,%r12,8),%xmm0
mov    0x1b0(%rsp),%rax
mulsd  (%rax,%r12,8),%xmm0
mov    0x1b8(%rsp),%rax
mulsd  (%rax,%r12,8),%xmm0
mov    0x180(%rsp),%rax
mov    0x90(%rsp),%rbp
addsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,(%rax,%rbp,8)
add    $0x1,%rbx
add    $0x1,%r13
mov    %r15,%rax
add    %r13,%rax
je     1580 <cu_compute_vector_stamps+0x470>
lea    (%r14,%r13,1),%eax
cmpb   $0x0,0x8(%rsp)
mov    0x44(%rsp),%edi
jne    1690 <cu_compute_vector_stamps+0x580>
cmpb   $0x0,0xc(%rsp)
je     1770 <cu_compute_vector_stamps+0x660>
mov    0x168(%rsp),%rcx
mov    0xf0(%rsp),%rdx
cmpl   $0x0,(%rcx,%rdx,4)
je     1786 <cu_compute_vector_stamps+0x676>
mov    0xe8(%rsp),%rcx
jmp    178e <cu_compute_vector_stamps+0x67e>
cs nopw 0x0(%rax,%rax,1)
mov    0x20(%rsp),%rcx
lea    (%rcx,%r13,1),%r12
movsd  0x887(%rip),%xmm0
jmp    16a4 <cu_compute_vector_stamps+0x594>
mov    0xe0(%rsp),%rcx
movsd  (%rcx,%r13,8),%xmm0
mov    0x18(%rsp),%rcx
lea    (%rbx,%rcx,1),%r12
mov    0x10(%rsp),%rcx
subsd  (%rcx,%r13,8),%xmm0
jmp    16a4 <cu_compute_vector_stamps+0x594>
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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