#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 /*<<< orphan*/ atou1(double *, int, int);
scalar_t__ calloc(int, int);
int fabs(double);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ ldvmat(double *, double *, int);
int /*<<< orphan*/ qrbdu1(double *, double *, double *, int, double *, int);
double sqrt(double);
int svdu1v(double *d, double *a, int m, double *v, int n) {
  double *p, *p1, *q, *pp, *w, *e;
  double s, h, r, t, sv;
  int i, j, k, mm, nm, ms;
  if (m < n)
    return -1;
  w = (double *)calloc(m + n, sizeof(double));
  e = w + m;
  for (i = 0, mm = m, nm = n - 1, p = a; i < n; ++i, --mm, --nm, p += n + 1) {
    if (mm > 1) {
      sv = h = 0.;
      for (j = 0, q = p, s = 0.; j < mm; ++j, q += n) {
        w[j] = *q;
        s += *q * *q;
      }
      if (s > 0.) {
        h = sqrt(s);
        if (*p < 0.)
          h = -h;
        s += *p * h;
        s = 1. / s;
        t = 1. / (w[0] += h);
        sv = 1. + fabs(*p / h);
        for (k = 1, ms = n - i; k < ms; ++k) {
          for (j = 0, q = p + k, r = 0.; j < mm; q += n)
            r += w[j++] * *q;
          r *= s;
          for (j = 0, q = p + k; j < mm; q += n)
            *q -= r * w[j++];
        }
        for (j = 1, q = p; j < mm;)
          *(q += n) = t * w[j++];
      }
      *p = sv;
      d[i] = -h;
    }
    if (mm == 1)
      d[i] = *p;
    p1 = p + 1;
    sv = h = 0.;
    if (nm > 1) {
      for (j = 0, q = p1, s = 0.; j < nm; ++j, ++q)
        s += *q * *q;
      if (s > 0.) {
        h = sqrt(s);
        if (*p1 < 0.)
          h = -h;
        sv = 1. + fabs(*p1 / h);
        s += *p1 * h;
        s = 1. / s;
        t = 1. / (*p1 += h);
        for (k = n, ms = n * (m - i); k < ms; k += n) {
          for (j = 0, q = p1, pp = p1 + k, r = 0.; j < nm; ++j)
            r += *q++ * *pp++;
          r *= s;
          for (j = 0, q = p1, pp = p1 + k; j < nm; ++j)
            *pp++ -= r * *q++;
        }
        for (j = 1, q = p1 + 1; j < nm; ++j)
          *q++ *= t;
      }
      *p1 = sv;
      e[i] = -h;
    }
    if (nm == 1)
      e[i] = *p1;
  }
  ldvmat(a, v, n);
  atou1(a, m, n);
  qrbdu1(d, e, a, m, v, n);
  for (i = 0; i < n; ++i) {
    if (d[i] < 0.) {
      d[i] = -d[i];
      for (j = 0, p = v + i; j < n; ++j, p += n)
        *p = -*p;
    }
  }
  free(w);
  return 0;
}


