#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);
void LowerTriangular(double **A, int jStrt, int jEnd, int na, int *rnk,
                     int *ib) {
  int i, j, k, itmp;
  double dtmp;
  ib[0] = jStrt;
  for (i = 1; i < na; i++)
    ib[i] = -1;
  *rnk = 1;
  k = jStrt + 1;
  while (*rnk < na && k <= jEnd) {
    dtmp = 1.0e-13;
    itmp = -1;
    for (i = *rnk; i < na; i++)
      if (fabs(A[k][i]) > dtmp) {
        dtmp = fabs(A[k][i]);
        itmp = i;
      }
    if (itmp >= 0) {
      for (i = 0; i < itmp; i++)
        A[k][i] /= A[k][itmp];
      for (i = itmp + 1; i < na; i++)
        A[k][i] /= A[k][itmp];
      for (j = k + 1; j <= jEnd; j++) {
        for (i = 0; i < itmp; i++)
          A[j][i] -= A[k][i] * A[j][itmp];
        for (i = itmp + 1; i < na; i++)
          A[j][i] -= A[k][i] * A[j][itmp];
        A[j][itmp] /= A[k][itmp];
      }
      if (itmp != *rnk) {
        for (j = jStrt; j <= jEnd; j++) {
          dtmp = A[j][*rnk];
          A[j][*rnk] = A[j][itmp];
          A[j][itmp] = dtmp;
        }
      }
      for (i = 0; i < *rnk; i++)
        A[k][i] = 0.0e0;
      A[k][*rnk] = 1.0e0;
      for (i = *rnk + 1; i < na; i++)
        A[k][i] = 0.0e0;
      ++(*rnk);
      ib[*rnk - 1] = k;
    }
    ++k;
  }
}


