void c_sgbco(double *abd, int lda, int n, int ml, int mu, int *ipvt,
             double *rcond, double *z) {
  int info;
  register int is, j, ju, k, kb, kp1, l, la, lm, lz, m, mm;
  double anorm, ek, s, sm, t, wk, wkm, ynorm;
  anorm = 0.;
  l = ml + 1;
  is = l + mu;
  for (j = 1; j <= n; j++) {
    anorm = ({
      const double _x = (double)(anorm);
      const double _y = (double)(c_sasum(l, &abd[is - 1 + (j - 1) * lda]));
      _x > _y ? _x : _y;
    });
    if (is > ml + 1) {
      is--;
    }
    if (j <= mu) {
      l++;
    }
    if (j >= n - ml) {
      l--;
    }
  }
  c_sgbfa(abd, lda, n, ml, mu, ipvt, &info);
  ek = 1.;
  memset(z, 0, n * sizeof(double));
  m = ml + mu + 1;
  ju = 0;
  for (k = 1; k <= n; k++) {
    if (z[k - 1] != 0.) {
      ek = ((-z[k - 1]) >= 0. ? fabs(ek) : -fabs(ek));
    }
    if (fabs(ek - z[k - 1]) > fabs(abd[m - 1 + (k - 1) * lda])) {
      s = fabs(abd[m - 1 + (k - 1) * lda]) / fabs(ek - z[k - 1]);
      c_sscal(n, s, z);
      ek *= s;
    }
    wk = ek - z[k - 1];
    wkm = -ek - z[k - 1];
    s = fabs(wk);
    sm = fabs(wkm);
    if (abd[m - 1 + (k - 1) * lda] != 0.) {
      wk /= abd[m - 1 + (k - 1) * lda];
      wkm /= abd[m - 1 + (k - 1) * lda];
    } else {
      wk = 1.;
      wkm = 1.;
    }
    kp1 = k + 1;
    ju = ({
      const int _i = (int)(({
        const int _i = (int)(ju);
        const int _j = (int)(mu + ipvt[k - 1]);
        _i > _j ? _i : _j;
      }));
      const int _j = (int)(n);
      _i < _j ? _i : _j;
    });
    mm = m;
    if (kp1 <= ju) {
      for (j = kp1; j <= ju; j++) {
        mm--;
        sm += fabs(z[j - 1] + wkm * abd[mm - 1 + (j - 1) * lda]);
        z[j - 1] += wk * abd[mm - 1 + (j - 1) * lda];
        s += fabs(z[j - 1]);
      }
      if (s < sm) {
        t = wkm - wk;
        wk = wkm;
        mm = m;
        for (j = kp1; j <= ju; j++) {
          mm--;
          z[j - 1] += t * abd[mm - 1 + (j - 1) * lda];
        }
      }
    }
    z[k - 1] = wk;
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  for (kb = 1; kb <= n; kb++) {
    k = n + 1 - kb;
    lm = ({
      const int _i = (int)(ml);
      const int _j = (int)(n - k);
      _i < _j ? _i : _j;
    });
    if (k < n) {
      z[k - 1] += c_sdot(lm, &abd[m + 1 - 1 + (k - 1) * lda], &z[k + 1 - 1]);
    }
    if (fabs(z[k - 1]) > 1.) {
      s = 1. / fabs(z[k - 1]);
      c_sscal(n, s, z);
    }
    l = ipvt[k - 1];
    t = z[l - 1];
    z[l - 1] = z[k - 1];
    z[k - 1] = t;
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  ynorm = 1.;
  for (k = 1; k <= n; k++) {
    l = ipvt[k - 1];
    t = z[l - 1];
    z[l - 1] = z[k - 1];
    z[k - 1] = t;
    lm = ({
      const int _i = (int)(ml);
      const int _j = (int)(n - k);
      _i < _j ? _i : _j;
    });
    if (k < n) {
      c_saxpy(lm, t, &abd[m + 1 - 1 + (k - 1) * lda], &z[k + 1 - 1]);
    }
    if (fabs(z[k - 1]) > 1.) {
      s = 1. / fabs(z[k - 1]);
      c_sscal(n, s, z);
      ynorm *= s;
    }
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  ynorm *= s;
  for (kb = 1; kb <= n; kb++) {
    k = n + 1 - kb;
    if (fabs(z[k - 1]) > fabs(abd[m - 1 + (k - 1) * lda])) {
      s = fabs(abd[m - 1 + (k - 1) * lda]) / fabs(z[k - 1]);
      c_sscal(n, s, z);
      ynorm *= s;
    }
    if (abd[m - 1 + (k - 1) * lda] != 0.) {
      z[k - 1] /= abd[m - 1 + (k - 1) * lda];
    } else {
      z[k - 1] = 1.;
    }
    lm = ({
           const int _i = (int)(k);
           const int _j = (int)(m);
           _i < _j ? _i : _j;
         }) -
         1;
    la = m - lm;
    lz = k - lm;
    t = -z[k - 1];
    c_saxpy(lm, t, &abd[la - 1 + (k - 1) * lda], &z[lz - 1]);
  }
  s = 1. / c_sasum(n, z);
  c_sscal(n, s, z);
  ynorm *= s;
  if (anorm != 0.) {
    *rcond = ynorm / anorm;
  } else {
    *rcond = 0.;
  }
  return;
}


x86 assembly (clang, -O0):
<c_sgbco>:
push   %rbp
mov    %rsp,%rbp
sub    $0x120,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
add    -0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     12ab <c_sgbco+0x14b>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x40(%rbp),%edi
mov    -0x8(%rbp),%rsi
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    $0x0,%al
call   1070 <c_sasum@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
ucomisd -0xa8(%rbp),%xmm0
ff
jbe    1222 <c_sgbco+0xc2>
movsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xf8(%rbp)
ff
jmp    1232 <c_sgbco+0xd2>
movsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0xf8(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0x60(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jle    1269 <c_sgbco+0x109>
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     127e <c_sgbco+0x11e>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x14(%rbp),%ecx
cmp    %ecx,%eax
jl     1298 <c_sgbco+0x138>
mov    -0x40(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x40(%rbp)
jmp    129d <c_sgbco+0x13d>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11a9 <c_sgbco+0x49>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%r8d
mov    -0x20(%rbp),%r9
lea    -0x24(%rbp),%rax
mov    %rax,(%rsp)
mov    $0x0,%al
call   1050 <c_sgbfa@plt>
movsd  0xd29(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    0x18(%rbp),%rdi
movslq -0x10(%rbp),%rdx
shl    $0x3,%rdx
xor    %esi,%esi
call   1040 <memset@plt>
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1874 <c_sgbco+0x714>
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    133f <c_sgbco+0x1df>
jp     133f <c_sgbco+0x1df>
jmp    13ce <c_sgbco+0x26e>
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1392 <c_sgbco+0x232>
movsd  -0x68(%rbp),%xmm0
movaps 0xc8f(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
jmp    13c1 <c_sgbco+0x261>
movsd  -0x68(%rbp),%xmm0
movaps 0xc72(%rip),%xmm1
pand   %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x100(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movaps 0xc24(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xbfb(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1499 <c_sgbco+0x339>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xbc8(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  -0x68(%rbp),%xmm1
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
movaps 0xba6(%rip),%xmm2
pand   %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  -0x70(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movaps 0xb13(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
movaps 0xafb(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1554 <c_sgbco+0x3f4>
jp     1554 <c_sgbco+0x3f4>
jmp    15bd <c_sgbco+0x45d>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x88(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x90(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
jmp    15dd <c_sgbco+0x47d>
movsd  0xa3b(%rip),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xa2b(%rip),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0xb8(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x20(%rbp),%rcx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
add    (%rcx,%rdx,4),%eax
mov    %eax,-0xbc(%rbp)
mov    -0xb8(%rbp),%eax
cmp    -0xbc(%rbp),%eax
jle    162b <c_sgbco+0x4cb>
mov    -0xb8(%rbp),%eax
mov    %eax,-0x104(%rbp)
jmp    1637 <c_sgbco+0x4d7>
mov    -0xbc(%rbp),%eax
mov    %eax,-0x104(%rbp)
mov    -0x104(%rbp),%eax
mov    %eax,-0xc0(%rbp)
mov    -0xc0(%rbp),%eax
mov    %eax,-0xb4(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0xc4(%rbp)
mov    -0xb4(%rbp),%eax
cmp    -0xc4(%rbp),%eax
jge    167b <c_sgbco+0x51b>
mov    -0xb4(%rbp),%eax
mov    %eax,-0x108(%rbp)
jmp    1687 <c_sgbco+0x527>
mov    -0xc4(%rbp),%eax
mov    %eax,-0x108(%rbp)
mov    -0x108(%rbp),%eax
mov    %eax,-0xc8(%rbp)
mov    -0xc8(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     184c <c_sgbco+0x6ec>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     1798 <c_sgbco+0x638>
mov    -0x54(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x54(%rbp)
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movaps 0x900(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0x894(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    16b4 <c_sgbco+0x554>
movsd  -0x70(%rbp),%xmm1
movsd  -0x78(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1847 <c_sgbco+0x6e7>
movsd  -0x90(%rbp),%xmm0
ff
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x80(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x50(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     1842 <c_sgbco+0x6e2>
mov    -0x54(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x54(%rbp)
movsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    17dd <c_sgbco+0x67d>
jmp    1847 <c_sgbco+0x6e7>
jmp    184c <c_sgbco+0x6ec>
movsd  -0x88(%rbp),%xmm0
ff
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1309 <c_sgbco+0x1a9>
mov    -0x10(%rbp),%edi
mov    0x18(%rbp),%rsi
mov    $0x0,%al
call   1070 <c_sasum@plt>
cvtsi2sd %eax,%xmm1
movsd  0x772(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1a6e <c_sgbco+0x90e>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
sub    -0x38(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0xcc(%rbp)
mov    -0x10(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0xd0(%rbp)
mov    -0xcc(%rbp),%eax
cmp    -0xd0(%rbp),%eax
jge    1901 <c_sgbco+0x7a1>
mov    -0xcc(%rbp),%eax
mov    %eax,-0x10c(%rbp)
jmp    190d <c_sgbco+0x7ad>
mov    -0xd0(%rbp),%eax
mov    %eax,-0x10c(%rbp)
mov    -0x10c(%rbp),%eax
mov    %eax,-0xd4(%rbp)
mov    -0xd4(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    198b <c_sgbco+0x82b>
mov    -0x48(%rbp),%edi
mov    -0x8(%rbp),%rsi
mov    -0x50(%rbp),%eax
add    $0x1,%eax
sub    $0x1,%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    0x18(%rbp),%rdx
mov    -0x34(%rbp),%eax
add    $0x1,%eax
sub    $0x1,%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
mov    $0x0,%al
call   1060 <c_sdot@plt>
cvtsi2sd %eax,%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0x66c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x650(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    19fb <c_sgbco+0x89b>
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0x63d(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0x621(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    18b1 <c_sgbco+0x751>
mov    -0x10(%rbp),%edi
mov    0x18(%rbp),%rsi
mov    $0x0,%al
call   1070 <c_sasum@plt>
cvtsi2sd %eax,%xmm1
movsd  0x578(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  0x554(%rip),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1c70 <c_sgbco+0xb10>
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
mov    %eax,-0xd8(%rbp)
mov    -0x10(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0xdc(%rbp)
mov    -0xd8(%rbp),%eax
cmp    -0xdc(%rbp),%eax
jge    1b64 <c_sgbco+0xa04>
mov    -0xd8(%rbp),%eax
mov    %eax,-0x110(%rbp)
jmp    1b70 <c_sgbco+0xa10>
mov    -0xdc(%rbp),%eax
mov    %eax,-0x110(%rbp)
mov    -0x110(%rbp),%eax
mov    %eax,-0xe0(%rbp)
mov    -0xe0(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1bd8 <c_sgbco+0xa78>
mov    -0x48(%rbp),%edi
movsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rsi
mov    -0x50(%rbp),%eax
add    $0x1,%eax
sub    $0x1,%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    0x18(%rbp),%rdx
mov    -0x34(%rbp),%eax
add    $0x1,%eax
sub    $0x1,%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
mov    $0x1,%al
call   1080 <c_saxpy@plt>
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0x41f(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0x403(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1c5d <c_sgbco+0xafd>
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0x3f0(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0x3d4(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  -0x70(%rbp),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1c62 <c_sgbco+0xb02>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1abb <c_sgbco+0x95b>
mov    -0x10(%rbp),%edi
mov    0x18(%rbp),%rsi
mov    $0x0,%al
call   1070 <c_sasum@plt>
cvtsi2sd %eax,%xmm1
movsd  0x376(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  -0x70(%rbp),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1f17 <c_sgbco+0xdb7>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
sub    -0x38(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0x31d(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0x2f4(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1da1 <c_sgbco+0xc41>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0x2c1(%rip),%xmm1
pand   %xmm1,%xmm0
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0x2a4(%rip),%xmm2
pand   %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  -0x70(%rbp),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1dd7 <c_sgbco+0xc77>
jp     1dd7 <c_sgbco+0xc77>
jmp    1e15 <c_sgbco+0xcb5>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1e2f <c_sgbco+0xccf>
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  0x1d6(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
mov    %eax,-0xe4(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0xe8(%rbp)
mov    -0xe4(%rbp),%eax
cmp    -0xe8(%rbp),%eax
jge    1e64 <c_sgbco+0xd04>
mov    -0xe4(%rbp),%eax
mov    %eax,-0x114(%rbp)
jmp    1e70 <c_sgbco+0xd10>
mov    -0xe8(%rbp),%eax
mov    %eax,-0x114(%rbp)
mov    -0x114(%rbp),%eax
mov    %eax,-0xec(%rbp)
mov    -0xec(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    -0x50(%rbp),%eax
sub    -0x48(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x34(%rbp),%eax
sub    -0x48(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x48(%rbp),%edi
movsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rsi
mov    -0x44(%rbp),%eax
sub    $0x1,%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
imul   -0xc(%rbp),%ecx
add    %ecx,%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    0x18(%rbp),%rdx
mov    -0x4c(%rbp),%eax
sub    $0x1,%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
mov    $0x1,%al
call   1080 <c_saxpy@plt>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1cc2 <c_sgbco+0xb62>
mov    -0x10(%rbp),%edi
mov    0x18(%rbp),%rsi
mov    $0x0,%al
call   1070 <c_sasum@plt>
cvtsi2sd %eax,%xmm1
movsd  0xcf(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%edi
movsd  -0x70(%rbp),%xmm0
mov    0x18(%rbp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  -0x70(%rbp),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1f7f <c_sgbco+0xe1f>
jp     1f7f <c_sgbco+0xe1f>
jmp    1f99 <c_sgbco+0xe39>
movsd  -0x98(%rbp),%xmm0
ff
divsd  -0x60(%rbp),%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    1fa4 <c_sgbco+0xe44>
mov    0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
add    $0x120,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_sgbco>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %r9,0x40(%rsp)
mov    %rdi,0x20(%rsp)
test   %edx,%edx
mov    %rdx,0x50(%rsp)
mov    %rcx,0x68(%rsp)
mov    %r8,0x18(%rsp)
jle    125a <c_sgbco+0xfa>
lea    0x1(%rcx),%eax
lea    (%rax,%r8,1),%r14d
mov    %edx,%edi
sub    %ecx,%edi
mov    %edi,0x30(%rsp)
mov    %rdx,%r12
neg    %r12d
mov    %r12,0x10(%rsp)
xorpd  %xmm0,%xmm0
mov    $0x1,%ebp
mov    $0xffffffff,%ebx
mov    %esi,%r12d
mov    %rax,0x28(%rsp)
mov    %eax,%r13d
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm0,0x48(%rsp)
lea    -0x1(%r14),%r15d
lea    (%r14,%rbx,1),%eax
cltq
mov    0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%rsi
mov    %r13d,%edi
xor    %eax,%eax
call   1070 <c_sasum@plt>
mov    0x20(%rsp),%rdi
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x48(%rsp),%xmm1
maxsd  %xmm0,%xmm1
movsd  %xmm1,0x48(%rsp)
movsd  0x48(%rsp),%xmm0
cmp    0x28(%rsp),%r14d
cmovle %r14d,%r15d
xor    %eax,%eax
cmp    0x18(%rsp),%ebp
setle  %al
add    %r13d,%eax
xor    %ecx,%ecx
cmp    0x30(%rsp),%ebp
setge  %cl
sub    %ecx,%eax
mov    0x10(%rsp),%rcx
add    %ebp,%ecx
add    $0x1,%ecx
mov    %ebp,%edx
add    $0x1,%edx
add    %r12d,%ebx
mov    %r15d,%r14d
mov    %edx,%ebp
mov    %eax,%r13d
cmp    $0x1,%ecx
jne    11d0 <c_sgbco+0x70>
jmp    1266 <c_sgbco+0x106>
xorpd  %xmm0,%xmm0
mov    %esi,%r12d
mov    0x20(%rsp),%rdi
mov    %r12d,%r15d
mov    %r12d,0x38(%rsp)
movsd  %xmm0,0x48(%rsp)
lea    0xec(%rsp),%rax
mov    %rax,(%rsp)
mov    %rdi,%r14
mov    %r12d,%esi
mov    0x50(%rsp),%r12
mov    %r12d,%edx
mov    0x68(%rsp),%r13
mov    %r13d,%ecx
mov    0x18(%rsp),%rbp
mov    %ebp,%r8d
mov    0x40(%rsp),%r9
xor    %eax,%eax
call   1050 <c_sgbfa@plt>
movslq %r12d,%rax
mov    %rax,0x80(%rsp)
lea    0x0(,%rax,8),%rdx
mov    0x148(%rsp),%rbx
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
lea    0x0(%rbp,%r13,1),%edx
lea    0x1(%rbp,%r13,1),%eax
mov    %r12d,%ecx
mov    %rcx,0x30(%rsp)
movslq %r15d,%rcx
test   %r12d,%r12d
mov    %eax,0x3c(%rsp)
mov    %rcx,0x70(%rsp)
mov    %rdx,0x78(%rsp)
jle    1839 <c_sgbco+0x6d9>
movslq %eax,%rsi
movslq %edx,%rax
mov    %rax,0xe0(%rsp)
lea    0x1(%r12),%eax
mov    %rax,0xd8(%rsp)
lea    -0x1(%rsi),%rax
mov    %rax,0xb8(%rsp)
lea    -0x2(%rsi),%rax
mov    %rax,0xb0(%rsp)
mov    %rsi,0xc0(%rsp)
lea    -0x10(%r14,%rsi,8),%rax
mov    %rax,0xd0(%rsp)
lea    -0x8(,%rcx,8),%r15
ff
lea    -0x8(%rbx),%rax
mov    %rax,0xc8(%rsp)
lea    -0x18(%r14),%rax
mov    %rax,0xa8(%rsp)
lea    -0x2(%rcx,%rcx,1),%rax
mov    %rax,0x88(%rsp)
lea    -0x10(%r14),%rax
mov    %rax,0xa0(%rsp)
movsd  0xc7c(%rip),%xmm7
mov    $0x2,%edi
mov    $0x1,%r9d
mov    $0xffffffff,%r13d
xorpd  %xmm8,%xmm8
movapd 0xc6e(%rip),%xmm6
movapd 0xc75(%rip),%xmm9
movapd 0xc5c(%rip),%xmm10
movapd 0xc73(%rip),%xmm11
xor    %r12d,%r12d
mov    %rbx,%rsi
mov    %rcx,%rdx
xor    %ebx,%ebx
jmp    13f3 <c_sgbco+0x293>
nopw   0x0(%rax,%rax,1)
mov    0x28(%rsp),%rax
movsd  %xmm0,-0x8(%rsi,%rax,8)
add    $0x1,%edi
add    $0x1,%ebx
add    $0xffffffff,%r13d
cmp    0xd8(%rsp),%r9
je     16d4 <c_sgbco+0x574>
mov    %r9,%rcx
lea    -0x1(%r9),%rbp
movsd  -0x8(%rsi,%r9,8),%xmm1
ucomisd %xmm8,%xmm1
mov    %rbx,0x10(%rsp)
jne    140f <c_sgbco+0x2af>
jnp    141f <c_sgbco+0x2bf>
ucomisd %xmm1,%xmm8
andpd  %xmm6,%xmm7
jae    141f <c_sgbco+0x2bf>
xorpd  %xmm9,%xmm7
movapd %xmm7,%xmm4
subsd  %xmm1,%xmm4
movapd %xmm4,%xmm3
andpd  %xmm10,%xmm3
imul   %rdx,%rbp
add    0xe0(%rsp),%rbp
mov    0x20(%rsp),%r14
movsd  (%r14,%rbp,8),%xmm5
movapd %xmm5,%xmm0
andpd  %xmm10,%xmm0
ucomisd %xmm0,%xmm3
mov    %rcx,0x28(%rsp)
jbe    14fb <c_sgbco+0x39b>
divsd  %xmm3,%xmm0
movapd %xmm0,0xf0(%rsp)
mov    %r12,%rbx
mov    %r13,%r12
mov    %edi,%r13d
mov    0x50(%rsp),%rdi
mov    $0x1,%al
movapd %xmm7,0x90(%rsp)
call   1030 <c_sscal@plt>
mov    0x28(%rsp),%rcx
movapd 0xb98(%rip),%xmm11
movapd 0xb6f(%rip),%xmm10
movapd 0xb76(%rip),%xmm9
xorpd  %xmm8,%xmm8
mov    %r13d,%edi
mov    %r12,%r13
mov    %rbx,%r12
movapd 0x90(%rsp),%xmm7
movapd 0xb47(%rip),%xmm6
mov    0x70(%rsp),%rdx
mov    0x148(%rsp),%rsi
mulsd  0xf0(%rsp),%xmm7
movsd  -0x8(%rsi,%rcx,8),%xmm1
movsd  (%r14,%rbp,8),%xmm5
movapd %xmm7,%xmm4
subsd  %xmm1,%xmm4
movapd %xmm4,%xmm3
andpd  %xmm6,%xmm3
mov    0x10(%rsp),%rbx
mov    0x88(%rsp),%r14
movapd %xmm7,%xmm2
xorpd  %xmm9,%xmm2
subsd  %xmm1,%xmm2
ucomisd %xmm8,%xmm5
movapd %xmm11,%xmm1
jne    1523 <c_sgbco+0x3c3>
jnp    1533 <c_sgbco+0x3d3>
movapd %xmm2,%xmm1
unpcklpd %xmm4,%xmm1
unpcklpd %xmm5,%xmm5
divpd  %xmm5,%xmm1
mov    0x40(%rsp),%rax
mov    -0x4(%rax,%rcx,4),%eax
add    0x18(%rsp),%eax
cmp    %eax,%r12d
cmovle %eax,%r12d
mov    0x50(%rsp),%rax
cmp    %eax,%r12d
cmovge %eax,%r12d
lea    0x1(%rcx),%r9
movslq %r12d,%rax
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
cmp    %rax,%rcx
jge    13d0 <c_sgbco+0x270>
mov    %edi,0x90(%rsp)
movslq %edi,%rcx
lea    -0x1(%rcx),%r10
imul   %rdx,%r10
mov    0xd0(%rsp),%rax
lea    (%rax,%r10,8),%rbp
mov    0xc8(%rsp),%rax
lea    (%rax,%rcx,8),%rdx
mov    %ebx,%r11d
not    %r11d
lea    0x2(%rbx),%r8d
andpd  %xmm6,%xmm2
unpcklpd %xmm3,%xmm2
lea    (%r12,%r13,1),%eax
xor    %edi,%edi
nopl   (%rax)
movsd  0x0(%rbp),%xmm3
unpcklpd %xmm3,%xmm3
movsd  (%rdx,%rdi,8),%xmm4
unpcklpd %xmm4,%xmm4
mulpd  %xmm1,%xmm3
addpd  %xmm4,%xmm3
movhpd %xmm3,(%rdx,%rdi,8)
andpd  %xmm6,%xmm3
addpd  %xmm3,%xmm2
add    %r15,%rbp
add    $0x1,%rdi
cmp    %edi,%eax
jne    15b0 <c_sgbco+0x450>
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
ucomisd %xmm3,%xmm2
mov    0x70(%rsp),%rdx
mov    0x90(%rsp),%edi
jbe    13d0 <c_sgbco+0x270>
movapd %xmm1,%xmm2
subsd  %xmm0,%xmm2
add    %r12d,%r11d
test   $0x1,%r11b
jne    161b <c_sgbco+0x4bb>
mov    0xc0(%rsp),%rbp
jmp    164a <c_sgbco+0x4ea>
add    0xb0(%rsp),%r10
mov    0x20(%rsp),%rax
movsd  (%rax,%r10,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  -0x8(%rsi,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rsi,%rcx,8)
add    $0x1,%rcx
mov    0xb8(%rsp),%rbp
movapd %xmm1,%xmm0
cmp    %r8d,%r12d
je     13d0 <c_sgbco+0x270>
lea    0x1(%r12),%r8d
mov    %rdx,%rax
imul   %rcx,%rax
mov    0xa8(%rsp),%rdi
lea    (%rdi,%rax,8),%rax
mov    %rdx,%rdi
lea    -0x1(%rcx),%rdx
imul   %rdi,%rdx
mov    0xa0(%rsp),%rdi
lea    (%rdi,%rdx,8),%rdx
mov    %r8d,%edi
nopl   0x0(%rax)
movsd  (%rdx,%rbp,8),%xmm1
mulsd  %xmm2,%xmm1
addsd  -0x8(%rsi,%rcx,8),%xmm1
movsd  %xmm1,-0x8(%rsi,%rcx,8)
movsd  (%rax,%rbp,8),%xmm1
mulsd  %xmm2,%xmm1
addsd  (%rsi,%rcx,8),%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x2,%rcx
add    %r14,%rbp
cmp    %ecx,%edi
jne    1690 <c_sgbco+0x530>
mov    0x70(%rsp),%rdx
mov    0x90(%rsp),%edi
jmp    13d0 <c_sgbco+0x270>
mov    0x50(%rsp),%r12
mov    %r12d,%edi
xor    %eax,%eax
mov    %rsi,%rbp
call   1070 <c_sasum@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x90b(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %r12d,%edi
mov    %rbp,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  0x8f2(%rip),%xmm3
test   %r12d,%r12d
movsd  %xmm3,0x18(%rsp)
jle    19da <c_sgbco+0x87a>
lea    -0x1(%r12),%eax
mov    0x38(%rsp),%r13d
imul   %r13d,%eax
mov    0x78(%rsp),%rcx
lea    (%rax,%rcx,1),%r14d
add    $0x1,%r14d
xor    %r15d,%r15d
mov    %r12d,%ebx
mov    %rbp,%rcx
mov    0x30(%rsp),%rdx
movapd 0x8c2(%rip),%xmm2
jmp    177d <c_sgbco+0x61d>
mov    0x40(%rsp),%rax
movslq (%rax,%rbp,4),%rax
movsd  -0x8(%rcx,%rax,8),%xmm1
movsd  %xmm0,-0x8(%rcx,%rax,8)
movsd  %xmm1,(%rcx,%rbp,8)
add    $0xffffffff,%ebx
sub    %r13d,%r14d
add    $0x1,%r15
cmp    %r15,%rdx
je     1869 <c_sgbco+0x709>
cmp    %r12d,%ebx
jge    17e0 <c_sgbco+0x680>
mov    0x68(%rsp),%rax
cmp    %eax,%r15d
mov    %r15d,%edi
cmovg  %eax,%edi
lea    -0x1(%rbx),%ebp
movslq %r14d,%rax
mov    0x20(%rsp),%rdx
lea    (%rdx,%rax,8),%rsi
movslq %ebx,%rax
lea    (%rcx,%rax,8),%rdx
xor    %eax,%eax
call   1060 <c_sdot@plt>
movsd  0x84b(%rip),%xmm3
movapd 0x853(%rip),%xmm2
mov    0x30(%rsp),%rdx
mov    0x148(%rsp),%rcx
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movslq %ebp,%rbp
addsd  (%rcx,%rbp,8),%xmm0
movsd  %xmm0,(%rcx,%rbp,8)
jmp    17e6 <c_sgbco+0x686>
lea    -0x1(%rbx),%eax
movslq %eax,%rbp
movsd  (%rcx,%rbp,8),%xmm0
movapd %xmm0,%xmm1
andpd  %xmm2,%xmm1
ucomisd %xmm3,%xmm1
jbe    1750 <c_sgbco+0x5f0>
movapd %xmm3,%xmm0
divsd  %xmm1,%xmm0
mov    %r12d,%edi
mov    %rcx,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  0x7e6(%rip),%xmm3
movapd 0x7ee(%rip),%xmm2
mov    0x30(%rsp),%rdx
mov    0x148(%rsp),%rcx
movsd  (%rcx,%rbp,8),%xmm0
jmp    1750 <c_sgbco+0x5f0>
mov    %r12d,%edi
mov    %rbx,%rsi
xor    %eax,%eax
call   1070 <c_sasum@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x7ab(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %r12d,%edi
mov    %rbx,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
mov    %rbx,%rcx
mov    %r12d,%edi
mov    %rcx,%rsi
xor    %eax,%eax
call   1070 <c_sasum@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x77b(%rip),%xmm0
movsd  %xmm0,0x18(%rsp)
divsd  %xmm1,%xmm0
mov    %r12d,%edi
mov    0x148(%rsp),%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
mov    0x148(%rsp),%rsi
test   %r12d,%r12d
mov    0x70(%rsp),%r14
jle    1a0a <c_sgbco+0x8aa>
movslq 0x3c(%rsp),%rax
movslq 0x68(%rsp),%rcx
mov    %rcx,0x10(%rsp)
mov    0x50(%rsp),%rcx
lea    0x1(%rcx),%r15d
mov    0x80(%rsp),%rcx
lea    -0x1(%rcx),%rbx
mov    0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%r12
shl    $0x3,%r14
movsd  0x710(%rip),%xmm2
mov    $0x1,%ebp
mov    $0x8,%r13d
movsd  %xmm2,0x18(%rsp)
jmp    192b <c_sgbco+0x7cb>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r14,%rax
add    $0x1,%rbp
add    $0xffffffffffffffff,%rbx
add    %r14,%r12
add    $0x8,%r13
cmp    %rbp,%r15
je     1a0a <c_sgbco+0x8aa>
mov    0x40(%rsp),%rax
movslq -0x4(%rax,%rbp,4),%rax
movsd  -0x8(%rsi,%rax,8),%xmm0
movsd  -0x8(%rsi,%rbp,8),%xmm1
movsd  %xmm1,-0x8(%rsi,%rax,8)
movsd  %xmm0,-0x8(%rsi,%rbp,8)
cmp    0x80(%rsp),%rbp
jge    1987 <c_sgbco+0x827>
cmp    0x10(%rsp),%rbx
mov    %ebx,%edi
cmovg  0x68(%rsp),%edi
lea    (%rsi,%r13,1),%rdx
mov    %r12,%rsi
mov    $0x1,%al
call   1080 <c_saxpy@plt>
movsd  0x687(%rip),%xmm2
mov    0x148(%rsp),%rsi
movsd  -0x8(%rsi,%rbp,8),%xmm0
andpd  0x681(%rip),%xmm0
ucomisd %xmm2,%xmm0
jbe    1910 <c_sgbco+0x7b0>
divsd  %xmm0,%xmm2
movsd  %xmm2,0x28(%rsp)
mov    0x50(%rsp),%rdi
movapd %xmm2,%xmm0
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  0x645(%rip),%xmm2
mov    0x148(%rsp),%rsi
movsd  0x18(%rsp),%xmm0
mulsd  0x28(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
jmp    1910 <c_sgbco+0x7b0>
mov    %r12d,%edi
mov    %rbp,%rsi
xor    %eax,%eax
call   1070 <c_sasum@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x60a(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %r12d,%edi
mov    %rbp,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
mov    %rbp,%rsi
xor    %eax,%eax
mov    %rax,0x28(%rsp)
mov    0x50(%rsp),%r14
mov    %r14d,%edi
xor    %eax,%eax
mov    %rsi,%rbx
call   1070 <c_sasum@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x5ce(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
mov    %r14d,%edi
movapd %xmm1,%xmm0
mov    %rbx,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
movsd  0x10(%rsp),%xmm0
mulsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
test   %r14d,%r14d
jle    1bbf <c_sgbco+0xa5f>
lea    -0x1(%r14),%eax
imul   0x38(%rsp),%eax
mov    0x78(%rsp),%rcx
mov    %ecx,%esi
add    %eax,%esi
lea    (%rax,%rcx,1),%r13d
add    $0x1,%r13d
mov    $0xffffffff,%edi
mov    0x148(%rsp),%rdx
mov    0x20(%rsp),%r15
mov    0x3c(%rsp),%ebp
mov    0x30(%rsp),%rcx
jmp    1b09 <c_sgbco+0x9a9>
nop
mov    0x3c(%rsp),%ebp
cmp    %ebp,%edi
cmovge %ebp,%edi
movsd  %xmm0,(%rdx,%rbx,8)
add    $0xffffffff,%edi
xorpd  0x567(%rip),%xmm0
movslq %ecx,%rcx
lea    (%r15,%rcx,8),%rsi
add    %r14d,%eax
cltq
lea    (%rdx,%rax,8),%rdx
mov    $0x1,%al
call   1080 <c_saxpy@plt>
mov    0x148(%rsp),%rdx
mov    0x28(%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,0x28(%rsp)
mov    0x18(%rsp),%rdi
add    $0xffffffff,%edi
mov    0x38(%rsp),%eax
sub    %eax,%r13d
mov    0x40(%rsp),%esi
sub    %eax,%esi
mov    0x30(%rsp),%rcx
add    $0xffffffffffffffff,%rcx
je     1bc7 <c_sgbco+0xa67>
cmp    %ebp,%ecx
cmovl  %ecx,%ebp
lea    (%r14,%rdi,1),%eax
movslq %eax,%rbx
movsd  (%rdx,%rbx,8),%xmm3
movapd 0x4ee(%rip),%xmm2
andpd  %xmm2,%xmm3
movslq %esi,%r12
movsd  (%r15,%r12,8),%xmm1
movapd %xmm1,%xmm0
andpd  %xmm2,%xmm0
ucomisd %xmm0,%xmm3
mov    %rcx,0x30(%rsp)
mov    %esi,0x40(%rsp)
mov    %rdi,0x18(%rsp)
jbe    1b87 <c_sgbco+0xa27>
divsd  %xmm3,%xmm0
movapd %xmm0,0x50(%rsp)
mov    %r14d,%edi
mov    %rdx,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
mov    0x20(%rsp),%r15
mov    0x148(%rsp),%rdx
movsd  0x10(%rsp),%xmm0
mulsd  0x50(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  (%r15,%r12,8),%xmm1
mov    0x28(%rsp),%rsi
mov    %esi,%eax
sub    %ebp,%eax
mov    %r13d,%ecx
sub    %ebp,%ecx
lea    (%r14,%rsi,1),%edi
ucomisd 0x467(%rip),%xmm1
movsd  0x457(%rip),%xmm0
jne    1bb1 <c_sgbco+0xa51>
jnp    1aa0 <c_sgbco+0x940>
movsd  (%rdx,%rbx,8),%xmm0
divsd  %xmm1,%xmm0
jmp    1aa0 <c_sgbco+0x940>
mov    0x148(%rsp),%rdx
mov    0x140(%rsp),%rbx
mov    %r14d,%edi
mov    %rdx,%rsi
xor    %eax,%eax
mov    %rdx,%rbp
call   1070 <c_sasum@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x412(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %r14d,%edi
movsd  %xmm0,0x20(%rsp)
mov    %rbp,%rsi
mov    $0x1,%al
call   1030 <c_sscal@plt>
xorpd  %xmm0,%xmm0
movsd  0x48(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jne    1c17 <c_sgbco+0xab7>
jnp    1c27 <c_sgbco+0xac7>
movsd  0x10(%rsp),%xmm0
mulsd  0x20(%rsp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx)
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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