#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);
int fabs(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, ...);
int /*<<< orphan*/ stderr;
void hc_ludcmp_3x3(double a[3][3], int n, int *indx) {
  int i, imax = 0, j, k;
  double big, dum, sum, temp;
  double vv[3];
  for (i = 0; i < n; i++) {
    big = 0.0;
    for (j = 0; j < n; j++)
      if ((temp = fabs(a[i][j])) > big)
        big = temp;
    if (fabs(big) < 3e-16) {
      fprintf(stderr, "hc_ludcmp_3x3: singular matrix in routine, big: %g\n",
              (double)big);
      for (j = 0; j < n; j++) {
        for (k = 0; k < n; k++)
          fprintf(stderr, "%g ", (double)a[j][k]);
        fprintf(stderr, "\n");
      }
      exit(-1);
    }
    vv[i] = 1.0 / big;
  }
  for (j = 0; j < n; j++) {
    for (i = 0; i < j; i++) {
      sum = a[i][j];
      for (k = 0; k < i; k++)
        sum -= a[i][k] * a[k][j];
      a[i][j] = sum;
    }
    big = 0.0;
    for (i = j; i < n; i++) {
      sum = a[i][j];
      for (k = 0; k < j; k++)
        sum -= a[i][k] * a[k][j];
      a[i][j] = sum;
      if ((dum = vv[i] * fabs(sum)) >= big) {
        big = dum;
        imax = i;
      }
    }
    if (j != imax) {
      for (k = 0; k < n; k++) {
        dum = a[imax][k];
        a[imax][k] = a[j][k];
        a[j][k] = dum;
      }
      vv[imax] = vv[j];
    }
    indx[j] = imax;
    if (fabs(a[j][j]) < 3e-16)
      a[j][j] = 1.0e-20;
    ;
    if (j != 2) {
      dum = 1.0 / (a[j][j]);
      for (i = j + 1; i < n; i++)
        a[i][j] *= dum;
    }
  }
}


