#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 sqrt(double);
int cholbnd_upper(double *Abnd, double *b, int n, int bw) {
  int i, j, k;
  double sum;
  for (i = 0; i < n; ++i) {
    for (j = ((0 > i - bw) ? 0 : i - bw); j < i; ++j) {
      sum = Abnd[i + n * (bw + j - i)];
      for (k = ((0 > i - bw) ? 0 : i - bw); k < j; ++k)
        sum -= Abnd[i + n * (bw + k - i)] * Abnd[j + n * (bw + k - j)];
      Abnd[i + n * (bw + j - i)] = sum / Abnd[j + n * bw];
    }
    sum = Abnd[i + n * bw];
    for (k = ((0 > i - bw) ? 0 : i - bw); k < i; ++k)
      sum -= Abnd[i + n * (bw + k - i)] * Abnd[i + n * (bw + k - i)];
    if (sum >= 0.0)
      Abnd[i + n * bw] = sqrt(sum);
    else
      return 1;
  }
  b[0] = b[0] / Abnd[0 + n * bw];
  for (i = 1; i < n; ++i) {
    sum = 0.;
    for (k = ((0 > i - bw) ? 0 : i - bw); k < i; ++k)
      sum += Abnd[i + n * (bw + k - i)] * b[k];
    b[i] = (b[i] - sum) / Abnd[i + n * bw];
  }
  b[n - 1] = b[n - 1] / Abnd[n - 1 + n * bw];
  for (i = n - 2; i >= 0; --i) {
    sum = 0.;
    for (k = i + 1; k < ((n < i + bw + 1) ? n : i + bw + 1); ++k)
      sum += Abnd[n * (bw + i - k) + k] * b[k];
    b[i] = (b[i] - sum) / Abnd[n * bw + i];
  }
  return 0;
}