x86 assembly (clang, -O0):
<svdu1v>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    11a6 <svdu1v+0x36>
movl   $0xffffffff,-0x4(%rbp)
jmp    1b14 <svdu1v+0x9a4>
mov    -0x1c(%rbp),%edi
add    -0x2c(%rbp),%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x8c(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x98(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x9c(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    19c3 <svdu1v+0x853>
cmpl   $0x1,-0x98(%rbp)
jle    15a5 <svdu1v+0x435>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movl   $0x0,-0x90(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x98(%rbp),%eax
jge    12b2 <svdu1v+0x142>
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
mov    -0x2c(%rbp),%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
jmp    123d <svdu1v+0xcd>
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1569 <svdu1v+0x3f9>
movsd  -0x68(%rbp),%xmm0
call   1080 <sqrt@plt>
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jbe    1305 <svdu1v+0x195>
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0xcd4(%rip),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x70(%rbp),%xmm1
mov    -0x58(%rbp),%rax
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
movsd  0xcb1(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
divsd  -0x70(%rbp),%xmm0
call   1050 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xc8a(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movl   $0x1,-0x94(%rbp)
mov    -0x2c(%rbp),%eax
sub    -0x8c(%rbp),%eax
mov    %eax,-0xa0(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xa0(%rbp),%eax
jge    14f9 <svdu1v+0x389>
movl   $0x0,-0x90(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x94(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x98(%rbp),%eax
jge    1438 <svdu1v+0x2c8>
mov    -0x58(%rbp),%rax
mov    -0x90(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x90(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x2c(%rbp),%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
jmp    13d5 <svdu1v+0x265>
movsd  -0x68(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x94(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x98(%rbp),%eax
jge    14e0 <svdu1v+0x370>
movsd  -0x78(%rbp),%xmm0
mov    -0x58(%rbp),%rax
mov    -0x90(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x90(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
jmp    1467 <svdu1v+0x2f7>
jmp    14e5 <svdu1v+0x375>
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    139b <svdu1v+0x22b>
movl   $0x1,-0x90(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x98(%rbp),%eax
jge    1564 <svdu1v+0x3f4>
movsd  -0x80(%rbp),%xmm0
mov    -0x58(%rbp),%rax
mov    -0x90(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x90(%rbp)
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x2c(%rbp),%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
mov    %rcx,%rsi
shl    $0x3,%rsi
mov    %rax,%rdx
add    %rsi,%rdx
mov    %rdx,-0x48(%rbp)
movsd  %xmm0,(%rax,%rcx,8)
jmp    150b <svdu1v+0x39b>
jmp    1569 <svdu1v+0x3f9>
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x1,-0x98(%rbp)
jne    15ca <svdu1v+0x45a>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
cmpl   $0x1,-0x9c(%rbp)
jle    194f <svdu1v+0x7df>
movl   $0x0,-0x90(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x9c(%rbp),%eax
jge    1664 <svdu1v+0x4f4>
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x48(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
mov    -0x48(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x48(%rbp)
jmp    1610 <svdu1v+0x4a0>
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1913 <svdu1v+0x7a3>
movsd  -0x68(%rbp),%xmm0
call   1080 <sqrt@plt>
movsd  %xmm0,-0x70(%rbp)
mov    -0x40(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jbe    16b7 <svdu1v+0x547>
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
divsd  -0x70(%rbp),%xmm0
call   1050 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0x92b(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  0x8f8(%rip),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x70(%rbp),%xmm1
mov    -0x40(%rbp),%rax
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
movsd  0x8d5(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x94(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    -0x8c(%rbp),%ecx
imul   %ecx,%eax
mov    %eax,-0xa0(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xa0(%rbp),%eax
jge    18b6 <svdu1v+0x746>
movl   $0x0,-0x90(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rax
movslq -0x94(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x9c(%rbp),%eax
jge    17f2 <svdu1v+0x682>
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  (%rax),%xmm0
mov    -0x50(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x50(%rbp)
movsd  (%rax),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1794 <svdu1v+0x624>
movsd  -0x68(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rax
movslq -0x94(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x9c(%rbp),%eax
jge    189d <svdu1v+0x72d>
movsd  -0x78(%rbp),%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  (%rax),%xmm2
mov    -0x50(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x50(%rbp)
movsd  (%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1829 <svdu1v+0x6b9>
jmp    18a2 <svdu1v+0x732>
mov    -0x2c(%rbp),%eax
add    -0x94(%rbp),%eax
mov    %eax,-0x94(%rbp)
jmp    1752 <svdu1v+0x5e2>
movl   $0x1,-0x90(%rbp)
mov    -0x40(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x48(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x9c(%rbp),%eax
jge    190e <svdu1v+0x79e>
movsd  -0x80(%rbp),%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    18cc <svdu1v+0x75c>
jmp    1913 <svdu1v+0x7a3>
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x40(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x60(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x1,-0x9c(%rbp)
jne    1974 <svdu1v+0x804>
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1979 <svdu1v+0x809>
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
mov    -0x98(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x98(%rbp)
mov    -0x9c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x9c(%rbp)
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    11f4 <svdu1v+0x84>
mov    -0x18(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    -0x2c(%rbp),%edx
call   1060 <ldvmat@plt>
mov    -0x18(%rbp),%rdi
mov    -0x1c(%rbp),%esi
mov    -0x2c(%rbp),%edx
call   1070 <atou1@plt>
mov    -0x10(%rbp),%rdi
mov    -0x60(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x28(%rbp),%r8
mov    -0x2c(%rbp),%r9d
call   1090 <qrbdu1@plt>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1b04 <svdu1v+0x994>
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1aeb <svdu1v+0x97b>
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x90(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x8c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1ae6 <svdu1v+0x976>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
mov    -0x2c(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    1a87 <svdu1v+0x917>
jmp    1aeb <svdu1v+0x97b>
jmp    1af0 <svdu1v+0x980>
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    1a08 <svdu1v+0x898>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<svdu1v>:
mov    $0xffffffff,%eax
cmp    %r8d,%edx
jl     215a <svdu1v+0xfea>
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x168,%rsp
mov    %r8d,%r14d
mov    %rcx,%r15
mov    %edx,%ebp
mov    %rsi,%r12
mov    %rdi,0x18(%rsp)
lea    (%r14,%rbp,1),%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,%rbx
movslq %ebp,%rax
mov    %rax,0x118(%rsp)
lea    (%rbx,%rax,8),%r10
test   %r14d,%r14d
jle    2109 <svdu1v+0xf99>
mov    %r15,0xe0(%rsp)
mov    %r14d,%r13d
lea    0x1(%r13),%rax
mov    %rax,0xf8(%rsp)
lea    -0x3(%r14),%r11d
lea    0x2(%r13),%rax
mov    %rax,0xf0(%rsp)
lea    -0x2(%r14),%edi
lea    -0x1(%r14),%eax
lea    0x8(,%r13,8),%r8
mov    %r13,%r9
shl    $0x4,%r9
mov    %rbp,%rcx
lea    0x8(%r12),%r15
mov    %r13,%rbp
shl    $0x5,%rbp
mov    %rbp,0xa8(%rsp)
lea    0x0(,%r13,8),%rbp
lea    0x0(,%rbp,2),%rsi
add    %rbp,%rsi
mov    %rsi,0x110(%rsp)
lea    0x3(%rcx),%esi
lea    0x8(%r12,%r13,8),%rdx
movapd 0x1db8(%rip),%xmm7
mov    %r11d,0xb8(%rsp)
mov    %r11d,0x5c(%rsp)
mov    %edi,0x64(%rsp)
mov    %edi,0x60(%rsp)
mov    %rdx,%rdi
mov    %eax,0x34(%rsp)
mov    %rax,0x50(%rsp)
mov    %r12,0xa0(%rsp)
mov    %r12,0x10(%rsp)
xor    %r12d,%r12d
mov    %r14,0xe8(%rsp)
mov    %r14,0x48(%rsp)
mov    %rcx,0xd8(%rsp)
mov    %ecx,%r14d
mov    0x18(%rsp),%r11
mov    %r10,0x70(%rsp)
mov    %r8,0x78(%rsp)
mov    %r9,0x88(%rsp)
jmp    12f3 <svdu1v+0x183>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r12
mov    0x28(%rsp),%rax
lea    -0x1(%rax),%r14d
add    %r8,0x10(%rsp)
add    %r8,%r15
add    $0x3,%sil
mov    0x50(%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,0x50(%rsp)
add    %r8,%rdi
addl   $0xffffffff,0x60(%rsp)
addl   $0xffffffff,0x5c(%rsp)
mov    %edx,%eax
mov    %rax,0x48(%rsp)
cmp    %r13,%r12
je     1f99 <svdu1v+0xe29>
movzbl %sil,%esi
and    $0x3,%esi
mov    %r14d,%eax
cmp    $0x2,%r14d
mov    %rax,0x28(%rsp)
jl     1792 <svdu1v+0x622>
mov    %r12,0x20(%rsp)
lea    -0x1(%rax),%rdx
mov    %rdx,%rcx
mov    %rdx,(%rsp)
and    $0xfffffffffffffffc,%rdx
mov    %rdx,0x80(%rsp)
and    $0x3,%r14d
mov    %eax,%r12d
and    $0xfffffffc,%r12d
mov    %eax,%ecx
and    $0xfffffffe,%ecx
add    $0xfffffffffffffffe,%rax
mov    %rax,0xb0(%rsp)
xorpd  %xmm0,%xmm0
mov    0x10(%rsp),%rax
xor    %edx,%edx
nop
movsd  (%rax),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
movsd  (%rax),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
movsd  (%rax,%r13,8),%xmm1
movsd  %xmm1,0x8(%rbx,%rdx,8)
movsd  (%rax,%r13,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
add    $0x2,%rdx
add    %r9,%rax
cmp    %rdx,%rcx
jne    1350 <svdu1v+0x1e0>
testb  $0x1,0x28(%rsp)
je     13a7 <svdu1v+0x237>
movsd  (%rax),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
movsd  (%rax),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
xorpd  %xmm5,%xmm5
ucomisd %xmm5,%xmm0
xorpd  %xmm1,%xmm1
jbe    1775 <svdu1v+0x605>
ucomisd 0x1c57(%rip),%xmm0
mov    %rsi,0x40(%rsp)
mov    %rdi,0x38(%rsp)
movsd  %xmm0,0x98(%rsp)
jb     13dc <svdu1v+0x26c>
sqrtsd %xmm0,%xmm0
jmp    13e9 <svdu1v+0x279>
call   1080 <sqrt@plt>
movapd 0x1c17(%rip),%xmm7
mov    0x10(%rsp),%rax
movsd  (%rax),%xmm2
movapd %xmm0,%xmm1
xorpd  %xmm7,%xmm1
movsd  %xmm2,0x90(%rsp)
cmpltsd 0x1c0c(%rip),%xmm2
andpd  %xmm2,%xmm1
andnpd %xmm0,%xmm2
orpd   %xmm1,%xmm2
movapd %xmm2,0xc0(%rsp)
movsd  (%rbx),%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  %xmm0,(%rbx)
movsd  (%rax),%xmm0
divsd  %xmm2,%xmm0
call   1050 <fabs@plt>
movapd 0xc0(%rsp),%xmm5
movapd 0x1baf(%rip),%xmm7
mov    %r13,%rcx
sub    0x20(%rsp),%rcx
cmp    $0x2,%rcx
mov    0x18(%rsp),%r11
movsd  0x1ba6(%rip),%xmm4
mov    0x28(%rsp),%rdi
jl     1680 <svdu1v+0x510>
movsd  0x90(%rsp),%xmm1
mulsd  %xmm5,%xmm1
addsd  0x98(%rsp),%xmm1
movapd %xmm4,%xmm0
divsd  %xmm1,%xmm0
mov    0x48(%rsp),%r8d
mov    $0x1,%r10d
mov    %r15,%r9
jmp    14c1 <svdu1v+0x351>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r10
add    $0x8,%r9
cmp    %r8,%r10
je     1680 <svdu1v+0x510>
test   %edi,%edi
jle    14b0 <svdu1v+0x340>
mov    0x10(%rsp),%rcx
lea    (%rcx,%r10,8),%rdx
cmpq   $0x3,(%rsp)
jae    14f0 <svdu1v+0x380>
xorpd  %xmm1,%xmm1
xor    %esi,%esi
mov    %rdx,%rcx
test   $0x3,%dil
jne    1567 <svdu1v+0x3f7>
jmp    1589 <svdu1v+0x419>
nopl   (%rax)
xorpd  %xmm1,%xmm1
mov    %r9,%rdi
xor    %esi,%esi
nopl   0x0(%rax)
mov    %rdi,%rcx
movsd  (%rbx,%rsi,8),%xmm2
movsd  0x8(%rbx,%rsi,8),%xmm3
mulsd  (%rdi),%xmm2
addsd  %xmm1,%xmm2
add    %rbp,%rdi
mulsd  (%rcx,%rbp,1),%xmm3
addsd  %xmm2,%xmm3
movsd  0x10(%rbx,%rsi,8),%xmm2
mulsd  0x0(%rbp,%rdi,1),%xmm2
add    %rbp,%rdi
addsd  %xmm3,%xmm2
movsd  0x18(%rbx,%rsi,8),%xmm1
add    $0x4,%rsi
mulsd  0x0(%rbp,%rdi,1),%xmm1
add    %rbp,%rdi
addsd  %xmm2,%xmm1
add    %rbp,%rdi
cmp    %rsi,%r12
jne    1500 <svdu1v+0x390>
add    0xa8(%rsp),%rcx
mov    0x28(%rsp),%rdi
test   $0x3,%dil
je     1589 <svdu1v+0x419>
lea    (%rbx,%rsi,8),%rsi
xor    %edi,%edi
nopl   (%rax)
movsd  (%rsi,%rdi,8),%xmm2
mulsd  (%rcx),%xmm2
addsd  %xmm2,%xmm1
add    $0x1,%rdi
add    %rbp,%rcx
cmp    %rdi,%r14
jne    1570 <svdu1v+0x400>
mov    0x28(%rsp),%rdi
test   %edi,%edi
jle    14b0 <svdu1v+0x340>
xorpd  %xmm7,%xmm1
mulsd  %xmm0,%xmm1
cmpq   $0x3,(%rsp)
jae    15b6 <svdu1v+0x446>
xor    %esi,%esi
test   $0x3,%dil
jne    1642 <svdu1v+0x4d2>
jmp    14b0 <svdu1v+0x340>
mov    %r9,%rcx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
mov    %rcx,%rdx
movsd  (%rbx,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rcx),%xmm2
movsd  %xmm2,(%rcx)
add    %rbp,%rcx
movsd  0x8(%rbx,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rdx,%rbp,1),%xmm2
movsd  %xmm2,(%rdx,%rbp,1)
lea    (%rcx,%rbp,1),%rdi
movsd  0x10(%rbx,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x0(%rbp,%rcx,1),%xmm2
movsd  %xmm2,0x0(%rbp,%rcx,1)
lea    (%rdi,%rbp,1),%rcx
movsd  0x18(%rbx,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x0(%rbp,%rdi,1),%xmm2
add    $0x4,%rsi
movsd  %xmm2,0x0(%rbp,%rdi,1)
add    %rbp,%rcx
cmp    %rsi,%r12
jne    15c0 <svdu1v+0x450>
add    0xa8(%rsp),%rdx
mov    0x28(%rsp),%rdi
test   $0x3,%dil
je     14b0 <svdu1v+0x340>
lea    (%rbx,%rsi,8),%rcx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rdx),%xmm2
movsd  %xmm2,(%rdx)
add    $0x1,%rsi
add    %rbp,%rdx
cmp    %rsi,%r14
jne    1650 <svdu1v+0x4e0>
jmp    14b0 <svdu1v+0x340>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm4,%xmm0
divsd  0x8(%rsp),%xmm0
mov    $0x1,%ecx
mov    0x10(%rsp),%r14
mov    %r14,%rdx
cmpq   $0x3,0xb0(%rsp)
mov    0x70(%rsp),%r10
mov    0x78(%rsp),%r8
mov    0x88(%rsp),%r9
mov    0xa8(%rsp),%rdi
mov    0x110(%rsp),%r12
jb     172e <svdu1v+0x5be>
mov    %r14,%rdx
xor    %ecx,%ecx
mov    0x80(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x8(%rbx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rbp,1)
movsd  0x10(%rbx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%r9,1)
movsd  0x18(%rbx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%r12,1)
movsd  0x20(%rbx,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rdi,1)
add    %rdi,%rdx
add    $0x4,%rcx
cmp    %rcx,%rsi
jne    16e0 <svdu1v+0x570>
add    $0x1,%rcx
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  %xmm4,%xmm1
testb  $0x3,(%rsp)
je     176b <svdu1v+0x5fb>
lea    (%rbx,%rcx,8),%rax
lea    (%rdx,%r13,8),%rcx
xor    %edx,%edx
mov    0x40(%rsp),%rsi
xchg   %ax,%ax
movsd  (%rax,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%rcx)
add    $0x1,%rdx
add    %rbp,%rcx
cmp    %rdx,%rsi
jne    1750 <svdu1v+0x5e0>
jmp    1770 <svdu1v+0x600>
mov    0x40(%rsp),%rsi
mov    0x38(%rsp),%rdi
mov    0x10(%rsp),%rax
movsd  %xmm1,(%rax)
xorpd  %xmm7,%xmm5
mov    0x20(%rsp),%r12
movlpd %xmm5,(%r11,%r12,8)
mov    0x28(%rsp),%rax
cmp    $0x1,%eax
mov    0x10(%rsp),%rax
jne    17a6 <svdu1v+0x636>
movsd  (%rax),%xmm0
movsd  %xmm0,(%r11,%r12,8)
mov    0x48(%rsp),%rcx
lea    -0x1(%rcx),%edx
add    $0x8,%rax
mov    %rax,(%rsp)
cmp    $0x3,%ecx
jl     1f6d <svdu1v+0xdfd>
mov    %r10,%r8
mov    %r11,%r10
mov    0x64(%rsp),%eax
mov    %eax,%ecx
sub    %r12d,%ecx
mov    0x50(%rsp),%rax
mov    %eax,%r14d
and    $0x3,%r14d
and    $0xfffffffc,%eax
mov    %eax,0xc0(%rsp)
mov    0x34(%rsp),%r11d
sub    %r12d,%r11d
xorpd  %xmm0,%xmm0
mov    (%rsp),%rax
mov    %rcx,0x8(%rsp)
cmp    $0x3,%ecx
mov    %r11d,0x80(%rsp)
jb     1864 <svdu1v+0x6f4>
xorpd  %xmm0,%xmm0
mov    0xc0(%rsp),%ecx
mov    (%rsp),%rax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rax),%xmm2
movsd  0x8(%rax),%xmm1
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0x18(%rax),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
add    $0x20,%rax
add    $0xfffffffc,%ecx
jne    1820 <svdu1v+0x6b0>
mov    0x80(%rsp),%r11d
test   $0x3,%r11b
je     1886 <svdu1v+0x716>
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rcx
cmp    %ecx,%r14d
jne    1870 <svdu1v+0x700>
xorpd  %xmm10,%xmm10
ucomisd %xmm10,%xmm0
xorpd  %xmm1,%xmm1
mov    %r10,%r11
mov    %r8,%r10
mov    0x78(%rsp),%r8
jbe    1f5a <svdu1v+0xdea>
mov    %rdi,0x38(%rsp)
ucomisd 0x1766(%rip),%xmm0
mov    %rsi,0x40(%rsp)
mov    %edx,0x58(%rsp)
movsd  %xmm0,0x68(%rsp)
jb     18c9 <svdu1v+0x759>
sqrtsd %xmm0,%xmm0
jmp    18d6 <svdu1v+0x766>
call   1080 <sqrt@plt>
movapd 0x172a(%rip),%xmm7
mov    0xb8(%rsp),%eax
sub    %r12d,%eax
mov    %rax,0x148(%rsp)
mov    (%rsp),%rax
movsd  (%rax),%xmm1
movapd %xmm0,%xmm2
xorpd  %xmm7,%xmm2
movapd %xmm1,%xmm3
cmpltsd 0x1713(%rip),%xmm3
andpd  %xmm3,%xmm2
andnpd %xmm0,%xmm3
orpd   %xmm2,%xmm3
movapd %xmm3,0x150(%rsp)
divsd  %xmm3,%xmm1
movapd %xmm1,%xmm0
call   1050 <fabs@plt>
movapd 0x150(%rsp),%xmm10
mov    %eax,0xbc(%rsp)
mov    (%rsp),%rax
movsd  (%rax),%xmm2
movapd %xmm2,%xmm8
addsd  %xmm10,%xmm8
movsd  %xmm8,(%rax)
mov    0x118(%rsp),%rax
mov    %r12,0x20(%rsp)
sub    %r12,%rax
imul   %r13,%rax
mov    %rax,0xb0(%rsp)
cmp    %r13,%rax
movsd  0x1699(%rip),%xmm9
mov    0x48(%rsp),%r10
mov    0x8(%rsp),%rax
jle    1dc0 <svdu1v+0xc50>
mov    0x60(%rsp),%esi
add    $0x1,%rsi
and    $0xfffffffffffffffc,%rsi
add    $0x1,%rax
add    $0xfffffffffffffffc,%rsi
mov    %rax,0x130(%rsp)
mov    %rax,%rdi
movabs $0x1fffffffc,%rax
and    %rax,%rdi
shr    $0x2,%rsi
lea    -0x4(%rdi),%rdx
mov    0x20(%rsp),%rax
mov    0xf8(%rsp),%r8
imul   %r8,%rax
add    $0x1,%rsi
mov    %rdx,0x128(%rsp)
mov    %rdx,%rcx
shr    $0x2,%rcx
mov    0xf0(%rsp),%rdx
add    %rax,%rdx
and    $0xfffffffffffffffe,%rsi
mov    %rsi,0x108(%rsp)
add    $0x1,%rcx
mov    %rcx,0x120(%rsp)
lea    (%r8,%rax,1),%rcx
mov    %rcx,0x98(%rsp)
add    0x8(%rsp),%rdx
mov    %rdx,0x140(%rsp)
mov    0xa0(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
add    $0x8,%rax
mov    0x8(%rsp),%rcx
mov    %rax,0x138(%rsp)
lea    0x8(%rax,%rcx,8),%rax
mov    %rax,0x90(%rsp)
mulsd  %xmm10,%xmm2
addsd  0x68(%rsp),%xmm2
movapd %xmm9,%xmm1
divsd  %xmm2,%xmm1
mov    (%rsp),%rax
mov    %rdi,0x68(%rsp)
lea    (%rax,%rdi,8),%rax
mov    %rax,0x100(%rsp)
mov    0x38(%rsp),%r9
xor    %r11d,%r11d
mov    %r13,%rax
jmp    1a98 <svdu1v+0x928>
nopl   0x0(%rax,%rax,1)
add    %r13,%rax
add    $0x1,%r11
add    %rbp,%r9
cmp    0xb0(%rsp),%rax
jge    1dc0 <svdu1v+0xc50>
cmp    $0x2,%r10d
jl     1a80 <svdu1v+0x910>
mov    (%rsp),%rdi
lea    (%rdi,%rax,8),%r8
xorpd  %xmm2,%xmm2
mov    %r8,%rsi
cmp    $0x3,%ecx
mov    0x80(%rsp),%edx
jb     1b26 <svdu1v+0x9b6>
xorpd  %xmm2,%xmm2
mov    0xc0(%rsp),%esi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movsd  (%r15,%rdi,1),%xmm0
movsd  0x8(%r15,%rdi,1),%xmm3
mulsd  (%r9,%rdi,1),%xmm0
mulsd  0x8(%r9,%rdi,1),%xmm3
addsd  %xmm2,%xmm0
movsd  0x10(%r15,%rdi,1),%xmm4
mulsd  0x10(%r9,%rdi,1),%xmm4
addsd  %xmm0,%xmm3
movsd  0x18(%r15,%rdi,1),%xmm2
mulsd  0x18(%r9,%rdi,1),%xmm2
addsd  %xmm3,%xmm4
addsd  %xmm4,%xmm2
add    $0x20,%rdi
add    $0xfffffffc,%esi
jne    1ad0 <svdu1v+0x960>
lea    (%r9,%rdi,1),%rsi
add    %r15,%rdi
test   $0x3,%dl
je     1b47 <svdu1v+0x9d7>
xor    %edx,%edx
nopl   (%rax)
movsd  (%rdi,%rdx,8),%xmm3
mulsd  (%rsi,%rdx,8),%xmm3
addsd  %xmm3,%xmm2
add    $0x1,%rdx
cmp    %edx,%r14d
jne    1b30 <svdu1v+0x9c0>
cmp    $0x2,%r10d
jl     1a80 <svdu1v+0x910>
xor    %edi,%edi
xorpd  0x14a5(%rip),%xmm2
mulsd  %xmm1,%xmm2
cmp    $0x3,%ecx
jb     1bab <svdu1v+0xa3b>
mov    %r11,%rdx
imul   %r13,%rdx
mov    0x98(%rsp),%rcx
lea    (%rcx,%rdx,1),%rsi
mov    0xa0(%rsp),%rcx
lea    (%rcx,%rsi,8),%rsi
cmp    0x90(%rsp),%rsi
jae    1c9e <svdu1v+0xb2e>
add    0x140(%rsp),%rdx
lea    (%rcx,%rdx,8),%rdx
cmp    %rdx,0x138(%rsp)
jae    1c9e <svdu1v+0xb2e>
mov    (%rsp),%r12
mov    0x20(%rsp),%rcx
lea    (%rdi,%rcx,1),%edx
mov    0x34(%rsp),%ecx
sub    %edx,%ecx
mov    0x64(%rsp),%esi
sub    %edx,%esi
test   $0x3,%cl
je     1c0b <svdu1v+0xa9b>
mov    %r14d,%ecx
sub    %dil,%cl
movzbl %cl,%r10d
and    $0x3,%r10d
neg    %r10d
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%r12),%xmm0
mulsd  %xmm2,%xmm0
addsd  (%r8),%xmm0
movsd  %xmm0,(%r8)
add    $0x8,%r8
add    $0x8,%r12
add    $0xffffffff,%edx
cmp    %edx,%r10d
jne    1be0 <svdu1v+0xa70>
sub    %edx,%edi
mov    0x48(%rsp),%r10
mov    0x8(%rsp),%rcx
cmp    $0x3,%esi
jb     1a80 <svdu1v+0x910>
mov    0x50(%rsp),%rsi
sub    %edi,%esi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r12,%rdi,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  (%r8,%rdi,8),%xmm0
movsd  %xmm0,(%r8,%rdi,8)
movsd  0x8(%r12,%rdi,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  0x8(%r8,%rdi,8),%xmm0
movsd  %xmm0,0x8(%r8,%rdi,8)
movsd  0x10(%r12,%rdi,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  0x10(%r8,%rdi,8),%xmm0
movsd  %xmm0,0x10(%r8,%rdi,8)
movsd  0x18(%r12,%rdi,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  0x18(%r8,%rdi,8),%xmm0
movsd  %xmm0,0x18(%r8,%rdi,8)
add    $0x4,%rdi
cmp    %edi,%esi
jne    1c30 <svdu1v+0xac0>
jmp    1a80 <svdu1v+0x910>
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
cmpq   $0x0,0x128(%rsp)
je     1dac <svdu1v+0xc3c>
mov    0x108(%rsp),%rsi
xor    %edi,%edi
nop
movupd (%r15,%rdi,8),%xmm4
movupd 0x10(%r15,%rdi,8),%xmm5
mulpd  %xmm3,%xmm4
movupd (%r9,%rdi,8),%xmm6
addpd  %xmm4,%xmm6
movupd 0x10(%r9,%rdi,8),%xmm4
movupd 0x20(%r9,%rdi,8),%xmm7
movupd 0x30(%r9,%rdi,8),%xmm0
mulpd  %xmm3,%xmm5
addpd  %xmm4,%xmm5
movupd %xmm6,(%r9,%rdi,8)
movupd %xmm5,0x10(%r9,%rdi,8)
movupd 0x20(%r15,%rdi,8),%xmm4
movupd 0x30(%r15,%rdi,8),%xmm5
mulpd  %xmm3,%xmm4
addpd  %xmm7,%xmm4
mulpd  %xmm3,%xmm5
addpd  %xmm0,%xmm5
movupd %xmm4,0x20(%r9,%rdi,8)
movupd %xmm5,0x30(%r9,%rdi,8)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1cc0 <svdu1v+0xb50>
testb  $0x1,0x120(%rsp)
01
je     1d7e <svdu1v+0xc0e>
mov    (%rsp),%rcx
movupd (%rcx,%rdi,8),%xmm0
movupd 0x10(%rcx,%rdi,8),%xmm4
movupd (%r8,%rdi,8),%xmm5
movupd 0x10(%r8,%rdi,8),%xmm6
mulpd  %xmm3,%xmm0
addpd  %xmm5,%xmm0
mulpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm0,(%r8,%rdi,8)
movupd %xmm4,0x10(%r8,%rdi,8)
mov    0x68(%rsp),%rcx
cmp    %rcx,0x130(%rsp)
mov    0x8(%rsp),%rcx
je     1a80 <svdu1v+0x910>
mov    0x68(%rsp),%rdi
lea    (%r8,%rdi,8),%r8
mov    0x100(%rsp),%r12
jmp    1baf <svdu1v+0xa3f>
xor    %edi,%edi
testb  $0x1,0x120(%rsp)
01
jne    1d45 <svdu1v+0xbd5>
jmp    1d7e <svdu1v+0xc0e>
nopw   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sdl 0xbc(%rsp),%xmm1
addsd  %xmm9,%xmm1
movapd %xmm9,%xmm2
divsd  %xmm8,%xmm2
mov    0x10(%rsp),%rcx
lea    0x10(%rcx),%rax
mov    $0x1,%edx
mov    0x148(%rsp),%r10
cmp    $0x3,%r10d
jae    1e25 <svdu1v+0xcb5>
mov    0x18(%rsp),%r11
movapd 0x11fc(%rip),%xmm7
mov    0x70(%rsp),%r10
mov    0x78(%rsp),%r8
mov    0x88(%rsp),%r9
mov    0x40(%rsp),%rsi
mov    0x20(%rsp),%r12
jmp    1f2d <svdu1v+0xdbd>
add    $0x1,%r10
mov    %r10,%r14
movabs $0x1fffffffc,%rdx
and    %rdx,%r14
lea    -0x4(%r14),%rsi
mov    %rsi,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
movapd %xmm2,%xmm0
unpcklpd %xmm2,%xmm0
test   %rsi,%rsi
mov    0x18(%rsp),%r11
movapd 0x11a0(%rip),%xmm7
mov    0x78(%rsp),%r8
mov    0x88(%rsp),%r9
mov    0x20(%rsp),%r12
je     1f89 <svdu1v+0xe19>
mov    0x5c(%rsp),%esi
add    $0x1,%rsi
and    $0xfffffffffffffffc,%rsi
add    $0xfffffffffffffffc,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movupd 0x10(%rcx,%rdi,8),%xmm3
movupd 0x20(%rcx,%rdi,8),%xmm4
movupd 0x30(%rcx,%rdi,8),%xmm5
movupd 0x40(%rcx,%rdi,8),%xmm6
mulpd  %xmm0,%xmm3
mulpd  %xmm0,%xmm4
movupd %xmm3,0x10(%rcx,%rdi,8)
movupd %xmm4,0x20(%rcx,%rdi,8)
mulpd  %xmm0,%xmm5
mulpd  %xmm0,%xmm6
movupd %xmm5,0x30(%rcx,%rdi,8)
movupd %xmm6,0x40(%rcx,%rdi,8)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1ea0 <svdu1v+0xd30>
test   $0x1,%dl
je     1f0d <svdu1v+0xd9d>
movupd (%rax,%rdi,8),%xmm3
movupd 0x10(%rax,%rdi,8),%xmm4
mulpd  %xmm0,%xmm3
mulpd  %xmm0,%xmm4
movupd %xmm3,(%rax,%rdi,8)
movupd %xmm4,0x10(%rax,%rdi,8)
cmp    %r14,%r10
mov    0x70(%rsp),%r10
mov    0x40(%rsp),%rsi
jne    1f23 <svdu1v+0xdb3>
mov    0x38(%rsp),%rdi
jmp    1f56 <svdu1v+0xde6>
mov    %r14d,%edx
or     $0x1,%edx
lea    (%rax,%r14,8),%rax
mov    0x50(%rsp),%rcx
sub    %edx,%ecx
xor    %edx,%edx
mov    0x38(%rsp),%rdi
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm0
mulsd  %xmm2,%xmm0
movsd  %xmm0,(%rax,%rdx,8)
add    $0x1,%rdx
cmp    %edx,%ecx
jne    1f40 <svdu1v+0xdd0>
mov    0x58(%rsp),%edx
mov    (%rsp),%rax
movsd  %xmm1,(%rax)
xorpd  %xmm7,%xmm10
movlpd %xmm10,(%r10,%r12,8)
cmp    $0x1,%edx
jne    12b0 <svdu1v+0x140>
mov    (%rsp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,(%r10,%r12,8)
jmp    12b0 <svdu1v+0x140>
xor    %edi,%edi
test   $0x1,%dl
jne    1eef <svdu1v+0xd7f>
jmp    1f0d <svdu1v+0xd9d>
mov    0xa0(%rsp),%rbp
mov    %rbp,%rdi
mov    0xe0(%rsp),%r15
mov    %r15,%rsi
mov    0xe8(%rsp),%r14
mov    %r14d,%edx
call   1060 <ldvmat@plt>
mov    %rbp,%rdi
mov    0xd8(%rsp),%r12
mov    %r12d,%esi
mov    %r14d,%edx
call   1070 <atou1@plt>
mov    0x18(%rsp),%rdi
mov    0x70(%rsp),%rsi
mov    %rbp,%rdx
mov    %r12d,%ecx
mov    %r15,%r8
mov    %r14d,%r9d
call   1090 <qrbdu1@plt>
movapd 0x1008(%rip),%xmm3
mov    0x18(%rsp),%r11
test   %r14d,%r14d
jle    213f <svdu1v+0xfcf>
mov    %r14d,%r9d
and    $0x3,%r9d
and    $0xfffffffc,%r14d
mov    %r13,%r8
shl    $0x5,%r8
lea    0x0(,%r13,8),%rdx
xor    %esi,%esi
xorpd  %xmm0,%xmm0
movapd 0xfd2(%rip),%xmm1
mov    %r15,%r10
jmp    2051 <svdu1v+0xee1>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rsi
add    $0x8,%r10
cmp    %r13,%rsi
je     213f <svdu1v+0xfcf>
movsd  (%r11,%rsi,8),%xmm2
ucomisd %xmm2,%xmm0
jbe    2040 <svdu1v+0xed0>
xorpd  %xmm3,%xmm2
movlpd %xmm2,(%r11,%rsi,8)
cmpl   $0x3,0x34(%rsp)
jae    2080 <svdu1v+0xf10>
lea    (%r15,%rsi,8),%rax
test   %r9d,%r9d
jne    20eb <svdu1v+0xf7b>
jmp    2040 <svdu1v+0xed0>
nopl   0x0(%rax)
mov    %r14d,%ecx
mov    %r10,%rbp
cs nopw 0x0(%rax,%rax,1)
mov    %rbp,%rax
movsd  0x0(%rbp),%xmm2
xorpd  %xmm1,%xmm2
movlpd %xmm2,0x0(%rbp)
add    %rdx,%rbp
movsd  (%rax,%rdx,1),%xmm2
xorpd  %xmm1,%xmm2
movlpd %xmm2,(%rax,%rdx,1)
lea    0x0(%rbp,%rdx,1),%rdi
movsd  (%rdx,%rbp,1),%xmm2
xorpd  %xmm1,%xmm2
movlpd %xmm2,(%rdx,%rbp,1)
lea    (%rdi,%rdx,1),%rbp
movsd  (%rdx,%rdi,1),%xmm2
xorpd  %xmm1,%xmm2
movlpd %xmm2,(%rdx,%rdi,1)
add    %rdx,%rbp
add    $0xfffffffc,%ecx
jne    2090 <svdu1v+0xf20>
add    %r8,%rax
test   %r9d,%r9d
je     2040 <svdu1v+0xed0>
mov    %r9d,%ecx
xchg   %ax,%ax
movsd  (%rax),%xmm2
xorpd  %xmm3,%xmm2
movlpd %xmm2,(%rax)
add    %rdx,%rax
add    $0xffffffff,%ecx
jne    20f0 <svdu1v+0xf80>
jmp    2040 <svdu1v+0xed0>
mov    %r12,%rdi
mov    %r10,%r13
mov    %r15,%rsi
mov    %r14d,%edx
call   1060 <ldvmat@plt>
mov    %r12,%rdi
mov    %ebp,%esi
mov    %r14d,%edx
call   1070 <atou1@plt>
mov    0x18(%rsp),%rdi
mov    %r13,%rsi
mov    %r12,%rdx
mov    %ebp,%ecx
mov    %r15,%r8
mov    %r14d,%r9d
call   1090 <qrbdu1@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
xor    %eax,%eax
add    $0x168,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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