x86 assembly (clang, -O0):
<hc_ludcmp_3x3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12b9 <hc_ludcmp_3x3+0x189>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11c6 <hc_ludcmp_3x3+0x96>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x48(%rbp)
ucomisd -0x30(%rbp),%xmm0
jbe    11b3 <hc_ludcmp_3x3+0x83>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    11b8 <hc_ludcmp_3x3+0x88>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    116c <hc_ludcmp_3x3+0x3c>
movsd  -0x30(%rbp),%xmm0
call   1030 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xe24(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1294 <hc_ludcmp_3x3+0x164>
mov    0x2df3(%rip),%rax
mov    (%rax),%edi
movsd  -0x30(%rbp),%xmm0
lea    0xe1d(%rip),%rsi
mov    $0x1,%al
call   1040 <fprintf@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    128a <hc_ludcmp_3x3+0x15a>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1265 <hc_ludcmp_3x3+0x135>
mov    0x2db1(%rip),%rax
mov    (%rax),%edi
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xdfc(%rip),%rsi
mov    $0x1,%al
call   1040 <fprintf@plt>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    121c <hc_ludcmp_3x3+0xec>
mov    0x2d74(%rip),%rax
mov    (%rax),%edi
lea    0xdd5(%rip),%rsi
mov    $0x0,%al
call   1040 <fprintf@plt>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1209 <hc_ludcmp_3x3+0xd9>
mov    $0xffffffff,%edi
call   1050 <exit@plt>
movsd  0xd74(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
movslq -0x1c(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1151 <hc_ludcmp_3x3+0x21>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1660 <hc_ludcmp_3x3+0x530>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13a1 <hc_ludcmp_3x3+0x271>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1376 <hc_ludcmp_3x3+0x246>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1303 <hc_ludcmp_3x3+0x1d3>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    12d3 <hc_ludcmp_3x3+0x1a3>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14c8 <hc_ludcmp_3x3+0x398>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1452 <hc_ludcmp_3x3+0x322>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    13df <hc_ludcmp_3x3+0x2af>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x1c(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1030 <fabs@plt>
movsd  -0x68(%rbp),%xmm0
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
ucomisd -0x30(%rbp),%xmm0
jb     14b5 <hc_ludcmp_3x3+0x385>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
jmp    14ba <hc_ludcmp_3x3+0x38a>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13af <hc_ludcmp_3x3+0x27f>
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
je     1573 <hc_ludcmp_3x3+0x443>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    155f <hc_ludcmp_3x3+0x42f>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    14db <hc_ludcmp_3x3+0x3ab>
movslq -0x24(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm0
movslq -0x20(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0x20(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xa56(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    15d4 <hc_ludcmp_3x3+0x4a4>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  0xa39(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x2,-0x24(%rbp)
je     164d <hc_ludcmp_3x3+0x51d>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  0xa17(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1648 <hc_ludcmp_3x3+0x518>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    160c <hc_ludcmp_3x3+0x4dc>
jmp    164d <hc_ludcmp_3x3+0x51d>
jmp    1652 <hc_ludcmp_3x3+0x522>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12c0 <hc_ludcmp_3x3+0x190>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<hc_ludcmp_3x3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xb8,%rsp
mov    %rdx,0x70(%rsp)
mov    %esi,0xc(%rsp)
test   %esi,%esi
jle    17b0 <hc_ludcmp_3x3+0x680>
mov    %rdi,%r13
mov    0xc(%rsp),%r12d
mov    %r12d,%r14d
and    $0xfffffffe,%r14d
lea    0x8(%rdi),%rbp
xor    %r15d,%r15d
mov    %r12,0x28(%rsp)
nopl   (%rax)
cmpl   $0x1,0xc(%rsp)
jne    1180 <hc_ludcmp_3x3+0x50>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
jmp    11d5 <hc_ludcmp_3x3+0xa5>
nop
xorps  %xmm0,%xmm0
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm0,0x10(%rsp)
movsd  -0x8(%rbp,%rbx,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
maxsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x0(%rbp,%rbx,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
maxsd  0x10(%rsp),%xmm0
add    $0x2,%rbx
cmp    %rbx,%r14
jne    1190 <hc_ludcmp_3x3+0x60>
test   $0x1,%r12b
je     1207 <hc_ludcmp_3x3+0xd7>
lea    (%r15,%r15,2),%rax
lea    0x0(,%rax,8),%rax
add    %r13,%rax
movsd  %xmm0,0x10(%rsp)
movsd  (%rax,%rbx,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
maxsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
call   1030 <fabs@plt>
test   %eax,%eax
jle    17c2 <hc_ludcmp_3x3+0x692>
movsd  0xdde(%rip),%xmm0
divsd  0x10(%rsp),%xmm0
movsd  %xmm0,0xa0(%rsp,%r15,8)
add    $0x1,%r15
add    $0x18,%rbp
cmp    %r12,%r15
jne    1170 <hc_ludcmp_3x3+0x40>
cmpl   $0x0,0xc(%rsp)
jle    17b0 <hc_ludcmp_3x3+0x680>
lea    0x0(,%r12,8),%rax
add    %r13,%rax
mov    %r12d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0x58(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x50(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rdx,%rcx
mov    %rdx,0x48(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x40(%rsp)
mov    %r12,%rcx
mov    %rax,%r12
neg    %rcx
mov    %rcx,0x60(%rsp)
lea    0x8(%r13),%rax
mov    %rax,0x20(%rsp)
lea    0x30(%r13),%rax
mov    %rax,0x38(%rsp)
lea    0x18(%r13),%rsi
xor    %r15d,%r15d
mov    $0x1,%edi
xor    %ebp,%ebp
mov    %r12,0x68(%rsp)
jmp    12d4 <hc_ludcmp_3x3+0x1a4>
xchg   %ax,%ax
add    $0x1,%rdi
add    $0x18,%rsi
mov    %rax,%r15
cmp    %rbx,%rax
je     17b0 <hc_ludcmp_3x3+0x680>
mov    %rdi,0x80(%rsp)
mov    %rsi,0x88(%rsp)
test   %r15,%r15
je     13c0 <hc_ludcmp_3x3+0x290>
lea    0x0(,%r15,8),%r8
add    %r13,%r8
mov    0x20(%rsp),%rcx
xor    %ebx,%ebx
jmp    1329 <hc_ludcmp_3x3+0x1f9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%rsi,%r15,8),%rax
movsd  %xmm0,(%rax)
add    $0x1,%rbx
add    $0x18,%rcx
cmp    %r15,%rbx
je     13c0 <hc_ludcmp_3x3+0x290>
lea    (%rbx,%rbx,2),%rax
lea    0x0(,%rax,8),%rsi
add    %r13,%rsi
movsd  (%rsi,%r15,8),%xmm0
test   %rbx,%rbx
je     1310 <hc_ludcmp_3x3+0x1e0>
cmp    $0x1,%rbx
jne    1350 <hc_ludcmp_3x3+0x220>
xor    %edi,%edi
jmp    1389 <hc_ludcmp_3x3+0x259>
nopl   (%rax)
mov    %rbx,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %r8,%rax
xor    %edi,%edi
nopl   0x0(%rax)
movsd  -0x8(%rcx,%rdi,8),%xmm1
mulsd  (%rax),%xmm1
movsd  (%rcx,%rdi,8),%xmm2
subsd  %xmm1,%xmm0
mulsd  0x18(%rax),%xmm2
subsd  %xmm2,%xmm0
add    $0x2,%rdi
add    $0x30,%rax
cmp    %rdi,%rdx
jne    1360 <hc_ludcmp_3x3+0x230>
test   $0x1,%bl
je     1310 <hc_ludcmp_3x3+0x1e0>
movsd  (%rsi,%rdi,8),%xmm1
lea    (%rdi,%rdi,2),%rax
lea    0x0(,%rax,8),%rax
add    %r13,%rax
mulsd  (%rax,%r15,8),%xmm1
subsd  %xmm1,%xmm0
jmp    1310 <hc_ludcmp_3x3+0x1e0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r15,%r12
and    $0xfffffffffffffffe,%r12
lea    (%r15,%r15,2),%rcx
mov    0x20(%rsp),%rax
mov    %rcx,0x78(%rsp)
lea    (%rax,%rcx,8),%rax
mov    %r13,%rbx
lea    0x0(%r13,%r15,8),%rcx
mov    %rcx,0x90(%rsp)
xorpd  %xmm1,%xmm1
mov    %rax,0x30(%rsp)
mov    %rax,%r13
mov    %r15,%r14
jmp    146a <hc_ludcmp_3x3+0x33a>
nopw   0x0(%rax,%rax,1)
lea    (%rax,%r15,8),%rax
movsd  %xmm0,(%rax)
movsd  0xa0(%rsp,%r14,8),%xmm1
movsd  %xmm1,0x98(%rsp)
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0x98(%rsp),%xmm0
movapd 0x10(%rsp),%xmm3
ucomisd %xmm3,%xmm0
cmovae %r14d,%ebp
movapd %xmm3,%xmm1
cmpnlesd %xmm0,%xmm1
movapd %xmm1,%xmm2
andnpd %xmm0,%xmm2
andpd  %xmm3,%xmm1
orpd   %xmm2,%xmm1
add    $0x1,%r14
add    $0x18,%r13
cmp    0x28(%rsp),%r14
je     1510 <hc_ludcmp_3x3+0x3e0>
lea    (%r14,%r14,2),%rax
lea    (%rbx,%rax,8),%rax
movsd  (%rax,%r15,8),%xmm0
test   %r15,%r15
movapd %xmm1,0x10(%rsp)
je     1400 <hc_ludcmp_3x3+0x2d0>
cmp    $0x1,%r15
jne    14a0 <hc_ludcmp_3x3+0x370>
xor    %ecx,%ecx
jmp    14dc <hc_ludcmp_3x3+0x3ac>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x90(%rsp),%rdx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%r13,%rcx,8),%xmm1
mulsd  (%rdx),%xmm1
movsd  0x0(%r13,%rcx,8),%xmm2
subsd  %xmm1,%xmm0
mulsd  0x18(%rdx),%xmm2
subsd  %xmm2,%xmm0
add    $0x2,%rcx
add    $0x30,%rdx
cmp    %rcx,%r12
jne    14b0 <hc_ludcmp_3x3+0x380>
test   $0x1,%r15b
je     1400 <hc_ludcmp_3x3+0x2d0>
movsd  (%rax,%rcx,8),%xmm1
lea    (%rcx,%rcx,2),%rcx
lea    (%rbx,%rcx,8),%rcx
mulsd  (%rcx,%r15,8),%xmm1
subsd  %xmm1,%xmm0
jmp    1400 <hc_ludcmp_3x3+0x2d0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %ebp,%eax
cmp    %rax,%r15
mov    %rbx,%r13
mov    0x28(%rsp),%rbx
mov    0x68(%rsp),%r12
mov    0x78(%rsp),%r8
je     1709 <hc_ludcmp_3x3+0x5d9>
movslq %ebp,%r9
cmpl   $0x4,0xc(%rsp)
jae    1550 <hc_ludcmp_3x3+0x420>
xor    %edx,%edx
mov    0x30(%rsp),%rax
jmp    1670 <hc_ludcmp_3x3+0x540>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    0x0(,%r9,8),%rcx
lea    (%rcx,%rcx,2),%rcx
lea    (%rcx,%r13,1),%rsi
lea    0x0(,%r15,8),%rdx
lea    (%rdx,%rdx,2),%rdx
lea    (%r12,%rdx,1),%rdi
cmp    %rdi,%rsi
mov    0x30(%rsp),%rax
jae    158c <hc_ludcmp_3x3+0x45c>
add    %r12,%rcx
add    %r13,%rdx
cmp    %rcx,%rdx
jae    158c <hc_ludcmp_3x3+0x45c>
xor    %edx,%edx
jmp    1670 <hc_ludcmp_3x3+0x540>
cmpq   $0x0,0x50(%rsp)
je     179e <hc_ludcmp_3x3+0x66e>
lea    (%r9,%r9,2),%rcx
mov    0x38(%rsp),%rsi
lea    (%rsi,%rcx,8),%rdx
lea    (%rsi,%r8,8),%rsi
mov    0x40(%rsp),%rdi
xor    %ecx,%ecx
movups -0x30(%rdx,%rcx,8),%xmm0
movups -0x20(%rdx,%rcx,8),%xmm1
movups -0x30(%rsi,%rcx,8),%xmm2
movups -0x20(%rsi,%rcx,8),%xmm3
movups %xmm2,-0x30(%rdx,%rcx,8)
movups %xmm3,-0x20(%rdx,%rcx,8)
movups %xmm0,-0x30(%rsi,%rcx,8)
movups %xmm1,-0x20(%rsi,%rcx,8)
movupd -0x10(%rdx,%rcx,8),%xmm0
movupd (%rdx,%rcx,8),%xmm1
movupd -0x10(%rsi,%rcx,8),%xmm2
movupd (%rsi,%rcx,8),%xmm3
movupd %xmm2,-0x10(%rdx,%rcx,8)
movupd %xmm3,(%rdx,%rcx,8)
movupd %xmm0,-0x10(%rsi,%rcx,8)
movupd %xmm1,(%rsi,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdi
jne    15b0 <hc_ludcmp_3x3+0x480>
testb  $0x1,0x48(%rsp)
je     1655 <hc_ludcmp_3x3+0x525>
lea    (%r9,%r9,2),%rdx
lea    0x0(,%rdx,8),%rdx
add    %r13,%rdx
movupd (%rdx,%rcx,8),%xmm0
movupd 0x10(%rdx,%rcx,8),%xmm1
lea    0x0(%r13,%r8,8),%rsi
movupd (%rsi,%rcx,8),%xmm2
movupd 0x10(%rsi,%rcx,8),%xmm3
movupd %xmm2,(%rdx,%rcx,8)
movupd %xmm3,0x10(%rdx,%rcx,8)
movupd %xmm0,(%rsi,%rcx,8)
movupd %xmm1,0x10(%rsi,%rcx,8)
mov    0x58(%rsp),%rcx
mov    %rcx,%rdx
cmp    %rbx,%rcx
je     16f5 <hc_ludcmp_3x3+0x5c5>
cs nopw 0x0(%rax,%rax,1)
mov    %rdx,%rcx
test   $0x1,%bl
je     16a7 <hc_ludcmp_3x3+0x577>
lea    (%r9,%r9,2),%rcx
lea    0x0(,%rcx,8),%rcx
add    %r13,%rcx
movsd  (%rcx,%rdx,8),%xmm0
lea    0x0(%r13,%r8,8),%rsi
movsd  (%rsi,%rdx,8),%xmm1
movsd  %xmm1,(%rcx,%rdx,8)
movsd  %xmm0,(%rsi,%rdx,8)
mov    %rdx,%rcx
or     $0x1,%rcx
not    %rdx
cmp    0x60(%rsp),%rdx
je     16f5 <hc_ludcmp_3x3+0x5c5>
lea    (%r9,%r9,2),%rdx
mov    0x20(%rsp),%rsi
lea    (%rsi,%rdx,8),%rdx
xchg   %ax,%ax
movsd  -0x8(%rdx,%rcx,8),%xmm0
movsd  -0x8(%rax,%rcx,8),%xmm1
movsd  %xmm1,-0x8(%rdx,%rcx,8)
movsd  %xmm0,-0x8(%rax,%rcx,8)
movsd  (%rdx,%rcx,8),%xmm0
movsd  (%rax,%rcx,8),%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
movsd  %xmm0,(%rax,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rbx
jne    16c0 <hc_ludcmp_3x3+0x590>
movsd  0xa0(%rsp,%r15,8),%xmm0
movsd  %xmm0,0xa0(%rsp,%r9,8)
mov    0x70(%rsp),%rax
mov    %ebp,(%rax,%r15,4)
lea    0x0(,%r8,8),%rax
add    %r13,%rax
lea    (%rax,%r15,8),%r14
movsd  (%rax,%r15,8),%xmm0
call   1030 <fabs@plt>
test   %eax,%eax
jg     173d <hc_ludcmp_3x3+0x60d>
movabs $0x3bc79ca10c924223,%rax
9c c7 3b
mov    %rax,(%r14)
mov    $0x3,%eax
cmp    $0x2,%r15
mov    0x88(%rsp),%rsi
mov    0x80(%rsp),%rdi
je     12c0 <hc_ludcmp_3x3+0x190>
lea    0x1(%r15),%rax
cmp    %rbx,%rax
jge    12c0 <hc_ludcmp_3x3+0x190>
movsd  0x88f(%rip),%xmm0
divsd  (%r14),%xmm0
lea    (%rsi,%r15,8),%rcx
mov    %rbx,%rdx
nopl   (%rax)
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
add    $0x18,%rcx
add    $0xffffffffffffffff,%rdx
cmp    %rdx,%rdi
jne    1780 <hc_ludcmp_3x3+0x650>
jmp    12c0 <hc_ludcmp_3x3+0x190>
xor    %ecx,%ecx
testb  $0x1,0x48(%rsp)
jne    1615 <hc_ludcmp_3x3+0x4e5>
jmp    1655 <hc_ludcmp_3x3+0x525>
add    $0xb8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x2817(%rip),%r12
mov    (%r12),%edi
lea    0x834(%rip),%rsi
movsd  0x10(%rsp),%xmm0
mov    $0x1,%al
call   1040 <fprintf@plt>
lea    0x854(%rip),%rbp
xor    %r15d,%r15d
mov    0x28(%rsp),%r14
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%r12),%edi
movsd  0x0(%r13,%rbx,8),%xmm0
mov    %rbp,%rsi
mov    $0x1,%al
call   1040 <fprintf@plt>
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1800 <hc_ludcmp_3x3+0x6d0>
mov    (%r12),%edi
lea    0x811(%rip),%rsi
xor    %eax,%eax
call   1040 <fprintf@plt>
add    $0x1,%r15
add    $0x18,%r13
cmp    %r14,%r15
jne    17f0 <hc_ludcmp_3x3+0x6c0>
mov    $0xffffffff,%edi
call   1050 <exit@plt>

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

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