void orthesx(n, low, igh, a, ort) int n, low, igh;
double *a, *ort;
{
  int i, j, m, ii, jj, la, mp, kp1;
  double f, g, h, scale;
  la = igh - 1;
  kp1 = low + 1;
  if (la < kp1)
    return;
  for (m = kp1; m <= la; m++) {
    h = 0.0;
    ort[m - 1] = 0.0;
    scale = 0.0;
    for (i = m; i <= igh; i++)
      scale = scale + fabs(a[i - 1 + (m - 2) * n]);
    if (scale == 0.0)
      continue;
    mp = m + igh;
    for (ii = m; ii <= igh; ii++) {
      i = mp - ii;
      ort[i - 1] = a[i - 1 + (m - 2) * n] / scale;
      h = h + ort[i - 1] * ort[i - 1];
    }
    g = -sign(sqrt(h), ort[m - 1]);
    h = h - ort[m - 1] * g;
    ort[m - 1] = ort[m - 1] - g;
    for (j = m; j <= n; j++) {
      f = 0.0;
      for (ii = m; ii <= igh; ii++) {
        i = mp - ii;
        f = f + ort[i - 1] * a[i - 1 + (j - 1) * n];
      }
      f = f / h;
      for (i = m; i <= igh; i++)
        a[i - 1 + (j - 1) * n] = a[i - 1 + (j - 1) * n] - f * ort[i - 1];
    }
    for (i = 1; i <= igh; i++) {
      f = 0.0;
      for (jj = m; jj <= igh; jj++) {
        j = mp - jj;
        f = f + ort[j - 1] * a[i - 1 + (j - 1) * n];
      }
      f = f / h;
      for (j = m; j <= igh; j++)
        a[i - 1 + (j - 1) * n] = a[i - 1 + (j - 1) * n] - f * ort[j - 1];
    }
    ort[m - 1] = scale * ort[m - 1];
    a[m - 1 + (m - 2) * n] = scale * g;
  }
}


