#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(int dp, int ds, int db, int nx, int ny, 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;
  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 (j = kernelRadius; j < ny - kernelRadius; j++) {
      y = (j - 0.5 * (ny - 1)) / (ny - 1);
      py = pow(y, m1);
      for (i = kernelRadius; i < nx - kernelRadius; 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>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x68(%rbp),%rax
mov    0x60(%rbp),%rax
mov    0x58(%rbp),%rax
mov    0x50(%rbp),%rax
mov    0x48(%rbp),%rax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9,-0x20(%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,-0x28(%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,-0x2c(%rbp)
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    0x40(%rbp),%eax
jge    1557 <cu_compute_vector+0x447>
mov    -0x84(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    11c4 <cu_compute_vector+0xb4>
xor    %eax,%eax
mov    %eax,-0x88(%rbp)
jmp    11da <cu_compute_vector+0xca>
mov    -0x84(%rbp),%eax
sub    -0x28(%rbp),%eax
cltd
idivl  -0x2c(%rbp)
add    $0x1,%eax
mov    %eax,-0x88(%rbp)
mov    -0x88(%rbp),%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x30(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1219 <cu_compute_vector+0x109>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x34(%rbp)
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x30(%rbp)
jne    1237 <cu_compute_vector+0x127>
mov    -0x4(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x84(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    1289 <cu_compute_vector+0x179>
mov    -0x30(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1266 <cu_compute_vector+0x156>
mov    -0x8(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x84(%rbp),%eax
sub    -0x28(%rbp),%eax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x24(%rbp)
jmp    1284 <cu_compute_vector+0x174>
mov    -0xc(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x84(%rbp),%eax
sub    -0x28(%rbp),%eax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x2c(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x24(%rbp)
jmp    1289 <cu_compute_vector+0x179>
movl   $0x0,-0x40(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     12fd <cu_compute_vector+0x1ed>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x3c(%rbp),%ecx
sub    -0x48(%rbp),%ecx
cmp    %ecx,%eax
jg     12ea <cu_compute_vector+0x1da>
mov    -0x40(%rbp),%eax
cmp    -0x24(%rbp),%eax
jne    12d3 <cu_compute_vector+0x1c3>
mov    -0x48(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12aa <cu_compute_vector+0x19a>
jmp    12ef <cu_compute_vector+0x1df>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1297 <cu_compute_vector+0x187>
movsd  0xcfb(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    0x30(%rbp),%rax
movslq -0x84(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x28(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    0x28(%rbp),%ecx
cmp    %ecx,%eax
jge    153e <cu_compute_vector+0x42e>
cvtsi2sdl -0x44(%rbp),%xmm1
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xcbd(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
mov    -0x54(%rbp),%edi
call   1030 <pow@plt>
movsd  %xmm0,-0x60(%rbp)
mov    0x28(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    0x28(%rbp),%ecx
cmp    %ecx,%eax
jge    152b <cu_compute_vector+0x41b>
cvtsi2sdl -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xc62(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
cmpl   $0x0,-0x30(%rbp)
jne    13ed <cu_compute_vector+0x2dd>
mov    0x48(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    14a1 <cu_compute_vector+0x391>
mov    -0x30(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    148f <cu_compute_vector+0x37f>
mov    0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     144d <cu_compute_vector+0x33d>
mov    0x50(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x34(%rbp),%ecx
mov    -0x10(%rbp),%edx
mov    -0x44(%rbp),%esi
add    -0x38(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    148a <cu_compute_vector+0x37a>
mov    0x48(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x34(%rbp),%ecx
mov    -0x10(%rbp),%edx
mov    -0x44(%rbp),%esi
add    -0x38(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    149c <cu_compute_vector+0x38c>
movsd  0xb69(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    14a1 <cu_compute_vector+0x391>
movsd  -0x68(%rbp),%xmm0
mov    -0x50(%rbp),%edi
call   1030 <pow@plt>
mulsd  -0x60(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
mov    0x58(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x60(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x68(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   -0x44(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x30(%rbp),%rax
movslq -0x84(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    137e <cu_compute_vector+0x26e>
jmp    1530 <cu_compute_vector+0x420>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1323 <cu_compute_vector+0x213>
jmp    1543 <cu_compute_vector+0x433>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    1199 <cu_compute_vector+0x89>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cu_compute_vector>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %r9,0x70(%rsp)
mov    %edx,0x60(%rsp)
mov    %rsi,0x68(%rsp)
mov    %rdi,0x40(%rsp)
mov    0x160(%rsp),%eax
test   %eax,%eax
jle    172d <cu_compute_vector+0x61d>
mov    0x168(%rsp),%rdi
mov    0x148(%rsp),%ebp
mov    0x40(%rsp),%rsi
lea    0x1(%rsi),%edx
add    $0x2,%esi
imul   %edx,%esi
mov    %esi,%ebx
shr    $0x1f,%ebx
add    %esi,%ebx
sar    %ebx
mov    0x68(%rsp),%rsi
lea    0x1(%rsi),%edx
add    $0x2,%esi
imul   %edx,%esi
mov    %esi,%edx
shr    $0x1f,%edx
add    %esi,%edx
sar    %edx
mov    %edx,0x28(%rsp)
mov    %r8d,%edx
sub    %ebp,%edx
mov    %edx,0x4c(%rsp)
add    $0xffffffff,%r8d
cvtsi2sd %r8d,%xmm0
mov    %ecx,%edx
sub    %ebp,%edx
mov    %edx,0x64(%rsp)
lea    -0x1(%rcx),%edx
cvtsi2sd %edx,%xmm1
movslq %ebp,%rdx
movslq %ecx,%rcx
mov    %eax,%eax
mov    %rax,0xb8(%rsp)
mov    %ecx,%esi
sub    %edx,%esi
sub    %edx,%esi
mov    %rsi,0xe8(%rsp)
lea    0x1(%rcx),%rax
imul   %rdx,%rax
movsd  0xe2d(%rip),%xmm2
movsd  %xmm0,0xd0(%rsp)
mulsd  %xmm0,%xmm2
movsd  %xmm2,0xc8(%rsp)
movsd  0xe0f(%rip),%xmm0
movsd  %xmm1,0xf0(%rsp)
mulsd  %xmm1,%xmm0
movsd  %xmm0,0xe0(%rsp)
mov    %ebx,0x2c(%rsp)
movslq %ebx,%rsi
mov    %rsi,0xb0(%rsp)
mov    %rax,0xa8(%rsp)
lea    (%rdi,%rax,8),%rax
mov    %rax,0xa0(%rsp)
mov    %rcx,0x50(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,0xc0(%rsp)
mov    %rdx,0x78(%rsp)
mov    %edx,%eax
mov    %rax,0xd8(%rsp)
xor    %ebp,%ebp
jmp    1272 <cu_compute_vector+0x162>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rbp
cmp    0xb8(%rsp),%rbp
je     172d <cu_compute_vector+0x61d>
cmp    0xb0(%rsp),%rbp
mov    %rbp,0x98(%rsp)
jge    12c0 <cu_compute_vector+0x1b0>
cmpl   $0x0,0x140(%rsp)
jle    1370 <cu_compute_vector+0x260>
movl   $0x0,0x20(%rsp)
mov    0x70(%rsp),%rax
mov    (%rax),%eax
mov    %eax,0x18(%rsp)
mov    0x130(%rsp),%rax
mov    (%rax),%eax
mov    %eax,0x10(%rsp)
mov    $0x1,%al
mov    %eax,0xc(%rsp)
jmp    1390 <cu_compute_vector+0x280>
xchg   %ax,%ax
mov    %ebp,%eax
sub    0x2c(%rsp),%eax
cltd
idivl  0x28(%rsp)
lea    0x1(%rax),%esi
cmp    0x140(%rsp),%esi
jge    1330 <cu_compute_vector+0x220>
movslq %esi,%rcx
mov    0x70(%rsp),%rdx
mov    (%rdx,%rcx,4),%edi
mov    0x130(%rsp),%rdx
mov    (%rdx,%rcx,4),%edx
mov    $0x1,%cl
mov    %ecx,0xc(%rsp)
test   %esi,%esi
mov    %edi,0x18(%rsp)
mov    %edx,0x10(%rsp)
je     1388 <cu_compute_vector+0x278>
mov    %esi,0x20(%rsp)
imul   0x28(%rsp),%eax
add    0x2c(%rsp),%eax
mov    %ebp,%r13d
sub    %eax,%r13d
movl   $0x0,0x8(%rsp)
mov    0x68(%rsp),%rax
mov    %eax,%r8d
jmp    13c9 <cu_compute_vector+0x2b9>
nopl   0x0(%rax,%rax,1)
test   %esi,%esi
je     1398 <cu_compute_vector+0x288>
mov    %esi,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    0x60(%rsp),%eax
mov    %eax,%r8d
jmp    13c9 <cu_compute_vector+0x2b9>
movl   $0x0,0x20(%rsp)
mov    $0x1,%al
mov    %eax,0x8(%rsp)
movl   $0x0,0xc(%rsp)
jmp    13ae <cu_compute_vector+0x29e>
movl   $0x0,0x20(%rsp)
mov    $0x1,%al
mov    %eax,0x8(%rsp)
jmp    13be <cu_compute_vector+0x2ae>
mov    $0x1,%al
mov    %eax,0x8(%rsp)
movl   $0x0,0xc(%rsp)
movl   $0x0,0x20(%rsp)
movl   $0x0,0x18(%rsp)
movl   $0x0,0x10(%rsp)
mov    0x40(%rsp),%rax
mov    %eax,%r8d
mov    %ebp,%r13d
test   %r8d,%r8d
js     14d0 <cu_compute_vector+0x3c0>
lea    0x1(%r8),%r15d
mov    %r13d,%r9d
neg    %r9d
xor    %ecx,%ecx
xor    %ebp,%ebp
mov    %r8,0x58(%rsp)
mov    %r9,0x30(%rsp)
jmp    1402 <cu_compute_vector+0x2f2>
nopl   0x0(%rax)
lea    0x1(%rcx),%eax
add    $0xffffffff,%r15d
cmp    %ecx,%r8d
mov    %eax,%ecx
je     14d0 <cu_compute_vector+0x3c0>
mov    %r8d,%edx
sub    %ecx,%edx
jl     13f0 <cu_compute_vector+0x2e0>
xor    %eax,%eax
cmp    $0x3,%edx
jb     148e <cu_compute_vector+0x37e>
mov    %r15d,%r12d
and    $0xfffffffc,%r12d
lea    (%r9,%rbp,1),%r11d
xor    %eax,%eax
mov    %rbp,0x38(%rsp)
mov    %ebp,%edi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
lea    0x1(%rdi),%r8d
lea    0x1(%rax),%r9d
lea    0x2(%rdi),%r10d
mov    %r11d,%esi
add    %eax,%esi
cmove  %eax,%r14d
lea    0x2(%rax),%esi
cmove  %ecx,%ebx
cmp    %r13d,%r8d
cmove  %r9d,%r14d
cmove  %ecx,%ebx
lea    0x3(%rdi),%edx
cmp    %r13d,%r10d
mov    %r14d,%ebp
cmove  %esi,%ebp
cmove  %ecx,%ebx
lea    0x3(%rax),%r14d
cmp    %r13d,%edx
cmove  %ecx,%ebx
cmovne %ebp,%r14d
add    $0x4,%edi
add    $0x4,%eax
cmp    %r12d,%eax
jne    1430 <cu_compute_vector+0x320>
mov    0x38(%rsp),%rbp
add    %eax,%ebp
mov    0x58(%rsp),%r8
mov    0x30(%rsp),%r9
test   $0x3,%r15b
je     13f0 <cu_compute_vector+0x2e0>
mov    %r15d,%edx
and    $0x3,%edx
neg    %edx
lea    (%r9,%rbp,1),%edi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
cmp    %esi,%edi
cmove  %ecx,%ebx
cmove  %eax,%r14d
add    $0x1,%eax
add    $0xffffffff,%esi
cmp    %esi,%edx
jne    14b0 <cu_compute_vector+0x3a0>
sub    %esi,%ebp
jmp    13f0 <cu_compute_vector+0x2e0>
nopw   0x0(%rax,%rax,1)
mov    0x150(%rsp),%rax
mov    0x98(%rsp),%rbp
movq   $0x0,(%rax,%rbp,8)
mov    0x4c(%rsp),%eax
cmp    0x148(%rsp),%eax
jle    1260 <cu_compute_vector+0x150>
movslq 0x20(%rsp),%rax
mov    %rax,0x20(%rsp)
movslq 0x18(%rsp),%rax
movslq 0x10(%rsp),%rcx
mov    0x78(%rsp),%rdx
add    %rdx,%rcx
imul   0x50(%rsp),%rcx
add    %rdx,%rcx
add    %rax,%rcx
mov    0x170(%rsp),%rax
lea    (%rax,%rcx,8),%rax
mov    %rax,0x90(%rsp)
mov    0x168(%rsp),%rax
lea    (%rax,%rcx,8),%rax
mov    %rax,0x88(%rsp)
mov    0xa0(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    0xa8(%rsp),%rax
mov    %rax,0x10(%rsp)
jmp    15b2 <cu_compute_vector+0x4a2>
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x80(%rsp),%rdx
add    $0x1,%rdx
mov    0x10(%rsp),%rax
add    0x50(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    0xc0(%rsp),%rax
add    %rax,0x30(%rsp)
add    %rax,0x90(%rsp)
add    %rax,0x88(%rsp)
cmp    %edx,0x4c(%rsp)
je     1260 <cu_compute_vector+0x150>
mov    %rdx,0x80(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
addsd  0xc8(%rsp),%xmm0
divsd  0xd0(%rsp),%xmm0
mov    %r14d,%edi
call   1030 <pow@plt>
movsd  %xmm0,0x58(%rsp)
mov    0x64(%rsp),%eax
cmp    0x148(%rsp),%eax
jle    1570 <cu_compute_vector+0x460>
mov    0x80(%rsp),%rax
imul   0x50(%rsp),%rax
mov    %rax,0x18(%rsp)
xor    %r15d,%r15d
mov    0x78(%rsp),%r12
jmp    16b0 <cu_compute_vector+0x5a0>
xchg   %ax,%ax
mov    0x18(%rsp),%rcx
lea    (%r12,%rcx,1),%r13
mov    0x30(%rsp),%rcx
movsd  (%rcx,%r15,8),%xmm0
movsd  %xmm0,0x38(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
addsd  0xe0(%rsp),%xmm0
divsd  0xf0(%rsp),%xmm0
mov    %ebx,%edi
call   1030 <pow@plt>
mulsd  0x58(%rsp),%xmm0
mulsd  0x38(%rsp),%xmm0
mov    0x178(%rsp),%rax
mulsd  (%rax,%r13,8),%xmm0
mov    0x180(%rsp),%rax
mulsd  (%rax,%r13,8),%xmm0
mov    0x188(%rsp),%rax
mulsd  (%rax,%r13,8),%xmm0
mov    0x150(%rsp),%rax
mov    0x98(%rsp),%rbp
addsd  (%rax,%rbp,8),%xmm0
movsd  %xmm0,(%rax,%rbp,8)
add    $0x1,%r12
add    $0x1,%r15
cmp    %r15d,0xe8(%rsp)
je     1570 <cu_compute_vector+0x460>
mov    0xd8(%rsp),%rax
add    %r15d,%eax
cmpb   $0x0,0x8(%rsp)
jne    1610 <cu_compute_vector+0x500>
cmpb   $0x0,0xc(%rsp)
je     16f0 <cu_compute_vector+0x5e0>
mov    0x138(%rsp),%rcx
mov    0x20(%rsp),%rdx
cmpl   $0x0,(%rcx,%rdx,4)
je     1706 <cu_compute_vector+0x5f6>
mov    0x90(%rsp),%rcx
jmp    170e <cu_compute_vector+0x5fe>
nopw   0x0(%rax,%rax,1)
mov    0x10(%rsp),%rcx
lea    (%rcx,%r15,1),%r13
movsd  0x907(%rip),%xmm0
jmp    1624 <cu_compute_vector+0x514>
mov    0x88(%rsp),%rcx
movsd  (%rcx,%r15,8),%xmm0
mov    0x18(%rsp),%rcx
lea    (%r12,%rcx,1),%r13
mov    0x30(%rsp),%rcx
subsd  (%rcx,%r15,8),%xmm0
jmp    1624 <cu_compute_vector+0x514>
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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