#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double fabs(double);
double sign(int /*<<< orphan*/, double);
int /*<<< orphan*/ sqrt(double);
void orthes(int n, int low, int igh, double *a, double *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):
<orthes>:
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 <orthes+0x3c>
jmp    161e <orthes+0x4fe>
mov    -0x40(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     161e <orthes+0x4fe>
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 <orthes+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 <orthes+0x79>
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1207 <orthes+0xe7>
jp     1207 <orthes+0xe7>
jmp    1610 <orthes+0x4f0>
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 <orthes+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 <orthes+0xf6>
movsd  -0x58(%rbp),%xmm0
call   1040 <sqrt@plt>
mov    %eax,%edi
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <sign@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%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     1485 <orthes+0x365>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     13d1 <orthes+0x2b1>
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    136c <orthes+0x24c>
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     1472 <orthes+0x352>
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    13e6 <orthes+0x2c6>
jmp    1477 <orthes+0x357>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1352 <orthes+0x232>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     15bf <orthes+0x49f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     150b <orthes+0x3eb>
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    14a6 <orthes+0x386>
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     15ac <orthes+0x48c>
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    1520 <orthes+0x400>
jmp    15b1 <orthes+0x491>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    148c <orthes+0x36c>
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 <orthes+0x42>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<orthes>:
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    1961 <orthes+0x841>
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),%r10
mov    %r10,%rax
imul   %r8,%rax
lea    (%rax,%rdx,1),%r9
lea    (%rax,%rdx,1),%rbp
add    $0xffffffffffffffff,%rbp
mov    %rbp,0x130(%rsp)
lea    -0x1(%rsi),%ebp
imul   %r11d,%ebp
add    %esi,%ebp
mov    %r13,%rcx
shl    $0x20,%rcx
add    %rdi,%rcx
mov    %rcx,0xb0(%rsp)
mov    $0x1,%r11d
sub    %rdx,%r11
lea    (%r14,%r9,8),%rcx
add    $0x8,%rcx
mov    %rcx,0x90(%rsp)
lea    (%r15,%rsi,1),%ecx
mov    %ecx,0x2c(%rsp)
lea    -0x1(%r15,%rsi,1),%ecx
mov    %ecx,0x28(%rsp)
mov    %ebp,%ecx
mov    %r8,%rsi
shl    $0x5,%rsi
mov    %rsi,0x110(%rsp)
movapd 0xe0f(%rip),%xmm1
xorpd  %xmm8,%xmm8
lea    0x0(,%r8,8),%r12
lea    0x8(,%r8,8),%rdi
lea    0x1(%r8),%ebp
mov    %ebp,0x80(%rsp)
lea    0x8(%rbx,%rdx,8),%rsi
lea    (%r14,%rax,8),%rbp
mov    %rbp,0x88(%rsp)
mov    %rsi,%rbp
lea    (%r8,%r8,1),%r15d
mov    %rax,0xf8(%rsp)
lea    (%r14,%rax,8),%rax
add    $0x10,%rax
mov    %rax,0x98(%rsp)
mov    %r11,%rax
xor    %r9d,%r9d
mov    %r10,0x100(%rsp)
mov    %rdx,0x108(%rsp)
mov    %rdx,%r11
mov    %r8,0x70(%rsp)
mov    %rdi,0xf0(%rsp)
jmp    1316 <orthes+0x1f6>
nopl   0x0(%rax)
xorpd  %xmm9,%xmm0
movsd  -0x8(%rbx,%r11,8),%xmm1
movsd  0xe0(%rsp),%xmm2
mulsd  %xmm2,%xmm1
movsd  %xmm1,-0x8(%rbx,%r11,8)
mulsd  %xmm2,%xmm0
mov    0xd8(%rsp),%rax
add    0xe8(%rsp),%rax
movsd  %xmm0,(%r14,%rax,8)
movapd 0xd4e(%rip),%xmm1
add    $0x1,%r11
add    $0x1,%r9
mov    0x30(%rsp),%r10
add    $0x1,%r10
mov    0x84(%rsp),%ecx
add    0x80(%rsp),%ecx
mov    0x118(%rsp),%rax
add    $0xffffffffffffffff,%rax
add    $0x8,%rbp
mov    0xf0(%rsp),%rdi
add    %rdi,0x90(%rsp)
add    %r12,0x88(%rsp)
addl   $0x1,0x2c(%rsp)
addl   $0x1,0x28(%rsp)
add    %r12,0x98(%rsp)
cmp    0x40(%rsp),%r11d
je     1961 <orthes+0x841>
cmp    %r13,%r11
mov    %r13,%rdx
cmovg  %r11,%rdx
mov    %rax,0x118(%rsp)
add    %rdx,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0xc8(%rsp)
mov    0x108(%rsp),%rax
add    %r9,%rax
mov    %rdx,%rsi
mov    %rax,0x140(%rsp)
sub    %rax,%rsi
imul   %r9,%rdi
mov    %rdi,0x68(%rsp)
mov    0xf8(%rsp),%rax
sub    %r9,%rax
mov    %rdx,0xd0(%rsp)
add    %rdx,%rax
mov    %rax,0x60(%rsp)
mov    0x100(%rsp),%rax
add    %r9,%rax
mov    %rax,0x128(%rsp)
lea    (%rbx,%rax,8),%rax
mov    %rax,0xa8(%rsp)
lea    (%rbx,%r11,8),%rax
lea    (%rax,%rsi,8),%rax
mov    %rax,0xb8(%rsp)
add    $0x1,%rsi
mov    %rsi,0x18(%rsp)
lea    -0x1(%r11),%rdx
movq   $0x0,-0x8(%rbx,%r11,8)
mov    %ecx,0x84(%rsp)
mov    %ecx,%eax
mov    %r10,0x30(%rsp)
mov    %r10,%rcx
xorpd  %xmm2,%xmm2
xchg   %ax,%ax
cltq
movsd  (%r14,%rax,8),%xmm0
andpd  %xmm1,%xmm0
addsd  %xmm0,%xmm2
add    $0x1,%rcx
add    $0x1,%eax
cmp    %r13,%rcx
jl     13d0 <orthes+0x2b0>
ucomisd %xmm8,%xmm2
jne    13f9 <orthes+0x2d9>
jnp    12b2 <orthes+0x192>
mov    %rdx,0xe8(%rsp)
mov    %rbp,0xa0(%rsp)
mov    %r9,0xc0(%rsp)
lea    (%r11,%r13,1),%rax
mov    %rax,0x48(%rsp)
mov    %r11,0x78(%rsp)
lea    -0x2(%r11),%rsi
imul   %r8,%rsi
xorpd  %xmm0,%xmm0
mov    0x30(%rsp),%rax
mov    0xb0(%rsp),%rcx
movabs $0xffffffff00000000,%rdi
ff ff ff
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rcx,%rdx
sar    $0x20,%rdx
add    %rsi,%rdx
movsd  (%r14,%rdx,8),%xmm1
divsd  %xmm2,%xmm1
mov    %rcx,%rdx
sar    $0x1d,%rdx
movsd  %xmm1,(%rbx,%rdx,1)
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    %rdi,%rcx
add    $0x1,%rax
cmp    %r13,%rax
jl     1450 <orthes+0x330>
mov    %rsi,0xd8(%rsp)
movsd  %xmm2,0xe0(%rsp)
movsd  %xmm0,0x10(%rsp)
call   1040 <sqrt@plt>
mov    0x78(%rsp),%rbp
movsd  -0x8(%rbx,%rbp,8),%xmm0
mov    %eax,%edi
call   1030 <sign@plt>
mov    %rbp,%r11
movsd  -0x8(%rbx,%rbp,8),%xmm2
movapd %xmm2,%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm2
movsd  %xmm2,-0x8(%rbx,%rbp,8)
mov    0x70(%rsp),%r8
cmp    %r8,%rbp
movapd 0xb2c(%rip),%xmm9
mov    0x30(%rsp),%r9
jle    17a2 <orthes+0x682>
cmpl   $0x0,0x40(%rsp)
xorpd  %xmm8,%xmm8
mov    0xc0(%rsp),%r9
mov    0xa0(%rsp),%rbp
jle    1270 <orthes+0x150>
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    0x98(%rsp),%rax
mov    %rax,0x38(%rsp)
xor    %eax,%eax
mov    %rax,0x20(%rsp)
jmp    15b2 <orthes+0x492>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rdi
addq   $0x1,0x20(%rsp)
addq   $0x8,0x38(%rsp)
add    $0x8,%rdx
cmp    0x138(%rsp),%rdi
je     1270 <orthes+0x150>
mov    %rdx,0x10(%rsp)
lea    -0x1(%rdi),%r9
xorpd  %xmm2,%xmm2
mov    %r11,%rdx
testb  $0x1,0x18(%rsp)
je     15eb <orthes+0x4cb>
mov    0x60(%rsp),%rax
movsd  (%rbx,%rax,8),%xmm2
mov    0x58(%rsp),%rax
add    %r9,%rax
mulsd  (%r14,%rax,8),%xmm2
addsd  %xmm8,%xmm2
mov    0x50(%rsp),%rdx
mov    0x140(%rsp),%rax
cmp    %rax,0xd0(%rsp)
je     1679 <orthes+0x559>
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
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movslq %r10d,%r10
movsd  (%rbx,%r10,8),%xmm3
movslq %r11d,%r11
lea    (%r9,%r11,1),%r8
mulsd  (%r14,%r8,8),%xmm3
addsd  %xmm2,%xmm3
movslq %ecx,%rcx
movsd  (%rbx,%rcx,8),%xmm2
movslq %edx,%rdx
lea    (%r9,%rdx,1),%rsi
mulsd  (%r14,%rsi,8),%xmm2
addsd  %xmm3,%xmm2
sub    %r15d,%r11d
add    $0xfffffffe,%r10d
add    $0x2,%rax
sub    %r15d,%edx
add    $0xfffffffe,%ecx
cmp    %r13,%rax
jl     1630 <orthes+0x510>
xorpd  %xmm9,%xmm2
divsd  %xmm1,%xmm2
mov    0x78(%rsp),%r11
mov    %r11,%rax
cmpb   $0x0,0x68(%rsp)
mov    0x70(%rsp),%r8
mov    0xc0(%rsp),%r9
mov    0xc8(%rsp),%rsi
mov    0x10(%rsp),%rdx
je     1770 <orthes+0x650>
mov    0x128(%rsp),%rax
mov    0x20(%rsp),%rcx
add    %rcx,%rax
lea    (%r14,%rax,8),%rax
cmp    0xb8(%rsp),%rax
jae    16f2 <orthes+0x5d2>
mov    0xd0(%rsp),%rax
mov    0x20(%rsp),%rcx
add    %rcx,%rax
lea    (%r14,%rax,8),%rcx
mov    %r11,%rax
cmp    %rcx,0xa8(%rsp)
jb     1770 <orthes+0x650>
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
mov    0x38(%rsp),%rax
xor    %ecx,%ecx
mov    0x110(%rsp),%rdx
nopl   0x0(%rax)
movupd -0x10(%rax),%xmm4
movupd (%rax),%xmm5
movupd -0x10(%rbp,%rcx,8),%xmm6
movupd 0x0(%rbp,%rcx,8),%xmm7
mulpd  %xmm3,%xmm6
addpd  %xmm4,%xmm6
mulpd  %xmm3,%xmm7
addpd  %xmm5,%xmm7
movupd %xmm6,-0x10(%rax)
movupd %xmm7,(%rax)
add    $0x4,%rcx
add    %rdx,%rax
cmp    %rcx,%rsi
jne    1710 <orthes+0x5f0>
mov    0x120(%rsp),%rax
mov    0x48(%rsp),%rcx
cmp    %rcx,0x18(%rsp)
mov    0x10(%rsp),%rdx
je     1590 <orthes+0x470>
nopw   0x0(%rax,%rax,1)
add    $0xffffffffffffffff,%rax
mov    %r8,%rcx
imul   %rax,%rcx
lea    (%rdx,%rcx,8),%rcx
nop
movsd  (%rbx,%rax,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rcx),%xmm3
movsd  %xmm3,(%rcx)
add    $0x1,%rax
add    %r12,%rcx
cmp    %r13,%rax
jl     1780 <orthes+0x660>
jmp    1590 <orthes+0x470>
mov    0x18(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x58(%rsp)
add    %r11,%rax
mov    %rax,0x50(%rsp)
mov    0x88(%rsp),%rdx
mov    0x90(%rsp),%rsi
xor    %ecx,%ecx
mov    %r11,%rbp
jmp    17f4 <orthes+0x6d4>
nop
lea    0x1(%rbp),%rax
mov    0x10(%rsp),%rcx
add    $0x1,%rcx
add    %r12,%rsi
add    %r12,%rdx
cmp    %r8,%rbp
mov    %rax,%rbp
mov    0x30(%rsp),%r9
jge    14ef <orthes+0x3cf>
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,0x38(%rsp)
mov    %rbp,0x20(%rsp)
lea    -0x1(%rbp),%rdi
imul   %r8,%rdi
xorpd  %xmm2,%xmm2
mov    %r9,%r8
mov    0xb0(%rsp),%r11
movabs $0xffffffff00000000,%rax
ff ff ff
xchg   %ax,%ax
mov    %r11,%r9
sar    $0x20,%r9
mov    %r11,%r10
sar    $0x1d,%r10
movsd  (%rbx,%r10,1),%xmm3
add    %rdi,%r9
mulsd  (%r14,%r9,8),%xmm3
addsd  %xmm3,%xmm2
add    %rax,%r11
add    $0x1,%r8
cmp    %r13,%r8
jl     1850 <orthes+0x730>
xorpd  %xmm9,%xmm2
divsd  %xmm1,%xmm2
mov    0x78(%rsp),%r11
mov    %r11,%rdi
cmpq   $0x4,0x18(%rsp)
mov    0x70(%rsp),%r8
mov    0x20(%rsp),%rbp
jb     1930 <orthes+0x810>
cmp    0xb8(%rsp),%rcx
jae    18c0 <orthes+0x7a0>
mov    %r11,%rdi
mov    0x38(%rsp),%rax
cmp    0xa8(%rsp),%rax
ja     1930 <orthes+0x810>
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
xor    %eax,%eax
mov    0xa0(%rsp),%rcx
mov    0xc8(%rsp),%rdi
nopw   0x0(%rax,%rax,1)
movupd -0x10(%rsi,%rax,8),%xmm4
movupd (%rsi,%rax,8),%xmm5
movupd -0x10(%rcx,%rax,8),%xmm6
movupd (%rcx,%rax,8),%xmm7
mulpd  %xmm3,%xmm6
addpd  %xmm4,%xmm6
mulpd  %xmm3,%xmm7
addpd  %xmm5,%xmm7
movupd %xmm6,-0x10(%rsi,%rax,8)
movupd %xmm7,(%rsi,%rax,8)
add    $0x4,%rax
cmp    %rax,%rdi
jne    18e0 <orthes+0x7c0>
mov    0x50(%rsp),%rdi
mov    0x58(%rsp),%rax
cmp    %rax,0x18(%rsp)
je     17d0 <orthes+0x6b0>
nop
add    $0xffffffffffffffff,%rdi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx,%rdi,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rdx,%rdi,8),%xmm3
movsd  %xmm3,(%rdx,%rdi,8)
add    $0x1,%rdi
cmp    %r13,%rdi
jl     1940 <orthes+0x820>
jmp    17d0 <orthes+0x6b0>
add    $0x148,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @orthes(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 %350
31: ; preds = %5
%32 = load i32, i32* %18
store i32 %32, i32* %13
br label %33
33: ; preds = %347, %31
%34 = load i32, i32* %13
%35 = load i32, i32* %16
%36 = icmp sle i32 %34, %35
br i1 %36, label %37, label %350
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 %347
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 i32 @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 double @sign(i32 %120, double %126)
%128 = fneg double %127
store double %128, double* %20
%129 = load double, double* %21
%130 = load double*, double** %10
%131 = load i32, i32* %13
%132 = sub nsw i32 %131, 1
%133 = sext i32 %132 to i64
%134 = getelementptr inbounds double, double* %130, i64 %133
%135 = load double, double* %134
%136 = load double, double* %20
%137 = fneg double %135
%138 = call double @llvm.fmuladd.f64(double %137, double %136, double %129)
store double %138, double* %21
%139 = load double*, double** %10
%140 = load i32, i32* %13
%141 = sub nsw i32 %140, 1
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %139, i64 %142
%144 = load double, double* %143
%145 = load double, double* %20
%146 = fsub double %144, %145
%147 = load double*, double** %10
%148 = load i32, i32* %13
%149 = sub nsw i32 %148, 1
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %147, i64 %150
store double %146, double* %151
%152 = load i32, i32* %13
store i32 %152, i32* %12
br label %153
153: ; preds = %233, %118
%154 = load i32, i32* %12
%155 = load i32, i32* %6
%156 = icmp sle i32 %154, %155
br i1 %156, label %157, label %236
157: ; preds = %153
store double 0.000000e+00, double* %19
%158 = load i32, i32* %13
store i32 %158, i32* %14
br label %159
159: ; preds = %186, %157
%160 = load i32, i32* %14
%161 = load i32, i32* %8
%162 = icmp sle i32 %160, %161
br i1 %162, label %163, label %189
163: ; preds = %159
%164 = load i32, i32* %17
%165 = load i32, i32* %14
%166 = sub nsw i32 %164, %165
store i32 %166, i32* %11
%167 = load double, double* %19
%168 = load double*, double** %10
%169 = load i32, i32* %11
%170 = sub nsw i32 %169, 1
%171 = sext i32 %170 to i64
%172 = getelementptr inbounds double, double* %168, i64 %171
%173 = load double, double* %172
%174 = load double*, double** %9
%175 = load i32, i32* %11
%176 = sub nsw i32 %175, 1
%177 = load i32, i32* %12
%178 = sub nsw i32 %177, 1
%179 = load i32, i32* %6
%180 = mul nsw i32 %178, %179
%181 = add nsw i32 %176, %180
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %174, i64 %182
%184 = load double, double* %183
%185 = call double @llvm.fmuladd.f64(double %173, double %184, double %167)
store double %185, double* %19
br label %186
186: ; preds = %163
%187 = load i32, i32* %14
%188 = add nsw i32 %187, 1
store i32 %188, i32* %14
br label %159
189: ; preds = %159
%190 = load double, double* %19
%191 = load double, double* %21
%192 = fdiv double %190, %191
store double %192, double* %19
%193 = load i32, i32* %13
store i32 %193, i32* %11
br label %194
194: ; preds = %229, %189
%195 = load i32, i32* %11
%196 = load i32, i32* %8
%197 = icmp sle i32 %195, %196
br i1 %197, label %198, label %232
198: ; preds = %194
%199 = load double*, double** %9
%200 = load i32, i32* %11
%201 = sub nsw i32 %200, 1
%202 = load i32, i32* %12
%203 = sub nsw i32 %202, 1
%204 = load i32, i32* %6
%205 = mul nsw i32 %203, %204
%206 = add nsw i32 %201, %205
%207 = sext i32 %206 to i64
%208 = getelementptr inbounds double, double* %199, i64 %207
%209 = load double, double* %208
%210 = load double, double* %19
%211 = load double*, double** %10
%212 = load i32, i32* %11
%213 = sub nsw i32 %212, 1
%214 = sext i32 %213 to i64
%215 = getelementptr inbounds double, double* %211, i64 %214
%216 = load double, double* %215
%217 = fneg double %210
%218 = call double @llvm.fmuladd.f64(double %217, double %216, double %209)
%219 = load double*, double** %9
%220 = load i32, i32* %11
%221 = sub nsw i32 %220, 1
%222 = load i32, i32* %12
%223 = sub nsw i32 %222, 1
%224 = load i32, i32* %6
%225 = mul nsw i32 %223, %224
%226 = add nsw i32 %221, %225
%227 = sext i32 %226 to i64
%228 = getelementptr inbounds double, double* %219, i64 %227
store double %218, double* %228
br label %229
229: ; preds = %198
%230 = load i32, i32* %11
%231 = add nsw i32 %230, 1
store i32 %231, i32* %11
br label %194
232: ; preds = %194
br label %233
233: ; preds = %232
%234 = load i32, i32* %12
%235 = add nsw i32 %234, 1
store i32 %235, i32* %12
br label %153
236: ; preds = %153
store i32 1, i32* %11
br label %237
237: ; preds = %317, %236
%238 = load i32, i32* %11
%239 = load i32, i32* %8
%240 = icmp sle i32 %238, %239
br i1 %240, label %241, label %320
241: ; preds = %237
store double 0.000000e+00, double* %19
%242 = load i32, i32* %13
store i32 %242, i32* %15
br label %243
243: ; preds = %270, %241
%244 = load i32, i32* %15
%245 = load i32, i32* %8
%246 = icmp sle i32 %244, %245
br i1 %246, label %247, label %273
247: ; preds = %243
%248 = load i32, i32* %17
%249 = load i32, i32* %15
%250 = sub nsw i32 %248, %249
store i32 %250, i32* %12
%251 = load double, double* %19
%252 = load double*, double** %10
%253 = load i32, i32* %12
%254 = sub nsw i32 %253, 1
%255 = sext i32 %254 to i64
%256 = getelementptr inbounds double, double* %252, i64 %255
%257 = load double, double* %256
%258 = load double*, double** %9
%259 = load i32, i32* %11
%260 = sub nsw i32 %259, 1
%261 = load i32, i32* %12
%262 = sub nsw i32 %261, 1
%263 = load i32, i32* %6
%264 = mul nsw i32 %262, %263
%265 = add nsw i32 %260, %264
%266 = sext i32 %265 to i64
%267 = getelementptr inbounds double, double* %258, i64 %266
%268 = load double, double* %267
%269 = call double @llvm.fmuladd.f64(double %257, double %268, double %251)
store double %269, double* %19
br label %270
270: ; preds = %247
%271 = load i32, i32* %15
%272 = add nsw i32 %271, 1
store i32 %272, i32* %15
br label %243
273: ; preds = %243
%274 = load double, double* %19
%275 = load double, double* %21
%276 = fdiv double %274, %275
store double %276, double* %19
%277 = load i32, i32* %13
store i32 %277, i32* %12
br label %278
278: ; preds = %313, %273
%279 = load i32, i32* %12
%280 = load i32, i32* %8
%281 = icmp sle i32 %279, %280
br i1 %281, label %282, label %316
282: ; preds = %278
%283 = load double*, double** %9
%284 = load i32, i32* %11
%285 = sub nsw i32 %284, 1
%286 = load i32, i32* %12
%287 = sub nsw i32 %286, 1
%288 = load i32, i32* %6
%289 = mul nsw i32 %287, %288
%290 = add nsw i32 %285, %289
%291 = sext i32 %290 to i64
%292 = getelementptr inbounds double, double* %283, i64 %291
%293 = load double, double* %292
%294 = load double, double* %19
%295 = load double*, double** %10
%296 = load i32, i32* %12
%297 = sub nsw i32 %296, 1
%298 = sext i32 %297 to i64
%299 = getelementptr inbounds double, double* %295, i64 %298
%300 = load double, double* %299
%301 = fneg double %294
%302 = call double @llvm.fmuladd.f64(double %301, double %300, double %293)
%303 = load double*, double** %9
%304 = load i32, i32* %11
%305 = sub nsw i32 %304, 1
%306 = load i32, i32* %12
%307 = sub nsw i32 %306, 1
%308 = load i32, i32* %6
%309 = mul nsw i32 %307, %308
%310 = add nsw i32 %305, %309
%311 = sext i32 %310 to i64
%312 = getelementptr inbounds double, double* %303, i64 %311
store double %302, double* %312
br label %313
313: ; preds = %282
%314 = load i32, i32* %12
%315 = add nsw i32 %314, 1
store i32 %315, i32* %12
br label %278
316: ; preds = %278
br label %317
317: ; preds = %316
%318 = load i32, i32* %11
%319 = add nsw i32 %318, 1
store i32 %319, i32* %11
br label %237
320: ; preds = %237
%321 = load double, double* %22
%322 = load double*, double** %10
%323 = load i32, i32* %13
%324 = sub nsw i32 %323, 1
%325 = sext i32 %324 to i64
%326 = getelementptr inbounds double, double* %322, i64 %325
%327 = load double, double* %326
%328 = fmul double %321, %327
%329 = load double*, double** %10
%330 = load i32, i32* %13
%331 = sub nsw i32 %330, 1
%332 = sext i32 %331 to i64
%333 = getelementptr inbounds double, double* %329, i64 %332
store double %328, double* %333
%334 = load double, double* %22
%335 = load double, double* %20
%336 = fmul double %334, %335
%337 = load double*, double** %9
%338 = load i32, i32* %13
%339 = sub nsw i32 %338, 1
%340 = load i32, i32* %13
%341 = sub nsw i32 %340, 2
%342 = load i32, i32* %6
%343 = mul nsw i32 %341, %342
%344 = add nsw i32 %339, %343
%345 = sext i32 %344 to i64
%346 = getelementptr inbounds double, double* %337, i64 %345
store double %336, double* %346
br label %347
347: ; preds = %320, %69
%348 = load i32, i32* %13
%349 = add nsw i32 %348, 1
store i32 %349, i32* %13
br label %33
350: ; preds = %30, %33
ret void
}

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