#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 */
int M_PI;
double acos(double);
int /*<<< orphan*/ approx(double *, double *, int, double *, double *, int, int,
                          int /*<<< orphan*/, int, int /*<<< orphan*/);
int /*<<< orphan*/ col_gradient(int, int, double *, double, double *);
int cos(double);
int /*<<< orphan*/ cumtrapz(int *, double *, double *, double *);
int exp(int);
int /*<<< orphan*/ group_action_by_gamma(int *, int *, double *, double *,
                                         double *);
double innerprod_q2(int *, double *, double *);
int /*<<< orphan*/ linspace(int /*<<< orphan*/, int, int, double *);
double log(int);
int /*<<< orphan*/ printf(char *, int, double);
int /*<<< orphan*/ product(int, int, int, double *, double *, double *);
int /*<<< orphan*/ pvecnorm(int *, double *, double *, double *);
int sin(double);
int sqrt(int);
void oclogit_warp_grad(int *n1, int *T1, double *alpha, double *nu, double *q,
                       int *y, int *max_itri, double *toli, double *deltaOi,
                       double *deltagi, int *displayi, double *gamout,
                       double *Oout) {
  int TT = *T1;
  int n = *n1;
  int max_itr = *max_itri;
  double tol = *toli, deltaO = *deltaOi, deltag = *deltagi;
  double display = *displayi;
  int j, k, l, jj;
  int p = 20;
  int p1 = 10;
  int itr = 0;
  double binsize1 = 1.0;
  double t[TT], gam1[TT], ones[TT], O1[4], f_basis[TT * p], q_tilde[TT * n];
  double O_tmp[4], q_tmp[TT * n], O2[4], q_tilde_diff[TT * n], cbar[TT];
  double ftmp[TT], c[TT], tmp5[TT * p], hpsi[TT], psi[TT], gam2[TT];
  double gam_tmp[TT], max_val[max_itr + 1], tmp7[TT];
  double binsize, A, theta, B, tmp1, tmp2, thetanew, tmpi;
  double max_val_change, res_cos, res_sin, hO;
  double *t_ptr, *gam1_ptr, *f_basis_ptr, *q_tilde_ptr, *A_ptr, *nu_ptr;
  double *O_tmp_ptr, *q_tmp_ptr, *alpha_ptr, *q_tilde_diff_ptr;
  double *ftmp_ptr, *c_ptr, *cbar_ptr, *tmp5_ptr, *hpsi_ptr, *psi_ptr;
  double *ones_ptr, *gam2_ptr, *gam_tmp_ptr, *tmp7_ptr, *O1_ptr, *q_ptr;
  int *y_ptr;
  t_ptr = t;
  linspace(0, 1, TT, t_ptr);
  binsize = 1.0 / (TT - 1);
  gam1_ptr = gam1;
  for (k = 0; k < TT; k++) {
    gam1[k] = t[k];
    ones[k] = 1.0;
  }
  O1[0] = 1;
  O1[1] = 0;
  O1[2] = 0;
  O1[3] = 1;
  f_basis_ptr = f_basis;
  for (k = 0; k < p1; k++) {
    for (j = 0; j < TT; j++) {
      f_basis[j + 2 * k * TT] = 1 / sqrt(M_PI) * sin(2 * M_PI * (k + 1) * t[j]);
      f_basis[j + (2 * k + 1) * TT] =
          1 / sqrt(M_PI) * cos(2 * M_PI * (k + 1) * t[j]);
    }
  }
  q_tilde_ptr = q_tilde;
  for (k = 0; k < TT * n; k++)
    q_tilde[k] = q[k];
  do {
    nu_ptr = nu;
    q_tilde_ptr = q_tilde;
    A = innerprod_q2(T1, q_tilde_ptr, nu_ptr);
    theta = acos(O1[0]);
    O_tmp[0] = -1 * sin(theta);
    O_tmp[1] = cos(theta);
    O_tmp[2] = -1 * cos(theta);
    O_tmp[3] = -1 * sin(theta);
    O_tmp_ptr = O_tmp;
    q_tilde_ptr = q_tilde;
    q_tmp_ptr = q_tmp;
    product(n, n, TT, O_tmp_ptr, q_tilde_ptr, q_tmp_ptr);
    nu_ptr = nu;
    B = innerprod_q2(T1, q_tmp_ptr, nu_ptr);
    y_ptr = y;
    alpha_ptr = alpha;
    tmp1 = exp((-1 * y_ptr[0]) * (alpha_ptr[0] + A));
    tmp2 = (y_ptr[0] * tmp1) / (1 + tmp1);
    hO = tmp2 * B;
    thetanew = theta + deltaO * hO;
    O2[0] = cos(thetanew);
    O2[1] = sin(thetanew);
    O2[2] = -1 * sin(thetanew);
    O2[3] = cos(thetanew);
    q_tilde_diff_ptr = q_tilde_diff;
    q_tilde_ptr = q_tilde;
    col_gradient(n, TT, q_tilde_ptr, binsize, q_tilde_diff_ptr);
    cbar_ptr = cbar;
    ftmp_ptr = ftmp;
    q_tmp_ptr = q_tmp;
    c_ptr = c;
    nu_ptr = nu;
    tmp5_ptr = tmp5;
    for (k = 0; k < p; k++) {
      for (l = 0; l < TT; l++)
        ftmp[l] = f_basis[l + k * TT];
      cumtrapz(T1, t_ptr, ftmp_ptr, cbar_ptr);
      for (jj = 0; jj < n; jj++) {
        for (l = 0; l < TT; l++)
          q_tmp[n * l + jj] = 2 * q_tilde_diff[n * l + jj] * cbar[l] +
                              q_tilde[n * l + jj] * f_basis[l + k * TT];
      }
      tmpi = innerprod_q2(T1, q_tmp_ptr, nu_ptr);
      for (l = 0; l < TT; l++)
        tmp5[l + k * TT] = tmpi * f_basis[l + k * TT];
    }
    for (jj = 0; jj < TT; jj++) {
      c[jj] = tmp5[jj];
      for (l = 1; l < p; l++)
        c[jj] = c[jj] + tmp5[jj + l * TT];
    }
    hpsi_ptr = hpsi;
    c_ptr = c;
    for (k = 0; k < TT; k++)
      hpsi_ptr[k] = tmp2 * c_ptr[k];
    psi_ptr = psi;
    ones_ptr = ones;
    gam2_ptr = gam2;
    gam_tmp_ptr = gam_tmp;
    hpsi_ptr = hpsi;
    pvecnorm(T1, hpsi_ptr, &binsize1, &tmpi);
    res_cos = cos(deltag * tmpi);
    res_sin = sin(deltag * tmpi);
    if (tmpi == 0) {
      for (j = 0; j < TT; j++)
        psi_ptr[j] = res_cos * ones_ptr[j];
    } else {
      for (j = 0; j < TT; j++)
        psi_ptr[j] = res_cos * ones_ptr[j] + res_sin * (hpsi_ptr[j] / tmpi);
    }
    for (j = 0; j < TT; j++)
      tmp7[j] = psi_ptr[j] * psi_ptr[j];
    tmp7_ptr = tmp7;
    cumtrapz(T1, t_ptr, tmp7_ptr, gam_tmp_ptr);
    for (j = 0; j < TT; j++)
      gam_tmp_ptr[j] = (gam_tmp_ptr[j] - gam_tmp_ptr[0]) /
                       (gam_tmp_ptr[TT - 1] - gam_tmp_ptr[0]);
    approx(t, gam1_ptr, TT, gam_tmp_ptr, gam2_ptr, TT, 1, 0, 1, 0);
    y_ptr = y;
    alpha_ptr = alpha;
    max_val[itr] = log(1 / (1 + exp((-1 * y_ptr[0]) * (alpha_ptr[0] + A))));
    if (display == 1)
      printf("Iteration %d : Cost %f\n", (itr + 1), max_val[itr]);
    gam1_ptr = gam2;
    O1_ptr = O2;
    O1[0] = O2[0];
    O1[1] = O2[1];
    O1[2] = O2[2];
    O1[3] = O2[3];
    q_tilde_ptr = q_tilde;
    q_tmp_ptr = q_tmp;
    q_ptr = q;
    product(n, n, TT, O1_ptr, q_ptr, q_tmp_ptr);
    group_action_by_gamma(n1, T1, q_tmp_ptr, gam1_ptr, q_tilde_ptr);
    if (itr >= 2) {
      max_val_change = max_val[itr] - max_val[itr - 1];
      if (hO < tol && tmpi < tol)
        break;
    }
    itr++;
  } while (max_itr >= itr);
  for (k = 0; k < TT; k++) {
    gamout[k] = gam1_ptr[k];
  }
  for (k = 0; k < n * n; k++) {
    Oout[k] = O1_ptr[k];
  }
}