x86 assembly (clang, -O0):
<orthesx>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    115c <orthesx+0x3c>
jmp    1611 <orthesx+0x4f1>
mov    -0x40(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     1611 <orthesx+0x4f1>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11ea <orthesx+0xca>
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xe31(%rip),%xmm2
pand   %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1199 <orthesx+0x79>
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1207 <orthesx+0xe7>
jp     1207 <orthesx+0xe7>
jmp    1603 <orthesx+0x4e3>
mov    -0x2c(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12a4 <orthesx+0x184>
mov    -0x3c(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x60(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1216 <orthesx+0xf6>
movsd  -0x58(%rbp),%xmm0
call   1040 <sqrt@plt>
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    $0x2,%al
call   1030 <sign@plt>
mov    %eax,%ecx
xor    %eax,%eax
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1478 <orthesx+0x358>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     13c4 <orthesx+0x2a4>
mov    -0x3c(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x24(%rbp)
movsd  -0x48(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    135f <orthesx+0x23f>
movsd  -0x48(%rbp),%xmm0
divsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1465 <orthesx+0x345>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    13d9 <orthesx+0x2b9>
jmp    146a <orthesx+0x34a>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1345 <orthesx+0x225>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     15b2 <orthesx+0x492>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     14fe <orthesx+0x3de>
mov    -0x3c(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x28(%rbp)
movsd  -0x48(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1499 <orthesx+0x379>
movsd  -0x48(%rbp),%xmm0
divsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     159f <orthesx+0x47f>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
sub    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1513 <orthesx+0x3f3>
jmp    15a4 <orthesx+0x484>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    147f <orthesx+0x35f>
movsd  -0x60(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%edx
sub    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1162 <orthesx+0x42>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<orthesx>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x148,%rsp
mov    %edi,0xc(%rsp)
lea    0x1(%rsi),%eax
mov    %rdx,0x40(%rsp)
cmp    %edx,%eax
jge    1981 <orthesx+0x861>
mov    %r8,%rbx
mov    %rcx,%r14
movabs $0xffffffff00000000,%rdi
ff ff ff
movslq %eax,%rdx
mov    0x40(%rsp),%r15
movslq %r15d,%r13
mov    0xc(%rsp),%r11d
movslq %r11d,%r8
lea    0x1(%r15),%eax
mov    %rax,0x138(%rsp)
lea    -0x1(%rdx),%rax
mov    %rax,0x88(%rsp)
imul   %r8,%rax
lea    (%rax,%rdx,1),%r10
lea    (%rax,%rdx,1),%rbp
add    $0xffffffffffffffff,%rbp
mov    %rbp,0x130(%rsp)
lea    -0x1(%rsi),%r9d
imul   %r11d,%r9d
add    %esi,%r9d
mov    %r13,%rcx
shl    $0x20,%rcx
add    %rdi,%rcx
mov    %rcx,0xc0(%rsp)
mov    $0x1,%r11d
sub    %rdx,%r11
lea    (%r14,%r10,8),%rcx
add    $0x8,%rcx
mov    %rcx,0x98(%rsp)
lea    (%r15,%rsi,1),%ecx
mov    %ecx,0x2c(%rsp)
lea    (%r15,%rsi,1),%ecx
add    $0xffffffff,%ecx
mov    %ecx,0x28(%rsp)
mov    %r8,%rsi
shl    $0x5,%rsi
mov    %rsi,0x110(%rsp)
movapd 0xe08(%rip),%xmm1
xorpd  %xmm9,%xmm9
lea    0x0(,%r8,8),%r12
lea    0x8(,%r8,8),%rdi
lea    0x1(%r8),%ebp
mov    %ebp,0x80(%rsp)
lea    0x8(%rbx,%rdx,8),%rcx
lea    (%r14,%rax,8),%rbp
mov    %rbp,0x90(%rsp)
mov    %rcx,%rbp
lea    (%r8,%r8,1),%r15d
mov    %rax,0x100(%rsp)
lea    (%r14,%rax,8),%rax
add    $0x10,%rax
mov    %rax,0xa0(%rsp)
mov    %r11,%rax
xor    %r10d,%r10d
mov    0x88(%rsp),%rcx
mov    %rdx,%r11
mov    %rdx,0x108(%rsp)
mov    %r8,0x78(%rsp)
mov    %rdi,0xf8(%rsp)
jmp    1324 <orthesx+0x204>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%rbx,%r11,8),%xmm1
movsd  0xe8(%rsp),%xmm0
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%rbx,%r11,8)
mulsd  %xmm0,%xmm8
mov    0xe0(%rsp),%rax
add    0xf0(%rsp),%rax
movsd  %xmm8,(%r14,%rax,8)
movapd 0xd42(%rip),%xmm1
add    $0x1,%r11
add    $0x1,%r10
mov    0x70(%rsp),%rcx
add    $0x1,%rcx
mov    0x84(%rsp),%r9d
add    0x80(%rsp),%r9d
mov    0x118(%rsp),%rax
add    $0xffffffffffffffff,%rax
add    $0x8,%rbp
mov    0xf8(%rsp),%rdi
add    %rdi,0x98(%rsp)
add    %r12,0x90(%rsp)
addl   $0x1,0x2c(%rsp)
addl   $0x1,0x28(%rsp)
add    %r12,0xa0(%rsp)
cmp    0x40(%rsp),%r11d
je     1981 <orthesx+0x861>
cmp    %r13,%r11
mov    %r13,%rdx
cmovg  %r11,%rdx
mov    %rax,0x118(%rsp)
add    %rdx,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0xb0(%rsp)
mov    0x108(%rsp),%rax
add    %r10,%rax
mov    %rdx,%rsi
mov    %rax,0x140(%rsp)
sub    %rax,%rsi
imul   %r10,%rdi
mov    %rdi,0x68(%rsp)
mov    0x100(%rsp),%rax
sub    %r10,%rax
mov    %rdx,0xd8(%rsp)
add    %rdx,%rax
mov    %rax,0x60(%rsp)
mov    0x88(%rsp),%rax
add    %r10,%rax
mov    %rax,0x128(%rsp)
lea    (%rbx,%rax,8),%rax
mov    %rax,0xb8(%rsp)
lea    (%rbx,%r11,8),%rax
lea    (%rax,%rsi,8),%rax
mov    %rax,0xc8(%rsp)
add    $0x1,%rsi
mov    %rsi,0x18(%rsp)
lea    -0x1(%r11),%rdx
mov    %r11,0x20(%rsp)
movq   $0x0,-0x8(%rbx,%r11,8)
mov    %r9d,0x84(%rsp)
mov    %r9d,%eax
mov    %rcx,0x70(%rsp)
xorpd  %xmm2,%xmm2
cltq
movsd  (%r14,%rax,8),%xmm0
andpd  %xmm1,%xmm0
addsd  %xmm0,%xmm2
add    $0x1,%rcx
add    $0x1,%eax
cmp    %r13,%rcx
jl     13e0 <orthesx+0x2c0>
ucomisd %xmm9,%xmm2
mov    0x20(%rsp),%r11
jne    140e <orthesx+0x2ee>
jnp    12be <orthesx+0x19e>
mov    %rdx,0xf0(%rsp)
lea    (%r11,%r13,1),%r9
lea    -0x2(%r11),%rsi
imul   %r8,%rsi
xorpd  %xmm1,%xmm1
mov    0x70(%rsp),%rax
mov    0xc0(%rsp),%rcx
movabs $0xffffffff00000000,%rdi
ff ff ff
nopl   (%rax)
mov    %rcx,%rdx
sar    $0x20,%rdx
add    %rsi,%rdx
movsd  (%r14,%rdx,8),%xmm0
divsd  %xmm2,%xmm0
mov    %rcx,%rdx
sar    $0x1d,%rdx
movsd  %xmm0,(%rbx,%rdx,1)
mulsd  %xmm0,%xmm0
addsd  %xmm0,%xmm1
add    %rdi,%rcx
add    $0x1,%rax
cmp    %r13,%rax
jl     1440 <orthesx+0x320>
ucomisd %xmm9,%xmm1
mov    %rbp,0xa8(%rsp)
mov    %r10,0xd0(%rsp)
movsd  %xmm2,0xe8(%rsp)
mov    %rsi,0xe0(%rsp)
mov    %r9,0x48(%rsp)
movsd  %xmm1,0x10(%rsp)
jb     14b0 <orthesx+0x390>
xorps  %xmm0,%xmm0
sqrtsd %xmm1,%xmm0
jmp    14be <orthesx+0x39e>
movapd %xmm1,%xmm0
call   1040 <sqrt@plt>
mov    0x20(%rsp),%r11
movsd  -0x8(%rbx,%r11,8),%xmm1
mov    $0x2,%al
call   1030 <sign@plt>
mov    0x20(%rsp),%r11
neg    %eax
cvtsi2sd %eax,%xmm8
movsd  -0x8(%rbx,%r11,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm8,%xmm2
movsd  0x10(%rsp),%xmm7
subsd  %xmm2,%xmm7
subsd  %xmm8,%xmm1
movsd  %xmm1,-0x8(%rbx,%r11,8)
mov    0x78(%rsp),%r8
cmp    %r8,%r11
xorpd  %xmm9,%xmm9
movapd 0xafd(%rip),%xmm0
jle    17c2 <orthesx+0x6a2>
cmpl   $0x0,0x40(%rsp)
mov    0xd0(%rsp),%r10
mov    0xa8(%rsp),%rbp
jle    1280 <orthesx+0x160>
mov    0xc(%rsp),%edx
cmp    $0x1,%edx
sete   %al
mov    %r11d,%ecx
not    %ecx
add    0x48(%rsp),%ecx
movslq %ecx,%rsi
mov    %rsi,%rcx
mov    %rsi,0x60(%rsp)
imul   %edx,%ecx
movslq %ecx,%rcx
mov    %rcx,0x58(%rsp)
lea    0x1(%r11),%rcx
mov    %rcx,0x50(%rsp)
mov    0x18(%rsp),%rcx
cmp    $0x4,%rcx
setae  %dl
and    %al,%dl
mov    %dl,0x68(%rsp)
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x48(%rsp)
lea    (%r11,%rcx,1),%rax
mov    %rax,0x120(%rsp)
mov    $0x1,%edi
mov    %r14,%rdx
mov    0xa0(%rsp),%rax
mov    %rax,0x30(%rsp)
xor    %esi,%esi
jmp    15d0 <orthesx+0x4b0>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rdi
add    $0x1,%rsi
addq   $0x8,0x30(%rsp)
add    $0x8,%rdx
cmp    0x138(%rsp),%rdi
je     1280 <orthesx+0x160>
mov    %rsi,0x38(%rsp)
mov    %rdx,0x10(%rsp)
lea    -0x1(%rdi),%r9
xorpd  %xmm1,%xmm1
mov    %r11,%rdx
testb  $0x1,0x18(%rsp)
je     160e <orthesx+0x4ee>
mov    0x60(%rsp),%rax
movsd  (%rbx,%rax,8),%xmm1
mov    0x58(%rsp),%rax
add    %r9,%rax
mulsd  (%r14,%rax,8),%xmm1
addsd  %xmm9,%xmm1
mov    0x50(%rsp),%rdx
mov    0x140(%rsp),%rax
cmp    %rax,0xd8(%rsp)
je     1699 <orthesx+0x579>
mov    0x2c(%rsp),%eax
mov    %eax,%r10d
sub    %edx,%r10d
mov    0xc(%rsp),%esi
mov    %esi,%r11d
imul   %r10d,%r11d
mov    %rdx,%rax
add    $0xffffffffffffffff,%rax
mov    0x28(%rsp),%ecx
sub    %edx,%ecx
mov    %esi,%edx
imul   %ecx,%edx
nopw   0x0(%rax,%rax,1)
movslq %r10d,%r10
movsd  (%rbx,%r10,8),%xmm2
movslq %r11d,%r11
lea    (%r9,%r11,1),%r8
mulsd  (%r14,%r8,8),%xmm2
addsd  %xmm1,%xmm2
movslq %ecx,%rcx
movsd  (%rbx,%rcx,8),%xmm1
movslq %edx,%rdx
lea    (%r9,%rdx,1),%rsi
mulsd  (%r14,%rsi,8),%xmm1
addsd  %xmm2,%xmm1
sub    %r15d,%r11d
add    $0xfffffffe,%r10d
add    $0x2,%rax
sub    %r15d,%edx
add    $0xfffffffe,%ecx
cmp    %r13,%rax
jl     1650 <orthesx+0x530>
xorpd  %xmm0,%xmm1
divsd  %xmm7,%xmm1
mov    0x20(%rsp),%r11
mov    %r11,%rax
cmpb   $0x0,0x68(%rsp)
mov    0x78(%rsp),%r8
mov    0xd0(%rsp),%r10
mov    0x10(%rsp),%rdx
mov    0x38(%rsp),%rsi
je     1790 <orthesx+0x670>
mov    0x128(%rsp),%rax
add    %rsi,%rax
lea    (%r14,%rax,8),%rax
cmp    0xc8(%rsp),%rax
jae    1704 <orthesx+0x5e4>
mov    0xd8(%rsp),%rax
add    %rsi,%rax
lea    (%r14,%rax,8),%rcx
mov    %r11,%rax
cmp    %rcx,0xb8(%rsp)
jb     1790 <orthesx+0x670>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
mov    0x30(%rsp),%rax
xor    %ecx,%ecx
mov    0x110(%rsp),%rdx
mov    0xb0(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x10(%rax),%xmm3
movupd (%rax),%xmm4
movupd -0x10(%rbp,%rcx,8),%xmm5
movupd 0x0(%rbp,%rcx,8),%xmm6
mulpd  %xmm2,%xmm5
addpd  %xmm3,%xmm5
mulpd  %xmm2,%xmm6
addpd  %xmm4,%xmm6
movupd %xmm5,-0x10(%rax)
movupd %xmm6,(%rax)
add    $0x4,%rcx
add    %rdx,%rax
cmp    %rcx,%rsi
jne    1730 <orthesx+0x610>
mov    0x120(%rsp),%rax
mov    0x48(%rsp),%rcx
cmp    %rcx,0x18(%rsp)
mov    0x10(%rsp),%rdx
mov    0x38(%rsp),%rsi
je     15b0 <orthesx+0x490>
nopl   0x0(%rax)
add    $0xffffffffffffffff,%rax
mov    %r8,%rcx
imul   %rax,%rcx
lea    (%rdx,%rcx,8),%rcx
nop
movsd  (%rbx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rcx),%xmm2
movsd  %xmm2,(%rcx)
add    $0x1,%rax
add    %r12,%rcx
cmp    %r13,%rax
jl     17a0 <orthesx+0x680>
jmp    15b0 <orthesx+0x490>
mov    0x18(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x58(%rsp)
add    %r11,%rax
mov    %rax,0x50(%rsp)
mov    0x90(%rsp),%rdx
mov    0x98(%rsp),%rsi
xor    %ecx,%ecx
mov    %r11,%rbp
jmp    180f <orthesx+0x6ef>
nop
lea    0x1(%rbp),%rax
mov    0x10(%rsp),%rcx
add    $0x1,%rcx
add    %r12,%rsi
add    %r12,%rdx
cmp    %r8,%rbp
mov    %rax,%rbp
jge    1519 <orthesx+0x3f9>
mov    %rcx,0x10(%rsp)
mov    %rcx,%rax
imul   %r8,%rax
mov    0x130(%rsp),%rcx
add    %rax,%rcx
lea    (%r14,%rcx,8),%rcx
mov    0x68(%rsp),%rdi
add    %rdi,%rcx
add    0x60(%rsp),%rax
lea    (%r14,%rax,8),%rax
add    %rdi,%rax
mov    %rax,0x30(%rsp)
mov    %rbp,0x38(%rsp)
lea    -0x1(%rbp),%rdi
imul   %r8,%rdi
xorpd  %xmm1,%xmm1
mov    0x70(%rsp),%r8
mov    0xc0(%rsp),%r11
movabs $0xffffffff00000000,%rax
ff ff ff
nopl   0x0(%rax,%rax,1)
mov    %r11,%r9
sar    $0x20,%r9
mov    %r11,%r10
sar    $0x1d,%r10
movsd  (%rbx,%r10,1),%xmm2
add    %rdi,%r9
mulsd  (%r14,%r9,8),%xmm2
addsd  %xmm2,%xmm1
add    %rax,%r11
add    $0x1,%r8
cmp    %r13,%r8
jl     1870 <orthesx+0x750>
xorpd  %xmm0,%xmm1
divsd  %xmm7,%xmm1
mov    0x20(%rsp),%r11
mov    %r11,%rdi
cmpq   $0x4,0x18(%rsp)
mov    0x78(%rsp),%r8
mov    0x38(%rsp),%rbp
jb     1950 <orthesx+0x830>
cmp    0xc8(%rsp),%rcx
jae    18df <orthesx+0x7bf>
mov    %r11,%rdi
mov    0x30(%rsp),%rax
cmp    0xb8(%rsp),%rax
ja     1950 <orthesx+0x830>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
xor    %eax,%eax
mov    0xb0(%rsp),%rcx
mov    0xa8(%rsp),%rdi
nopl   0x0(%rax)
movupd -0x10(%rsi,%rax,8),%xmm3
movupd (%rsi,%rax,8),%xmm4
movupd -0x10(%rdi,%rax,8),%xmm5
movupd (%rdi,%rax,8),%xmm6
mulpd  %xmm2,%xmm5
addpd  %xmm3,%xmm5
mulpd  %xmm2,%xmm6
addpd  %xmm4,%xmm6
movupd %xmm5,-0x10(%rsi,%rax,8)
movupd %xmm6,(%rsi,%rax,8)
add    $0x4,%rax
cmp    %rax,%rcx
jne    1900 <orthesx+0x7e0>
mov    0x50(%rsp),%rdi
mov    0x58(%rsp),%rax
cmp    %rax,0x18(%rsp)
je     17f0 <orthesx+0x6d0>
nop
add    $0xffffffffffffffff,%rdi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx,%rdi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rdx,%rdi,8),%xmm2
movsd  %xmm2,(%rdx,%rdi,8)
add    $0x1,%rdi
cmp    %r13,%rdi
jl     1960 <orthesx+0x840>
jmp    17f0 <orthesx+0x6d0>
add    $0x148,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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