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

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ exit(int);
double fabs(double);
int /*<<< orphan*/ printf(char *);
void gaussj(int n, double a[], int m, double b[]) {
  int i, j, k, l, ll;
  int icol, irow;
  int ipiv[50], indxr[50], indxc[50];
  double big, dum, pivinv;
  for (j = 0; j < n; ++j)
    ipiv[j] = 0;
  for (i = 0; i < n; ++i) {
    big = 0.0;
    for (j = 0; j < n; ++j) {
      if (ipiv[j] != 1) {
        for (k = 0; k < n; ++k) {
          if (ipiv[k] == 0) {
            if (fabs(a[j + k * n]) >= big) {
              big = fabs(a[j + k * n]);
              irow = j;
              icol = k;
            }
          } else if (ipiv[k] > 1) {
            printf("Error: singular matrix in gaussj\n");
            exit(99);
          }
        }
      }
    }
    ++ipiv[icol];
    if (irow != icol) {
      for (l = 0; l < n; ++l) {
        dum = a[irow + l * n];
        a[irow + l * n] = a[icol + l * n];
        a[icol + l * n] = dum;
      }
      for (l = 0; l < m; ++l) {
        dum = b[irow + l * m];
        b[irow + l * m] = b[icol + l * m];
        b[icol + l * m] = dum;
      }
    }
    indxr[i] = irow;
    indxc[i] = icol;
    if (a[icol + n * icol] == 0.0) {
      printf("Error: singular matrix in gaussj\n");
      exit(99);
    }
    pivinv = 1.0 / a[icol + icol * n];
    a[icol + icol * n] = 1.0;
    for (l = 0; l < n; ++l)
      a[icol + l * n] *= pivinv;
    for (l = 0; l < m; ++l)
      b[icol + l * m] *= pivinv;
    for (ll = 0; ll < n; ++ll) {
      if (ll != icol) {
        dum = a[ll + icol * n];
        a[ll + icol * n] = 0.0;
        for (l = 0; l < n; ++l)
          a[ll + l * n] -= a[icol + l * n] * dum;
        for (l = 0; l < m; ++l)
          b[ll + l * m] -= b[icol + l * m] * dum;
      }
    }
  }
  for (l = (n - 1); l >= 0; --l) {
    if (indxr[l] != indxc[l]) {
      for (k = 0; k < n; ++k) {
        dum = a[k + indxr[l] * n];
        a[k + indxr[l] * n] = a[k + indxc[l] * n];
        a[k + indxc[l] * n] = dum;
      }
    }
  }
}