x86 assembly (clang, -O0):
<oclogit_warp_grad>:
push   %rbp
mov    %rsp,%rbp
sub    $0x330,%rsp
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x38(%rbp)
mov    0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x3c(%rbp)
mov    0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    0x30(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x14,-0x74(%rbp)
movl   $0xa,-0x78(%rbp)
movl   $0x0,-0x7c(%rbp)
movabs $0x3ff0000000000000,%rax
mov    %rax,-0x88(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,%ecx
mov    %rsp,%rax
mov    %rax,-0x90(%rbp)
lea    0xf(,%rcx,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rax
sub    %rdx,%rax
mov    %rax,-0x328(%rbp)
mov    %rax,%rsp
mov    %rcx,-0x98(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2a8(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0xa0(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x320(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0xa8(%rbp)
mov    -0x34(%rbp),%ecx
mov    -0x74(%rbp),%edx
imul   %edx,%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x318(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0xd8(%rbp)
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   %edx,%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x310(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0xe0(%rbp)
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   %edx,%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x308(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x108(%rbp)
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   %edx,%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x300(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x138(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2f8(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x140(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2f0(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x148(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2e8(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x150(%rbp)
mov    -0x34(%rbp),%ecx
mov    -0x74(%rbp),%edx
imul   %edx,%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2e0(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x158(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2d8(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x160(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2d0(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x168(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2c8(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x170(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2c0(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x178(%rbp)
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2b8(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x180(%rbp)
mov    -0x34(%rbp),%ecx
lea    0xf(,%rcx,8),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x2b0(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x188(%rbp)
mov    %rax,-0x1f0(%rbp)
mov    -0x34(%rbp),%edx
mov    -0x1f0(%rbp),%rcx
xor    %edi,%edi
mov    $0x1,%esi
call   10e0 <linspace@plt>
mov    -0x2a8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cvtsi2sd %ecx,%xmm1
movsd  0x1a73(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
mov    %rax,-0x1f8(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    15f9 <oclogit_warp_grad+0x409>
mov    -0x320(%rbp),%rax
mov    -0x2a8(%rbp),%rcx
mov    -0x328(%rbp),%rdx
movslq -0x68(%rbp),%rsi
movsd  (%rdx,%rsi,8),%xmm0
movslq -0x68(%rbp),%rdx
movsd  %xmm0,(%rcx,%rdx,8)
movslq -0x68(%rbp),%rcx
movsd  0x1a1a(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    15a7 <oclogit_warp_grad+0x3b7>
mov    -0x318(%rbp),%rax
movsd  0x19f8(%rip),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  0x19d2(%rip),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    %rax,-0x200(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x78(%rbp),%eax
jge    177d <oclogit_warp_grad+0x58d>
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    176a <oclogit_warp_grad+0x57a>
mov    0x397e(%rip),%rax
mov    (%rax),%edi
call   1110 <sqrt@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x328(%rbp),%rax
mov    %ecx,-0x330(%rbp)
mov    0x3957(%rip),%rcx
mov    (%rcx),%ecx
shl    $0x1,%ecx
mov    -0x68(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
cvtsi2sd %ecx,%xmm0
movslq -0x64(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
call   1090 <sin@plt>
mov    -0x330(%rbp),%ecx
mov    %eax,%edx
mov    -0x318(%rbp),%rax
imul   %edx,%ecx
cvtsi2sd %ecx,%xmm0
mov    -0x64(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x3903(%rip),%rax
mov    (%rax),%edi
call   1110 <sqrt@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x328(%rbp),%rax
mov    %ecx,-0x32c(%rbp)
mov    0x38dc(%rip),%rcx
mov    (%rcx),%ecx
shl    $0x1,%ecx
mov    -0x68(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
cvtsi2sd %ecx,%xmm0
movslq -0x64(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
call   10f0 <cos@plt>
mov    -0x32c(%rbp),%ecx
mov    %eax,%edx
mov    -0x318(%rbp),%rax
imul   %edx,%ecx
cvtsi2sd %ecx,%xmm0
mov    -0x64(%rbp),%ecx
mov    -0x68(%rbp),%edx
shl    $0x1,%edx
add    $0x1,%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1657 <oclogit_warp_grad+0x467>
jmp    176f <oclogit_warp_grad+0x57f>
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1644 <oclogit_warp_grad+0x454>
mov    -0x310(%rbp),%rax
mov    %rax,-0x208(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    -0x34(%rbp),%ecx
imul   -0x38(%rbp),%ecx
cmp    %ecx,%eax
jge    17cf <oclogit_warp_grad+0x5df>
mov    -0x310(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x68(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1792 <oclogit_warp_grad+0x5a2>
jmp    17d4 <oclogit_warp_grad+0x5e4>
mov    -0x310(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    %rcx,-0x218(%rbp)
mov    %rax,-0x208(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x208(%rbp),%rsi
mov    -0x218(%rbp),%rdx
call   1080 <innerprod_q2@plt>
movsd  %xmm0,-0x198(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
call   1060 <acos@plt>
movsd  %xmm0,-0x1a0(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
call   1090 <sin@plt>
imul   $0xffffffff,%eax,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
call   10f0 <cos@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
call   10f0 <cos@plt>
imul   $0xffffffff,%eax,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
call   1090 <sin@plt>
mov    -0x310(%rbp),%rcx
mov    %eax,%edx
mov    -0x308(%rbp),%rax
imul   $0xffffffff,%edx,%edx
cvtsi2sd %edx,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
lea    -0x100(%rbp),%rdx
mov    %rdx,-0x220(%rbp)
mov    %rcx,-0x208(%rbp)
mov    %rax,-0x228(%rbp)
mov    -0x38(%rbp),%edi
mov    -0x38(%rbp),%esi
mov    -0x34(%rbp),%edx
mov    -0x220(%rbp),%rcx
mov    -0x208(%rbp),%r8
mov    -0x228(%rbp),%r9
call   10c0 <product@plt>
mov    -0x20(%rbp),%rax
mov    %rax,-0x218(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x228(%rbp),%rsi
mov    -0x218(%rbp),%rdx
call   1080 <innerprod_q2@plt>
movsd  %xmm0,-0x1a8(%rbp)
ff
mov    -0x30(%rbp),%rax
mov    %rax,-0x2a0(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x230(%rbp)
mov    -0x2a0(%rbp),%rax
imul   $0xffffffff,(%rax),%eax
cvtsi2sd %eax,%xmm0
mov    -0x230(%rbp),%rax
movsd  (%rax),%xmm1
addsd  -0x198(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   10d0 <exp@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1b0(%rbp)
ff
mov    -0x2a0(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
mulsd  -0x1b0(%rbp),%xmm0
ff
movsd  0x1682(%rip),%xmm1
addsd  -0x1b0(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0x1b8(%rbp),%xmm0
ff
mulsd  -0x1a8(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e8(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm1
ff
movsd  -0x50(%rbp),%xmm0
movsd  -0x1e8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x1c0(%rbp),%xmm0
ff
call   10f0 <cos@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x1c0(%rbp),%xmm0
ff
call   1090 <sin@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x128(%rbp)
movsd  -0x1c0(%rbp),%xmm0
ff
call   1090 <sin@plt>
imul   $0xffffffff,%eax,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0x1c0(%rbp),%xmm0
ff
call   10f0 <cos@plt>
mov    -0x300(%rbp),%rcx
mov    %eax,%edx
mov    -0x310(%rbp),%rax
cvtsi2sd %edx,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
mov    %rcx,-0x238(%rbp)
mov    %rax,-0x208(%rbp)
mov    -0x38(%rbp),%edi
mov    -0x34(%rbp),%esi
mov    -0x208(%rbp),%rdx
movsd  -0x190(%rbp),%xmm0
ff
mov    -0x238(%rbp),%rcx
call   10a0 <col_gradient@plt>
mov    -0x2f8(%rbp),%rdi
mov    -0x2f0(%rbp),%rsi
mov    -0x308(%rbp),%rdx
mov    -0x2e8(%rbp),%rcx
mov    -0x2e0(%rbp),%rax
mov    %rdi,-0x250(%rbp)
mov    %rsi,-0x240(%rbp)
mov    %rdx,-0x228(%rbp)
mov    %rcx,-0x248(%rbp)
mov    -0x20(%rbp),%rcx
mov    %rcx,-0x218(%rbp)
mov    %rax,-0x258(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x74(%rbp),%eax
jge    1ca8 <oclogit_warp_grad+0xab8>
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1b28 <oclogit_warp_grad+0x938>
mov    -0x2f0(%rbp),%rax
mov    -0x318(%rbp),%rcx
mov    -0x6c(%rbp),%edx
mov    -0x68(%rbp),%esi
imul   -0x34(%rbp),%esi
add    %esi,%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x6c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1ae3 <oclogit_warp_grad+0x8f3>
mov    -0x10(%rbp),%rdi
mov    -0x1f0(%rbp),%rsi
mov    -0x240(%rbp),%rdx
mov    -0x250(%rbp),%rcx
call   1070 <cumtrapz@plt>
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1c17 <oclogit_warp_grad+0xa27>
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1c04 <oclogit_warp_grad+0xa14>
mov    -0x308(%rbp),%rax
mov    -0x318(%rbp),%rcx
mov    -0x310(%rbp),%rdx
mov    -0x2f8(%rbp),%rsi
mov    -0x300(%rbp),%rdi
mov    -0x38(%rbp),%r8d
imul   -0x6c(%rbp),%r8d
add    -0x70(%rbp),%r8d
movslq %r8d,%r8
movsd  0x1461(%rip),%xmm0
mulsd  (%rdi,%r8,8),%xmm0
movslq -0x6c(%rbp),%rdi
movsd  (%rsi,%rdi,8),%xmm2
mov    -0x38(%rbp),%esi
imul   -0x6c(%rbp),%esi
add    -0x70(%rbp),%esi
movslq %esi,%rsi
movsd  (%rdx,%rsi,8),%xmm1
mov    -0x6c(%rbp),%edx
mov    -0x68(%rbp),%esi
imul   -0x34(%rbp),%esi
add    %esi,%edx
movslq %edx,%rdx
mulsd  (%rcx,%rdx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x38(%rbp),%ecx
imul   -0x6c(%rbp),%ecx
add    -0x70(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1b60 <oclogit_warp_grad+0x970>
jmp    1c09 <oclogit_warp_grad+0xa19>
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    1b4d <oclogit_warp_grad+0x95d>
mov    -0x10(%rbp),%rdi
mov    -0x228(%rbp),%rsi
mov    -0x218(%rbp),%rdx
call   1080 <innerprod_q2@plt>
movsd  %xmm0,-0x1c8(%rbp)
ff
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1c95 <oclogit_warp_grad+0xaa5>
mov    -0x2e0(%rbp),%rax
mov    -0x318(%rbp),%rcx
movsd  -0x1c8(%rbp),%xmm0
ff
mov    -0x6c(%rbp),%edx
mov    -0x68(%rbp),%esi
imul   -0x34(%rbp),%esi
add    %esi,%edx
movslq %edx,%rdx
mulsd  (%rcx,%rdx,8),%xmm0
mov    -0x6c(%rbp),%ecx
mov    -0x68(%rbp),%edx
imul   -0x34(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1c3d <oclogit_warp_grad+0xa4d>
jmp    1c9a <oclogit_warp_grad+0xaaa>
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1ad0 <oclogit_warp_grad+0x8e0>
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1d43 <oclogit_warp_grad+0xb53>
mov    -0x2e8(%rbp),%rax
mov    -0x2e0(%rbp),%rcx
movslq -0x70(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x70(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x74(%rbp),%eax
jge    1d30 <oclogit_warp_grad+0xb40>
mov    -0x2e8(%rbp),%rax
mov    -0x2e0(%rbp),%rcx
movslq -0x70(%rbp),%rdx
movsd  (%rax,%rdx,8),%xmm0
mov    -0x70(%rbp),%edx
mov    -0x6c(%rbp),%esi
imul   -0x34(%rbp),%esi
add    %esi,%edx
movslq %edx,%rdx
addsd  (%rcx,%rdx,8),%xmm0
movslq -0x70(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1ce2 <oclogit_warp_grad+0xaf2>
jmp    1d35 <oclogit_warp_grad+0xb45>
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    1caf <oclogit_warp_grad+0xabf>
mov    -0x2e8(%rbp),%rax
mov    -0x2d8(%rbp),%rcx
mov    %rcx,-0x260(%rbp)
mov    %rax,-0x248(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1da8 <oclogit_warp_grad+0xbb8>
movsd  -0x1b8(%rbp),%xmm0
ff
mov    -0x248(%rbp),%rax
movslq -0x68(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x260(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    1d66 <oclogit_warp_grad+0xb76>
mov    -0x2d8(%rbp),%rax
mov    -0x2c0(%rbp),%rcx
mov    -0x2c8(%rbp),%rdx
mov    -0x320(%rbp),%rsi
mov    -0x2d0(%rbp),%rdi
mov    %rdi,-0x268(%rbp)
mov    %rsi,-0x270(%rbp)
mov    %rdx,-0x278(%rbp)
mov    %rcx,-0x280(%rbp)
mov    %rax,-0x260(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x260(%rbp),%rsi
lea    -0x88(%rbp),%rdx
lea    -0x1c8(%rbp),%rcx
call   10b0 <pvecnorm@plt>
movsd  -0x58(%rbp),%xmm0
mulsd  -0x1c8(%rbp),%xmm0
ff
call   10f0 <cos@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
mulsd  -0x1c8(%rbp),%xmm0
ff
call   1090 <sin@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  -0x1c8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1eb1 <oclogit_warp_grad+0xcc1>
jp     1eb1 <oclogit_warp_grad+0xcc1>
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1eac <oclogit_warp_grad+0xcbc>
movsd  -0x1d8(%rbp),%xmm0
ff
mov    -0x270(%rbp),%rax
movslq -0x64(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x268(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1e6a <oclogit_warp_grad+0xc7a>
jmp    1f2b <oclogit_warp_grad+0xd3b>
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1f26 <oclogit_warp_grad+0xd36>
movsd  -0x1d8(%rbp),%xmm0
ff
mov    -0x270(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x1e0(%rbp),%xmm1
ff
mov    -0x260(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm3
divsd  -0x1c8(%rbp),%xmm3
ff
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x268(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1eb8 <oclogit_warp_grad+0xcc8>
jmp    1f2b <oclogit_warp_grad+0xd3b>
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1f7c <oclogit_warp_grad+0xd8c>
mov    -0x2b0(%rbp),%rax
mov    -0x268(%rbp),%rcx
movslq -0x64(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x268(%rbp),%rcx
movslq -0x64(%rbp),%rdx
mulsd  (%rcx,%rdx,8),%xmm0
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1f32 <oclogit_warp_grad+0xd42>
mov    -0x2b0(%rbp),%rax
mov    %rax,-0x288(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x1f0(%rbp),%rsi
mov    -0x288(%rbp),%rdx
mov    -0x280(%rbp),%rcx
call   1070 <cumtrapz@plt>
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    2018 <oclogit_warp_grad+0xe28>
mov    -0x280(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x280(%rbp),%rax
subsd  (%rax),%xmm0
mov    -0x280(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x280(%rbp),%rax
subsd  (%rax),%xmm1
divsd  %xmm1,%xmm0
mov    -0x280(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1faf <oclogit_warp_grad+0xdbf>
mov    -0x328(%rbp),%rdi
mov    -0x1f8(%rbp),%rsi
mov    -0x34(%rbp),%edx
mov    -0x280(%rbp),%rcx
mov    -0x278(%rbp),%r8
mov    -0x34(%rbp),%r9d
xor    %eax,%eax
sub    $0x20,%rsp
movl   $0x1,(%rsp)
movl   $0x0,0x8(%rsp)
movl   $0x1,0x10(%rsp)
movl   $0x0,0x18(%rsp)
call   1030 <approx@plt>
add    $0x20,%rsp
mov    -0x30(%rbp),%rax
mov    %rax,-0x2a0(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x230(%rbp)
mov    -0x2a0(%rbp),%rax
imul   $0xffffffff,(%rax),%eax
cvtsi2sd %eax,%xmm0
mov    -0x230(%rbp),%rax
movsd  (%rax),%xmm1
addsd  -0x198(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   10d0 <exp@plt>
mov    %eax,%ecx
add    $0x1,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
mov    %eax,%edi
call   1100 <log@plt>
mov    -0x2b8(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
movsd  0xf1f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    2113 <oclogit_warp_grad+0xf23>
jp     2113 <oclogit_warp_grad+0xf23>
mov    -0x2b8(%rbp),%rax
mov    -0x7c(%rbp),%esi
add    $0x1,%esi
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xf02(%rip),%rdi
call   1050 <printf@plt>
mov    -0x308(%rbp),%rax
mov    -0x310(%rbp),%rcx
mov    -0x2c8(%rbp),%rdx
mov    %rdx,-0x1f8(%rbp)
lea    -0x130(%rbp),%rdx
mov    %rdx,-0x290(%rbp)
movsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x120(%rbp),%xmm0
ff
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x118(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
mov    %rcx,-0x208(%rbp)
mov    %rax,-0x228(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x298(%rbp)
mov    -0x38(%rbp),%edi
mov    -0x38(%rbp),%esi
mov    -0x34(%rbp),%edx
mov    -0x290(%rbp),%rcx
mov    -0x298(%rbp),%r8
mov    -0x228(%rbp),%r9
call   10c0 <product@plt>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x228(%rbp),%rdx
mov    -0x1f8(%rbp),%rcx
mov    -0x208(%rbp),%r8
call   1040 <group_action_by_gamma@plt>
cmpl   $0x2,-0x7c(%rbp)
jl     2243 <oclogit_warp_grad+0x1053>
mov    -0x2b8(%rbp),%rax
movslq -0x7c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x7c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x1e8(%rbp),%xmm1
ff
movsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    223e <oclogit_warp_grad+0x104e>
movsd  -0x1c8(%rbp),%xmm1
ff
movsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    223e <oclogit_warp_grad+0x104e>
jmp    2258 <oclogit_warp_grad+0x1068>
jmp    2243 <oclogit_warp_grad+0x1053>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x7c(%rbp),%eax
jge    17d4 <oclogit_warp_grad+0x5e4>
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    2296 <oclogit_warp_grad+0x10a6>
mov    -0x1f8(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x38(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    225f <oclogit_warp_grad+0x106f>
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x38(%rbp),%ecx
cmp    %ecx,%eax
jge    22da <oclogit_warp_grad+0x10ea>
mov    -0x290(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x40(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    229d <oclogit_warp_grad+0x10ad>
mov    -0x90(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<oclogit_warp_grad>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1d8,%rsp
mov    %r9,-0x140(%rbp)
mov    %r8,-0x120(%rbp)
mov    %rcx,-0x100(%rbp)
mov    %rdx,-0x138(%rbp)
mov    %rsp,-0x160(%rbp)
mov    0x30(%rbp),%r8
mov    0x28(%rbp),%rcx
mov    0x20(%rbp),%rdx
mov    0x18(%rbp),%rbx
mov    0x10(%rbp),%rax
movsd  (%rbx),%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
mov    %rsi,-0x48(%rbp)
movslq (%rsi),%rbx
movsd  (%rdx),%xmm0
movsd  %xmm0,-0x188(%rbp)
ff
mov    %rdi,-0x198(%rbp)
movslq (%rdi),%r14
mov    (%rax),%eax
mov    %eax,-0x104(%rbp)
movsd  (%rcx),%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
mov    (%r8),%eax
mov    %eax,-0x108(%rbp)
movabs $0x3ff0000000000000,%r12
mov    %r12,-0x1b8(%rbp)
mov    %ebx,%r15d
mov    %rsp,%rcx
lea    0xf(,%r15,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rcx
mov    %rcx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0xf0(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0x98(%rbp)
mov    %rdx,%rsp
lea    0x0(,%rbx,4),%edx
lea    (%rdx,%rdx,4),%esi
shl    $0x3,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x70(%rbp)
mov    %rdx,%rsp
mov    %r14d,%edx
imul   %ebx,%edx
mov    %rdx,-0x60(%rbp)
lea    0xf(,%rdx,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rdi
sub    %rdx,%rdi
mov    %rdi,-0xa8(%rbp)
mov    %rdi,%rsp
mov    %rsp,%rdi
sub    %rdx,%rdi
mov    %rdi,-0xd8(%rbp)
mov    %rdi,%rsp
mov    %rsp,%rdi
sub    %rdx,%rdi
mov    %rdi,-0x158(%rbp)
mov    %rdi,%rsp
mov    %rsp,%r13
sub    %rax,%r13
mov    %r13,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0x1b0(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0xa0(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0xc8(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0xc0(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0x118(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0x128(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0xe8(%rbp)
mov    %rdx,%rsp
mov    %rsp,%rdx
sub    %rax,%rdx
mov    %rdx,-0xb8(%rbp)
mov    %rdx,%rsp
xor    %edi,%edi
mov    $0x1,%esi
mov    %ebx,%edx
mov    %rcx,-0x78(%rbp)
call   10f0 <linspace@plt>
lea    -0x1(%rbx),%eax
mov    %r14d,%ecx
mov    %rcx,-0xd0(%rbp)
cvtsi2sd %eax,%xmm1
test   %rbx,%rbx
mov    %rbx,-0x30(%rbp)
jle    154d <oclogit_warp_grad+0x34d>
movsd  %xmm1,-0x38(%rbp)
lea    0x0(,%r15,8),%rdx
mov    -0xf0(%rbp),%rdi
mov    -0x78(%rbp),%rsi
call   10a0 <memcpy@plt>
mov    -0x30(%rbp),%rbx
cmp    $0x4,%ebx
jae    142d <oclogit_warp_grad+0x22d>
xor    %eax,%eax
mov    -0x98(%rbp),%rdi
movsd  -0x38(%rbp),%xmm1
jmp    1540 <oclogit_warp_grad+0x340>
mov    %r15d,%eax
and    $0xfffffffc,%eax
lea    -0x4(%rax),%rsi
mov    %rsi,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %edx,%ecx
and    $0x7,%ecx
cmp    $0x1c,%rsi
jae    1460 <oclogit_warp_grad+0x260>
xor    %esi,%esi
mov    -0x98(%rbp),%rdi
movsd  -0x38(%rbp),%xmm1
jmp    1501 <oclogit_warp_grad+0x301>
movabs $0x7ffffffffffffff8,%rsi
ff ff 7f
and    %rsi,%rdx
xor    %esi,%esi
movaps 0x1b9a(%rip),%xmm0
mov    -0x98(%rbp),%rdi
movsd  -0x38(%rbp),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movaps %xmm0,(%rdi,%rsi,8)
movaps %xmm0,0x10(%rdi,%rsi,8)
movaps %xmm0,0x20(%rdi,%rsi,8)
movaps %xmm0,0x30(%rdi,%rsi,8)
movaps %xmm0,0x40(%rdi,%rsi,8)
movaps %xmm0,0x50(%rdi,%rsi,8)
movaps %xmm0,0x60(%rdi,%rsi,8)
movaps %xmm0,0x70(%rdi,%rsi,8)
movaps %xmm0,0x80(%rdi,%rsi,8)
movaps %xmm0,0x90(%rdi,%rsi,8)
movaps %xmm0,0xa0(%rdi,%rsi,8)
movaps %xmm0,0xb0(%rdi,%rsi,8)
movaps %xmm0,0xc0(%rdi,%rsi,8)
movaps %xmm0,0xd0(%rdi,%rsi,8)
movaps %xmm0,0xe0(%rdi,%rsi,8)
movaps %xmm0,0xf0(%rdi,%rsi,8)
add    $0x20,%rsi
add    $0xfffffffffffffff8,%rdx
jne    1490 <oclogit_warp_grad+0x290>
test   %rcx,%rcx
je     1532 <oclogit_warp_grad+0x332>
lea    (%rdi,%rsi,8),%rdx
add    $0x10,%rdx
shl    $0x5,%rcx
xor    %esi,%esi
movaps 0x1af5(%rip),%xmm0
nopl   0x0(%rax,%rax,1)
movaps %xmm0,-0x10(%rdx,%rsi,1)
movaps %xmm0,(%rdx,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%rcx
jne    1520 <oclogit_warp_grad+0x320>
cmp    %r15,%rax
je     154d <oclogit_warp_grad+0x34d>
nopw   0x0(%rax,%rax,1)
mov    %r12,(%rdi,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    1540 <oclogit_warp_grad+0x340>
movsd  0x1aab(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
mov    -0x70(%rbp),%rax
lea    (%rax,%rbx,8),%rcx
mov    %rbx,%rdx
shl    $0x4,%rdx
mov    %rdx,-0x90(%rbp)
mov    0x3a6a(%rip),%rdx
xor    %esi,%esi
jmp    15b6 <oclogit_warp_grad+0x3b6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x90(%rbp),%rax
add    %rax,%rcx
mov    -0x58(%rbp),%rdi
add    %rax,%rdi
mov    -0x50(%rbp),%rsi
cmp    $0xa,%rsi
mov    -0x30(%rbp),%rbx
mov    %rdi,%rax
je     169c <oclogit_warp_grad+0x49c>
mov    %rax,-0x58(%rbp)
add    $0x1,%rsi
mov    %rsi,-0x50(%rbp)
test   %ebx,%ebx
jle    1590 <oclogit_warp_grad+0x390>
mov    -0x50(%rbp),%rax
add    %eax,%eax
mov    %eax,-0x38(%rbp)
xor    %r12d,%r12d
mov    %rcx,-0x68(%rbp)
cs nopw 0x0(%rax,%rax,1)
mov    (%rdx),%edi
call   1120 <sqrt@plt>
mov    %eax,%ebx
lea    0x1(%rbx),%eax
cmp    $0x3,%eax
mov    $0x0,%eax
cmovae %eax,%ebx
mov    0x39ea(%rip),%rax
mov    (%rax),%eax
imul   -0x38(%rbp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    -0x78(%rbp),%rax
mulsd  (%rax,%r12,8),%xmm0
call   1090 <sin@plt>
imul   %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    -0x58(%rbp),%rax
movsd  %xmm0,(%rax,%r12,8)
mov    0x39b3(%rip),%rax
mov    (%rax),%edi
call   1120 <sqrt@plt>
mov    %eax,%ebx
lea    0x1(%rbx),%eax
cmp    $0x3,%eax
mov    $0x0,%eax
cmovae %eax,%ebx
mov    0x3995(%rip),%rax
mov    (%rax),%eax
imul   -0x38(%rbp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    -0x78(%rbp),%rax
mulsd  (%rax,%r12,8),%xmm0
call   1100 <cos@plt>
mov    0x3972(%rip),%rdx
mov    -0x68(%rbp),%rcx
imul   %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rcx,%r12,8)
add    $0x1,%r12
cmp    %r12,%r15
jne    15e0 <oclogit_warp_grad+0x3e0>
jmp    1590 <oclogit_warp_grad+0x390>
mov    -0x60(%rbp),%rdx
test   %edx,%edx
jle    16bf <oclogit_warp_grad+0x4bf>
shl    $0x3,%rdx
mov    -0xa8(%rbp),%rdi
mov    -0x120(%rbp),%rsi
call   10a0 <memcpy@plt>
mov    -0x30(%rbp),%rbx
lea    0x0(,%r15,8),%rax
mov    %rax,-0x38(%rbp)
lea    0x0(,%rbx,8),%rax
mov    %rax,-0x110(%rbp)
mov    %r15d,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x58(%rbp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x68(%rbp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %r15d,%eax
and    $0xfffffffe,%eax
mov    %rax,-0xe0(%rbp)
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x178(%rbp)
mov    %rax,%rsi
shr    %rsi
add    $0x1,%rsi
movabs $0x7ffffffffffffff8,%rax
ff ff 7f
add    $0x6,%rax
mov    %rcx,%rdx
mov    %rcx,-0x50(%rbp)
and    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    %rsi,%rax
mov    %rsi,-0x170(%rbp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x168(%rbp)
mov    -0xc8(%rbp),%rax
add    $0x30,%rax
mov    %rax,-0x1a8(%rbp)
mov    -0x70(%rbp),%rax
add    $0x30,%rax
mov    %rax,-0x1a0(%rbp)
movsd  0x188d(%rip),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
xor    %eax,%eax
nopl   (%rax)
mov    %rax,-0x150(%rbp)
mov    -0x48(%rbp),%rdi
mov    -0xa8(%rbp),%rbx
mov    %rbx,%rsi
mov    -0x100(%rbp),%rdx
call   1080 <innerprod_q2@plt>
movsd  %xmm0,-0x148(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
call   1060 <acos@plt>
movsd  %xmm0,-0x40(%rbp)
call   1090 <sin@plt>
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x200(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
call   1100 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1f8(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
call   1100 <cos@plt>
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1f0(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
call   1090 <sin@plt>
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1e8(%rbp)
ff
mov    -0xd0(%rbp),%rdi
mov    %edi,%esi
mov    -0x30(%rbp),%r12
mov    %r12d,%edx
lea    -0x200(%rbp),%rcx
mov    %rbx,%r8
mov    -0xd8(%rbp),%rbx
mov    %rbx,%r9
call   10d0 <product@plt>
mov    -0x48(%rbp),%rdi
mov    %rbx,%rsi
mov    -0x100(%rbp),%rdx
call   1080 <innerprod_q2@plt>
xor    %eax,%eax
mov    -0x140(%rbp),%rbx
sub    (%rbx),%eax
cvtsi2sd %eax,%xmm1
movsd  -0x148(%rbp),%xmm2
ff
mov    -0x138(%rbp),%rax
addsd  (%rax),%xmm2
movsd  %xmm0,-0xf8(%rbp)
ff
mulsd  %xmm1,%xmm2
cvttsd2si %xmm2,%edi
call   10e0 <exp@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sdl (%rbx),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x1753(%rip),%xmm0
divsd  %xmm0,%xmm1
movapd %xmm1,-0x90(%rbp)
ff
movsd  -0xf8(%rbp),%xmm2
ff
mulsd  %xmm1,%xmm2
movsd  -0x188(%rbp),%xmm0
ff
movsd  %xmm2,-0xf8(%rbp)
ff
mulsd  %xmm2,%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
call   1100 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
call   1090 <sin@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
call   1090 <sin@plt>
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
call   1100 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
mov    %r14d,%edi
mov    %r12d,%esi
mov    -0xa8(%rbp),%rdx
movsd  -0x180(%rbp),%xmm0
ff
mov    -0x158(%rbp),%rcx
call   10b0 <col_gradient@plt>
mov    %r12,%r8
xor    %r12d,%r12d
jmp    1982 <oclogit_warp_grad+0x782>
nopl   (%rax)
add    $0x1,%r12
cmp    $0x14,%r12
mov    -0x30(%rbp),%r8
je     1e90 <oclogit_warp_grad+0xc90>
test   %r8d,%r8d
mov    -0x1b0(%rbp),%rbx
jle    19ab <oclogit_warp_grad+0x7ab>
mov    %r12d,%eax
imul   %r15d,%eax
cltq
mov    -0x70(%rbp),%rcx
lea    (%rcx,%rax,8),%rsi
mov    %rbx,%rdi
mov    -0x38(%rbp),%rdx
call   10a0 <memcpy@plt>
mov    -0x48(%rbp),%rdi
mov    -0x78(%rbp),%rsi
mov    %rbx,%rdx
mov    %r13,%rcx
call   1070 <cumtrapz@plt>
test   %r14d,%r14d
mov    -0x30(%rbp),%r9
mov    -0xd0(%rbp),%r10
jle    1a60 <oclogit_warp_grad+0x860>
mov    %r12,%rax
imul   %r9,%rax
mov    -0x70(%rbp),%rcx
lea    (%rcx,%rax,8),%rax
mov    -0xd8(%rbp),%rcx
mov    -0xa8(%rbp),%rdx
mov    -0x158(%rbp),%rsi
xor    %r8d,%r8d
jmp    1a15 <oclogit_warp_grad+0x815>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
add    $0x8,%rsi
add    $0x8,%rdx
add    $0x8,%rcx
cmp    %r10,%r8
je     1a60 <oclogit_warp_grad+0x860>
test   %r9d,%r9d
jle    1a00 <oclogit_warp_grad+0x800>
xor    %ebx,%ebx
xor    %edi,%edi
xchg   %ax,%ax
movsd  (%rsi,%rbx,8),%xmm0
movsd  (%rdx,%rbx,8),%xmm1
mulsd  (%rax,%rdi,8),%xmm1
addsd  %xmm0,%xmm0
mulsd  0x0(%r13,%rdi,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x1,%rdi
add    %r14,%rbx
cmp    %rdi,%r15
jne    1a20 <oclogit_warp_grad+0x820>
jmp    1a00 <oclogit_warp_grad+0x800>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    -0x48(%rbp),%rdi
mov    -0xd8(%rbp),%rsi
mov    -0x100(%rbp),%rdx
call   1080 <innerprod_q2@plt>
movsd  %xmm0,-0xb0(%rbp)
ff
test   %r15d,%r15d
jle    1970 <oclogit_warp_grad+0x770>
mov    %r12,%rax
mov    -0x30(%rbp),%r8
imul   %r8,%rax
cmp    $0x4,%r15d
mov    -0xc8(%rbp),%rbx
jae    1ab0 <oclogit_warp_grad+0x8b0>
xor    %ecx,%ecx
jmp    1b5f <oclogit_warp_grad+0x95f>
nopw   0x0(%rax,%rax,1)
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,-0x68(%rbp)
je     1b97 <oclogit_warp_grad+0x997>
mov    -0x1a8(%rbp),%rcx
lea    (%rcx,%rax,8),%rdx
mov    -0x1a0(%rbp),%rcx
lea    (%rcx,%rax,8),%rsi
mov    -0x60(%rbp),%rdi
xor    %ecx,%ecx
nop
movupd -0x30(%rsi,%rcx,8),%xmm2
movupd -0x20(%rsi,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,-0x30(%rdx,%rcx,8)
movupd %xmm3,-0x20(%rdx,%rcx,8)
movupd -0x10(%rsi,%rcx,8),%xmm2
movupd (%rsi,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,-0x10(%rdx,%rcx,8)
movupd %xmm3,(%rdx,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdi
jne    1ae0 <oclogit_warp_grad+0x8e0>
testb  $0x1,-0x50(%rbp)
je     1b53 <oclogit_warp_grad+0x953>
add    %rax,%rcx
mov    -0x70(%rbp),%rdx
movupd (%rdx,%rcx,8),%xmm2
movupd 0x10(%rdx,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rbx,%rcx,8)
movupd %xmm3,0x10(%rbx,%rcx,8)
mov    -0x58(%rbp),%rdx
mov    %rdx,%rcx
cmp    %r15,%rdx
je     1b87 <oclogit_warp_grad+0x987>
mov    -0x70(%rbp),%rdx
lea    (%rdx,%rax,8),%rdx
lea    (%rbx,%rax,8),%rax
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r15
jne    1b70 <oclogit_warp_grad+0x970>
add    $0x1,%r12
cmp    $0x14,%r12
jne    1982 <oclogit_warp_grad+0x782>
jmp    1bb0 <oclogit_warp_grad+0x9b0>
xor    %ecx,%ecx
testb  $0x1,-0x50(%rbp)
jne    1b2e <oclogit_warp_grad+0x92e>
jmp    1b53 <oclogit_warp_grad+0x953>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
test   %r8d,%r8d
jle    1e90 <oclogit_warp_grad+0xc90>
xor    %eax,%eax
cmp    $0x1,%r8d
jne    1bdb <oclogit_warp_grad+0x9db>
mov    -0xa0(%rbp),%rdx
mov    -0xc8(%rbp),%rsi
mov    -0x110(%rbp),%rdi
jmp    1ce0 <oclogit_warp_grad+0xae0>
mov    -0xa0(%rbp),%rdx
mov    -0xc8(%rbp),%rsi
mov    -0x110(%rbp),%rdi
mov    -0xe0(%rbp),%rbx
nopw   0x0(%rax,%rax,1)
lea    (%rsi,%rax,8),%rcx
movupd (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addpd  (%rsi,%rax,8),%xmm0
addpd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
addpd  %xmm0,%xmm1
addpd  (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addpd  %xmm1,%xmm0
addpd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
addpd  %xmm0,%xmm1
addpd  (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addpd  %xmm1,%xmm0
addpd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
addpd  %xmm0,%xmm1
addpd  (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addpd  %xmm1,%xmm0
addpd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
addpd  %xmm0,%xmm1
addpd  (%rdi,%rcx,1),%xmm1
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addpd  %xmm1,%xmm0
addpd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
movupd (%rdi,%rcx,1),%xmm1
addpd  %xmm0,%xmm1
movapd %xmm1,(%rdx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rbx
jne    1c00 <oclogit_warp_grad+0xa00>
mov    %rbx,%rax
cmp    %r15,%rbx
je     1d90 <oclogit_warp_grad+0xb90>
lea    (%rsi,%rax,8),%rcx
movsd  (%rsi,%rax,8),%xmm0
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
add    %rdi,%rcx
addsd  (%rdi,%rcx,1),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    1ce0 <oclogit_warp_grad+0xae0>
test   %r8d,%r8d
jle    1e90 <oclogit_warp_grad+0xc90>
cmp    $0x4,%r8d
jae    1da6 <oclogit_warp_grad+0xba6>
xor    %eax,%eax
jmp    1e60 <oclogit_warp_grad+0xc60>
movapd -0x90(%rbp),%xmm0
ff
unpcklpd %xmm0,%xmm0
cmpq   $0x0,-0x68(%rbp)
je     243a <oclogit_warp_grad+0x123a>
mov    -0x60(%rbp),%rcx
xor    %eax,%eax
mov    -0xa0(%rbp),%rdx
mov    -0xc0(%rbp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd (%rdx,%rax,8),%xmm1
mulpd  %xmm0,%xmm1
movapd 0x10(%rdx,%rax,8),%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm1,(%rsi,%rax,8)
movapd %xmm2,0x10(%rsi,%rax,8)
movapd 0x20(%rdx,%rax,8),%xmm1
mulpd  %xmm0,%xmm1
movapd 0x30(%rdx,%rax,8),%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm1,0x20(%rsi,%rax,8)
movapd %xmm2,0x30(%rsi,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    1de0 <oclogit_warp_grad+0xbe0>
testb  $0x1,-0x50(%rbp)
je     1e48 <oclogit_warp_grad+0xc48>
movapd (%rdx,%rax,8),%xmm1
mulpd  %xmm0,%xmm1
mulpd  0x10(%rdx,%rax,8),%xmm0
movapd %xmm1,(%rsi,%rax,8)
movapd %xmm0,0x10(%rsi,%rax,8)
mov    -0x58(%rbp),%rcx
mov    %rcx,%rax
cmp    %r15,%rcx
je     1e90 <oclogit_warp_grad+0xc90>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    -0xa0(%rbp),%rcx
movsd  (%rcx,%rax,8),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
mov    -0xc0(%rbp),%rcx
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    1e60 <oclogit_warp_grad+0xc60>
nopl   0x0(%rax)
mov    -0x48(%rbp),%rdi
mov    -0xc0(%rbp),%rbx
mov    %rbx,%rsi
lea    -0x1b8(%rbp),%rdx
lea    -0xb0(%rbp),%rcx
call   10c0 <pvecnorm@plt>
movsd  -0xb0(%rbp),%xmm0
ff
mulsd  -0x130(%rbp),%xmm0
ff
call   1100 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm0,-0x90(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
mulsd  -0x130(%rbp),%xmm0
ff
call   1090 <sin@plt>
movsd  -0xb0(%rbp),%xmm0
ff
ucomisd 0x110e(%rip),%xmm0
jne    1f30 <oclogit_warp_grad+0xd30>
jp     1f30 <oclogit_warp_grad+0xd30>
mov    -0x30(%rbp),%rdx
test   %edx,%edx
mov    -0x78(%rbp),%r12
mov    -0x98(%rbp),%rsi
mov    -0x118(%rbp),%rdi
jle    2150 <oclogit_warp_grad+0xf50>
cmp    $0x4,%edx
jae    1f76 <oclogit_warp_grad+0xd76>
xor    %eax,%eax
movapd -0x90(%rbp),%xmm3
ff
jmp    2020 <oclogit_warp_grad+0xe20>
test   %r15d,%r15d
mov    -0x30(%rbp),%rdx
mov    -0x78(%rbp),%r12
mov    -0x98(%rbp),%rsi
mov    -0x118(%rbp),%rdi
mov    -0xe0(%rbp),%r8
jle    2150 <oclogit_warp_grad+0xf50>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
cmp    $0x1,%r15d
jne    203c <oclogit_warp_grad+0xe3c>
xor    %eax,%eax
movapd -0x90(%rbp),%xmm7
ff
jmp    2100 <oclogit_warp_grad+0xf00>
movapd -0x90(%rbp),%xmm3
ff
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
cmpq   $0x0,-0x68(%rbp)
je     231a <oclogit_warp_grad+0x111a>
mov    -0x60(%rbp),%rcx
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movapd (%rsi,%rax,8),%xmm1
mulpd  %xmm0,%xmm1
movapd 0x10(%rsi,%rax,8),%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm1,(%rdi,%rax,8)
movapd %xmm2,0x10(%rdi,%rax,8)
movapd 0x20(%rsi,%rax,8),%xmm1
mulpd  %xmm0,%xmm1
movapd 0x30(%rsi,%rax,8),%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm1,0x20(%rdi,%rax,8)
movapd %xmm2,0x30(%rdi,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    1fa0 <oclogit_warp_grad+0xda0>
testb  $0x1,-0x50(%rbp)
je     2008 <oclogit_warp_grad+0xe08>
movapd (%rsi,%rax,8),%xmm1
mulpd  %xmm0,%xmm1
mulpd  0x10(%rsi,%rax,8),%xmm0
movapd %xmm1,(%rdi,%rax,8)
movapd %xmm0,0x10(%rdi,%rax,8)
mov    -0x58(%rbp),%rcx
mov    %rcx,%rax
cmp    %r15,%rcx
je     2128 <oclogit_warp_grad+0xf28>
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%rdi,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    2020 <oclogit_warp_grad+0xe20>
jmp    2128 <oclogit_warp_grad+0xf28>
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
movapd %xmm1,%xmm4
unpcklpd %xmm1,%xmm4
movapd -0x90(%rbp),%xmm7
ff
movapd %xmm7,%xmm2
unpcklpd %xmm7,%xmm2
cmpq   $0x0,-0x178(%rbp)
je     2459 <oclogit_warp_grad+0x1259>
mov    -0x168(%rbp),%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd (%rbx,%rax,8),%xmm5
divpd  %xmm3,%xmm5
mulpd  %xmm4,%xmm5
movapd (%rsi,%rax,8),%xmm6
mulpd  %xmm2,%xmm6
addpd  %xmm5,%xmm6
movapd %xmm6,(%rdi,%rax,8)
movapd 0x10(%rbx,%rax,8),%xmm5
divpd  %xmm3,%xmm5
mulpd  %xmm4,%xmm5
movapd 0x10(%rsi,%rax,8),%xmm6
mulpd  %xmm2,%xmm6
addpd  %xmm5,%xmm6
movapd %xmm6,0x10(%rdi,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rcx
jne    2080 <oclogit_warp_grad+0xe80>
testb  $0x1,-0x170(%rbp)
je     20ef <oclogit_warp_grad+0xeef>
movapd (%rbx,%rax,8),%xmm5
divpd  %xmm3,%xmm5
mulpd  %xmm4,%xmm5
mulpd  (%rsi,%rax,8),%xmm2
addpd  %xmm5,%xmm2
movapd %xmm2,(%rdi,%rax,8)
mov    %r8,%rax
cmp    %r15,%r8
je     2128 <oclogit_warp_grad+0xf28>
nopw   0x0(%rax,%rax,1)
movsd  (%rbx,%rax,8),%xmm2
divsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rsi,%rax,8),%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    2100 <oclogit_warp_grad+0xf00>
test   %edx,%edx
jle    2150 <oclogit_warp_grad+0xf50>
cmp    $0x4,%edx
jae    22b4 <oclogit_warp_grad+0x10b4>
xor    %eax,%eax
mov    -0xe8(%rbp),%rbx
mov    -0xb8(%rbp),%rdx
jmp    2370 <oclogit_warp_grad+0x1170>
nopw   0x0(%rax,%rax,1)
mov    -0x48(%rbp),%rdi
mov    %r12,%rsi
mov    -0xb8(%rbp),%rdx
mov    -0xe8(%rbp),%rbx
mov    %rbx,%rcx
call   1070 <cumtrapz@plt>
mov    -0x30(%rbp),%rdx
mov    %r12,%rdi
mov    -0xf0(%rbp),%rsi
mov    %rbx,%rcx
mov    -0x128(%rbp),%r8
mov    %rdx,%r9
push   $0x0
push   $0x1
push   $0x0
push   $0x1
call   1030 <approx@plt>
add    $0x20,%rsp
xor    %eax,%eax
mov    -0x140(%rbp),%rcx
sub    (%rcx),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    -0x138(%rbp),%rax
movsd  (%rax),%xmm1
addsd  -0x148(%rbp),%xmm1
ff
mulsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
call   10e0 <exp@plt>
lea    0x1(%rax),%edi
add    $0x2,%eax
cmp    $0x3,%eax
mov    $0x0,%eax
cmovae %eax,%edi
call   1110 <log@plt>
mov    -0x150(%rbp),%rax
lea    0x1(%rax),%r12d
cmpl   $0x1,-0x108(%rbp)
jne    2204 <oclogit_warp_grad+0x1004>
lea    0xe24(%rip),%rdi
mov    %r12d,%esi
call   1050 <printf@plt>
movsd  -0x1e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0xd0(%rbp),%rdi
mov    %r14d,%esi
mov    %r15d,%edx
lea    -0x1e0(%rbp),%rcx
mov    -0x120(%rbp),%r8
mov    -0xd8(%rbp),%rbx
mov    %rbx,%r9
call   10d0 <product@plt>
mov    -0x198(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    %rbx,%rdx
mov    -0x128(%rbp),%rbx
mov    %rbx,%rcx
mov    -0xa8(%rbp),%r8
call   1040 <group_action_by_gamma@plt>
mov    -0x150(%rbp),%rdx
cmp    $0x2,%edx
setb   %al
movsd  -0x190(%rbp),%xmm0
ff
ucomisd -0xf8(%rbp),%xmm0
ff
setbe  %cl
or     %al,%cl
ucomisd -0xb0(%rbp),%xmm0
ff
setbe  %al
or     %cl,%al
cmp    $0x1,%al
jne    246d <oclogit_warp_grad+0x126d>
mov    %rbx,-0xf0(%rbp)
cmp    %edx,-0x104(%rbp)
mov    %r12d,%eax
jg     1780 <oclogit_warp_grad+0x580>
jmp    246d <oclogit_warp_grad+0x126d>
cmpq   $0x0,-0x68(%rbp)
je     232b <oclogit_warp_grad+0x112b>
mov    -0x60(%rbp),%rcx
xor    %eax,%eax
mov    -0xb8(%rbp),%rdx
nopl   0x0(%rax,%rax,1)
movapd (%rdi,%rax,8),%xmm0
movapd 0x10(%rdi,%rax,8),%xmm1
mulpd  %xmm0,%xmm0
mulpd  %xmm1,%xmm1
movapd %xmm0,(%rdx,%rax,8)
movapd %xmm1,0x10(%rdx,%rax,8)
movapd 0x20(%rdi,%rax,8),%xmm0
movapd 0x30(%rdi,%rax,8),%xmm1
mulpd  %xmm0,%xmm0
mulpd  %xmm1,%xmm1
movapd %xmm0,0x20(%rdx,%rax,8)
movapd %xmm1,0x30(%rdx,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    22d0 <oclogit_warp_grad+0x10d0>
jmp    2334 <oclogit_warp_grad+0x1134>
xor    %eax,%eax
testb  $0x1,-0x50(%rbp)
jne    1fee <oclogit_warp_grad+0xdee>
jmp    2008 <oclogit_warp_grad+0xe08>
xor    %eax,%eax
mov    -0xb8(%rbp),%rdx
testb  $0x1,-0x50(%rbp)
mov    -0xe8(%rbp),%rbx
je     235f <oclogit_warp_grad+0x115f>
movapd (%rdi,%rax,8),%xmm0
movapd 0x10(%rdi,%rax,8),%xmm1
mulpd  %xmm0,%xmm0
mulpd  %xmm1,%xmm1
movapd %xmm0,(%rdx,%rax,8)
movapd %xmm1,0x10(%rdx,%rax,8)
mov    -0x58(%rbp),%rcx
mov    %rcx,%rax
cmp    %r15,%rcx
je     2387 <oclogit_warp_grad+0x1187>
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    2370 <oclogit_warp_grad+0x1170>
mov    -0x48(%rbp),%rdi
mov    %r12,%rsi
mov    %rbx,%rcx
call   1070 <cumtrapz@plt>
test   %r15d,%r15d
mov    -0x30(%rbp),%rdx
jle    2171 <oclogit_warp_grad+0xf71>
xor    %eax,%eax
cmp    $0x1,%r15d
je     240b <oclogit_warp_grad+0x120b>
mov    -0xe0(%rbp),%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rbx),%xmm0
movsd  (%rbx,%rax,8),%xmm1
movsd  0x8(%rbx,%rax,8),%xmm2
subsd  %xmm0,%xmm1
movsd  -0x8(%rbx,%rdx,8),%xmm3
subsd  %xmm0,%xmm3
divsd  %xmm3,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
movsd  (%rbx),%xmm0
subsd  %xmm0,%xmm2
movsd  -0x8(%rbx,%rdx,8),%xmm1
subsd  %xmm0,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%rbx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rcx
jne    23c0 <oclogit_warp_grad+0x11c0>
test   $0x1,%r15b
je     2171 <oclogit_warp_grad+0xf71>
movsd  (%rbx,%rax,8),%xmm0
movsd  (%rbx),%xmm1
subsd  %xmm1,%xmm0
movsd  -0x8(%rbx,%rdx,8),%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm2,%xmm0
movsd  %xmm0,(%rbx,%rax,8)
jmp    2171 <oclogit_warp_grad+0xf71>
xor    %eax,%eax
mov    -0xa0(%rbp),%rdx
mov    -0xc0(%rbp),%rsi
testb  $0x1,-0x50(%rbp)
jne    1e2e <oclogit_warp_grad+0xc2e>
jmp    1e48 <oclogit_warp_grad+0xc48>
xor    %eax,%eax
testb  $0x1,-0x170(%rbp)
jne    20d4 <oclogit_warp_grad+0xed4>
jmp    20ef <oclogit_warp_grad+0xeef>
cmpl   $0x0,-0x30(%rbp)
jle    2483 <oclogit_warp_grad+0x1283>
mov    0x38(%rbp),%rdi
mov    %rbx,%rsi
mov    -0x38(%rbp),%rdx
call   10a0 <memcpy@plt>
mov    -0xd0(%rbp),%rax
imul   %r14d,%eax
test   %eax,%eax
je     24a8 <oclogit_warp_grad+0x12a8>
mov    %eax,%edx
shl    $0x3,%rdx
lea    -0x1e0(%rbp),%rsi
mov    0x40(%rbp),%rdi
call   10a0 <memcpy@plt>
mov    -0x160(%rbp),%rsp
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@M_PI = global i32 0
@.str = [24 x i8] c"Iteration %d : Cost %f\0A\00"
define void @oclogit_warp_grad(i32* %0, i32* %1, double* %2, double* %3, double* %4, i32* %5, i32* %6, double* %7, double* %8, double* %9, i32* %10, double* %11, double* %12) {
%14 = alloca i32*
%15 = alloca i32*
%16 = alloca double*
%17 = alloca double*
%18 = alloca double*
%19 = alloca i32*
%20 = alloca i32*
%21 = alloca double*
%22 = alloca double*
%23 = alloca double*
%24 = alloca i32*
%25 = alloca double*
%26 = alloca double*
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca i32
%38 = alloca i32
%39 = alloca i32
%40 = alloca i32
%41 = alloca double
%42 = alloca i8*
%43 = alloca i64
%44 = alloca i64
%45 = alloca i64
%46 = alloca [4 x double]
%47 = alloca i64
%48 = alloca i64
%49 = alloca [4 x double]
%50 = alloca i64
%51 = alloca [4 x double]
%52 = alloca i64
%53 = alloca i64
%54 = alloca i64
%55 = alloca i64
%56 = alloca i64
%57 = alloca i64
%58 = alloca i64
%59 = alloca i64
%60 = alloca i64
%61 = alloca i64
%62 = alloca i64
%63 = alloca double
%64 = alloca double
%65 = alloca double
%66 = alloca double
%67 = alloca double
%68 = alloca double
%69 = alloca double
%70 = alloca double
%71 = alloca double
%72 = alloca double
%73 = alloca double
%74 = alloca double
%75 = alloca double*
%76 = alloca double*
%77 = alloca double*
%78 = alloca double*
%79 = alloca double*
%80 = alloca double*
%81 = alloca double*
%82 = alloca double*
%83 = alloca double*
%84 = alloca double*
%85 = alloca double*
%86 = alloca double*
%87 = alloca double*
%88 = alloca double*
%89 = alloca double*
%90 = alloca double*
%91 = alloca double*
%92 = alloca double*
%93 = alloca double*
%94 = alloca double*
%95 = alloca double*
%96 = alloca double*
%97 = alloca i32*
store i32* %0, i32** %14
store i32* %1, i32** %15
store double* %2, double** %16
store double* %3, double** %17
store double* %4, double** %18
store i32* %5, i32** %19
store i32* %6, i32** %20
store double* %7, double** %21
store double* %8, double** %22
store double* %9, double** %23
store i32* %10, i32** %24
store double* %11, double** %25
store double* %12, double** %26
%98 = load i32*, i32** %15
%99 = load i32, i32* %98
store i32 %99, i32* %27
%100 = load i32*, i32** %14
%101 = load i32, i32* %100
store i32 %101, i32* %28
%102 = load i32*, i32** %20
%103 = load i32, i32* %102
store i32 %103, i32* %29
%104 = load double*, double** %21
%105 = load double, double* %104
store double %105, double* %30
%106 = load double*, double** %22
%107 = load double, double* %106
store double %107, double* %31
%108 = load double*, double** %23
%109 = load double, double* %108
store double %109, double* %32
%110 = load i32*, i32** %24
%111 = load i32, i32* %110
%112 = sitofp i32 %111 to double
store double %112, double* %33
store i32 20, i32* %38
store i32 10, i32* %39
store i32 0, i32* %40
store double 1.000000e+00, double* %41
%113 = load i32, i32* %27
%114 = zext i32 %113 to i64
%115 = call i8* @llvm.stacksave()
store i8* %115, i8** %42
%116 = alloca double, i64 %114
store i64 %114, i64* %43
%117 = load i32, i32* %27
%118 = zext i32 %117 to i64
%119 = alloca double, i64 %118
store i64 %118, i64* %44
%120 = load i32, i32* %27
%121 = zext i32 %120 to i64
%122 = alloca double, i64 %121
store i64 %121, i64* %45
%123 = load i32, i32* %27
%124 = load i32, i32* %38
%125 = mul nsw i32 %123, %124
%126 = zext i32 %125 to i64
%127 = alloca double, i64 %126
store i64 %126, i64* %47
%128 = load i32, i32* %27
%129 = load i32, i32* %28
%130 = mul nsw i32 %128, %129
%131 = zext i32 %130 to i64
%132 = alloca double, i64 %131
store i64 %131, i64* %48
%133 = load i32, i32* %27
%134 = load i32, i32* %28
%135 = mul nsw i32 %133, %134
%136 = zext i32 %135 to i64
%137 = alloca double, i64 %136
store i64 %136, i64* %50
%138 = load i32, i32* %27
%139 = load i32, i32* %28
%140 = mul nsw i32 %138, %139
%141 = zext i32 %140 to i64
%142 = alloca double, i64 %141
store i64 %141, i64* %52
%143 = load i32, i32* %27
%144 = zext i32 %143 to i64
%145 = alloca double, i64 %144
store i64 %144, i64* %53
%146 = load i32, i32* %27
%147 = zext i32 %146 to i64
%148 = alloca double, i64 %147
store i64 %147, i64* %54
%149 = load i32, i32* %27
%150 = zext i32 %149 to i64
%151 = alloca double, i64 %150
store i64 %150, i64* %55
%152 = load i32, i32* %27
%153 = load i32, i32* %38
%154 = mul nsw i32 %152, %153
%155 = zext i32 %154 to i64
%156 = alloca double, i64 %155
store i64 %155, i64* %56
%157 = load i32, i32* %27
%158 = zext i32 %157 to i64
%159 = alloca double, i64 %158
store i64 %158, i64* %57
%160 = load i32, i32* %27
%161 = zext i32 %160 to i64
%162 = alloca double, i64 %161
store i64 %161, i64* %58
%163 = load i32, i32* %27
%164 = zext i32 %163 to i64
%165 = alloca double, i64 %164
store i64 %164, i64* %59
%166 = load i32, i32* %27
%167 = zext i32 %166 to i64
%168 = alloca double, i64 %167
store i64 %167, i64* %60
%169 = load i32, i32* %29
%170 = add nsw i32 %169, 1
%171 = zext i32 %170 to i64
%172 = alloca double, i64 %171
store i64 %171, i64* %61
%173 = load i32, i32* %27
%174 = zext i32 %173 to i64
%175 = alloca double, i64 %174
store i64 %174, i64* %62
store double* %116, double** %75
%176 = load i32, i32* %27
%177 = load double*, double** %75
%178 = call i32 @linspace(i32 0, i32 1, i32 %176, double* %177)
%179 = load i32, i32* %27
%180 = sub nsw i32 %179, 1
%181 = sitofp i32 %180 to double
%182 = fdiv double 1.000000e+00, %181
store double %182, double* %63
store double* %119, double** %76
store i32 0, i32* %35
br label %183
183: ; preds = %198, %13
%184 = load i32, i32* %35
%185 = load i32, i32* %27
%186 = icmp slt i32 %184, %185
br i1 %186, label %187, label %201
187: ; preds = %183
%188 = load i32, i32* %35
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds double, double* %116, i64 %189
%191 = load double, double* %190
%192 = load i32, i32* %35
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds double, double* %119, i64 %193
store double %191, double* %194
%195 = load i32, i32* %35
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds double, double* %122, i64 %196
store double 1.000000e+00, double* %197
br label %198
198: ; preds = %187
%199 = load i32, i32* %35
%200 = add nsw i32 %199, 1
store i32 %200, i32* %35
br label %183
201: ; preds = %183
%202 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 0
store double 1.000000e+00, double* %202
%203 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 1
store double 0.000000e+00, double* %203
%204 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 2
store double 0.000000e+00, double* %204
%205 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 3
store double 1.000000e+00, double* %205
store double* %127, double** %77
store i32 0, i32* %35
br label %206
206: ; preds = %271, %201
%207 = load i32, i32* %35
%208 = load i32, i32* %39
%209 = icmp slt i32 %207, %208
br i1 %209, label %210, label %274
210: ; preds = %206
store i32 0, i32* %34
br label %211
211: ; preds = %267, %210
%212 = load i32, i32* %34
%213 = load i32, i32* %27
%214 = icmp slt i32 %212, %213
br i1 %214, label %215, label %270
215: ; preds = %211
%216 = load i32, i32* @M_PI
%217 = call i32 @sqrt(i32 %216)
%218 = sdiv i32 1, %217
%219 = load i32, i32* @M_PI
%220 = mul nsw i32 2, %219
%221 = load i32, i32* %35
%222 = add nsw i32 %221, 1
%223 = mul nsw i32 %220, %222
%224 = sitofp i32 %223 to double
%225 = load i32, i32* %34
%226 = sext i32 %225 to i64
%227 = getelementptr inbounds double, double* %116, i64 %226
%228 = load double, double* %227
%229 = fmul double %224, %228
%230 = call i32 @sin(double %229)
%231 = mul nsw i32 %218, %230
%232 = sitofp i32 %231 to double
%233 = load i32, i32* %34
%234 = load i32, i32* %35
%235 = mul nsw i32 2, %234
%236 = load i32, i32* %27
%237 = mul nsw i32 %235, %236
%238 = add nsw i32 %233, %237
%239 = sext i32 %238 to i64
%240 = getelementptr inbounds double, double* %127, i64 %239
store double %232, double* %240
%241 = load i32, i32* @M_PI
%242 = call i32 @sqrt(i32 %241)
%243 = sdiv i32 1, %242
%244 = load i32, i32* @M_PI
%245 = mul nsw i32 2, %244
%246 = load i32, i32* %35
%247 = add nsw i32 %246, 1
%248 = mul nsw i32 %245, %247
%249 = sitofp i32 %248 to double
%250 = load i32, i32* %34
%251 = sext i32 %250 to i64
%252 = getelementptr inbounds double, double* %116, i64 %251
%253 = load double, double* %252
%254 = fmul double %249, %253
%255 = call i32 @cos(double %254)
%256 = mul nsw i32 %243, %255
%257 = sitofp i32 %256 to double
%258 = load i32, i32* %34
%259 = load i32, i32* %35
%260 = mul nsw i32 2, %259
%261 = add nsw i32 %260, 1
%262 = load i32, i32* %27
%263 = mul nsw i32 %261, %262
%264 = add nsw i32 %258, %263
%265 = sext i32 %264 to i64
%266 = getelementptr inbounds double, double* %127, i64 %265
store double %257, double* %266
br label %267
267: ; preds = %215
%268 = load i32, i32* %34
%269 = add nsw i32 %268, 1
store i32 %269, i32* %34
br label %211
270: ; preds = %211
br label %271
271: ; preds = %270
%272 = load i32, i32* %35
%273 = add nsw i32 %272, 1
store i32 %273, i32* %35
br label %206
274: ; preds = %206
store double* %132, double** %78
store i32 0, i32* %35
br label %275
275: ; preds = %290, %274
%276 = load i32, i32* %35
%277 = load i32, i32* %27
%278 = load i32, i32* %28
%279 = mul nsw i32 %277, %278
%280 = icmp slt i32 %276, %279
br i1 %280, label %281, label %293
281: ; preds = %275
%282 = load double*, double** %18
%283 = load i32, i32* %35
%284 = sext i32 %283 to i64
%285 = getelementptr inbounds double, double* %282, i64 %284
%286 = load double, double* %285
%287 = load i32, i32* %35
%288 = sext i32 %287 to i64
%289 = getelementptr inbounds double, double* %132, i64 %288
store double %286, double* %289
br label %290
290: ; preds = %281
%291 = load i32, i32* %35
%292 = add nsw i32 %291, 1
store i32 %292, i32* %35
br label %275
293: ; preds = %275
br label %294
294: ; preds = %790, %293
%295 = load double*, double** %17
store double* %295, double** %80
store double* %132, double** %78
%296 = load i32*, i32** %15
%297 = load double*, double** %78
%298 = load double*, double** %80
%299 = call double @innerprod_q2(i32* %296, double* %297, double* %298)
store double %299, double* %64
%300 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 0
%301 = load double, double* %300
%302 = call double @acos(double %301)
store double %302, double* %65
%303 = load double, double* %65
%304 = call i32 @sin(double %303)
%305 = mul nsw i32 -1, %304
%306 = sitofp i32 %305 to double
%307 = getelementptr inbounds [4 x double], [4 x double]* %49, i64 0, i64 0
store double %306, double* %307
%308 = load double, double* %65
%309 = call i32 @cos(double %308)
%310 = sitofp i32 %309 to double
%311 = getelementptr inbounds [4 x double], [4 x double]* %49, i64 0, i64 1
store double %310, double* %311
%312 = load double, double* %65
%313 = call i32 @cos(double %312)
%314 = mul nsw i32 -1, %313
%315 = sitofp i32 %314 to double
%316 = getelementptr inbounds [4 x double], [4 x double]* %49, i64 0, i64 2
store double %315, double* %316
%317 = load double, double* %65
%318 = call i32 @sin(double %317)
%319 = mul nsw i32 -1, %318
%320 = sitofp i32 %319 to double
%321 = getelementptr inbounds [4 x double], [4 x double]* %49, i64 0, i64 3
store double %320, double* %321
%322 = getelementptr inbounds [4 x double], [4 x double]* %49, i64 0, i64 0
store double* %322, double** %81
store double* %132, double** %78
store double* %137, double** %82
%323 = load i32, i32* %28
%324 = load i32, i32* %28
%325 = load i32, i32* %27
%326 = load double*, double** %81
%327 = load double*, double** %78
%328 = load double*, double** %82
%329 = call i32 @product(i32 %323, i32 %324, i32 %325, double* %326, double* %327, double* %328)
%330 = load double*, double** %17
store double* %330, double** %80
%331 = load i32*, i32** %15
%332 = load double*, double** %82
%333 = load double*, double** %80
%334 = call double @innerprod_q2(i32* %331, double* %332, double* %333)
store double %334, double* %66
%335 = load i32*, i32** %19
store i32* %335, i32** %97
%336 = load double*, double** %16
store double* %336, double** %83
%337 = load i32*, i32** %97
%338 = getelementptr inbounds i32, i32* %337, i64 0
%339 = load i32, i32* %338
%340 = mul nsw i32 -1, %339
%341 = sitofp i32 %340 to double
%342 = load double*, double** %83
%343 = getelementptr inbounds double, double* %342, i64 0
%344 = load double, double* %343
%345 = load double, double* %64
%346 = fadd double %344, %345
%347 = fmul double %341, %346
%348 = fptosi double %347 to i32
%349 = call i32 @exp(i32 %348)
%350 = sitofp i32 %349 to double
store double %350, double* %67
%351 = load i32*, i32** %97
%352 = getelementptr inbounds i32, i32* %351, i64 0
%353 = load i32, i32* %352
%354 = sitofp i32 %353 to double
%355 = load double, double* %67
%356 = fmul double %354, %355
%357 = load double, double* %67
%358 = fadd double 1.000000e+00, %357
%359 = fdiv double %356, %358
store double %359, double* %68
%360 = load double, double* %68
%361 = load double, double* %66
%362 = fmul double %360, %361
store double %362, double* %74
%363 = load double, double* %65
%364 = load double, double* %31
%365 = load double, double* %74
%366 = call double @llvm.fmuladd.f64(double %364, double %365, double %363)
store double %366, double* %69
%367 = load double, double* %69
%368 = call i32 @cos(double %367)
%369 = sitofp i32 %368 to double
%370 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 0
store double %369, double* %370
%371 = load double, double* %69
%372 = call i32 @sin(double %371)
%373 = sitofp i32 %372 to double
%374 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 1
store double %373, double* %374
%375 = load double, double* %69
%376 = call i32 @sin(double %375)
%377 = mul nsw i32 -1, %376
%378 = sitofp i32 %377 to double
%379 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 2
store double %378, double* %379
%380 = load double, double* %69
%381 = call i32 @cos(double %380)
%382 = sitofp i32 %381 to double
%383 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 3
store double %382, double* %383
store double* %142, double** %84
store double* %132, double** %78
%384 = load i32, i32* %28
%385 = load i32, i32* %27
%386 = load double*, double** %78
%387 = load double, double* %63
%388 = load double*, double** %84
%389 = call i32 @col_gradient(i32 %384, i32 %385, double* %386, double %387, double* %388)
store double* %145, double** %87
store double* %148, double** %85
store double* %137, double** %82
store double* %151, double** %86
%390 = load double*, double** %17
store double* %390, double** %80
store double* %156, double** %88
store i32 0, i32* %35
br label %391
391: ; preds = %507, %294
%392 = load i32, i32* %35
%393 = load i32, i32* %38
%394 = icmp slt i32 %392, %393
br i1 %394, label %395, label %510
395: ; preds = %391
store i32 0, i32* %36
br label %396
396: ; preds = %412, %395
%397 = load i32, i32* %36
%398 = load i32, i32* %27
%399 = icmp slt i32 %397, %398
br i1 %399, label %400, label %415
400: ; preds = %396
%401 = load i32, i32* %36
%402 = load i32, i32* %35
%403 = load i32, i32* %27
%404 = mul nsw i32 %402, %403
%405 = add nsw i32 %401, %404
%406 = sext i32 %405 to i64
%407 = getelementptr inbounds double, double* %127, i64 %406
%408 = load double, double* %407
%409 = load i32, i32* %36
%410 = sext i32 %409 to i64
%411 = getelementptr inbounds double, double* %148, i64 %410
store double %408, double* %411
br label %412
412: ; preds = %400
%413 = load i32, i32* %36
%414 = add nsw i32 %413, 1
store i32 %414, i32* %36
br label %396
415: ; preds = %396
%416 = load i32*, i32** %15
%417 = load double*, double** %75
%418 = load double*, double** %85
%419 = load double*, double** %87
%420 = call i32 @cumtrapz(i32* %416, double* %417, double* %418, double* %419)
store i32 0, i32* %37
br label %421
421: ; preds = %473, %415
%422 = load i32, i32* %37
%423 = load i32, i32* %28
%424 = icmp slt i32 %422, %423
br i1 %424, label %425, label %476
425: ; preds = %421
store i32 0, i32* %36
br label %426
426: ; preds = %469, %425
%427 = load i32, i32* %36
%428 = load i32, i32* %27
%429 = icmp slt i32 %427, %428
br i1 %429, label %430, label %472
430: ; preds = %426
%431 = load i32, i32* %28
%432 = load i32, i32* %36
%433 = mul nsw i32 %431, %432
%434 = load i32, i32* %37
%435 = add nsw i32 %433, %434
%436 = sext i32 %435 to i64
%437 = getelementptr inbounds double, double* %142, i64 %436
%438 = load double, double* %437
%439 = fmul double 2.000000e+00, %438
%440 = load i32, i32* %36
%441 = sext i32 %440 to i64
%442 = getelementptr inbounds double, double* %145, i64 %441
%443 = load double, double* %442
%444 = load i32, i32* %28
%445 = load i32, i32* %36
%446 = mul nsw i32 %444, %445
%447 = load i32, i32* %37
%448 = add nsw i32 %446, %447
%449 = sext i32 %448 to i64
%450 = getelementptr inbounds double, double* %132, i64 %449
%451 = load double, double* %450
%452 = load i32, i32* %36
%453 = load i32, i32* %35
%454 = load i32, i32* %27
%455 = mul nsw i32 %453, %454
%456 = add nsw i32 %452, %455
%457 = sext i32 %456 to i64
%458 = getelementptr inbounds double, double* %127, i64 %457
%459 = load double, double* %458
%460 = fmul double %451, %459
%461 = call double @llvm.fmuladd.f64(double %439, double %443, double %460)
%462 = load i32, i32* %28
%463 = load i32, i32* %36
%464 = mul nsw i32 %462, %463
%465 = load i32, i32* %37
%466 = add nsw i32 %464, %465
%467 = sext i32 %466 to i64
%468 = getelementptr inbounds double, double* %137, i64 %467
store double %461, double* %468
br label %469
469: ; preds = %430
%470 = load i32, i32* %36
%471 = add nsw i32 %470, 1
store i32 %471, i32* %36
br label %426
472: ; preds = %426
br label %473
473: ; preds = %472
%474 = load i32, i32* %37
%475 = add nsw i32 %474, 1
store i32 %475, i32* %37
br label %421
476: ; preds = %421
%477 = load i32*, i32** %15
%478 = load double*, double** %82
%479 = load double*, double** %80
%480 = call double @innerprod_q2(i32* %477, double* %478, double* %479)
store double %480, double* %70
store i32 0, i32* %36
br label %481
481: ; preds = %503, %476
%482 = load i32, i32* %36
%483 = load i32, i32* %27
%484 = icmp slt i32 %482, %483
br i1 %484, label %485, label %506
485: ; preds = %481
%486 = load double, double* %70
%487 = load i32, i32* %36
%488 = load i32, i32* %35
%489 = load i32, i32* %27
%490 = mul nsw i32 %488, %489
%491 = add nsw i32 %487, %490
%492 = sext i32 %491 to i64
%493 = getelementptr inbounds double, double* %127, i64 %492
%494 = load double, double* %493
%495 = fmul double %486, %494
%496 = load i32, i32* %36
%497 = load i32, i32* %35
%498 = load i32, i32* %27
%499 = mul nsw i32 %497, %498
%500 = add nsw i32 %496, %499
%501 = sext i32 %500 to i64
%502 = getelementptr inbounds double, double* %156, i64 %501
store double %495, double* %502
br label %503
503: ; preds = %485
%504 = load i32, i32* %36
%505 = add nsw i32 %504, 1
store i32 %505, i32* %36
br label %481
506: ; preds = %481
br label %507
507: ; preds = %506
%508 = load i32, i32* %35
%509 = add nsw i32 %508, 1
store i32 %509, i32* %35
br label %391
510: ; preds = %391
store i32 0, i32* %37
br label %511
511: ; preds = %548, %510
%512 = load i32, i32* %37
%513 = load i32, i32* %27
%514 = icmp slt i32 %512, %513
br i1 %514, label %515, label %551
515: ; preds = %511
%516 = load i32, i32* %37
%517 = sext i32 %516 to i64
%518 = getelementptr inbounds double, double* %156, i64 %517
%519 = load double, double* %518
%520 = load i32, i32* %37
%521 = sext i32 %520 to i64
%522 = getelementptr inbounds double, double* %151, i64 %521
store double %519, double* %522
store i32 1, i32* %36
br label %523
523: ; preds = %544, %515
%524 = load i32, i32* %36
%525 = load i32, i32* %38
%526 = icmp slt i32 %524, %525
br i1 %526, label %527, label %547
527: ; preds = %523
%528 = load i32, i32* %37
%529 = sext i32 %528 to i64
%530 = getelementptr inbounds double, double* %151, i64 %529
%531 = load double, double* %530
%532 = load i32, i32* %37
%533 = load i32, i32* %36
%534 = load i32, i32* %27
%535 = mul nsw i32 %533, %534
%536 = add nsw i32 %532, %535
%537 = sext i32 %536 to i64
%538 = getelementptr inbounds double, double* %156, i64 %537
%539 = load double, double* %538
%540 = fadd double %531, %539
%541 = load i32, i32* %37
%542 = sext i32 %541 to i64
%543 = getelementptr inbounds double, double* %151, i64 %542
store double %540, double* %543
br label %544
544: ; preds = %527
%545 = load i32, i32* %36
%546 = add nsw i32 %545, 1
store i32 %546, i32* %36
br label %523
547: ; preds = %523
br label %548
548: ; preds = %547
%549 = load i32, i32* %37
%550 = add nsw i32 %549, 1
store i32 %550, i32* %37
br label %511
551: ; preds = %511
store double* %159, double** %89
store double* %151, double** %86
store i32 0, i32* %35
br label %552
552: ; preds = %568, %551
%553 = load i32, i32* %35
%554 = load i32, i32* %27
%555 = icmp slt i32 %553, %554
br i1 %555, label %556, label %571
556: ; preds = %552
%557 = load double, double* %68
%558 = load double*, double** %86
%559 = load i32, i32* %35
%560 = sext i32 %559 to i64
%561 = getelementptr inbounds double, double* %558, i64 %560
%562 = load double, double* %561
%563 = fmul double %557, %562
%564 = load double*, double** %89
%565 = load i32, i32* %35
%566 = sext i32 %565 to i64
%567 = getelementptr inbounds double, double* %564, i64 %566
store double %563, double* %567
br label %568
568: ; preds = %556
%569 = load i32, i32* %35
%570 = add nsw i32 %569, 1
store i32 %570, i32* %35
br label %552
571: ; preds = %552
store double* %162, double** %90
store double* %122, double** %91
store double* %165, double** %92
store double* %168, double** %93
store double* %159, double** %89
%572 = load i32*, i32** %15
%573 = load double*, double** %89
%574 = call i32 @pvecnorm(i32* %572, double* %573, double* %41, double* %70)
%575 = load double, double* %32
%576 = load double, double* %70
%577 = fmul double %575, %576
%578 = call i32 @cos(double %577)
%579 = sitofp i32 %578 to double
store double %579, double* %72
%580 = load double, double* %32
%581 = load double, double* %70
%582 = fmul double %580, %581
%583 = call i32 @sin(double %582)
%584 = sitofp i32 %583 to double
store double %584, double* %73
%585 = load double, double* %70
%586 = fcmp oeq double %585, 0.000000e+00
br i1 %586, label %587, label %608
587: ; preds = %571
store i32 0, i32* %34
br label %588
588: ; preds = %604, %587
%589 = load i32, i32* %34
%590 = load i32, i32* %27
%591 = icmp slt i32 %589, %590
br i1 %591, label %592, label %607
592: ; preds = %588
%593 = load double, double* %72
%594 = load double*, double** %91
%595 = load i32, i32* %34
%596 = sext i32 %595 to i64
%597 = getelementptr inbounds double, double* %594, i64 %596
%598 = load double, double* %597
%599 = fmul double %593, %598
%600 = load double*, double** %90
%601 = load i32, i32* %34
%602 = sext i32 %601 to i64
%603 = getelementptr inbounds double, double* %600, i64 %602
store double %599, double* %603
br label %604
604: ; preds = %592
%605 = load i32, i32* %34
%606 = add nsw i32 %605, 1
store i32 %606, i32* %34
br label %588
607: ; preds = %588
br label %638
608: ; preds = %571
store i32 0, i32* %34
br label %609
609: ; preds = %634, %608
%610 = load i32, i32* %34
%611 = load i32, i32* %27
%612 = icmp slt i32 %610, %611
br i1 %612, label %613, label %637
613: ; preds = %609
%614 = load double, double* %72
%615 = load double*, double** %91
%616 = load i32, i32* %34
%617 = sext i32 %616 to i64
%618 = getelementptr inbounds double, double* %615, i64 %617
%619 = load double, double* %618
%620 = load double, double* %73
%621 = load double*, double** %89
%622 = load i32, i32* %34
%623 = sext i32 %622 to i64
%624 = getelementptr inbounds double, double* %621, i64 %623
%625 = load double, double* %624
%626 = load double, double* %70
%627 = fdiv double %625, %626
%628 = fmul double %620, %627
%629 = call double @llvm.fmuladd.f64(double %614, double %619, double %628)
%630 = load double*, double** %90
%631 = load i32, i32* %34
%632 = sext i32 %631 to i64
%633 = getelementptr inbounds double, double* %630, i64 %632
store double %629, double* %633
br label %634
634: ; preds = %613
%635 = load i32, i32* %34
%636 = add nsw i32 %635, 1
store i32 %636, i32* %34
br label %609
637: ; preds = %609
br label %638
638: ; preds = %637, %607
store i32 0, i32* %34
br label %639
639: ; preds = %658, %638
%640 = load i32, i32* %34
%641 = load i32, i32* %27
%642 = icmp slt i32 %640, %641
br i1 %642, label %643, label %661
643: ; preds = %639
%644 = load double*, double** %90
%645 = load i32, i32* %34
%646 = sext i32 %645 to i64
%647 = getelementptr inbounds double, double* %644, i64 %646
%648 = load double, double* %647
%649 = load double*, double** %90
%650 = load i32, i32* %34
%651 = sext i32 %650 to i64
%652 = getelementptr inbounds double, double* %649, i64 %651
%653 = load double, double* %652
%654 = fmul double %648, %653
%655 = load i32, i32* %34
%656 = sext i32 %655 to i64
%657 = getelementptr inbounds double, double* %175, i64 %656
store double %654, double* %657
br label %658
658: ; preds = %643
%659 = load i32, i32* %34
%660 = add nsw i32 %659, 1
store i32 %660, i32* %34
br label %639
661: ; preds = %639
store double* %175, double** %94
%662 = load i32*, i32** %15
%663 = load double*, double** %75
%664 = load double*, double** %94
%665 = load double*, double** %93
%666 = call i32 @cumtrapz(i32* %662, double* %663, double* %664, double* %665)
store i32 0, i32* %34
br label %667
667: ; preds = %696, %661
%668 = load i32, i32* %34
%669 = load i32, i32* %27
%670 = icmp slt i32 %668, %669
br i1 %670, label %671, label %699
671: ; preds = %667
%672 = load double*, double** %93
%673 = load i32, i32* %34
%674 = sext i32 %673 to i64
%675 = getelementptr inbounds double, double* %672, i64 %674
%676 = load double, double* %675
%677 = load double*, double** %93
%678 = getelementptr inbounds double, double* %677, i64 0
%679 = load double, double* %678
%680 = fsub double %676, %679
%681 = load double*, double** %93
%682 = load i32, i32* %27
%683 = sub nsw i32 %682, 1
%684 = sext i32 %683 to i64
%685 = getelementptr inbounds double, double* %681, i64 %684
%686 = load double, double* %685
%687 = load double*, double** %93
%688 = getelementptr inbounds double, double* %687, i64 0
%689 = load double, double* %688
%690 = fsub double %686, %689
%691 = fdiv double %680, %690
%692 = load double*, double** %93
%693 = load i32, i32* %34
%694 = sext i32 %693 to i64
%695 = getelementptr inbounds double, double* %692, i64 %694
store double %691, double* %695
br label %696
696: ; preds = %671
%697 = load i32, i32* %34
%698 = add nsw i32 %697, 1
store i32 %698, i32* %34
br label %667
699: ; preds = %667
%700 = load double*, double** %76
%701 = load i32, i32* %27
%702 = load double*, double** %93
%703 = load double*, double** %92
%704 = load i32, i32* %27
%705 = call i32 @approx(double* %116, double* %700, i32 %701, double* %702, double* %703, i32 %704, i32 1, i32 0, i32 1, i32 0)
%706 = load i32*, i32** %19
store i32* %706, i32** %97
%707 = load double*, double** %16
store double* %707, double** %83
%708 = load i32*, i32** %97
%709 = getelementptr inbounds i32, i32* %708, i64 0
%710 = load i32, i32* %709
%711 = mul nsw i32 -1, %710
%712 = sitofp i32 %711 to double
%713 = load double*, double** %83
%714 = getelementptr inbounds double, double* %713, i64 0
%715 = load double, double* %714
%716 = load double, double* %64
%717 = fadd double %715, %716
%718 = fmul double %712, %717
%719 = fptosi double %718 to i32
%720 = call i32 @exp(i32 %719)
%721 = add nsw i32 1, %720
%722 = sdiv i32 1, %721
%723 = call double @log(i32 %722)
%724 = load i32, i32* %40
%725 = sext i32 %724 to i64
%726 = getelementptr inbounds double, double* %172, i64 %725
store double %723, double* %726
%727 = load double, double* %33
%728 = fcmp oeq double %727, 1.000000e+00
br i1 %728, label %729, label %737
729: ; preds = %699
%730 = load i32, i32* %40
%731 = add nsw i32 %730, 1
%732 = load i32, i32* %40
%733 = sext i32 %732 to i64
%734 = getelementptr inbounds double, double* %172, i64 %733
%735 = load double, double* %734
%736 = call i32 @printf(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i32 %731, double %735)
br label %737
737: ; preds = %729, %699
store double* %165, double** %76
%738 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 0
store double* %738, double** %95
%739 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 0
%740 = load double, double* %739
%741 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 0
store double %740, double* %741
%742 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 1
%743 = load double, double* %742
%744 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 1
store double %743, double* %744
%745 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 2
%746 = load double, double* %745
%747 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 2
store double %746, double* %747
%748 = getelementptr inbounds [4 x double], [4 x double]* %51, i64 0, i64 3
%749 = load double, double* %748
%750 = getelementptr inbounds [4 x double], [4 x double]* %46, i64 0, i64 3
store double %749, double* %750
store double* %132, double** %78
store double* %137, double** %82
%751 = load double*, double** %18
store double* %751, double** %96
%752 = load i32, i32* %28
%753 = load i32, i32* %28
%754 = load i32, i32* %27
%755 = load double*, double** %95
%756 = load double*, double** %96
%757 = load double*, double** %82
%758 = call i32 @product(i32 %752, i32 %753, i32 %754, double* %755, double* %756, double* %757)
%759 = load i32*, i32** %14
%760 = load i32*, i32** %15
%761 = load double*, double** %82
%762 = load double*, double** %76
%763 = load double*, double** %78
%764 = call i32 @group_action_by_gamma(i32* %759, i32* %760, double* %761, double* %762, double* %763)
%765 = load i32, i32* %40
%766 = icmp sge i32 %765, 2
br i1 %766, label %767, label %787
767: ; preds = %737
%768 = load i32, i32* %40
%769 = sext i32 %768 to i64
%770 = getelementptr inbounds double, double* %172, i64 %769
%771 = load double, double* %770
%772 = load i32, i32* %40
%773 = sub nsw i32 %772, 1
%774 = sext i32 %773 to i64
%775 = getelementptr inbounds double, double* %172, i64 %774
%776 = load double, double* %775
%777 = fsub double %771, %776
store double %777, double* %71
%778 = load double, double* %74
%779 = load double, double* %30
%780 = fcmp olt double %778, %779
br i1 %780, label %781, label %786
781: ; preds = %767
%782 = load double, double* %70
%783 = load double, double* %30
%784 = fcmp olt double %782, %783
br i1 %784, label %785, label %786
785: ; preds = %781
br label %794
786: ; preds = %781, %767
br label %787
787: ; preds = %786, %737
%788 = load i32, i32* %40
%789 = add nsw i32 %788, 1
store i32 %789, i32* %40
br label %790
790: ; preds = %787
%791 = load i32, i32* %29
%792 = load i32, i32* %40
%793 = icmp sge i32 %791, %792
br i1 %793, label %294, label %794
794: ; preds = %790, %785
store i32 0, i32* %35
br label %795
795: ; preds = %809, %794
%796 = load i32, i32* %35
%797 = load i32, i32* %27
%798 = icmp slt i32 %796, %797
br i1 %798, label %799, label %812
799: ; preds = %795
%800 = load double*, double** %76
%801 = load i32, i32* %35
%802 = sext i32 %801 to i64
%803 = getelementptr inbounds double, double* %800, i64 %802
%804 = load double, double* %803
%805 = load double*, double** %25
%806 = load i32, i32* %35
%807 = sext i32 %806 to i64
%808 = getelementptr inbounds double, double* %805, i64 %807
store double %804, double* %808
br label %809
809: ; preds = %799
%810 = load i32, i32* %35
%811 = add nsw i32 %810, 1
store i32 %811, i32* %35
br label %795
812: ; preds = %795
store i32 0, i32* %35
br label %813
813: ; preds = %829, %812
%814 = load i32, i32* %35
%815 = load i32, i32* %28
%816 = load i32, i32* %28
%817 = mul nsw i32 %815, %816
%818 = icmp slt i32 %814, %817
br i1 %818, label %819, label %832
819: ; preds = %813
%820 = load double*, double** %95
%821 = load i32, i32* %35
%822 = sext i32 %821 to i64
%823 = getelementptr inbounds double, double* %820, i64 %822
%824 = load double, double* %823
%825 = load double*, double** %26
%826 = load i32, i32* %35
%827 = sext i32 %826 to i64
%828 = getelementptr inbounds double, double* %825, i64 %827
store double %824, double* %828
br label %829
829: ; preds = %819
%830 = load i32, i32* %35
%831 = add nsw i32 %830, 1
store i32 %831, i32* %35
br label %813
832: ; preds = %813
%833 = load i8*, i8** %42
ret void
}

LLVM IR (-O2):
@M_PI = global i32 0
@.str = [24 x i8] c"Iteration %d : Cost %f\0A\00"
define void @oclogit_warp_grad(i32* %0, i32* %1, double* nocapture %2, double* %3, double* %4, i32* nocapture %5, i32* nocapture %6, double* nocapture %7, double* nocapture %8, double* nocapture %9, i32* nocapture %10, double* nocapture writeonly %11, double* nocapture writeonly %12) {
%14 = bitcast double* %4 to i8*
%15 = bitcast double* %11 to i8*
%16 = bitcast double* %12 to i8*
%17 = alloca double
%18 = alloca [4 x double]
%19 = alloca [4 x double]
%20 = bitcast [4 x double]* %19 to i8*
%21 = alloca double
%22 = load i32, i32* %1
%23 = load i32, i32* %0
%24 = load i32, i32* %6
%25 = load double, double* %7
%26 = load double, double* %8
%27 = load double, double* %9
%28 = load i32, i32* %10
%29 = bitcast double* %17 to i8*
store double 1.000000e+00, double* %17
%30 = zext i32 %22 to i64
%31 = tail call i8* @llvm.stacksave()
%32 = alloca double, i64 %30
%33 = alloca double, i64 %30
%34 = alloca double, i64 %30
%35 = mul nsw i32 %22, 20
%36 = zext i32 %35 to i64
%37 = alloca double, i64 %36
%38 = mul i32 %23, %22
%39 = zext i32 %38 to i64
%40 = alloca double, i64 %39
%41 = bitcast double* %40 to i8*
%42 = bitcast [4 x double]* %18 to i8*
%43 = alloca double, i64 %39
%44 = alloca double, i64 %39
%45 = alloca double, i64 %30
%46 = alloca double, i64 %30
%47 = bitcast double* %46 to i8*
%48 = alloca double, i64 %30
%49 = alloca double, i64 %36
%50 = alloca double, i64 %30
%51 = alloca double, i64 %30
%52 = alloca double, i64 %30
%53 = bitcast double* %52 to i8*
%54 = alloca double, i64 %30
%55 = alloca double, i64 %30
%56 = bitcast double* %21 to i8*
%57 = call i32 @linspace(i32 0, i32 1, i32 %22, double* %32)
%58 = add nsw i32 %22, -1
%59 = sitofp i32 %58 to double
%60 = fdiv double 1.000000e+00, %59
%61 = icmp sgt i32 %22, 0
br i1 %61, label %62, label %138
62: ; preds = %13
%63 = bitcast double* %33 to i8*
%64 = bitcast double* %32 to i8*
%65 = shl nuw nsw i64 %30, 3
%66 = icmp ult i32 %22, 4
br i1 %66, label %136, label %67
67: ; preds = %62
%68 = and i64 %30, 4294967292
%69 = add nsw i64 %68, -4
%70 = lshr exact i64 %69, 2
%71 = add nuw nsw i64 %70, 1
%72 = and i64 %71, 7
%73 = icmp ult i64 %69, 28
br i1 %73, label %121, label %74
74: ; preds = %67
%75 = and i64 %71, 9223372036854775800
br label %76
76: ; preds = %76, %74
%77 = phi i64 [ 0, %74 ], [ %118, %76 ]
%78 = phi i64 [ 0, %74 ], [ %119, %76 ]
%79 = getelementptr inbounds double, double* %34, i64 %77
%80 = bitcast double* %79 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %80
%81 = getelementptr inbounds double, double* %79, i64 2
%82 = bitcast double* %81 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %82
%83 = or i64 %77, 4
%84 = getelementptr inbounds double, double* %34, i64 %83
%85 = bitcast double* %84 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %85
%86 = getelementptr inbounds double, double* %84, i64 2
%87 = bitcast double* %86 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %87
%88 = or i64 %77, 8
%89 = getelementptr inbounds double, double* %34, i64 %88
%90 = bitcast double* %89 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %90
%91 = getelementptr inbounds double, double* %89, i64 2
%92 = bitcast double* %91 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %92
%93 = or i64 %77, 12
%94 = getelementptr inbounds double, double* %34, i64 %93
%95 = bitcast double* %94 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %95
%96 = getelementptr inbounds double, double* %94, i64 2
%97 = bitcast double* %96 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %97
%98 = or i64 %77, 16
%99 = getelementptr inbounds double, double* %34, i64 %98
%100 = bitcast double* %99 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %100
%101 = getelementptr inbounds double, double* %99, i64 2
%102 = bitcast double* %101 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %102
%103 = or i64 %77, 20
%104 = getelementptr inbounds double, double* %34, i64 %103
%105 = bitcast double* %104 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %105
%106 = getelementptr inbounds double, double* %104, i64 2
%107 = bitcast double* %106 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %107
%108 = or i64 %77, 24
%109 = getelementptr inbounds double, double* %34, i64 %108
%110 = bitcast double* %109 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %110
%111 = getelementptr inbounds double, double* %109, i64 2
%112 = bitcast double* %111 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %112
%113 = or i64 %77, 28
%114 = getelementptr inbounds double, double* %34, i64 %113
%115 = bitcast double* %114 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %115
%116 = getelementptr inbounds double, double* %114, i64 2
%117 = bitcast double* %116 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %117
%118 = add nuw i64 %77, 32
%119 = add i64 %78, 8
%120 = icmp eq i64 %119, %75
br i1 %120, label %121, label %76
121: ; preds = %76, %67
%122 = phi i64 [ 0, %67 ], [ %118, %76 ]
%123 = icmp eq i64 %72, 0
br i1 %123, label %134, label %124
124: ; preds = %121, %124
%125 = phi i64 [ %131, %124 ], [ %122, %121 ]
%126 = phi i64 [ %132, %124 ], [ 0, %121 ]
%127 = getelementptr inbounds double, double* %34, i64 %125
%128 = bitcast double* %127 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %128
%129 = getelementptr inbounds double, double* %127, i64 2
%130 = bitcast double* %129 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %130
%131 = add nuw i64 %125, 4
%132 = add i64 %126, 1
%133 = icmp eq i64 %132, %72
br i1 %133, label %134, label %124
134: ; preds = %124, %121
%135 = icmp eq i64 %68, %30
br i1 %135, label %138, label %136
136: ; preds = %62, %134
%137 = phi i64 [ 0, %62 ], [ %68, %134 ]
br label %140
138: ; preds = %140, %134, %13
%139 = sext i32 %22 to i64
br label %145
140: ; preds = %136, %140
%141 = phi i64 [ %143, %140 ], [ %137, %136 ]
%142 = getelementptr inbounds double, double* %34, i64 %141
store double 1.000000e+00, double* %142
%143 = add nuw nsw i64 %141, 1
%144 = icmp eq i64 %143, %30
br i1 %144, label %138, label %140
145: ; preds = %138, %194
%146 = phi i64 [ 0, %138 ], [ %147, %194 ]
%147 = add nuw nsw i64 %146, 1
br i1 %61, label %148, label %194
148: ; preds = %145
%149 = shl nuw nsw i64 %146, 1
%150 = mul nsw i64 %149, %139
%151 = or i64 %149, 1
%152 = mul nsw i64 %151, %139
%153 = trunc i64 %147 to i32
%154 = shl i32 %153, 1
br label %159
155: ; preds = %194
%156 = icmp sgt i32 %38, 0
br i1 %156, label %157, label %196
157: ; preds = %155
%158 = shl nuw nsw i64 %39, 3
br label %196
159: ; preds = %148, %159
%160 = phi i64 [ 0, %148 ], [ %192, %159 ]
%161 = load i32, i32* @M_PI
%162 = call i32 @sqrt(i32 %161)
%163 = add i32 %162, 1
%164 = icmp ult i32 %163, 3
%165 = select i1 %164, i32 %162, i32 0
%166 = load i32, i32* @M_PI
%167 = mul i32 %154, %166
%168 = sitofp i32 %167 to double
%169 = getelementptr inbounds double, double* %32, i64 %160
%170 = load double, double* %169
%171 = fmul double %170, %168
%172 = call i32 @sin(double %171)
%173 = mul nsw i32 %165, %172
%174 = sitofp i32 %173 to double
%175 = add nsw i64 %160, %150
%176 = getelementptr inbounds double, double* %37, i64 %175
store double %174, double* %176
%177 = load i32, i32* @M_PI
%178 = call i32 @sqrt(i32 %177)
%179 = add i32 %178, 1
%180 = icmp ult i32 %179, 3
%181 = select i1 %180, i32 %178, i32 0
%182 = load i32, i32* @M_PI
%183 = mul i32 %154, %182
%184 = sitofp i32 %183 to double
%185 = load double, double* %169
%186 = fmul double %185, %184
%187 = call i32 @cos(double %186)
%188 = mul nsw i32 %181, %187
%189 = sitofp i32 %188 to double
%190 = add nsw i64 %160, %152
%191 = getelementptr inbounds double, double* %37, i64 %190
store double %189, double* %191
%192 = add nuw nsw i64 %160, 1
%193 = icmp eq i64 %192, %30
br i1 %193, label %194, label %159
194: ; preds = %159, %145
%195 = icmp eq i64 %147, 10
br i1 %195, label %155, label %145
196: ; preds = %157, %155
%197 = getelementptr inbounds [4 x double], [4 x double]* %18, i64 0, i64 0
%198 = getelementptr inbounds [4 x double], [4 x double]* %18, i64 0, i64 1
%199 = getelementptr inbounds [4 x double], [4 x double]* %18, i64 0, i64 2
%200 = getelementptr inbounds [4 x double], [4 x double]* %18, i64 0, i64 3
%201 = getelementptr inbounds [4 x double], [4 x double]* %19, i64 0, i64 0
%202 = getelementptr inbounds [4 x double], [4 x double]* %19, i64 0, i64 1
%203 = getelementptr inbounds [4 x double], [4 x double]* %19, i64 0, i64 2
%204 = getelementptr inbounds [4 x double], [4 x double]* %19, i64 0, i64 3
%205 = icmp sgt i32 %23, 0
%206 = sext i32 %58 to i64
%207 = getelementptr inbounds double, double* %54, i64 %206
%208 = icmp eq i32 %28, 1
%209 = shl nuw nsw i64 %30, 3
%210 = sext i32 %23 to i64
%211 = zext i32 %23 to i64
%212 = shl nsw i64 %139, 1
%213 = mul nsw i64 %139, 3
%214 = shl nsw i64 %139, 2
%215 = mul nsw i64 %139, 5
%216 = mul nsw i64 %139, 6
%217 = mul nsw i64 %139, 7
%218 = shl nsw i64 %139, 3
%219 = mul nsw i64 %139, 9
%220 = mul nsw i64 %139, 10
%221 = mul nsw i64 %139, 11
%222 = mul nsw i64 %139, 12
%223 = mul nsw i64 %139, 13
%224 = mul nsw i64 %139, 14
%225 = mul nsw i64 %139, 15
%226 = shl nsw i64 %139, 4
%227 = mul nsw i64 %139, 17
%228 = mul nsw i64 %139, 18
%229 = mul nsw i64 %139, 19
%230 = and i64 %30, 4294967292
%231 = add nsw i64 %230, -4
%232 = lshr exact i64 %231, 2
%233 = add nuw nsw i64 %232, 1
%234 = and i64 %30, 4294967294
%235 = add nsw i64 %234, -2
%236 = lshr exact i64 %235, 1
%237 = add nuw i64 %236, 1
%238 = icmp ult i32 %22, 4
%239 = and i64 %30, 4294967292
%240 = and i64 %233, 1
%241 = icmp eq i64 %231, 0
%242 = and i64 %233, 9223372036854775806
%243 = icmp eq i64 %240, 0
%244 = icmp eq i64 %239, %30
%245 = icmp eq i32 %22, 1
%246 = and i64 %30, 4294967294
%247 = icmp eq i64 %246, %30
%248 = icmp ult i32 %22, 4
%249 = and i64 %30, 4294967292
%250 = and i64 %233, 1
%251 = icmp eq i64 %231, 0
%252 = and i64 %233, 9223372036854775806
%253 = icmp eq i64 %250, 0
%254 = icmp eq i64 %249, %30
%255 = icmp eq i32 %22, 1
%256 = and i64 %30, 4294967294
%257 = and i64 %237, 1
%258 = icmp eq i64 %235, 0
%259 = and i64 %237, -2
%260 = icmp eq i64 %257, 0
%261 = icmp eq i64 %256, %30
%262 = icmp ult i32 %22, 4
%263 = and i64 %30, 4294967292
%264 = and i64 %233, 1
%265 = icmp eq i64 %231, 0
%266 = and i64 %233, 9223372036854775806
%267 = icmp eq i64 %264, 0
%268 = icmp eq i64 %263, %30
%269 = icmp ult i32 %22, 4
%270 = and i64 %30, 4294967292
%271 = and i64 %233, 1
%272 = icmp eq i64 %231, 0
%273 = and i64 %233, 9223372036854775806
%274 = icmp eq i64 %271, 0
%275 = icmp eq i64 %270, %30
%276 = and i64 %30, 1
%277 = icmp eq i32 %22, 1
%278 = and i64 %30, 4294967294
%279 = icmp eq i64 %276, 0
br label %280
280: ; preds = %196, %952
%281 = phi double* [ %52, %952 ], [ %33, %196 ]
%282 = phi double [ %953, %952 ], [ 1.000000e+00, %196 ]
%283 = phi i32 [ %949, %952 ], [ 0, %196 ]
%284 = call double @innerprod_q2(i32* %1, double* %40, double* %3)
%285 = call double @acos(double %282)
%286 = call i32 @sin(double %285)
%287 = sub nsw i32 0, %286
%288 = sitofp i32 %287 to double
store double %288, double* %197
%289 = call i32 @cos(double %285)
%290 = sitofp i32 %289 to double
store double %290, double* %198
%291 = call i32 @cos(double %285)
%292 = sub nsw i32 0, %291
%293 = sitofp i32 %292 to double
store double %293, double* %199
%294 = call i32 @sin(double %285)
%295 = sub nsw i32 0, %294
%296 = sitofp i32 %295 to double
store double %296, double* %200
%297 = call i32 @product(i32 %23, i32 %23, i32 %22, double* %197, double* %40, double* %43)
%298 = call double @innerprod_q2(i32* %1, double* %43, double* %3)
%299 = load i32, i32* %5
%300 = sub nsw i32 0, %299
%301 = sitofp i32 %300 to double
%302 = load double, double* %2
%303 = fadd double %284, %302
%304 = fmul double %303, %301
%305 = fptosi double %304 to i32
%306 = call i32 @exp(i32 %305)
%307 = sitofp i32 %306 to double
%308 = load i32, i32* %5
%309 = sitofp i32 %308 to double
%310 = fmul double %307, %309
%311 = fadd double %307, 1.000000e+00
%312 = fdiv double %310, %311
%313 = fmul double %298, %312
%314 = call double @llvm.fmuladd.f64(double %26, double %313, double %285)
%315 = call i32 @cos(double %314)
%316 = sitofp i32 %315 to double
store double %316, double* %201
%317 = call i32 @sin(double %314)
%318 = sitofp i32 %317 to double
store double %318, double* %202
%319 = call i32 @sin(double %314)
%320 = sub nsw i32 0, %319
%321 = sitofp i32 %320 to double
store double %321, double* %203
%322 = call i32 @cos(double %314)
%323 = sitofp i32 %322 to double
store double %323, double* %204
%324 = call i32 @col_gradient(i32 %23, i32 %22, double* %40, double %60, double* %44)
br label %434
325: ; preds = %540
br i1 %61, label %326, label %695
326: ; preds = %325
br i1 %245, label %432, label %327
327: ; preds = %326, %327
%328 = phi i64 [ %429, %327 ], [ 0, %326 ]
%329 = getelementptr inbounds double, double* %49, i64 %328
%330 = bitcast double* %329 to <2 x double>*
%331 = load <2 x double>, <2 x double>* %330
%332 = getelementptr inbounds double, double* %48, i64 %328
%333 = add nsw i64 %328, %139
%334 = getelementptr inbounds double, double* %49, i64 %333
%335 = bitcast double* %334 to <2 x double>*
%336 = load <2 x double>, <2 x double>* %335
%337 = fadd <2 x double> %331, %336
%338 = add nsw i64 %212, %328
%339 = getelementptr inbounds double, double* %49, i64 %338
%340 = bitcast double* %339 to <2 x double>*
%341 = load <2 x double>, <2 x double>* %340
%342 = fadd <2 x double> %337, %341
%343 = add nsw i64 %213, %328
%344 = getelementptr inbounds double, double* %49, i64 %343
%345 = bitcast double* %344 to <2 x double>*
%346 = load <2 x double>, <2 x double>* %345
%347 = fadd <2 x double> %342, %346
%348 = add nsw i64 %214, %328
%349 = getelementptr inbounds double, double* %49, i64 %348
%350 = bitcast double* %349 to <2 x double>*
%351 = load <2 x double>, <2 x double>* %350
%352 = fadd <2 x double> %347, %351
%353 = add nsw i64 %215, %328
%354 = getelementptr inbounds double, double* %49, i64 %353
%355 = bitcast double* %354 to <2 x double>*
%356 = load <2 x double>, <2 x double>* %355
%357 = fadd <2 x double> %352, %356
%358 = add nsw i64 %216, %328
%359 = getelementptr inbounds double, double* %49, i64 %358
%360 = bitcast double* %359 to <2 x double>*
%361 = load <2 x double>, <2 x double>* %360
%362 = fadd <2 x double> %357, %361
%363 = add nsw i64 %217, %328
%364 = getelementptr inbounds double, double* %49, i64 %363
%365 = bitcast double* %364 to <2 x double>*
%366 = load <2 x double>, <2 x double>* %365
%367 = fadd <2 x double> %362, %366
%368 = add nsw i64 %218, %328
%369 = getelementptr inbounds double, double* %49, i64 %368
%370 = bitcast double* %369 to <2 x double>*
%371 = load <2 x double>, <2 x double>* %370
%372 = fadd <2 x double> %367, %371
%373 = add nsw i64 %219, %328
%374 = getelementptr inbounds double, double* %49, i64 %373
%375 = bitcast double* %374 to <2 x double>*
%376 = load <2 x double>, <2 x double>* %375
%377 = fadd <2 x double> %372, %376
%378 = add nsw i64 %220, %328
%379 = getelementptr inbounds double, double* %49, i64 %378
%380 = bitcast double* %379 to <2 x double>*
%381 = load <2 x double>, <2 x double>* %380
%382 = fadd <2 x double> %377, %381
%383 = add nsw i64 %221, %328
%384 = getelementptr inbounds double, double* %49, i64 %383
%385 = bitcast double* %384 to <2 x double>*
%386 = load <2 x double>, <2 x double>* %385
%387 = fadd <2 x double> %382, %386
%388 = add nsw i64 %222, %328
%389 = getelementptr inbounds double, double* %49, i64 %388
%390 = bitcast double* %389 to <2 x double>*
%391 = load <2 x double>, <2 x double>* %390
%392 = fadd <2 x double> %387, %391
%393 = add nsw i64 %223, %328
%394 = getelementptr inbounds double, double* %49, i64 %393
%395 = bitcast double* %394 to <2 x double>*
%396 = load <2 x double>, <2 x double>* %395
%397 = fadd <2 x double> %392, %396
%398 = add nsw i64 %224, %328
%399 = getelementptr inbounds double, double* %49, i64 %398
%400 = bitcast double* %399 to <2 x double>*
%401 = load <2 x double>, <2 x double>* %400
%402 = fadd <2 x double> %397, %401
%403 = add nsw i64 %225, %328
%404 = getelementptr inbounds double, double* %49, i64 %403
%405 = bitcast double* %404 to <2 x double>*
%406 = load <2 x double>, <2 x double>* %405
%407 = fadd <2 x double> %402, %406
%408 = add nsw i64 %226, %328
%409 = getelementptr inbounds double, double* %49, i64 %408
%410 = bitcast double* %409 to <2 x double>*
%411 = load <2 x double>, <2 x double>* %410
%412 = fadd <2 x double> %407, %411
%413 = add nsw i64 %227, %328
%414 = getelementptr inbounds double, double* %49, i64 %413
%415 = bitcast double* %414 to <2 x double>*
%416 = load <2 x double>, <2 x double>* %415
%417 = fadd <2 x double> %412, %416
%418 = add nsw i64 %228, %328
%419 = getelementptr inbounds double, double* %49, i64 %418
%420 = bitcast double* %419 to <2 x double>*
%421 = load <2 x double>, <2 x double>* %420
%422 = fadd <2 x double> %417, %421
%423 = add nsw i64 %229, %328
%424 = getelementptr inbounds double, double* %49, i64 %423
%425 = bitcast double* %424 to <2 x double>*
%426 = load <2 x double>, <2 x double>* %425
%427 = fadd <2 x double> %422, %426
%428 = bitcast double* %332 to <2 x double>*
store <2 x double> %427, <2 x double>* %428
%429 = add nuw i64 %328, 2
%430 = icmp eq i64 %429, %246
br i1 %430, label %431, label %327
431: ; preds = %327
br i1 %247, label %548, label %432
432: ; preds = %326, %431
%433 = phi i64 [ 0, %326 ], [ %246, %431 ]
br label %604
434: ; preds = %543, %280
%435 = phi i64 [ 0, %280 ], [ %544, %543 ]
br i1 %61, label %436, label %442
436: ; preds = %434
%437 = trunc i64 %435 to i32
%438 = mul i32 %22, %437
%439 = sext i32 %438 to i64
%440 = getelementptr double, double* %37, i64 %439
%441 = bitcast double* %440 to i8*
br label %442
442: ; preds = %436, %434
%443 = call i32 @cumtrapz(i32* %1, double* %32, double* %46, double* %45)
br i1 %205, label %444, label %470
444: ; preds = %442
%445 = mul nsw i64 %435, %139
br label %446
446: ; preds = %444, %467
%447 = phi i64 [ 0, %444 ], [ %468, %467 ]
br i1 %61, label %448, label %467
448: ; preds = %446, %448
%449 = phi i64 [ %465, %448 ], [ 0, %446 ]
%450 = mul nsw i64 %449, %210
%451 = add nsw i64 %450, %447
%452 = getelementptr inbounds double, double* %44, i64 %451
%453 = load double, double* %452
%454 = fmul double %453, 2.000000e+00
%455 = getelementptr inbounds double, double* %45, i64 %449
%456 = load double, double* %455
%457 = getelementptr inbounds double, double* %40, i64 %451
%458 = load double, double* %457
%459 = add nsw i64 %449, %445
%460 = getelementptr inbounds double, double* %37, i64 %459
%461 = load double, double* %460
%462 = fmul double %458, %461
%463 = call double @llvm.fmuladd.f64(double %454, double %456, double %462)
%464 = getelementptr inbounds double, double* %43, i64 %451
store double %463, double* %464
%465 = add nuw nsw i64 %449, 1
%466 = icmp eq i64 %465, %30
br i1 %466, label %467, label %448
467: ; preds = %448, %446
%468 = add nuw nsw i64 %447, 1
%469 = icmp eq i64 %468, %211
br i1 %469, label %470, label %446
470: ; preds = %467, %442
%471 = call double @innerprod_q2(i32* %1, double* %43, double* %3)
store double %471, double* %21
br i1 %61, label %472, label %545
472: ; preds = %470
%473 = mul nsw i64 %435, %139
br i1 %238, label %529, label %474
474: ; preds = %472
%475 = insertelement <2 x double> poison, double %471, i64 0
%476 = shufflevector <2 x double> %475, <2 x double> poison, <2 x i32> zeroinitializer
%477 = insertelement <2 x double> poison, double %471, i64 0
%478 = shufflevector <2 x double> %477, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %241, label %512, label %479
479: ; preds = %474, %479
%480 = phi i64 [ %509, %479 ], [ 0, %474 ]
%481 = phi i64 [ %510, %479 ], [ 0, %474 ]
%482 = add nsw i64 %480, %473
%483 = getelementptr inbounds double, double* %37, i64 %482
%484 = bitcast double* %483 to <2 x double>*
%485 = load <2 x double>, <2 x double>* %484
%486 = getelementptr inbounds double, double* %483, i64 2
%487 = bitcast double* %486 to <2 x double>*
%488 = load <2 x double>, <2 x double>* %487
%489 = fmul <2 x double> %476, %485
%490 = fmul <2 x double> %478, %488
%491 = getelementptr inbounds double, double* %49, i64 %482
%492 = bitcast double* %491 to <2 x double>*
store <2 x double> %489, <2 x double>* %492
%493 = getelementptr inbounds double, double* %491, i64 2
%494 = bitcast double* %493 to <2 x double>*
store <2 x double> %490, <2 x double>* %494
%495 = or i64 %480, 4
%496 = add nsw i64 %495, %473
%497 = getelementptr inbounds double, double* %37, i64 %496
%498 = bitcast double* %497 to <2 x double>*
%499 = load <2 x double>, <2 x double>* %498
%500 = getelementptr inbounds double, double* %497, i64 2
%501 = bitcast double* %500 to <2 x double>*
%502 = load <2 x double>, <2 x double>* %501
%503 = fmul <2 x double> %476, %499
%504 = fmul <2 x double> %478, %502
%505 = getelementptr inbounds double, double* %49, i64 %496
%506 = bitcast double* %505 to <2 x double>*
store <2 x double> %503, <2 x double>* %506
%507 = getelementptr inbounds double, double* %505, i64 2
%508 = bitcast double* %507 to <2 x double>*
store <2 x double> %504, <2 x double>* %508
%509 = add nuw i64 %480, 8
%510 = add i64 %481, 2
%511 = icmp eq i64 %510, %242
br i1 %511, label %512, label %479
512: ; preds = %479, %474
%513 = phi i64 [ 0, %474 ], [ %509, %479 ]
br i1 %243, label %528, label %514
514: ; preds = %512
%515 = add nsw i64 %513, %473
%516 = getelementptr inbounds double, double* %37, i64 %515
%517 = bitcast double* %516 to <2 x double>*
%518 = load <2 x double>, <2 x double>* %517
%519 = getelementptr inbounds double, double* %516, i64 2
%520 = bitcast double* %519 to <2 x double>*
%521 = load <2 x double>, <2 x double>* %520
%522 = fmul <2 x double> %476, %518
%523 = fmul <2 x double> %478, %521
%524 = getelementptr inbounds double, double* %49, i64 %515
%525 = bitcast double* %524 to <2 x double>*
store <2 x double> %522, <2 x double>* %525
%526 = getelementptr inbounds double, double* %524, i64 2
%527 = bitcast double* %526 to <2 x double>*
store <2 x double> %523, <2 x double>* %527
br label %528
528: ; preds = %512, %514
br i1 %244, label %540, label %529
529: ; preds = %472, %528
%530 = phi i64 [ 0, %472 ], [ %239, %528 ]
br label %531
531: ; preds = %529, %531
%532 = phi i64 [ %538, %531 ], [ %530, %529 ]
%533 = add nsw i64 %532, %473
%534 = getelementptr inbounds double, double* %37, i64 %533
%535 = load double, double* %534
%536 = fmul double %471, %535
%537 = getelementptr inbounds double, double* %49, i64 %533
store double %536, double* %537
%538 = add nuw nsw i64 %532, 1
%539 = icmp eq i64 %538, %30
br i1 %539, label %540, label %531
540: ; preds = %531, %528
%541 = add nuw nsw i64 %435, 1
%542 = icmp eq i64 %541, 20
br i1 %542, label %325, label %543
543: ; preds = %540, %545
%544 = phi i64 [ %541, %540 ], [ %546, %545 ]
br label %434
545: ; preds = %470
%546 = add nuw nsw i64 %435, 1
%547 = icmp eq i64 %546, 20
br i1 %547, label %695, label %543
548: ; preds = %604, %431
br i1 %61, label %549, label %695
549: ; preds = %548
br i1 %248, label %602, label %550
550: ; preds = %549
%551 = insertelement <2 x double> poison, double %312, i64 0
%552 = shufflevector <2 x double> %551, <2 x double> poison, <2 x i32> zeroinitializer
%553 = insertelement <2 x double> poison, double %312, i64 0
%554 = shufflevector <2 x double> %553, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %251, label %586, label %555
555: ; preds = %550, %555
%556 = phi i64 [ %583, %555 ], [ 0, %550 ]
%557 = phi i64 [ %584, %555 ], [ 0, %550 ]
%558 = getelementptr inbounds double, double* %48, i64 %556
%559 = bitcast double* %558 to <2 x double>*
%560 = load <2 x double>, <2 x double>* %559
%561 = getelementptr inbounds double, double* %558, i64 2
%562 = bitcast double* %561 to <2 x double>*
%563 = load <2 x double>, <2 x double>* %562
%564 = fmul <2 x double> %552, %560
%565 = fmul <2 x double> %554, %563
%566 = getelementptr inbounds double, double* %50, i64 %556
%567 = bitcast double* %566 to <2 x double>*
store <2 x double> %564, <2 x double>* %567
%568 = getelementptr inbounds double, double* %566, i64 2
%569 = bitcast double* %568 to <2 x double>*
store <2 x double> %565, <2 x double>* %569
%570 = or i64 %556, 4
%571 = getelementptr inbounds double, double* %48, i64 %570
%572 = bitcast double* %571 to <2 x double>*
%573 = load <2 x double>, <2 x double>* %572
%574 = getelementptr inbounds double, double* %571, i64 2
%575 = bitcast double* %574 to <2 x double>*
%576 = load <2 x double>, <2 x double>* %575
%577 = fmul <2 x double> %552, %573
%578 = fmul <2 x double> %554, %576
%579 = getelementptr inbounds double, double* %50, i64 %570
%580 = bitcast double* %579 to <2 x double>*
store <2 x double> %577, <2 x double>* %580
%581 = getelementptr inbounds double, double* %579, i64 2
%582 = bitcast double* %581 to <2 x double>*
store <2 x double> %578, <2 x double>* %582
%583 = add nuw i64 %556, 8
%584 = add i64 %557, 2
%585 = icmp eq i64 %584, %252
br i1 %585, label %586, label %555
586: ; preds = %555, %550
%587 = phi i64 [ 0, %550 ], [ %583, %555 ]
br i1 %253, label %601, label %588
588: ; preds = %586
%589 = getelementptr inbounds double, double* %48, i64 %587
%590 = bitcast double* %589 to <2 x double>*
%591 = load <2 x double>, <2 x double>* %590
%592 = getelementptr inbounds double, double* %589, i64 2
%593 = bitcast double* %592 to <2 x double>*
%594 = load <2 x double>, <2 x double>* %593
%595 = fmul <2 x double> %552, %591
%596 = fmul <2 x double> %554, %594
%597 = getelementptr inbounds double, double* %50, i64 %587
%598 = bitcast double* %597 to <2 x double>*
store <2 x double> %595, <2 x double>* %598
%599 = getelementptr inbounds double, double* %597, i64 2
%600 = bitcast double* %599 to <2 x double>*
store <2 x double> %596, <2 x double>* %600
br label %601
601: ; preds = %586, %588
br i1 %254, label %695, label %602
602: ; preds = %549, %601
%603 = phi i64 [ 0, %549 ], [ %249, %601 ]
br label %687
604: ; preds = %432, %604
%605 = phi i64 [ %685, %604 ], [ %433, %432 ]
%606 = getelementptr inbounds double, double* %49, i64 %605
%607 = load double, double* %606
%608 = getelementptr inbounds double, double* %48, i64 %605
%609 = add nsw i64 %605, %139
%610 = getelementptr inbounds double, double* %49, i64 %609
%611 = load double, double* %610
%612 = fadd double %607, %611
%613 = add nsw i64 %212, %605
%614 = getelementptr inbounds double, double* %49, i64 %613
%615 = load double, double* %614
%616 = fadd double %612, %615
%617 = add nsw i64 %213, %605
%618 = getelementptr inbounds double, double* %49, i64 %617
%619 = load double, double* %618
%620 = fadd double %616, %619
%621 = add nsw i64 %214, %605
%622 = getelementptr inbounds double, double* %49, i64 %621
%623 = load double, double* %622
%624 = fadd double %620, %623
%625 = add nsw i64 %215, %605
%626 = getelementptr inbounds double, double* %49, i64 %625
%627 = load double, double* %626
%628 = fadd double %624, %627
%629 = add nsw i64 %216, %605
%630 = getelementptr inbounds double, double* %49, i64 %629
%631 = load double, double* %630
%632 = fadd double %628, %631
%633 = add nsw i64 %217, %605
%634 = getelementptr inbounds double, double* %49, i64 %633
%635 = load double, double* %634
%636 = fadd double %632, %635
%637 = add nsw i64 %218, %605
%638 = getelementptr inbounds double, double* %49, i64 %637
%639 = load double, double* %638
%640 = fadd double %636, %639
%641 = add nsw i64 %219, %605
%642 = getelementptr inbounds double, double* %49, i64 %641
%643 = load double, double* %642
%644 = fadd double %640, %643
%645 = add nsw i64 %220, %605
%646 = getelementptr inbounds double, double* %49, i64 %645
%647 = load double, double* %646
%648 = fadd double %644, %647
%649 = add nsw i64 %221, %605
%650 = getelementptr inbounds double, double* %49, i64 %649
%651 = load double, double* %650
%652 = fadd double %648, %651
%653 = add nsw i64 %222, %605
%654 = getelementptr inbounds double, double* %49, i64 %653
%655 = load double, double* %654
%656 = fadd double %652, %655
%657 = add nsw i64 %223, %605
%658 = getelementptr inbounds double, double* %49, i64 %657
%659 = load double, double* %658
%660 = fadd double %656, %659
%661 = add nsw i64 %224, %605
%662 = getelementptr inbounds double, double* %49, i64 %661
%663 = load double, double* %662
%664 = fadd double %660, %663
%665 = add nsw i64 %225, %605
%666 = getelementptr inbounds double, double* %49, i64 %665
%667 = load double, double* %666
%668 = fadd double %664, %667
%669 = add nsw i64 %226, %605
%670 = getelementptr inbounds double, double* %49, i64 %669
%671 = load double, double* %670
%672 = fadd double %668, %671
%673 = add nsw i64 %227, %605
%674 = getelementptr inbounds double, double* %49, i64 %673
%675 = load double, double* %674
%676 = fadd double %672, %675
%677 = add nsw i64 %228, %605
%678 = getelementptr inbounds double, double* %49, i64 %677
%679 = load double, double* %678
%680 = fadd double %676, %679
%681 = add nsw i64 %229, %605
%682 = getelementptr inbounds double, double* %49, i64 %681
%683 = load double, double* %682
%684 = fadd double %680, %683
store double %684, double* %608
%685 = add nuw nsw i64 %605, 1
%686 = icmp eq i64 %685, %30
br i1 %686, label %548, label %604
687: ; preds = %602, %687
%688 = phi i64 [ %693, %687 ], [ %603, %602 ]
%689 = getelementptr inbounds double, double* %48, i64 %688
%690 = load double, double* %689
%691 = fmul double %312, %690
%692 = getelementptr inbounds double, double* %50, i64 %688
store double %691, double* %692
%693 = add nuw nsw i64 %688, 1
%694 = icmp eq i64 %693, %30
br i1 %694, label %695, label %687
695: ; preds = %545, %687, %601, %325, %548
%696 = call i32 @pvecnorm(i32* %1, double* %50, double* %17, double* %21)
%697 = load double, double* %21
%698 = fmul double %27, %697
%699 = call i32 @cos(double %698)
%700 = sitofp i32 %699 to double
%701 = load double, double* %21
%702 = fmul double %27, %701
%703 = call i32 @sin(double %702)
%704 = sitofp i32 %703 to double
%705 = load double, double* %21
%706 = fcmp oeq double %705, 0.000000e+00
br i1 %706, label %762, label %707
707: ; preds = %695
br i1 %61, label %708, label %898
708: ; preds = %707
br i1 %255, label %760, label %709
709: ; preds = %708
%710 = insertelement <2 x double> poison, double %705, i64 0
%711 = shufflevector <2 x double> %710, <2 x double> poison, <2 x i32> zeroinitializer
%712 = insertelement <2 x double> poison, double %704, i64 0
%713 = shufflevector <2 x double> %712, <2 x double> poison, <2 x i32> zeroinitializer
%714 = insertelement <2 x double> poison, double %700, i64 0
%715 = shufflevector <2 x double> %714, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %258, label %745, label %716
716: ; preds = %709, %716
%717 = phi i64 [ %742, %716 ], [ 0, %709 ]
%718 = phi i64 [ %743, %716 ], [ 0, %709 ]
%719 = getelementptr inbounds double, double* %34, i64 %717
%720 = bitcast double* %719 to <2 x double>*
%721 = load <2 x double>, <2 x double>* %720
%722 = getelementptr inbounds double, double* %50, i64 %717
%723 = bitcast double* %722 to <2 x double>*
%724 = load <2 x double>, <2 x double>* %723
%725 = fdiv <2 x double> %724, %711
%726 = fmul <2 x double> %725, %713
%727 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %715, <2 x double> %721, <2 x double> %726)
%728 = getelementptr inbounds double, double* %51, i64 %717
%729 = bitcast double* %728 to <2 x double>*
store <2 x double> %727, <2 x double>* %729
%730 = or i64 %717, 2
%731 = getelementptr inbounds double, double* %34, i64 %730
%732 = bitcast double* %731 to <2 x double>*
%733 = load <2 x double>, <2 x double>* %732
%734 = getelementptr inbounds double, double* %50, i64 %730
%735 = bitcast double* %734 to <2 x double>*
%736 = load <2 x double>, <2 x double>* %735
%737 = fdiv <2 x double> %736, %711
%738 = fmul <2 x double> %737, %713
%739 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %715, <2 x double> %733, <2 x double> %738)
%740 = getelementptr inbounds double, double* %51, i64 %730
%741 = bitcast double* %740 to <2 x double>*
store <2 x double> %739, <2 x double>* %741
%742 = add nuw i64 %717, 4
%743 = add nuw i64 %718, 2
%744 = icmp eq i64 %743, %259
br i1 %744, label %745, label %716
745: ; preds = %716, %709
%746 = phi i64 [ 0, %709 ], [ %742, %716 ]
br i1 %260, label %759, label %747
747: ; preds = %745
%748 = getelementptr inbounds double, double* %34, i64 %746
%749 = bitcast double* %748 to <2 x double>*
%750 = load <2 x double>, <2 x double>* %749
%751 = getelementptr inbounds double, double* %50, i64 %746
%752 = bitcast double* %751 to <2 x double>*
%753 = load <2 x double>, <2 x double>* %752
%754 = fdiv <2 x double> %753, %711
%755 = fmul <2 x double> %754, %713
%756 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %715, <2 x double> %750, <2 x double> %755)
%757 = getelementptr inbounds double, double* %51, i64 %746
%758 = bitcast double* %757 to <2 x double>*
store <2 x double> %756, <2 x double>* %758
br label %759
759: ; preds = %745, %747
br i1 %261, label %838, label %760
760: ; preds = %708, %759
%761 = phi i64 [ 0, %708 ], [ %256, %759 ]
br label %826
762: ; preds = %695
br i1 %61, label %763, label %898
763: ; preds = %762
br i1 %262, label %816, label %764
764: ; preds = %763
%765 = insertelement <2 x double> poison, double %700, i64 0
%766 = shufflevector <2 x double> %765, <2 x double> poison, <2 x i32> zeroinitializer
%767 = insertelement <2 x double> poison, double %700, i64 0
%768 = shufflevector <2 x double> %767, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %265, label %800, label %769
769: ; preds = %764, %769
%770 = phi i64 [ %797, %769 ], [ 0, %764 ]
%771 = phi i64 [ %798, %769 ], [ 0, %764 ]
%772 = getelementptr inbounds double, double* %34, i64 %770
%773 = bitcast double* %772 to <2 x double>*
%774 = load <2 x double>, <2 x double>* %773
%775 = getelementptr inbounds double, double* %772, i64 2
%776 = bitcast double* %775 to <2 x double>*
%777 = load <2 x double>, <2 x double>* %776
%778 = fmul <2 x double> %774, %766
%779 = fmul <2 x double> %777, %768
%780 = getelementptr inbounds double, double* %51, i64 %770
%781 = bitcast double* %780 to <2 x double>*
store <2 x double> %778, <2 x double>* %781
%782 = getelementptr inbounds double, double* %780, i64 2
%783 = bitcast double* %782 to <2 x double>*
store <2 x double> %779, <2 x double>* %783
%784 = or i64 %770, 4
%785 = getelementptr inbounds double, double* %34, i64 %784
%786 = bitcast double* %785 to <2 x double>*
%787 = load <2 x double>, <2 x double>* %786
%788 = getelementptr inbounds double, double* %785, i64 2
%789 = bitcast double* %788 to <2 x double>*
%790 = load <2 x double>, <2 x double>* %789
%791 = fmul <2 x double> %787, %766
%792 = fmul <2 x double> %790, %768
%793 = getelementptr inbounds double, double* %51, i64 %784
%794 = bitcast double* %793 to <2 x double>*
store <2 x double> %791, <2 x double>* %794
%795 = getelementptr inbounds double, double* %793, i64 2
%796 = bitcast double* %795 to <2 x double>*
store <2 x double> %792, <2 x double>* %796
%797 = add nuw i64 %770, 8
%798 = add i64 %771, 2
%799 = icmp eq i64 %798, %266
br i1 %799, label %800, label %769
800: ; preds = %769, %764
%801 = phi i64 [ 0, %764 ], [ %797, %769 ]
br i1 %267, label %815, label %802
802: ; preds = %800
%803 = getelementptr inbounds double, double* %34, i64 %801
%804 = bitcast double* %803 to <2 x double>*
%805 = load <2 x double>, <2 x double>* %804
%806 = getelementptr inbounds double, double* %803, i64 2
%807 = bitcast double* %806 to <2 x double>*
%808 = load <2 x double>, <2 x double>* %807
%809 = fmul <2 x double> %805, %766
%810 = fmul <2 x double> %808, %768
%811 = getelementptr inbounds double, double* %51, i64 %801
%812 = bitcast double* %811 to <2 x double>*
store <2 x double> %809, <2 x double>* %812
%813 = getelementptr inbounds double, double* %811, i64 2
%814 = bitcast double* %813 to <2 x double>*
store <2 x double> %810, <2 x double>* %814
br label %815
815: ; preds = %800, %802
br i1 %268, label %838, label %816
816: ; preds = %763, %815
%817 = phi i64 [ 0, %763 ], [ %263, %815 ]
br label %818
818: ; preds = %816, %818
%819 = phi i64 [ %824, %818 ], [ %817, %816 ]
%820 = getelementptr inbounds double, double* %34, i64 %819
%821 = load double, double* %820
%822 = fmul double %821, %700
%823 = getelementptr inbounds double, double* %51, i64 %819
store double %822, double* %823
%824 = add nuw nsw i64 %819, 1
%825 = icmp eq i64 %824, %30
br i1 %825, label %838, label %818
826: ; preds = %760, %826
%827 = phi i64 [ %836, %826 ], [ %761, %760 ]
%828 = getelementptr inbounds double, double* %34, i64 %827
%829 = load double, double* %828
%830 = getelementptr inbounds double, double* %50, i64 %827
%831 = load double, double* %830
%832 = fdiv double %831, %705
%833 = fmul double %832, %704
%834 = call double @llvm.fmuladd.f64(double %700, double %829, double %833)
%835 = getelementptr inbounds double, double* %51, i64 %827
store double %834, double* %835
%836 = add nuw nsw i64 %827, 1
%837 = icmp eq i64 %836, %30
br i1 %837, label %838, label %826
838: ; preds = %826, %818, %759, %815
br i1 %61, label %839, label %898
839: ; preds = %838
br i1 %269, label %888, label %840
840: ; preds = %839
br i1 %272, label %872, label %841
841: ; preds = %840, %841
%842 = phi i64 [ %869, %841 ], [ 0, %840 ]
%843 = phi i64 [ %870, %841 ], [ 0, %840 ]
%844 = getelementptr inbounds double, double* %51, i64 %842
%845 = bitcast double* %844 to <2 x double>*
%846 = load <2 x double>, <2 x double>* %845
%847 = getelementptr inbounds double, double* %844, i64 2
%848 = bitcast double* %847 to <2 x double>*
%849 = load <2 x double>, <2 x double>* %848
%850 = fmul <2 x double> %846, %846
%851 = fmul <2 x double> %849, %849
%852 = getelementptr inbounds double, double* %55, i64 %842
%853 = bitcast double* %852 to <2 x double>*
store <2 x double> %850, <2 x double>* %853
%854 = getelementptr inbounds double, double* %852, i64 2
%855 = bitcast double* %854 to <2 x double>*
store <2 x double> %851, <2 x double>* %855
%856 = or i64 %842, 4
%857 = getelementptr inbounds double, double* %51, i64 %856
%858 = bitcast double* %857 to <2 x double>*
%859 = load <2 x double>, <2 x double>* %858
%860 = getelementptr inbounds double, double* %857, i64 2
%861 = bitcast double* %860 to <2 x double>*
%862 = load <2 x double>, <2 x double>* %861
%863 = fmul <2 x double> %859, %859
%864 = fmul <2 x double> %862, %862
%865 = getelementptr inbounds double, double* %55, i64 %856
%866 = bitcast double* %865 to <2 x double>*
store <2 x double> %863, <2 x double>* %866
%867 = getelementptr inbounds double, double* %865, i64 2
%868 = bitcast double* %867 to <2 x double>*
store <2 x double> %864, <2 x double>* %868
%869 = add nuw i64 %842, 8
%870 = add i64 %843, 2
%871 = icmp eq i64 %870, %273
br i1 %871, label %872, label %841
872: ; preds = %841, %840
%873 = phi i64 [ 0, %840 ], [ %869, %841 ]
br i1 %274, label %887, label %874
874: ; preds = %872
%875 = getelementptr inbounds double, double* %51, i64 %873
%876 = bitcast double* %875 to <2 x double>*
%877 = load <2 x double>, <2 x double>* %876
%878 = getelementptr inbounds double, double* %875, i64 2
%879 = bitcast double* %878 to <2 x double>*
%880 = load <2 x double>, <2 x double>* %879
%881 = fmul <2 x double> %877, %877
%882 = fmul <2 x double> %880, %880
%883 = getelementptr inbounds double, double* %55, i64 %873
%884 = bitcast double* %883 to <2 x double>*
store <2 x double> %881, <2 x double>* %884
%885 = getelementptr inbounds double, double* %883, i64 2
%886 = bitcast double* %885 to <2 x double>*
store <2 x double> %882, <2 x double>* %886
br label %887
887: ; preds = %872, %874
br i1 %275, label %900, label %888
888: ; preds = %839, %887
%889 = phi i64 [ 0, %839 ], [ %270, %887 ]
br label %890
890: ; preds = %888, %890
%891 = phi i64 [ %896, %890 ], [ %889, %888 ]
%892 = getelementptr inbounds double, double* %51, i64 %891
%893 = load double, double* %892
%894 = fmul double %893, %893
%895 = getelementptr inbounds double, double* %55, i64 %891
store double %894, double* %895
%896 = add nuw nsw i64 %891, 1
%897 = icmp eq i64 %896, %30
br i1 %897, label %900, label %890
898: ; preds = %838, %707, %762
%899 = call i32 @cumtrapz(i32* %1, double* %32, double* %55, double* %54)
br label %934
900: ; preds = %890, %887
%901 = call i32 @cumtrapz(i32* %1, double* %32, double* %55, double* %54)
br i1 %61, label %902, label %934
902: ; preds = %900
br i1 %277, label %924, label %903
903: ; preds = %902, %903
%904 = phi i64 [ %921, %903 ], [ 0, %902 ]
%905 = phi i64 [ %922, %903 ], [ 0, %902 ]
%906 = getelementptr inbounds double, double* %54, i64 %904
%907 = load double, double* %906
%908 = load double, double* %54
%909 = fsub double %907, %908
%910 = load double, double* %207
%911 = fsub double %910, %908
%912 = fdiv double %909, %911
store double %912, double* %906
%913 = or i64 %904, 1
%914 = getelementptr inbounds double, double* %54, i64 %913
%915 = load double, double* %914
%916 = load double, double* %54
%917 = fsub double %915, %916
%918 = load double, double* %207
%919 = fsub double %918, %916
%920 = fdiv double %917, %919
store double %920, double* %914
%921 = add nuw nsw i64 %904, 2
%922 = add i64 %905, 2
%923 = icmp eq i64 %922, %278
br i1 %923, label %924, label %903
924: ; preds = %903, %902
%925 = phi i64 [ 0, %902 ], [ %921, %903 ]
br i1 %279, label %934, label %926
926: ; preds = %924
%927 = getelementptr inbounds double, double* %54, i64 %925
%928 = load double, double* %927
%929 = load double, double* %54
%930 = fsub double %928, %929
%931 = load double, double* %207
%932 = fsub double %931, %929
%933 = fdiv double %930, %932
store double %933, double* %927
br label %934
934: ; preds = %926, %924, %898, %900
%935 = call i32 @approx(double* %32, double* %281, i32 %22, double* %54, double* %52, i32 %22, i32 1, i32 0, i32 1, i32 0)
%936 = load i32, i32* %5
%937 = sub nsw i32 0, %936
%938 = sitofp i32 %937 to double
%939 = load double, double* %2
%940 = fadd double %284, %939
%941 = fmul double %940, %938
%942 = fptosi double %941 to i32
%943 = call i32 @exp(i32 %942)
%944 = add nsw i32 %943, 1
%945 = add i32 %943, 2
%946 = icmp ult i32 %945, 3
%947 = select i1 %946, i32 %944, i32 0
%948 = call double @log(i32 %947)
%949 = add nuw nsw i32 %283, 1
br i1 %208, label %950, label %952
950: ; preds = %934
%951 = call i32 @printf(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0), i32 %949, double %948)
br label %952
952: ; preds = %934, %950
%953 = load double, double* %201
%954 = call i32 @product(i32 %23, i32 %23, i32 %22, double* %201, double* %4, double* %43)
%955 = call i32 @group_action_by_gamma(i32* %0, i32* %1, double* %43, double* %52, double* %40)
%956 = icmp ult i32 %283, 2
%957 = fcmp uge double %313, %25
%958 = select i1 %956, i1 true, i1 %957
%959 = load double, double* %21
%960 = fcmp uge double %959, %25
%961 = select i1 %958, i1 true, i1 %960
%962 = icmp sgt i32 %24, %283
%963 = select i1 %961, i1 %962, i1 false
br i1 %963, label %280, label %964
964: ; preds = %952
br i1 %61, label %965, label %966
965: ; preds = %964
br label %966
966: ; preds = %965, %964
%967 = mul i32 %23, %23
%968 = icmp eq i32 %967, 0
br i1 %968, label %972, label %969
969: ; preds = %966
%970 = zext i32 %967 to i64
%971 = shl nuw nsw i64 %970, 3
br label %972
972: ; preds = %969, %966
ret void
}