void matinv(a, n) double a[];
unsigned int n;
{
  int indxc[100], indxr[100], ipiv[100];
  int i, icol, irow, j, k, l, ll;
  double big, dum, pivinv, temp;
  if (n > 100)
    mexErrMsgTxt("Matinv: matrix too big");
  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)
            mexErrMsgTxt("matinv: singular-1");
        }
    ++(ipiv[icol]);
    if (irow != icol) {
      for (l = 0; l < n; l++) {
        temp = (a[irow + l * n]);
        (a[irow + l * n]) = (a[icol + l * n]);
        (a[icol + l * n]) = temp;
      }
    }
    indxr[i] = irow;
    indxc[i] = icol;
    if (a[icol + icol * n] == 0.0)
      mexErrMsgTxt("matinv: singular-2");
    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 (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 = n - 1; l >= 0; l--) {
    if (indxr[l] != indxc[l])
      for (k = 0; k < n; k++) {
        temp = (a[k + indxr[l] * n]);
        (a[k + indxr[l] * n]) = (a[k + indxc[l] * n]);
        (a[k + indxc[l] * n]) = temp;
      };
  }
}


x86 assembly (clang, -O0):
<matinv>:
push   %rbp
mov    %rsp,%rbp
sub    $0x500,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
cmpl   $0x64,-0xc(%rbp)
jbe    113a <matinv+0x2a>
lea    0xeed(%rip),%rdi
mov    $0x0,%al
call   1030 <mexErrMsgTxt@plt>
movl   $0x0,-0x4d0(%rbp)
mov    -0x4d0(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    1179 <matinv+0x69>
movslq -0x4d0(%rbp),%rax
movl   $0x0,-0x4c0(%rbp,%rax,4)
mov    -0x4d0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d0(%rbp)
jmp    1144 <matinv+0x34>
movl   $0x0,-0x4c4(%rbp)
mov    -0x4c4(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    1604 <matinv+0x4f4>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x4e8(%rbp)
ff
movl   $0x0,-0x4d0(%rbp)
mov    -0x4d0(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    12db <matinv+0x1cb>
movslq -0x4d0(%rbp),%rax
cmpl   $0x1,-0x4c0(%rbp,%rax,4)
01
je     12c2 <matinv+0x1b2>
movl   $0x0,-0x4d4(%rbp)
mov    -0x4d4(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    12bd <matinv+0x1ad>
movslq -0x4d4(%rbp),%rax
cmpl   $0x0,-0x4c0(%rbp,%rax,4)
jne    127c <matinv+0x16c>
mov    -0x8(%rbp),%rax
mov    -0x4d0(%rbp),%ecx
mov    -0x4d4(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xdf3(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x4e8(%rbp),%xmm0
ff
jb     1277 <matinv+0x167>
mov    -0x8(%rbp),%rax
mov    -0x4d0(%rbp),%ecx
mov    -0x4d4(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xdbd(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x4e8(%rbp)
ff
mov    -0x4d0(%rbp),%eax
mov    %eax,-0x4cc(%rbp)
mov    -0x4d4(%rbp),%eax
mov    %eax,-0x4c8(%rbp)
jmp    12a4 <matinv+0x194>
movslq -0x4d4(%rbp),%rax
cmpl   $0x1,-0x4c0(%rbp,%rax,4)
01
jle    129f <matinv+0x18f>
lea    0xd9f(%rip),%rdi
mov    $0x0,%al
call   1030 <mexErrMsgTxt@plt>
jmp    12a4 <matinv+0x194>
jmp    12a9 <matinv+0x199>
mov    -0x4d4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d4(%rbp)
jmp    11d5 <matinv+0xc5>
jmp    12c2 <matinv+0x1b2>
jmp    12c7 <matinv+0x1b7>
mov    -0x4d0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d0(%rbp)
jmp    11a7 <matinv+0x97>
movslq -0x4c8(%rbp),%rax
mov    -0x4c0(%rbp,%rax,4),%ecx
add    $0x1,%ecx
mov    %ecx,-0x4c0(%rbp,%rax,4)
mov    -0x4cc(%rbp),%eax
cmp    -0x4c8(%rbp),%eax
je     13bb <matinv+0x2ab>
movl   $0x0,-0x4d8(%rbp)
mov    -0x4d8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    13b6 <matinv+0x2a6>
mov    -0x8(%rbp),%rax
mov    -0x4cc(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x500(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x4cc(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x500(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4d8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d8(%rbp)
jmp    130f <matinv+0x1ff>
jmp    13bb <matinv+0x2ab>
mov    -0x4cc(%rbp),%ecx
movslq -0x4c4(%rbp),%rax
mov    %ecx,-0x330(%rbp,%rax,4)
mov    -0x4c8(%rbp),%ecx
movslq -0x4c4(%rbp),%rax
mov    %ecx,-0x1a0(%rbp,%rax,4)
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4c8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1421 <matinv+0x311>
jp     1421 <matinv+0x311>
lea    0xc30(%rip),%rdi
mov    $0x0,%al
call   1030 <mexErrMsgTxt@plt>
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4c8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  0xbbf(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x4f8(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4c8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  0xb92(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x4d8(%rbp)
mov    -0x4d8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    14ca <matinv+0x3ba>
movsd  -0x4f8(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4d8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d8(%rbp)
jmp    147d <matinv+0x36d>
movl   $0x0,-0x4dc(%rbp)
mov    -0x4dc(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    15eb <matinv+0x4db>
mov    -0x4dc(%rbp),%eax
cmp    -0x4c8(%rbp),%eax
je     15d2 <matinv+0x4c2>
mov    -0x8(%rbp),%rax
mov    -0x4dc(%rbp),%ecx
mov    -0x4c8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x4f0(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x4dc(%rbp),%ecx
mov    -0x4c8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x4d8(%rbp)
mov    -0x4d8(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    15cd <matinv+0x4bd>
mov    -0x8(%rbp),%rax
mov    -0x4c8(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x4f0(%rbp),%xmm2
ff
mov    -0x8(%rbp),%rax
mov    -0x4dc(%rbp),%ecx
mov    -0x4d8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
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    -0x4d8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d8(%rbp)
jmp    1544 <matinv+0x434>
jmp    15d2 <matinv+0x4c2>
jmp    15d7 <matinv+0x4c7>
mov    -0x4dc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4dc(%rbp)
jmp    14d4 <matinv+0x3c4>
jmp    15f0 <matinv+0x4e0>
mov    -0x4c4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c4(%rbp)
jmp    1183 <matinv+0x73>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4d8(%rbp)
cmpl   $0x0,-0x4d8(%rbp)
jl     172e <matinv+0x61e>
movslq -0x4d8(%rbp),%rax
mov    -0x330(%rbp,%rax,4),%eax
movslq -0x4d8(%rbp),%rcx
cmp    -0x1a0(%rbp,%rcx,4),%eax
je     1715 <matinv+0x605>
movl   $0x0,-0x4d4(%rbp)
mov    -0x4d4(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    1710 <matinv+0x600>
mov    -0x8(%rbp),%rax
mov    -0x4d4(%rbp),%ecx
movslq -0x4d8(%rbp),%rdx
mov    -0x330(%rbp,%rdx,4),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x500(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x4d4(%rbp),%ecx
movslq -0x4d8(%rbp),%rdx
mov    -0x1a0(%rbp,%rdx,4),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x4d4(%rbp),%ecx
movslq -0x4d8(%rbp),%rdx
mov    -0x330(%rbp,%rdx,4),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x500(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x4d4(%rbp),%ecx
movslq -0x4d8(%rbp),%rdx
mov    -0x1a0(%rbp,%rdx,4),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4d4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4d4(%rbp)
jmp    1649 <matinv+0x539>
jmp    1715 <matinv+0x605>
jmp    171a <matinv+0x60a>
mov    -0x4d8(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4d8(%rbp)
jmp    1610 <matinv+0x500>
add    $0x500,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matinv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x4f8,%rsp
mov    %esi,%r14d
mov    %rdi,%rbx
cmp    $0x64,%esi
jbe    114c <matinv+0x2c>
lea    0xedd(%rip),%rdi
xor    %eax,%eax
call   1040 <mexErrMsgTxt@plt>
jmp    1155 <matinv+0x35>
test   %r14d,%r14d
je     1732 <matinv+0x612>
mov    %r14d,%ebp
lea    0x0(,%rbp,4),%rdx
lea    0x40(%rsp),%rdi
xor    %esi,%esi
call   1030 <memset@plt>
lea    -0x1(%rbp),%rax
mov    %rax,0x30(%rsp)
mov    %ebp,%r13d
and    $0xfffffffe,%r13d
lea    (%r14,%r14,1),%ecx
xor    %eax,%eax
mov    %rax,0x18(%rsp)
mov    %rax,0x8(%rsp)
mov    %rax,0x28(%rsp)
mov    %rbp,0x20(%rsp)
mov    %r13,0x38(%rsp)
mov    %ecx,0x4(%rsp)
jmp    11ca <matinv+0xaa>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x18(%rsp),%rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,0x18(%rsp)
cmp    %rbp,%rdx
je     14bb <matinv+0x39b>
xor    %r12d,%r12d
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
jmp    11e9 <matinv+0xc9>
nopl   0x0(%rax)
add    $0x1,%r12
cmp    %rbp,%r12
je     1260 <matinv+0x140>
cmpl   $0x1,0x40(%rsp,%r12,4)
je     11e0 <matinv+0xc0>
mov    %r12d,%r15d
xor    %r13d,%r13d
jmp    120c <matinv+0xec>
nopl   0x0(%rax)
add    $0x1,%r13
add    %r14d,%r15d
cmp    %r13,%rbp
je     11e0 <matinv+0xc0>
mov    0x40(%rsp,%r13,4),%eax
test   %eax,%eax
je     1230 <matinv+0x110>
cmp    $0x2,%eax
jl     1200 <matinv+0xe0>
lea    0xe16(%rip),%rdi
xor    %eax,%eax
call   1040 <mexErrMsgTxt@plt>
jmp    1200 <matinv+0xe0>
nopw   0x0(%rax,%rax,1)
mov    %r15d,%eax
movsd  (%rbx,%rax,8),%xmm0
andpd  0xdc0(%rip),%xmm0
ucomisd 0x10(%rsp),%xmm0
jb     1200 <matinv+0xe0>
mov    %r12d,%eax
mov    %rax,0x28(%rsp)
movsd  %xmm0,0x10(%rsp)
mov    %r13d,%eax
mov    %rax,0x8(%rsp)
jmp    1200 <matinv+0xe0>
mov    0x8(%rsp),%r10
movslq %r10d,%rax
addl   $0x1,0x40(%rsp,%rax,4)
mov    0x28(%rsp),%r9
cmp    %eax,%r9d
mov    0x30(%rsp),%r12
mov    0x38(%rsp),%r13
je     1308 <matinv+0x1e8>
test   %r12,%r12
je     14aa <matinv+0x38a>
mov    %r9d,%ecx
mov    %r10d,%edx
xor    %eax,%eax
mov    0x4(%rsp),%r8d
nopl   0x0(%rax,%rax,1)
mov    %ecx,%esi
movsd  (%rbx,%rsi,8),%xmm0
mov    %edx,%edi
movsd  (%rbx,%rdi,8),%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
movsd  %xmm0,(%rbx,%rdi,8)
lea    (%r14,%rcx,1),%esi
movsd  (%rbx,%rsi,8),%xmm0
lea    (%r14,%rdx,1),%edi
movsd  (%rbx,%rdi,8),%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
movsd  %xmm0,(%rbx,%rdi,8)
add    $0x2,%rax
add    %r8d,%edx
add    %r8d,%ecx
cmp    %rax,%r13
jne    12a0 <matinv+0x180>
test   $0x1,%bpl
je     1308 <matinv+0x1e8>
imul   %r14d,%eax
lea    (%rax,%r9,1),%ecx
movsd  (%rbx,%rcx,8),%xmm0
add    %r10d,%eax
movsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
movsd  %xmm0,(%rbx,%rax,8)
mov    0x18(%rsp),%rax
mov    %r9d,0x1d0(%rsp,%rax,4)
mov    %r10d,0x360(%rsp,%rax,4)
mov    %r10d,%r11d
imul   %r14d,%r11d
lea    (%r11,%r10,1),%r15d
movsd  (%rbx,%r15,8),%xmm1
ucomisd 0xce2(%rip),%xmm1
jne    1358 <matinv+0x238>
jp     1358 <matinv+0x238>
lea    0xd09(%rip),%rdi
xor    %eax,%eax
mov    %r11,0x10(%rsp)
call   1040 <mexErrMsgTxt@plt>
mov    0x10(%rsp),%r11
movsd  (%rbx,%r15,8),%xmm1
movsd  0xcb0(%rip),%xmm0
divsd  %xmm1,%xmm0
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rbx,%r15,8)
test   %r12,%r12
je     1490 <matinv+0x370>
mov    0x8(%rsp),%r10
mov    %r10d,%esi
xor    %eax,%eax
mov    0x4(%rsp),%ecx
nopl   0x0(%rax)
mov    %esi,%edx
movsd  (%rbx,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
lea    (%r14,%rsi,1),%edx
movsd  (%rbx,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
add    $0x2,%rax
add    %ecx,%esi
cmp    %rax,%r13
jne    1390 <matinv+0x270>
test   $0x1,%bpl
je     13d4 <matinv+0x2b4>
imul   %r14d,%eax
add    %r10d,%eax
mulsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
mov    %r10d,%r8d
xor    %r9d,%r9d
jmp    13ed <matinv+0x2cd>
nopl   0x0(%rax)
add    $0x1,%r9
cmp    %rbp,%r9
je     11b0 <matinv+0x90>
cmp    %r8,%r9
je     13e0 <matinv+0x2c0>
lea    (%r11,%r9,1),%edx
movsd  (%rbx,%rdx,8),%xmm0
movq   $0x0,(%rbx,%rdx,8)
test   %r12,%r12
je     1459 <matinv+0x339>
mov    %r10d,%esi
mov    %r9d,%edi
xor    %edx,%edx
mov    %esi,%eax
mov    %edi,%ebp
movsd  (%rbx,%rbp,8),%xmm1
movsd  (%rbx,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rbx,%rbp,8)
lea    (%r14,%rsi,1),%eax
lea    (%r14,%rdi,1),%ebp
movsd  (%rbx,%rbp,8),%xmm1
movsd  (%rbx,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rbx,%rbp,8)
add    $0x2,%rdx
add    %ecx,%edi
add    %ecx,%esi
cmp    %rdx,%r13
jne    1410 <matinv+0x2f0>
jmp    145b <matinv+0x33b>
xor    %edx,%edx
mov    0x20(%rsp),%rbp
test   $0x1,%bpl
je     13e0 <matinv+0x2c0>
imul   %r14d,%edx
lea    (%rdx,%r10,1),%eax
add    %r9d,%edx
movsd  (%rbx,%rdx,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
jmp    13e0 <matinv+0x2c0>
nopl   (%rax)
xor    %eax,%eax
mov    0x4(%rsp),%ecx
mov    0x8(%rsp),%r10
test   $0x1,%bpl
jne    13c3 <matinv+0x2a3>
jmp    13d4 <matinv+0x2b4>
xor    %eax,%eax
test   $0x1,%bpl
jne    12e9 <matinv+0x1c9>
jmp    1308 <matinv+0x1e8>
mov    %r14d,%r13d
add    $0xffffffff,%r13d
js     1732 <matinv+0x612>
mov    0x20(%rsp),%r12
lea    (%rbx,%r12,8),%r15
mov    %r12d,%eax
and    $0xfffffffc,%eax
mov    %rax,0x10(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,0x8(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,0x28(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x18(%rsp)
neg    %r12
jmp    151a <matinv+0x3fa>
cs nopw 0x0(%rax,%rax,1)
add    $0xffffffff,%r13d
js     1732 <matinv+0x612>
mov    %r13d,%eax
mov    0x1d0(%rsp,%rax,4),%ebp
mov    0x360(%rsp,%rax,4),%ecx
cmp    %ecx,%ebp
je     1510 <matinv+0x3f0>
imul   %r14d,%ebp
imul   %r14d,%ecx
cmp    $0x4,%r14d
jae    15e0 <matinv+0x4c0>
xor    %eax,%eax
mov    %rax,%rdx
testb  $0x1,0x20(%rsp)
je     156e <matinv+0x44e>
lea    (%rax,%rbp,1),%edx
movsd  (%rbx,%rdx,8),%xmm0
lea    (%rcx,%rax,1),%esi
movsd  (%rbx,%rsi,8),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
movsd  %xmm0,(%rbx,%rsi,8)
mov    %rax,%rdx
or     $0x1,%rdx
not    %rax
cmp    %r12,%rax
je     1510 <matinv+0x3f0>
add    %edx,%ebp
add    %edx,%ecx
mov    0x20(%rsp),%rax
sub    %rdx,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %ebp,%edx
movsd  (%rbx,%rdx,8),%xmm0
mov    %ecx,%esi
movsd  (%rbx,%rsi,8),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
movsd  %xmm0,(%rbx,%rsi,8)
lea    0x1(%rbp),%edx
movsd  (%rbx,%rdx,8),%xmm0
lea    0x1(%rcx),%esi
movsd  (%rbx,%rsi,8),%xmm1
movsd  %xmm1,(%rbx,%rdx,8)
movsd  %xmm0,(%rbx,%rsi,8)
add    $0x2,%ebp
add    $0x2,%ecx
add    $0xfffffffffffffffe,%rax
jne    1590 <matinv+0x470>
jmp    1510 <matinv+0x3f0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %ecx,%eax
mov    0x30(%rsp),%rsi
add    %esi,%eax
setb   %al
mov    %ebp,%edx
add    %esi,%edx
jb     1541 <matinv+0x421>
mov    %rsi,%rdx
shr    $0x20,%rdx
jne    1541 <matinv+0x421>
test   %al,%al
jne    1541 <matinv+0x421>
mov    $0x0,%eax
test   %rdx,%rdx
jne    1543 <matinv+0x423>
mov    %ebp,%eax
lea    (%rbx,%rax,8),%rsi
mov    %ecx,%edx
lea    (%r15,%rdx,8),%rdi
cmp    %rdi,%rsi
jae    163b <matinv+0x51b>
lea    (%r15,%rax,8),%rax
lea    (%rbx,%rdx,8),%rdx
cmp    %rax,%rdx
jb     1541 <matinv+0x421>
cmpq   $0x0,0x8(%rsp)
je     1726 <matinv+0x606>
mov    0x18(%rsp),%r11
mov    %ecx,%esi
mov    %ebp,%eax
xor    %r10d,%r10d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %eax,%r8d
movups (%rbx,%r8,8),%xmm0
movups 0x10(%rbx,%r8,8),%xmm1
mov    %esi,%r9d
movups (%rbx,%r9,8),%xmm2
movups 0x10(%rbx,%r9,8),%xmm3
movups %xmm2,(%rbx,%r8,8)
movups %xmm3,0x10(%rbx,%r8,8)
movups %xmm0,(%rbx,%r9,8)
movups %xmm1,0x10(%rbx,%r9,8)
lea    0x4(%rax),%edx
movupd (%rbx,%rdx,8),%xmm0
movupd 0x10(%rbx,%rdx,8),%xmm1
lea    0x4(%rsi),%edi
movupd (%rbx,%rdi,8),%xmm2
movups 0x10(%rbx,%rdi,8),%xmm3
movupd %xmm2,(%rbx,%rdx,8)
movups %xmm3,0x10(%rbx,%rdx,8)
movupd %xmm0,(%rbx,%rdi,8)
movupd %xmm1,0x10(%rbx,%rdi,8)
add    $0x8,%r10
add    $0x8,%eax
add    $0x8,%esi
add    $0xfffffffffffffffe,%r11
jne    1660 <matinv+0x540>
testb  $0x1,0x28(%rsp)
je     170e <matinv+0x5ee>
lea    (%r10,%rbp,1),%eax
movupd (%rbx,%rax,8),%xmm0
movupd 0x10(%rbx,%rax,8),%xmm1
add    %ecx,%r10d
movupd (%rbx,%r10,8),%xmm2
movups 0x10(%rbx,%r10,8),%xmm3
movupd %xmm2,(%rbx,%rax,8)
movups %xmm3,0x10(%rbx,%rax,8)
movupd %xmm0,(%rbx,%r10,8)
movupd %xmm1,0x10(%rbx,%r10,8)
mov    0x10(%rsp),%rdx
mov    %rdx,%rax
cmp    0x20(%rsp),%rdx
je     1510 <matinv+0x3f0>
jmp    1543 <matinv+0x423>
xor    %r10d,%r10d
testb  $0x1,0x28(%rsp)
jne    16d9 <matinv+0x5b9>
jmp    170e <matinv+0x5ee>
add    $0x4f8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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