x86 assembly (clang, -O0):
<cholbnd_upper>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13ab <cholbnd_upper+0x29b>
mov    -0x24(%rbp),%ecx
sub    -0x20(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    1153 <cholbnd_upper+0x43>
xor    %eax,%eax
mov    %eax,-0x3c(%rbp)
jmp    115c <cholbnd_upper+0x4c>
mov    -0x24(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    128a <cholbnd_upper+0x17a>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x28(%rbp),%esi
sub    -0x24(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%ecx
sub    -0x20(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    11ad <cholbnd_upper+0x9d>
xor    %eax,%eax
mov    %eax,-0x40(%rbp)
jmp    11b6 <cholbnd_upper+0xa6>
mov    -0x24(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    123f <cholbnd_upper+0x12f>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x2c(%rbp),%esi
sub    -0x24(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x2c(%rbp),%esi
sub    -0x28(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11bc <cholbnd_upper+0xac>
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x28(%rbp),%esi
sub    -0x24(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1162 <cholbnd_upper+0x52>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%ecx
sub    -0x20(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    12c1 <cholbnd_upper+0x1b1>
xor    %eax,%eax
mov    %eax,-0x44(%rbp)
jmp    12ca <cholbnd_upper+0x1ba>
mov    -0x24(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1353 <cholbnd_upper+0x243>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x2c(%rbp),%esi
sub    -0x24(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x2c(%rbp),%esi
sub    -0x24(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12d0 <cholbnd_upper+0x1c0>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     138c <cholbnd_upper+0x27c>
movsd  -0x38(%rbp),%xmm0
call   1030 <sqrt@plt>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1398 <cholbnd_upper+0x288>
movl   $0x1,-0x4(%rbp)
jmp    15ea <cholbnd_upper+0x4da>
jmp    139d <cholbnd_upper+0x28d>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112d <cholbnd_upper+0x1d>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    14b3 <cholbnd_upper+0x3a3>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%ecx
sub    -0x20(%rbp),%ecx
xor    %eax,%eax
cmp    %ecx,%eax
jle    1406 <cholbnd_upper+0x2f6>
xor    %eax,%eax
mov    %eax,-0x48(%rbp)
jmp    140f <cholbnd_upper+0x2ff>
mov    -0x24(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    146e <cholbnd_upper+0x35e>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    -0x20(%rbp),%esi
add    -0x2c(%rbp),%esi
sub    -0x24(%rbp),%esi
imul   %esi,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1415 <cholbnd_upper+0x305>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    13d8 <cholbnd_upper+0x2c8>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x1c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     15e3 <cholbnd_upper+0x4d3>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    -0x20(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    153b <cholbnd_upper+0x42b>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1547 <cholbnd_upper+0x437>
mov    -0x24(%rbp),%eax
add    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x50(%rbp),%ecx
cmp    %ecx,%eax
jge    15a0 <cholbnd_upper+0x490>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%edx
add    -0x24(%rbp),%edx
sub    -0x2c(%rbp),%edx
imul   %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1516 <cholbnd_upper+0x406>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    14fb <cholbnd_upper+0x3eb>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cholbnd_upper>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %edx,%r14d
mov    %rdi,%r15
mov    %ecx,%r12d
imul   %edx,%r12d
test   %edx,%edx
mov    %rdi,0x38(%rsp)
jle    144f <cholbnd_upper+0x33f>
mov    %rsi,0x68(%rsp)
mov    %ecx,%edx
neg    %edx
mov    %rcx,0x40(%rsp)
movslq %ecx,%rcx
mov    %r14d,%ebp
mov    %r12,0x60(%rsp)
movslq %r12d,%rax
mov    %rax,0x8(%rsp)
mov    %rbp,0x10(%rsp)
shl    $0x4,%rbp
lea    0x1(%rcx),%rax
xorpd  %xmm2,%xmm2
mov    %rax,0x70(%rsp)
mov    %rax,0x50(%rsp)
mov    %rcx,0x30(%rsp)
mov    %r15,0x48(%rsp)
xor    %r13d,%r13d
xor    %eax,%eax
mov    %rax,0x20(%rsp)
mov    %r14,0xa0(%rsp)
mov    %rcx,0x28(%rsp)
jmp    11ec <cholbnd_upper+0xdc>
nopw   0x0(%rax,%rax,1)
mov    %edx,%ebx
call   1030 <sqrt@plt>
xorpd  %xmm2,%xmm2
mov    0x28(%rsp),%rcx
mov    %ebx,%edx
movsd  %xmm0,(%r12)
mov    %r13d,%esi
mov    %r13,%rax
add    $0x1,%rax
not    %esi
mov    %rsi,0x20(%rsp)
add    $0x1,%edx
addq   $0x8,0x48(%rsp)
addq   $0xffffffffffffffff,0x30(%rsp)
addq   $0xffffffffffffffff,0x50(%rsp)
mov    %rax,%r13
cmp    0x10(%rsp),%rax
je     1474 <cholbnd_upper+0x364>
test   %edx,%edx
mov    $0x0,%r10d
cmovg  %edx,%r10d
mov    %r13,%rax
sub    %rcx,%rax
mov    $0x0,%esi
cmovs  %rsi,%rax
mov    %eax,%eax
mov    %rax,0x58(%rsp)
cmp    %r13,%rax
jae    13d0 <cholbnd_upper+0x2c0>
mov    %edx,0x1c(%rsp)
mov    0x30(%rsp),%rax
lea    (%r10,%rax,1),%edx
imul   %r14d,%edx
mov    %rcx,%rdi
sub    %r10,%rdi
mov    0x70(%rsp),%rsi
sub    %r10,%rsi
mov    0x20(%rsp),%rax
add    0x40(%rsp),%eax
mov    %rax,0x20(%rsp)
lea    (%r10,%rcx,1),%rax
mov    %rax,0x88(%rsp)
sub    %r13,%rax
imul   0x10(%rsp),%rax
add    %r13,%rax
mov    %rax,0x80(%rsp)
lea    0x1(%r10),%rax
mov    %rax,0x78(%rsp)
xor    %r11d,%r11d
mov    %r10,%rax
mov    %r13,0x98(%rsp)
mov    %r10,0x90(%rsp)
jmp    12ca <cholbnd_upper+0x1ba>
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x8(%rsp),%rcx
add    %rax,%rcx
mov    0x38(%rsp),%r15
divsd  (%r15,%rcx,8),%xmm0
movsd  %xmm0,(%r15,%r12,8)
add    $0x1,%rax
add    $0x1,%r11
add    $0xffffffffffffffff,%rdi
add    $0xffffffffffffffff,%rsi
cmp    %r13,%rax
mov    0xa0(%rsp),%r14
jae    13f0 <cholbnd_upper+0x2e0>
mov    0x20(%rsp),%rcx
add    %eax,%ecx
imul   %r14d,%ecx
movslq %ecx,%r12
add    %r13,%r12
movsd  (%r15,%r12,8),%xmm0
cmp    %rax,0x58(%rsp)
jae    1290 <cholbnd_upper+0x180>
test   $0x1,%r11b
jne    1300 <cholbnd_upper+0x1f0>
mov    %r10,%r9
cmp    $0x1,%r11
je     1290 <cholbnd_upper+0x180>
jmp    133b <cholbnd_upper+0x22b>
nopl   0x0(%rax)
mov    0x80(%rsp),%rcx
movsd  (%r15,%rcx,8),%xmm1
mov    0x88(%rsp),%rcx
sub    %rax,%rcx
imul   0x10(%rsp),%rcx
add    %rax,%rcx
mulsd  (%r15,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
mov    0x78(%rsp),%r9
cmp    $0x1,%r11
je     1290 <cholbnd_upper+0x180>
lea    (%r9,%rdi,1),%rcx
mov    0x10(%rsp),%r13
imul   %r13,%rcx
add    %rax,%rcx
mov    0x38(%rsp),%r10
lea    (%r10,%rcx,8),%r14
mov    0x30(%rsp),%rcx
add    %r9,%rcx
imul   %r13,%rcx
mov    0x48(%rsp),%rbx
lea    (%rbx,%rcx,8),%r15
lea    (%r9,%rsi,1),%rcx
imul   %r13,%rcx
add    %rax,%rcx
lea    (%r10,%rcx,8),%r8
mov    0x50(%rsp),%rcx
add    %r9,%rcx
imul   %r13,%rcx
mov    0x98(%rsp),%r13
add    %r13,%rcx
lea    (%r10,%rcx,8),%rcx
mov    0x90(%rsp),%r10
xor    %ebx,%ebx
nopl   (%rax)
movsd  (%r15,%rbx,1),%xmm1
mulsd  (%r14,%rbx,1),%xmm1
subsd  %xmm1,%xmm0
movsd  (%rcx,%rbx,1),%xmm1
mulsd  (%r8,%rbx,1),%xmm1
subsd  %xmm1,%xmm0
add    $0x2,%r9
add    %rbp,%rbx
cmp    %r9,%rax
jne    13a0 <cholbnd_upper+0x290>
jmp    1290 <cholbnd_upper+0x180>
mov    0x8(%rsp),%rax
add    %r13,%rax
lea    (%r15,%rax,8),%r12
movsd  (%r15,%rax,8),%xmm0
ucomisd %xmm2,%xmm0
jae    1440 <cholbnd_upper+0x330>
jmp    146a <cholbnd_upper+0x35a>
nopl   (%rax)
mov    0x8(%rsp),%rax
add    %r13,%rax
lea    (%r15,%rax,8),%r12
movsd  (%r15,%rax,8),%xmm0
cmp    %r13,0x58(%rsp)
jae    1431 <cholbnd_upper+0x321>
nopl   0x0(%rax)
movslq %edx,%rdx
lea    (%rdx,%r13,1),%rax
movsd  (%r15,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
add    $0x1,%r10
add    %r14d,%edx
cmp    %r13,%r10
jb     1410 <cholbnd_upper+0x300>
mov    0x1c(%rsp),%edx
mov    0x28(%rsp),%rcx
ucomisd %xmm2,%xmm0
jb     146a <cholbnd_upper+0x35a>
jb     11a0 <cholbnd_upper+0x90>
sqrtsd %xmm0,%xmm0
jmp    11b2 <cholbnd_upper+0xa2>
movslq %r12d,%rax
movsd  (%rsi),%xmm0
mov    %rax,0x8(%rsp)
divsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,(%rsi)
jmp    1557 <cholbnd_upper+0x447>
mov    $0x1,%eax
jmp    165d <cholbnd_upper+0x54d>
mov    0x68(%rsp),%rsi
movsd  (%rsi),%xmm0
mov    0x8(%rsp),%rax
divsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,(%rsi)
xor    %r9d,%r9d
cmp    $0x2,%r14d
setge  %r8b
mov    0x40(%rsp),%rcx
mov    0x60(%rsp),%r12
jl     1557 <cholbnd_upper+0x447>
mov    $0x1,%edx
sub    %ecx,%edx
movslq %ecx,%r10
lea    -0x1(%rcx),%r11d
mov    $0x1,%ebp
jmp    14ee <cholbnd_upper+0x3de>
nopl   0x0(%rax)
movsd  (%rsi,%rbp,8),%xmm1
subsd  %xmm0,%xmm1
mov    0x8(%rsp),%rax
add    %rbp,%rax
divsd  (%r15,%rax,8),%xmm1
movsd  %xmm1,(%rsi,%rbp,8)
add    $0x1,%rbp
add    $0x1,%edx
add    $0xffffffff,%r11d
cmp    0x10(%rsp),%rbp
je     155a <cholbnd_upper+0x44a>
test   %edx,%edx
mov    $0x0,%edi
cmovg  %edx,%edi
mov    %rbp,%rax
sub    %r10,%rax
cmovs  %r9,%rax
mov    %eax,%eax
xorpd  %xmm0,%xmm0
cmp    %rbp,%rax
jae    14c0 <cholbnd_upper+0x3b0>
add    %r11d,%edi
imul   %r14d,%edi
test   %edx,%edx
mov    $0x0,%eax
cmovg  %edx,%eax
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %edi,%rdi
lea    (%rdi,%rbp,1),%rbx
movsd  (%r15,%rbx,8),%xmm1
mulsd  (%rsi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
add    %r14d,%edi
cmp    %rbp,%rax
jb     1530 <cholbnd_upper+0x420>
jmp    14c0 <cholbnd_upper+0x3b0>
xor    %r8d,%r8d
movslq %r14d,%r9
movsd  -0x8(%rsi,%r9,8),%xmm0
lea    (%r14,%r12,1),%eax
add    $0xffffffff,%eax
cltq
divsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rsi,%r9,8)
xor    %eax,%eax
test   %r8b,%r8b
je     165d <cholbnd_upper+0x54d>
lea    -0x2(%r14),%edi
mov    %r15,%rdx
movslq %edi,%r15
movslq %ecx,%r8
lea    (%r8,%r15,1),%rcx
sub    %r9,%rcx
add    $0x2,%rcx
imul   %r9,%rcx
lea    (%rdx,%rcx,8),%r10
add    $0xfffffffffffffff8,%r10
lea    0x0(,%r9,8),%rcx
mov    $0x8,%ebp
sub    %rcx,%rbp
lea    (%rsi,%r9,8),%rcx
mov    %r14d,%r11d
neg    %r11d
jmp    160a <cholbnd_upper+0x4fa>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi,%r15,8),%xmm1
subsd  %xmm0,%xmm1
mov    0x8(%rsp),%rdx
add    %r15,%rdx
mov    0x38(%rsp),%rdi
divsd  (%rdi,%rdx,8),%xmm1
movsd  %xmm1,(%rsi,%r15,8)
lea    -0x1(%r15),%rdx
add    $0xfffffffffffffff8,%r10
add    $0xfffffffffffffff8,%rcx
add    $0x1,%r11d
test   %r15,%r15
mov    %rdx,%r15
jle    165d <cholbnd_upper+0x54d>
lea    (%r15,%r8,1),%rdx
lea    0x1(%rdx),%edi
cmp    %r9,%rdx
cmovge %r14d,%edi
lea    0x1(%r15),%rdx
movslq %edi,%rbx
xorpd  %xmm0,%xmm0
cmp    %rbx,%rdx
jge    15d0 <cholbnd_upper+0x4c0>
add    %r11d,%edi
xorpd  %xmm0,%xmm0
mov    $0xffffffffffffffff,%rbx
mov    %r10,%rdx
nopl   0x0(%rax)
movsd  (%rdx),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
addsd  %xmm1,%xmm0
add    %rbp,%rdx
add    $0x1,%rbx
cmp    %ebx,%edi
jne    1640 <cholbnd_upper+0x530>
jmp    15d0 <cholbnd_upper+0x4c0>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @cholbnd_upper(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 0
%6 = mul nsw i32 %3, %2
br i1 %5, label %13, label %7
7: ; preds = %4
%8 = sext i32 %6 to i64
%9 = load double, double* %1
%10 = getelementptr inbounds double, double* %0, i64 %8
%11 = load double, double* %10
%12 = fdiv double %9, %11
store double %12, double* %1
br label %192
13: ; preds = %4
%14 = sub i32 0, %3
%15 = sext i32 %3 to i64
%16 = zext i32 %2 to i64
%17 = sext i32 %6 to i64
%18 = zext i32 %2 to i64
br label %19
19: ; preds = %13, %136
%20 = phi i64 [ 0, %13 ], [ %138, %136 ]
%21 = phi i32 [ %14, %13 ], [ %140, %136 ]
%22 = phi i32 [ 0, %13 ], [ %139, %136 ]
%23 = trunc i64 %20 to i32
%24 = call i32 @llvm.smax.i32(i32 %21, i32 0)
%25 = zext i32 %24 to i64
%26 = sub nsw i64 %20, %15
%27 = icmp sgt i64 %26, 0
%28 = select i1 %27, i64 %26, i64 0
%29 = and i64 %28, 4294967295
%30 = icmp ult i64 %29, %20
br i1 %30, label %35, label %31
31: ; preds = %19
%32 = add nsw i64 %20, %17
%33 = getelementptr inbounds double, double* %0, i64 %32
%34 = load double, double* %33
br label %132
35: ; preds = %19
%36 = add i32 %22, %3
%37 = add nsw i64 %25, %15
%38 = sub nsw i64 %37, %20
%39 = mul nsw i64 %38, %16
%40 = add nsw i64 %39, %20
%41 = getelementptr inbounds double, double* %0, i64 %40
%42 = add nuw nsw i64 %25, 1
br label %43
43: ; preds = %35, %103
%44 = phi i64 [ 0, %35 ], [ %111, %103 ]
%45 = phi i64 [ %25, %35 ], [ %109, %103 ]
%46 = trunc i64 %45 to i32
%47 = add i32 %36, %46
%48 = mul nsw i32 %47, %2
%49 = sext i32 %48 to i64
%50 = add nsw i64 %20, %49
%51 = getelementptr inbounds double, double* %0, i64 %50
%52 = load double, double* %51
%53 = icmp ult i64 %29, %45
br i1 %53, label %54, label %103
54: ; preds = %43
%55 = and i64 %44, 1
%56 = icmp eq i64 %55, 0
br i1 %56, label %66, label %57
57: ; preds = %54
%58 = load double, double* %41
%59 = sub nsw i64 %37, %45
%60 = mul nsw i64 %59, %16
%61 = add nsw i64 %60, %45
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = load double, double* %62
%64 = fneg double %58
%65 = tail call double @llvm.fmuladd.f64(double %64, double %63, double %52)
br label %66
66: ; preds = %57, %54
%67 = phi double [ undef, %54 ], [ %65, %57 ]
%68 = phi i64 [ %25, %54 ], [ %42, %57 ]
%69 = phi double [ %52, %54 ], [ %65, %57 ]
%70 = icmp eq i64 %44, 1
br i1 %70, label %103, label %71
71: ; preds = %66, %71
%72 = phi i64 [ %101, %71 ], [ %68, %66 ]
%73 = phi double [ %100, %71 ], [ %69, %66 ]
%74 = add nsw i64 %72, %15
%75 = sub nsw i64 %74, %20
%76 = mul nsw i64 %75, %16
%77 = add nsw i64 %76, %20
%78 = getelementptr inbounds double, double* %0, i64 %77
%79 = load double, double* %78
%80 = sub nsw i64 %74, %45
%81 = mul nsw i64 %80, %16
%82 = add nsw i64 %81, %45
%83 = getelementptr inbounds double, double* %0, i64 %82
%84 = load double, double* %83
%85 = fneg double %79
%86 = tail call double @llvm.fmuladd.f64(double %85, double %84, double %73)
%87 = add nuw nsw i64 %72, 1
%88 = add nsw i64 %87, %15
%89 = sub nsw i64 %88, %20
%90 = mul nsw i64 %89, %16
%91 = add nsw i64 %90, %20
%92 = getelementptr inbounds double, double* %0, i64 %91
%93 = load double, double* %92
%94 = sub nsw i64 %88, %45
%95 = mul nsw i64 %94, %16
%96 = add nsw i64 %95, %45
%97 = getelementptr inbounds double, double* %0, i64 %96
%98 = load double, double* %97
%99 = fneg double %93
%100 = tail call double @llvm.fmuladd.f64(double %99, double %98, double %86)
%101 = add nuw nsw i64 %72, 2
%102 = icmp eq i64 %101, %45
br i1 %102, label %103, label %71
103: ; preds = %66, %71, %43
%104 = phi double [ %52, %43 ], [ %67, %66 ], [ %100, %71 ]
%105 = add nsw i64 %45, %17
%106 = getelementptr inbounds double, double* %0, i64 %105
%107 = load double, double* %106
%108 = fdiv double %104, %107
store double %108, double* %51
%109 = add nuw nsw i64 %45, 1
%110 = icmp ult i64 %109, %20
%111 = add i64 %44, 1
br i1 %110, label %43, label %112
112: ; preds = %103
%113 = add nsw i64 %20, %17
%114 = getelementptr inbounds double, double* %0, i64 %113
%115 = load double, double* %114
br i1 %30, label %116, label %132
116: ; preds = %112
%117 = add i32 %22, %3
br label %118
118: ; preds = %116, %118
%119 = phi i64 [ %25, %116 ], [ %130, %118 ]
%120 = phi double [ %115, %116 ], [ %129, %118 ]
%121 = trunc i64 %119 to i32
%122 = add i32 %117, %121
%123 = mul nsw i32 %122, %2
%124 = sext i32 %123 to i64
%125 = add nsw i64 %20, %124
%126 = getelementptr inbounds double, double* %0, i64 %125
%127 = load double, double* %126
%128 = fneg double %127
%129 = tail call double @llvm.fmuladd.f64(double %128, double %127, double %120)
%130 = add nuw nsw i64 %119, 1
%131 = icmp ult i64 %130, %20
br i1 %131, label %118, label %132
132: ; preds = %118, %31, %112
%133 = phi double* [ %114, %112 ], [ %33, %31 ], [ %114, %118 ]
%134 = phi double [ %115, %112 ], [ %34, %31 ], [ %129, %118 ]
%135 = fcmp ult double %134, 0.000000e+00
br i1 %135, label %246, label %136
136: ; preds = %132
%137 = tail call double @sqrt(double %134)
store double %137, double* %133
%138 = add nuw nsw i64 %20, 1
%139 = xor i32 %23, -1
%140 = add i32 %21, 1
%141 = icmp eq i64 %138, %18
br i1 %141, label %142, label %19
142: ; preds = %136
%143 = load double, double* %1
%144 = getelementptr inbounds double, double* %0, i64 %17
%145 = load double, double* %144
%146 = fdiv double %143, %145
store double %146, double* %1
%147 = icmp sgt i32 %2, 1
br i1 %147, label %148, label %192
148: ; preds = %142
%149 = sub i32 1, %3
%150 = sext i32 %3 to i64
%151 = zext i32 %2 to i64
br label %152
152: ; preds = %148, %180
%153 = phi i64 [ 1, %148 ], [ %189, %180 ]
%154 = phi i32 [ %149, %148 ], [ %190, %180 ]
%155 = sub nsw i64 %153, %150
%156 = icmp sgt i64 %155, 0
%157 = select i1 %156, i64 %155, i64 0
%158 = and i64 %157, 4294967295
%159 = icmp ult i64 %158, %153
br i1 %159, label %160, label %180
160: ; preds = %152
%161 = call i32 @llvm.smax.i32(i32 %154, i32 0)
%162 = zext i32 %161 to i64
%163 = trunc i64 %153 to i32
%164 = sub i32 %3, %163
br label %165
165: ; preds = %160, %165
%166 = phi i64 [ %162, %160 ], [ %178, %165 ]
%167 = phi double [ 0.000000e+00, %160 ], [ %177, %165 ]
%168 = trunc i64 %166 to i32
%169 = add i32 %164, %168
%170 = mul nsw i32 %169, %2
%171 = sext i32 %170 to i64
%172 = add nsw i64 %153, %171
%173 = getelementptr inbounds double, double* %0, i64 %172
%174 = load double, double* %173
%175 = getelementptr inbounds double, double* %1, i64 %166
%176 = load double, double* %175
%177 = tail call double @llvm.fmuladd.f64(double %174, double %176, double %167)
%178 = add nuw nsw i64 %166, 1
%179 = icmp ult i64 %178, %153
br i1 %179, label %165, label %180
180: ; preds = %165, %152
%181 = phi double [ 0.000000e+00, %152 ], [ %177, %165 ]
%182 = getelementptr inbounds double, double* %1, i64 %153
%183 = load double, double* %182
%184 = fsub double %183, %181
%185 = add nsw i64 %153, %17
%186 = getelementptr inbounds double, double* %0, i64 %185
%187 = load double, double* %186
%188 = fdiv double %184, %187
store double %188, double* %182
%189 = add nuw nsw i64 %153, 1
%190 = add i32 %154, 1
%191 = icmp eq i64 %189, %151
br i1 %191, label %192, label %152
192: ; preds = %180, %7, %142
%193 = phi i1 [ false, %7 ], [ false, %142 ], [ %147, %180 ]
%194 = phi i64 [ %8, %7 ], [ %17, %142 ], [ %17, %180 ]
%195 = add nsw i32 %2, -1
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds double, double* %1, i64 %196
%198 = load double, double* %197
%199 = add nsw i32 %6, %195
%200 = sext i32 %199 to i64
%201 = getelementptr inbounds double, double* %0, i64 %200
%202 = load double, double* %201
%203 = fdiv double %198, %202
store double %203, double* %197
br i1 %193, label %204, label %246
204: ; preds = %192
%205 = add i32 %2, -2
%206 = sext i32 %2 to i64
%207 = sext i32 %205 to i64
%208 = sext i32 %3 to i64
br label %209
209: ; preds = %204, %235
%210 = phi i64 [ %207, %204 ], [ %244, %235 ]
%211 = phi i64 [ %206, %204 ], [ %212, %235 ]
%212 = add nsw i64 %211, -1
%213 = add nsw i64 %210, %208
%214 = icmp slt i64 %213, %206
%215 = trunc i64 %213 to i32
%216 = add i32 %215, 1
%217 = select i1 %214, i32 %216, i32 %2
%218 = add nsw i64 %210, 1
%219 = sext i32 %217 to i64
%220 = icmp slt i64 %218, %219
br i1 %220, label %221, label %235
221: ; preds = %209, %221
%222 = phi i64 [ %232, %221 ], [ %212, %209 ]
%223 = phi double [ %231, %221 ], [ 0.000000e+00, %209 ]
%224 = sub nsw i64 %213, %222
%225 = mul nsw i64 %224, %206
%226 = add nsw i64 %225, %222
%227 = getelementptr inbounds double, double* %0, i64 %226
%228 = load double, double* %227
%229 = getelementptr inbounds double, double* %1, i64 %222
%230 = load double, double* %229
%231 = tail call double @llvm.fmuladd.f64(double %228, double %230, double %223)
%232 = add nsw i64 %222, 1
%233 = trunc i64 %232 to i32
%234 = icmp eq i32 %217, %233
br i1 %234, label %235, label %221
235: ; preds = %221, %209
%236 = phi double [ 0.000000e+00, %209 ], [ %231, %221 ]
%237 = getelementptr inbounds double, double* %1, i64 %210
%238 = load double, double* %237
%239 = fsub double %238, %236
%240 = add nsw i64 %210, %194
%241 = getelementptr inbounds double, double* %0, i64 %240
%242 = load double, double* %241
%243 = fdiv double %239, %242
store double %243, double* %237
%244 = add nsw i64 %210, -1
%245 = icmp sgt i64 %210, 0
br i1 %245, label %209, label %246
246: ; preds = %132, %235, %192
%247 = phi i32 [ 0, %192 ], [ 0, %235 ], [ 1, %132 ]
ret i32 %247
}