x86 assembly (clang, -O0):
<gaussj>:
push   %rbp
mov    %rsp,%rbp
sub    $0x2d0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1169 <gaussj+0x49>
movslq -0x28(%rbp),%rax
movl   $0x0,-0x110(%rbp,%rax,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1140 <gaussj+0x20>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    168e <gaussj+0x56e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x2b8(%rbp)
ff
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1296 <gaussj+0x176>
movslq -0x28(%rbp),%rax
cmpl   $0x1,-0x110(%rbp,%rax,4)
01
je     1283 <gaussj+0x163>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    127e <gaussj+0x15e>
movslq -0x2c(%rbp),%rax
cmpl   $0x0,-0x110(%rbp,%rax,4)
jne    123e <gaussj+0x11e>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe20(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x2b8(%rbp),%xmm0
ff
jb     1239 <gaussj+0x119>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xdef(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x2b8(%rbp)
ff
mov    -0x28(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    126b <gaussj+0x14b>
movslq -0x2c(%rbp),%rax
cmpl   $0x1,-0x110(%rbp,%rax,4)
01
jle    1266 <gaussj+0x146>
lea    0xdc9(%rip),%rdi
call   1030 <printf@plt>
mov    $0x63,%edi
call   1040 <exit@plt>
jmp    126b <gaussj+0x14b>
jmp    1270 <gaussj+0x150>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b3 <gaussj+0x93>
jmp    1283 <gaussj+0x163>
jmp    1288 <gaussj+0x168>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    118e <gaussj+0x6e>
movslq -0x38(%rbp),%rax
mov    -0x110(%rbp,%rax,4),%ecx
add    $0x1,%ecx
mov    %ecx,-0x110(%rbp,%rax,4)
mov    -0x3c(%rbp),%eax
cmp    -0x38(%rbp),%eax
je     13de <gaussj+0x2be>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1348 <gaussj+0x228>
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2c0(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x2c0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12be <gaussj+0x19e>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13d9 <gaussj+0x2b9>
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2c0(%rbp)
ff
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x2c0(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    134f <gaussj+0x22f>
jmp    13de <gaussj+0x2be>
mov    -0x3c(%rbp),%ecx
movslq -0x24(%rbp),%rax
mov    %ecx,-0x1e0(%rbp,%rax,4)
mov    -0x38(%rbp),%ecx
movslq -0x24(%rbp),%rax
mov    %ecx,-0x2b0(%rbp,%rax,4)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x4(%rbp),%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    143b <gaussj+0x31b>
jp     143b <gaussj+0x31b>
lea    0xbf4(%rip),%rdi
call   1030 <printf@plt>
mov    $0x63,%edi
call   1040 <exit@plt>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xbaa(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2c8(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xb82(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14c9 <gaussj+0x3a9>
movsd  -0x2c8(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    148a <gaussj+0x36a>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    150f <gaussj+0x3ef>
movsd  -0x2c8(%rbp),%xmm0
ff
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    14d0 <gaussj+0x3b0>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    167b <gaussj+0x55b>
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
je     1668 <gaussj+0x548>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2c0(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15e6 <gaussj+0x4c6>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x2c0(%rbp),%xmm2
ff
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%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    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1570 <gaussj+0x450>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1663 <gaussj+0x543>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x2c0(%rbp),%xmm2
ff
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%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    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    15ed <gaussj+0x4cd>
jmp    1668 <gaussj+0x548>
jmp    166d <gaussj+0x54d>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1516 <gaussj+0x3f6>
jmp    1680 <gaussj+0x560>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1170 <gaussj+0x50>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
jl     1786 <gaussj+0x666>
movslq -0x30(%rbp),%rax
mov    -0x1e0(%rbp,%rax,4),%eax
movslq -0x30(%rbp),%rcx
cmp    -0x2b0(%rbp,%rcx,4),%eax
je     1773 <gaussj+0x653>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    176e <gaussj+0x64e>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
movslq -0x30(%rbp),%rdx
mov    -0x1e0(%rbp,%rdx,4),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x2c0(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
movslq -0x30(%rbp),%rdx
mov    -0x2b0(%rbp,%rdx,4),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
movslq -0x30(%rbp),%rdx
mov    -0x1e0(%rbp,%rdx,4),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x2c0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
movslq -0x30(%rbp),%rdx
mov    -0x2b0(%rbp,%rdx,4),%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    16c4 <gaussj+0x5a4>
jmp    1773 <gaussj+0x653>
jmp    1778 <gaussj+0x658>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    1697 <gaussj+0x577>
add    $0x2d0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gaussj>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x2e8,%rsp
mov    %rcx,0x18(%rsp)
mov    %rsi,(%rsp)
mov    %edi,0xc(%rsp)
test   %edi,%edi
jle    18e0 <gaussj+0x7b0>
mov    %edx,%r15d
mov    0xc(%rsp),%ebp
lea    0x0(,%rbp,4),%rdx
lea    0x80(%rsp),%rdi
xor    %r12d,%r12d
xor    %esi,%esi
call   1040 <memset@plt>
movslq %r15d,%rax
mov    %eax,%ecx
lea    -0x1(%rbp),%rdx
mov    %rdx,0x20(%rsp)
lea    -0x1(%rcx),%rdx
mov    %rdx,0x30(%rsp)
mov    %ebp,%r14d
and    $0xfffffffe,%r14d
mov    %rcx,0x38(%rsp)
and    $0xfffffffe,%ecx
mov    %rcx,0x40(%rsp)
lea    0x0(,%rbp,8),%rbx
mov    %rbp,%r8
shl    $0x4,%r8
mov    %rax,%r10
shl    $0x4,%r10
mov    %rax,0x10(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,0x70(%rsp)
xorpd  %xmm0,%xmm0
movsd  0xe3a(%rip),%xmm1
movapd 0xe22(%rip),%xmm2
mov    (%rsp),%rdi
mov    %r15d,0x28(%rsp)
mov    %r10,0x68(%rsp)
jmp    1216 <gaussj+0xe6>
xchg   %ax,%ax
mov    0x78(%rsp),%r12
add    $0x1,%r12
cmp    %rbp,%r12
mov    0x28(%rsp),%r15d
mov    0x68(%rsp),%r10
mov    0x5c(%rsp),%r11d
mov    0x60(%rsp),%r9
je     16ab <gaussj+0x57b>
mov    %rdi,%rdx
xor    %ecx,%ecx
xorpd  %xmm3,%xmm3
jmp    123d <gaussj+0x10d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%rcx
add    $0x8,%rdx
cmp    %rbp,%rcx
je     1290 <gaussj+0x160>
cmpl   $0x1,0x80(%rsp,%rcx,4)
01
je     1230 <gaussj+0x100>
mov    %rdx,%rsi
xor    %edi,%edi
jmp    1265 <gaussj+0x135>
xchg   %ax,%ax
cmp    $0x2,%eax
jge    18f2 <gaussj+0x7c2>
add    $0x1,%rdi
add    %rbx,%rsi
cmp    %rdi,%rbp
je     1230 <gaussj+0x100>
mov    0x80(%rsp,%rdi,4),%eax
test   %eax,%eax
jne    1250 <gaussj+0x120>
movsd  (%rsi),%xmm4
andpd  %xmm2,%xmm4
ucomisd %xmm3,%xmm4
jb     1259 <gaussj+0x129>
mov    %edi,%r9d
mov    %ecx,%r11d
movapd %xmm4,%xmm3
jmp    1259 <gaussj+0x129>
nopw   0x0(%rax,%rax,1)
mov    %r9,0x60(%rsp)
movslq %r9d,%r13
addl   $0x1,0x80(%rsp,%r13,4)
cmp    %r13d,%r11d
mov    (%rsp),%rdi
je     13d0 <gaussj+0x2a0>
movslq %r11d,%r9
cmpq   $0x0,0x20(%rsp)
je     130a <gaussj+0x1da>
mov    %rbp,%rdx
add    %r13,%rdx
lea    (%r9,%rbp,1),%rdi
mov    (%rsp),%rsi
xor    %ecx,%ecx
nopl   0x0(%rax)
movsd  (%rsi,%r9,8),%xmm3
movsd  (%rsi,%r13,8),%xmm4
movsd  %xmm4,(%rsi,%r9,8)
movsd  %xmm3,(%rsi,%r13,8)
movsd  (%rsi,%rdi,8),%xmm3
movsd  (%rsi,%rdx,8),%xmm4
movsd  %xmm4,(%rsi,%rdi,8)
movsd  %xmm3,(%rsi,%rdx,8)
add    $0x2,%rcx
add    %r8,%rsi
cmp    %rcx,%r14
jne    12d0 <gaussj+0x1a0>
jmp    130c <gaussj+0x1dc>
xor    %ecx,%ecx
test   $0x1,%bpl
mov    (%rsp),%rdi
je     1335 <gaussj+0x205>
imul   %rbp,%rcx
lea    (%rcx,%r9,1),%rax
movsd  (%rdi,%rax,8),%xmm3
add    %r13,%rcx
movsd  (%rdi,%rcx,8),%xmm4
movsd  %xmm4,(%rdi,%rax,8)
movsd  %xmm3,(%rdi,%rcx,8)
test   %r15d,%r15d
jle    13d0 <gaussj+0x2a0>
cmpq   $0x0,0x30(%rsp)
je     139a <gaussj+0x26a>
mov    0x10(%rsp),%rax
lea    (%rax,%r13,1),%rdx
lea    (%rax,%r9,1),%rdi
mov    0x18(%rsp),%rsi
xor    %ecx,%ecx
mov    0x40(%rsp),%rax
nop
movsd  (%rsi,%r9,8),%xmm3
movsd  (%rsi,%r13,8),%xmm4
movsd  %xmm4,(%rsi,%r9,8)
movsd  %xmm3,(%rsi,%r13,8)
movsd  (%rsi,%rdi,8),%xmm3
movsd  (%rsi,%rdx,8),%xmm4
movsd  %xmm4,(%rsi,%rdi,8)
movsd  %xmm3,(%rsi,%rdx,8)
add    $0x2,%rcx
add    %r10,%rsi
cmp    %rcx,%rax
jne    1360 <gaussj+0x230>
jmp    139c <gaussj+0x26c>
xor    %ecx,%ecx
testb  $0x1,0x38(%rsp)
mov    (%rsp),%rdi
je     13d0 <gaussj+0x2a0>
imul   0x10(%rsp),%rcx
add    %rcx,%r9
mov    0x18(%rsp),%rax
movsd  (%rax,%r9,8),%xmm3
add    %r13,%rcx
movsd  (%rax,%rcx,8),%xmm4
movsd  %xmm4,(%rax,%r9,8)
movsd  %xmm3,(%rax,%rcx,8)
xchg   %ax,%ax
mov    %r11d,0x220(%rsp,%r12,4)
mov    0x60(%rsp),%r9
mov    %r9d,0x150(%rsp,%r12,4)
mov    %r9d,%esi
imul   0xc(%rsp),%esi
lea    (%rsi,%r9,1),%eax
movslq %eax,%rcx
movsd  (%rdi,%rcx,8),%xmm4
ucomisd %xmm0,%xmm4
jne    1405 <gaussj+0x2d5>
jnp    18f2 <gaussj+0x7c2>
movapd %xmm1,%xmm3
divsd  %xmm4,%xmm3
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rdi,%rcx,8)
cmpq   $0x0,0x20(%rsp)
je     1460 <gaussj+0x330>
mov    %rbp,%rdx
add    %r13,%rdx
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%r13,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%r13,8)
movsd  (%rdi,%rdx,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%rdx,8)
add    $0x2,%rcx
add    %r8,%rdi
cmp    %rcx,%r14
jne    1430 <gaussj+0x300>
jmp    1462 <gaussj+0x332>
nopl   0x0(%rax)
xor    %ecx,%ecx
test   $0x1,%bpl
mov    (%rsp),%rdi
je     1481 <gaussj+0x351>
imul   %rbp,%rcx
add    %r13,%rcx
movsd  (%rdi,%rcx,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%rcx,8)
mov    %r11d,0x5c(%rsp)
mov    %r12,0x78(%rsp)
test   %r15d,%r15d
jle    1505 <gaussj+0x3d5>
cmpq   $0x0,0x30(%rsp)
je     14e0 <gaussj+0x3b0>
mov    0x10(%rsp),%rax
lea    (%rax,%r13,1),%rdx
mov    0x18(%rsp),%rdi
xor    %ecx,%ecx
mov    0x40(%rsp),%rax
nopl   (%rax)
movsd  (%rdi,%r13,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%r13,8)
movsd  (%rdi,%rdx,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%rdx,8)
add    $0x2,%rcx
add    %r10,%rdi
cmp    %rcx,%rax
jne    14b0 <gaussj+0x380>
jmp    14e2 <gaussj+0x3b2>
nopl   0x0(%rax)
xor    %ecx,%ecx
testb  $0x1,0x38(%rsp)
mov    (%rsp),%rdi
je     1505 <gaussj+0x3d5>
imul   0x10(%rsp),%rcx
add    %r13,%rcx
mov    0x18(%rsp),%rax
mulsd  (%rax,%rcx,8),%xmm3
movsd  %xmm3,(%rax,%rcx,8)
movslq %esi,%rax
mov    %rax,0x48(%rsp)
mov    %r9d,%eax
mov    %rax,0x50(%rsp)
mov    %rbp,%r15
add    %r13,%r15
mov    0x10(%rsp),%rax
lea    (%rax,%r13,1),%r11
mov    0x70(%rsp),%r12
xor    %r9d,%r9d
mov    %rbx,%rsi
xor    %r10d,%r10d
jmp    1559 <gaussj+0x429>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r10
add    $0x8,%rsi
add    $0x8,%r9
add    $0x8,%r12
cmp    %rbp,%r10
je     11f0 <gaussj+0xc0>
cmp    0x50(%rsp),%r10
je     1540 <gaussj+0x410>
mov    0x48(%rsp),%rax
add    %r10,%rax
movsd  (%rdi,%rax,8),%xmm3
movq   $0x0,(%rdi,%rax,8)
cmpq   $0x0,0x20(%rsp)
je     1662 <gaussj+0x532>
mov    %rdi,%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rdx,%r9,1),%xmm4
movsd  (%rdx,%r13,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,(%rdx,%r9,1)
movsd  (%rdx,%rsi,1),%xmm4
movsd  (%rdx,%r15,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,(%rdx,%rsi,1)
add    $0x2,%rcx
add    %r8,%rdx
cmp    %rcx,%r14
jne    1590 <gaussj+0x460>
test   $0x1,%bpl
je     15f6 <gaussj+0x4c6>
imul   %rbp,%rcx
lea    (%rcx,%r13,1),%rax
add    %r10,%rcx
movsd  (%rdi,%rcx,8),%xmm4
movsd  (%rdi,%rax,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,(%rdi,%rcx,8)
cmpl   $0x0,0x28(%rsp)
mov    0x40(%rsp),%rax
jle    1540 <gaussj+0x410>
cmpq   $0x0,0x30(%rsp)
je     1670 <gaussj+0x540>
mov    0x18(%rsp),%rdx
xor    %ecx,%ecx
mov    0x68(%rsp),%rdi
nopw   0x0(%rax,%rax,1)
movsd  (%rdx,%r9,1),%xmm4
movsd  (%rdx,%r13,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,(%rdx,%r9,1)
movsd  (%rdx,%r12,1),%xmm4
movsd  (%rdx,%r11,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm5,%xmm4
movsd  %xmm4,(%rdx,%r12,1)
add    $0x2,%rcx
add    %rdi,%rdx
cmp    %rcx,%rax
jne    1620 <gaussj+0x4f0>
jmp    1672 <gaussj+0x542>
xor    %ecx,%ecx
test   $0x1,%bpl
jne    15d4 <gaussj+0x4a4>
jmp    15f6 <gaussj+0x4c6>
xor    %ecx,%ecx
testb  $0x1,0x38(%rsp)
mov    (%rsp),%rdi
je     1540 <gaussj+0x410>
imul   0x10(%rsp),%rcx
lea    (%rcx,%r13,1),%rax
add    %r10,%rcx
mov    0x18(%rsp),%rdx
movsd  (%rdx,%rcx,8),%xmm4
mulsd  (%rdx,%rax,8),%xmm3
subsd  %xmm3,%xmm4
movsd  %xmm4,(%rdx,%rcx,8)
jmp    1540 <gaussj+0x410>
cmpl   $0x0,0xc(%rsp)
jle    18e0 <gaussj+0x7b0>
mov    (%rsp),%rax
lea    (%rax,%rbp,8),%r15
mov    %ebp,%ecx
and    $0xfffffffc,%ecx
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x50(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %ebp,%r14d
and    $0xfffffffc,%r14d
mov    %rdx,%rcx
mov    %rdx,0x20(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x28(%rsp)
mov    %rbp,%r13
neg    %r13
lea    0x30(%rax),%rcx
mov    %rcx,0x48(%rsp)
lea    0x8(%rax),%r12
mov    %rbp,%rdx
jmp    171a <gaussj+0x5ea>
nopw   0x0(%rax,%rax,1)
cmp    $0x1,%rbx
jle    18e0 <gaussj+0x7b0>
mov    %rdx,%rbx
add    $0xffffffffffffffff,%rdx
mov    %edx,%ecx
mov    0x220(%rsp,%rcx,4),%eax
mov    0x150(%rsp,%rcx,4),%ecx
cmp    %ecx,%eax
je     1710 <gaussj+0x5e0>
mov    0xc(%rsp),%esi
imul   %esi,%eax
imul   %esi,%ecx
movslq %eax,%r10
movslq %ecx,%r11
cmp    $0x4,%esi
jb     1770 <gaussj+0x640>
mov    (%rsp),%rsi
lea    (%rsi,%r10,8),%rcx
lea    (%r15,%r11,8),%rdi
cmp    %rdi,%rcx
jae    17fa <gaussj+0x6ca>
lea    (%r15,%r10,8),%rcx
lea    (%rsi,%r11,8),%rdi
cmp    %rcx,%rdi
jae    17fa <gaussj+0x6ca>
xor    %edi,%edi
mov    %rdi,%rcx
test   $0x1,%bpl
je     17a2 <gaussj+0x672>
lea    (%rdi,%r10,1),%rcx
mov    (%rsp),%rax
movsd  (%rax,%rcx,8),%xmm0
lea    (%rdi,%r11,1),%rsi
movsd  (%rax,%rsi,8),%xmm1
movsd  %xmm1,(%rax,%rcx,8)
movsd  %xmm0,(%rax,%rsi,8)
mov    %rdi,%rcx
or     $0x1,%rcx
not    %rdi
cmp    %r13,%rdi
je     1710 <gaussj+0x5e0>
lea    (%r12,%r11,8),%rax
lea    (%r12,%r10,8),%rsi
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%rsi,%rcx,8),%xmm0
movsd  -0x8(%rax,%rcx,8),%xmm1
movsd  %xmm1,-0x8(%rsi,%rcx,8)
movsd  %xmm0,-0x8(%rax,%rcx,8)
movsd  (%rsi,%rcx,8),%xmm0
movsd  (%rax,%rcx,8),%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
movsd  %xmm0,(%rax,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rbp
jne    17c0 <gaussj+0x690>
jmp    1710 <gaussj+0x5e0>
cmpq   $0x0,0x50(%rsp)
je     18d5 <gaussj+0x7a5>
mov    0x48(%rsp),%rax
lea    (%rax,%r10,8),%rdi
lea    (%rax,%r11,8),%r8
mov    0x28(%rsp),%r9
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movups -0x30(%rdi,%rcx,8),%xmm0
movups -0x20(%rdi,%rcx,8),%xmm1
movups -0x30(%r8,%rcx,8),%xmm2
movups -0x20(%r8,%rcx,8),%xmm3
movups %xmm2,-0x30(%rdi,%rcx,8)
movups %xmm3,-0x20(%rdi,%rcx,8)
movups %xmm0,-0x30(%r8,%rcx,8)
movups %xmm1,-0x20(%r8,%rcx,8)
movupd -0x10(%rdi,%rcx,8),%xmm0
movupd (%rdi,%rcx,8),%xmm1
movupd -0x10(%r8,%rcx,8),%xmm2
movupd (%r8,%rcx,8),%xmm3
movupd %xmm2,-0x10(%rdi,%rcx,8)
movupd %xmm3,(%rdi,%rcx,8)
movupd %xmm0,-0x10(%r8,%rcx,8)
movupd %xmm1,(%r8,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%r9
jne    1820 <gaussj+0x6f0>
testb  $0x1,0x20(%rsp)
je     18c4 <gaussj+0x794>
lea    (%rcx,%r10,1),%rdi
mov    (%rsp),%rsi
movupd (%rsi,%rdi,8),%xmm0
movupd 0x10(%rsi,%rdi,8),%xmm1
add    %r11,%rcx
movupd (%rsi,%rcx,8),%xmm2
movupd 0x10(%rsi,%rcx,8),%xmm3
movupd %xmm2,(%rsi,%rdi,8)
movupd %xmm3,0x10(%rsi,%rdi,8)
movupd %xmm0,(%rsi,%rcx,8)
movupd %xmm1,0x10(%rsi,%rcx,8)
mov    %r14,%rdi
cmp    %rbp,%r14
je     1710 <gaussj+0x5e0>
jmp    1772 <gaussj+0x642>
xor    %ecx,%ecx
testb  $0x1,0x20(%rsp)
jne    188d <gaussj+0x75d>
jmp    18c4 <gaussj+0x794>
add    $0x2e8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0x71f(%rip),%rdi
call   1030 <puts@plt>
mov    $0x63,%edi
call   1050 <exit@plt>

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

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