x86 assembly (clang, -O0):
<LowerTriangular>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1152 <LowerTriangular+0x52>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0xffffffff,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1129 <LowerTriangular+0x29>
mov    -0x20(%rbp),%rax
movl   $0x1,(%rax)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%ecx
xor    %eax,%eax
cmp    -0x14(%rbp),%ecx
mov    %al,-0x41(%rbp)
jge    1185 <LowerTriangular+0x85>
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
setle  %al
mov    %al,-0x41(%rbp)
mov    -0x41(%rbp),%al
test   $0x1,%al
jne    1195 <LowerTriangular+0x95>
jmp    159c <LowerTriangular+0x49c>
movsd  0xe63(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0xffffffff,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1227 <LowerTriangular+0x127>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe36(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1214 <LowerTriangular+0x114>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe0b(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    1219 <LowerTriangular+0x119>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b2 <LowerTriangular+0xb2>
cmpl   $0x0,-0x38(%rbp)
jl     158e <LowerTriangular+0x48e>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1285 <LowerTriangular+0x185>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1238 <LowerTriangular+0x138>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12db <LowerTriangular+0x1db>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    128e <LowerTriangular+0x18e>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     143b <LowerTriangular+0x33b>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1374 <LowerTriangular+0x274>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12f7 <LowerTriangular+0x1f7>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13fa <LowerTriangular+0x2fa>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    137d <LowerTriangular+0x27d>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12e4 <LowerTriangular+0x1e4>
mov    -0x38(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
je     14d3 <LowerTriangular+0x3d3>
mov    -0xc(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     14ce <LowerTriangular+0x3ce>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movslq (%rcx),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movslq (%rcx),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1450 <LowerTriangular+0x350>
jmp    14d3 <LowerTriangular+0x3d3>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    150f <LowerTriangular+0x40f>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    14da <LowerTriangular+0x3da>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x20(%rbp),%rcx
movslq (%rcx),%rcx
movsd  0xade(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    156d <LowerTriangular+0x46d>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    153b <LowerTriangular+0x43b>
mov    -0x20(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    %ecx,(%rax)
mov    -0x34(%rbp),%edx
mov    -0x28(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rcx),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1165 <LowerTriangular+0x65>
pop    %rbp
ret

x86 assembly (clang, -O2):
<LowerTriangular>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r8,%r12
mov    %rdx,0x28(%rsp)
mov    %rsi,0x20(%rsp)
mov    %esi,(%r9)
cmp    $0x1,%ecx
jle    1635 <LowerTriangular+0x525>
mov    %ecx,%r13d
mov    %rdi,%r14
mov    %r9,0x70(%rsp)
lea    0x4(%r9),%rdi
lea    -0x1(%r13),%edx
shl    $0x2,%rdx
mov    $0xff,%esi
call   1030 <memset@plt>
movl   $0x1,(%r12)
cmp    $0x2,%r13d
jl     163d <LowerTriangular+0x52d>
mov    0x20(%rsp),%rax
cmp    %eax,0x28(%rsp)
jle    163d <LowerTriangular+0x52d>
mov    0x20(%rsp),%rcx
lea    0x2(%rcx),%esi
mov    0x28(%rsp),%rax
movslq %eax,%rdx
mov    %rdx,0x30(%rsp)
movslq %ecx,%r10
lea    -0x2(%r13),%edx
mov    %edx,0xc(%rsp)
movslq %r13d,%rdx
add    $0x1,%eax
mov    %rax,0x60(%rsp)
sub    %ecx,%eax
mov    %eax,0x3c(%rsp)
mov    %rdx,0x80(%rsp)
neg    %rdx
mov    %rdx,0x78(%rsp)
lea    0x1(%r10),%rax
mov    %rax,0x58(%rsp)
mov    %r13d,%r15d
mov    $0x1,%r8d
movapd 0xe39(%rip),%xmm3
mov    %r10,%rax
mov    %r10,0x40(%rsp)
mov    %r12,0x68(%rsp)
mov    %r13,0x48(%rsp)
nopl   0x0(%rax)
lea    0x1(%r10),%r9
cmp    %r13d,%r8d
jge    14b5 <LowerTriangular+0x3a5>
mov    (%r14,%r9,8),%rdx
movslq %r8d,%rcx
mov    0x80(%rsp),%rbp
sub    %rcx,%rbp
mov    %rcx,%rdi
not    %rdi
mov    $0xffffffff,%eax
cmp    0x78(%rsp),%rdi
jne    1230 <LowerTriangular+0x120>
movsd  0xdd7(%rip),%xmm0
jmp    128a <LowerTriangular+0x17a>
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0xffffffff,%eax
movsd  0xdbc(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rdx,%rcx,8),%xmm2
andpd  %xmm3,%xmm2
ucomisd %xmm1,%xmm2
movsd  0x8(%rdx,%rcx,8),%xmm0
cmova  %ecx,%eax
maxsd  %xmm1,%xmm2
andpd  %xmm3,%xmm0
ucomisd %xmm2,%xmm0
lea    0x1(%rcx),%ebx
cmova  %ebx,%eax
maxsd  %xmm2,%xmm0
add    $0x2,%rcx
movapd %xmm0,%xmm1
add    $0xfffffffffffffffe,%rdi
jne    1250 <LowerTriangular+0x140>
test   $0x1,%bpl
je     12a4 <LowerTriangular+0x194>
movsd  (%rdx,%rcx,8),%xmm1
andpd  0xd73(%rip),%xmm1
ucomisd %xmm0,%xmm1
cmova  %ecx,%eax
test   %eax,%eax
js     14b5 <LowerTriangular+0x3a5>
mov    %eax,%r11d
je     1311 <LowerTriangular+0x201>
mov    (%r14,%r9,8),%rcx
cmp    $0x1,%eax
jne    12be <LowerTriangular+0x1ae>
xor    %edx,%edx
jmp    12fb <LowerTriangular+0x1eb>
mov    %r11d,%edi
and    $0xfffffffe,%edi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rcx,%rdx,8),%xmm0
divsd  (%rcx,%r11,8),%xmm0
movsd  0x8(%rcx,%rdx,8),%xmm1
movsd  %xmm0,(%rcx,%rdx,8)
divsd  (%rcx,%r11,8),%xmm1
movsd  %xmm1,0x8(%rcx,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%rdi
jne    12d0 <LowerTriangular+0x1c0>
test   $0x1,%r11b
je     1311 <LowerTriangular+0x201>
movsd  (%rcx,%rdx,8),%xmm0
divsd  (%rcx,%r11,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
mov    %r9,0x10(%rsp)
mov    %r8,0x50(%rsp)
lea    0x1(%rax),%r9d
movslq %eax,%rdx
mov    %eax,%r8d
not    %r8d
cmp    %r13d,%r9d
jge    1389 <LowerTriangular+0x279>
mov    0x10(%rsp),%rcx
mov    (%r14,%rcx,8),%rcx
movslq %r9d,%rdi
lea    (%r8,%r13,1),%ebp
test   $0x1,%bpl
je     1356 <LowerTriangular+0x246>
movsd  (%rcx,%rdi,8),%xmm0
divsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdi,8)
add    $0x1,%rdi
cmp    %eax,0xc(%rsp)
je     1389 <LowerTriangular+0x279>
nopl   0x0(%rax)
movsd  (%rcx,%rdi,8),%xmm0
divsd  (%rcx,%rdx,8),%xmm0
movsd  0x8(%rcx,%rdi,8),%xmm1
movsd  %xmm0,(%rcx,%rdi,8)
divsd  (%rcx,%rdx,8),%xmm1
movsd  %xmm1,0x8(%rcx,%rdi,8)
add    $0x2,%rdi
cmp    %edi,%r15d
jne    1360 <LowerTriangular+0x250>
mov    %esi,0x1c(%rsp)
add    $0x2,%r10
cmp    0x30(%rsp),%r10
jle    14d4 <LowerTriangular+0x3c4>
mov    0x50(%rsp),%rbx
cmp    %ebx,%eax
mov    0x60(%rsp),%rdi
mov    0x10(%rsp),%r9
je     1429 <LowerTriangular+0x319>
movslq %ebx,%rax
mov    0x40(%rsp),%rcx
testb  $0x1,0x3c(%rsp)
je     13e0 <LowerTriangular+0x2d0>
mov    0x40(%rsp),%rcx
mov    (%r14,%rcx,8),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  (%rcx,%rdx,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  %xmm0,(%rcx,%rdx,8)
mov    0x58(%rsp),%rcx
mov    0x20(%rsp),%rsi
cmp    %esi,0x28(%rsp)
je     1429 <LowerTriangular+0x319>
nopl   0x0(%rax,%rax,1)
mov    (%r14,%rcx,8),%rsi
movsd  (%rsi,%rax,8),%xmm0
movsd  (%rsi,%rdx,8),%xmm1
movsd  %xmm1,(%rsi,%rax,8)
movsd  %xmm0,(%rsi,%rdx,8)
mov    0x8(%r14,%rcx,8),%rsi
movsd  (%rsi,%rax,8),%xmm0
movsd  (%rsi,%rdx,8),%xmm1
movsd  %xmm1,(%rsi,%rax,8)
movsd  %xmm0,(%rsi,%rdx,8)
add    $0x2,%rcx
cmp    %ecx,%edi
jne    13f0 <LowerTriangular+0x2e0>
mov    (%r14,%r9,8),%r12
test   %ebx,%ebx
jle    1453 <LowerTriangular+0x343>
mov    %ebx,%edx
shl    $0x3,%rdx
mov    %r12,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x10(%rsp),%r9
mov    0x50(%rsp),%rbx
movapd 0xbbd(%rip),%xmm3
movslq %ebx,%r13
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r12,%r13,8)
lea    0x1(%rbx),%ebp
cmp    0x48(%rsp),%ebp
jge    1496 <LowerTriangular+0x386>
movslq %ebp,%rax
lea    (%r12,%rax,8),%rdi
mov    0xc(%rsp),%eax
sub    %ebx,%eax
lea    0x8(,%rax,8),%rdx
xor    %esi,%esi
call   1030 <memset@plt>
mov    0x10(%rsp),%r9
movapd 0xb7a(%rip),%xmm3
mov    0x68(%rsp),%r12
mov    %ebp,(%r12)
mov    0x70(%rsp),%rax
mov    %r9d,(%rax,%r13,4)
mov    (%r12),%r8d
mov    0x1c(%rsp),%esi
mov    0x48(%rsp),%r13
cmp    %r13d,%r8d
jge    163d <LowerTriangular+0x52d>
add    $0x1,%esi
mov    %r9,%r10
cmp    0x30(%rsp),%r9
jl     11f0 <LowerTriangular+0xe0>
jmp    163d <LowerTriangular+0x52d>
movslq 0x1c(%rsp),%rbx
movslq %r9d,%r12
mov    0x10(%rsp),%rcx
mov    (%r14,%rcx,8),%rdi
add    %r13d,%r8d
mov    %r11d,%ecx
and    $0xfffffffe,%ecx
lea    0x1(%r12),%r10
jmp    1523 <LowerTriangular+0x413>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  0x0(%rbp,%rdx,8),%xmm0
divsd  (%rdi,%rdx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rdx,8)
lea    0x1(%rbx),%rsi
cmp    0x30(%rsp),%rbx
mov    %rsi,%rbx
jge    139c <LowerTriangular+0x28c>
test   %eax,%eax
jle    15a1 <LowerTriangular+0x491>
mov    (%r14,%rbx,8),%rbp
xor    %esi,%esi
cmp    $0x1,%eax
je     1580 <LowerTriangular+0x470>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdi,%rsi,8),%xmm0
mulsd  0x0(%rbp,%rdx,8),%xmm0
movsd  0x0(%rbp,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  0x8(%rbp,%rsi,8),%xmm0
movsd  %xmm1,0x0(%rbp,%rsi,8)
movsd  0x8(%rdi,%rsi,8),%xmm1
mulsd  0x0(%rbp,%rdx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rbp,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%rcx
jne    1540 <LowerTriangular+0x430>
test   $0x1,%r11b
je     15a1 <LowerTriangular+0x491>
movsd  (%rdi,%rsi,8),%xmm0
movsd  0x0(%rbp,%rsi,8),%xmm1
mulsd  0x0(%rbp,%rdx,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp,%rsi,8)
mov    (%r14,%rbx,8),%rbp
cmp    %r13d,%r9d
jge    1500 <LowerTriangular+0x3f0>
mov    %r12,%rsi
test   $0x1,%r8b
je     15d8 <LowerTriangular+0x4c8>
movsd  (%rdi,%r12,8),%xmm0
movsd  0x0(%rbp,%r12,8),%xmm1
mulsd  0x0(%rbp,%rdx,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp,%r12,8)
mov    %r10,%rsi
cmp    %eax,0xc(%rsp)
je     1500 <LowerTriangular+0x3f0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdi,%rsi,8),%xmm0
mulsd  0x0(%rbp,%rdx,8),%xmm0
movsd  0x0(%rbp,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  0x8(%rbp,%rsi,8),%xmm0
movsd  %xmm1,0x0(%rbp,%rsi,8)
movsd  0x8(%rdi,%rsi,8),%xmm1
mulsd  0x0(%rbp,%rdx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rbp,%rsi,8)
add    $0x2,%rsi
cmp    %esi,%r15d
jne    15f0 <LowerTriangular+0x4e0>
jmp    1500 <LowerTriangular+0x3f0>
movl   $0x1,(%r12)
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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