#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*/ dswap(double *, double *);
double *dvector(int, long);
double fabs(double);
int /*<<< orphan*/ free_dvector(double *, int, long);
int /*<<< orphan*/ nrerror(char *);
void dludcmp(double **a, long n, long *indx, double *d) {
  double big, dum, sum, temp;
  double *vv;
  long i, j, k;
  long imax = 0;
  vv = dvector(1, n);
  *d = 1.0;
  for (i = 1; i <= n; i++) {
    big = 0.0;
    for (j = 1; j <= n; j++)
      if ((temp = fabs(a[i][j])) > big)
        big = temp;
    if (big == 0.0)
      nrerror("singular matrix in routine dludcmp");
    vv[i] = 1.0 / big;
  }
  for (j = 1; j <= n; j++) {
    for (i = 1; i < j; i++) {
      sum = a[i][j];
      for (k = 1; 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 = 1; 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 = 1; k <= n; k++)
        dswap(&a[imax][k], &a[j][k]);
      *d = -(*d);
      vv[imax] = vv[j];
    }
    indx[j] = imax;
    if (a[j][j] == 0.0)
      a[j][j] = 1.0e-7;
    if (j != n) {
      dum = 1.0 / a[j][j];
      for (i = j + 1; i <= n; i++)
        a[i][j] *= dum;
    }
  }
  free_dvector(vv, 1, n);
}


x86 assembly (clang, -O0):
<dludcmp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movq   $0x0,-0x68(%rbp)
mov    -0x10(%rbp),%rsi
mov    $0x1,%edi
call   1050 <dvector@plt>
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movsd  0xe82(%rip),%xmm0
movsd  %xmm0,(%rax)
movq   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     1255 <dludcmp+0x115>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movq   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     1206 <dludcmp+0xc6>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe3e(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
ucomisd -0x28(%rbp),%xmm0
jbe    11f0 <dludcmp+0xb0>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11f5 <dludcmp+0xb5>
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    11a8 <dludcmp+0x68>
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    122a <dludcmp+0xea>
jp     122a <dludcmp+0xea>
lea    0xdfb(%rip),%rdi
call   1030 <nrerror@plt>
movsd  0xdce(%rip),%xmm0
divsd  -0x28(%rbp),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    118a <dludcmp+0x4a>
movq   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     15f8 <dludcmp+0x4b8>
movq   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x58(%rbp),%rax
jge    1340 <dludcmp+0x200>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movq   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%rax
cmp    -0x50(%rbp),%rax
jge    1315 <dludcmp+0x1d5>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    12a3 <dludcmp+0x163>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    1273 <dludcmp+0x133>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x58(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     1465 <dludcmp+0x325>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movq   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%rax
cmp    -0x58(%rbp),%rax
jge    13f2 <dludcmp+0x2b2>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    1380 <dludcmp+0x240>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x38(%rbp),%xmm1
movaps 0xbeb(%rip),%xmm2
pand   %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
ucomisd -0x28(%rbp),%xmm0
jb     144f <dludcmp+0x30f>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x50(%rbp),%rax
mov    %rax,-0x68(%rbp)
jmp    1454 <dludcmp+0x314>
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    1350 <dludcmp+0x210>
mov    -0x58(%rbp),%rax
cmp    -0x68(%rbp),%rax
je     150e <dludcmp+0x3ce>
movq   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     14cd <dludcmp+0x38d>
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x60(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
mov    -0x60(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
call   1040 <dswap@plt>
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    147b <dludcmp+0x33b>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x48(%rbp),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x68(%rbp),%rdx
mov    -0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1563 <dludcmp+0x423>
jp     1563 <dludcmp+0x423>
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  0xaaa(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
cmp    -0x10(%rbp),%rax
je     15e2 <dludcmp+0x4a2>
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
movsd  0xa77(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x10(%rbp),%rax
jg     15dd <dludcmp+0x49d>
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x58(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    159f <dludcmp+0x45f>
jmp    15e2 <dludcmp+0x4a2>
jmp    15e7 <dludcmp+0x4a7>
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    125d <dludcmp+0x11d>
mov    -0x48(%rbp),%rdi
mov    -0x10(%rbp),%rdx
mov    $0x1,%esi
call   1060 <free_dvector@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dludcmp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%rbp
mov    %rdx,0x18(%rsp)
mov    %rsi,%r13
mov    %rdi,%rbx
mov    $0x1,%edi
call   1050 <dvector@plt>
mov    %rax,%r9
movabs $0x3ff0000000000000,%rax
mov    %rbp,0x10(%rsp)
mov    %rax,0x0(%rbp)
test   %r13,%r13
jle    1573 <dludcmp+0x433>
lea    -0x1(%r13),%r15
mov    %r13d,%ebp
and    $0x3,%ebp
mov    %r13,%r12
and    $0xfffffffffffffffc,%r12
mov    $0x1,%r14d
movapd 0xe5c(%rip),%xmm4
movapd 0xe54(%rip),%xmm5
mov    %r9,(%rsp)
jmp    11e2 <dludcmp+0xa2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0xe58(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r9,%r14,8)
lea    0x1(%r14),%rax
cmp    %r13,%r14
mov    %rax,%r14
je     12bf <dludcmp+0x17f>
mov    (%rbx,%r14,8),%rax
mov    $0x1,%ecx
xorpd  %xmm0,%xmm0
cmp    $0x3,%r15
jb     1245 <dludcmp+0x105>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  0x8(%rax,%rcx,8),%xmm1
movsd  0x10(%rax,%rcx,8),%xmm2
movsd  0x18(%rax,%rcx,8),%xmm3
andpd  %xmm5,%xmm1
maxsd  %xmm0,%xmm1
movsd  0x20(%rax,%rcx,8),%xmm0
andpd  %xmm5,%xmm2
maxsd  %xmm1,%xmm2
andpd  %xmm5,%xmm3
maxsd  %xmm2,%xmm3
andpd  %xmm5,%xmm0
maxsd  %xmm3,%xmm0
add    $0x4,%rcx
cmp    %rcx,%r12
jne    1200 <dludcmp+0xc0>
add    $0x1,%rcx
movapd %xmm0,%xmm1
test   %rbp,%rbp
je     127a <dludcmp+0x13a>
lea    (%rax,%rcx,8),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rax,%rcx,8),%xmm1
andpd  %xmm4,%xmm1
maxsd  %xmm0,%xmm1
add    $0x1,%rcx
movapd %xmm1,%xmm0
cmp    %rcx,%rbp
jne    1260 <dludcmp+0x120>
ucomisd 0xdae(%rip),%xmm1
jne    11c0 <dludcmp+0x80>
jp     11c0 <dludcmp+0x80>
lea    0xda3(%rip),%rdi
movapd %xmm1,0x20(%rsp)
call   1030 <nrerror@plt>
movapd 0x20(%rsp),%xmm1
movapd 0xd52(%rip),%xmm5
movapd 0xd4a(%rip),%xmm4
mov    (%rsp),%r9
jmp    11c0 <dludcmp+0x80>
test   %r13,%r13
jle    1573 <dludcmp+0x433>
mov    $0x1,%r15d
xor    %r12d,%r12d
xorpd  %xmm5,%xmm5
xor    %r14d,%r14d
mov    %r13,0x8(%rsp)
jmp    12f4 <dludcmp+0x1b4>
nop
lea    0x1(%r15),%rax
add    $0x1,%r12
cmp    %r13,%r15
mov    %rax,%r15
je     1573 <dludcmp+0x433>
mov    %r12,%rax
and    $0xfffffffffffffffe,%rax
cmp    $0x2,%r15
jb     13c0 <dludcmp+0x280>
mov    $0x1,%ecx
xor    %r8d,%r8d
jmp    1327 <dludcmp+0x1e7>
nop
movsd  %xmm0,(%rsi,%r15,8)
add    $0x1,%rcx
add    $0x1,%r8
cmp    %r15,%rcx
je     13c0 <dludcmp+0x280>
mov    (%rbx,%rcx,8),%rsi
movsd  (%rsi,%r15,8),%xmm0
cmp    $0x2,%rcx
jb     1310 <dludcmp+0x1d0>
mov    $0x1,%edi
cmp    $0x1,%r8
je     139c <dludcmp+0x25c>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
neg    %rdx
mov    $0x1,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rdi,8),%xmm1
movsd  0x8(%rsi,%rdi,8),%xmm2
mov    (%rbx,%rdi,8),%rbp
mulsd  0x0(%rbp,%r15,8),%xmm1
mov    0x8(%rbx,%rdi,8),%rbp
subsd  %xmm1,%xmm0
mulsd  0x0(%rbp,%r15,8),%xmm2
subsd  %xmm2,%xmm0
lea    (%rdx,%rdi,1),%rbp
add    $0x2,%rbp
add    $0x2,%rdi
cmp    $0x1,%rbp
jne    1360 <dludcmp+0x220>
test   $0x1,%r8b
je     1310 <dludcmp+0x1d0>
movsd  (%rsi,%rdi,8),%xmm1
mov    (%rbx,%rdi,8),%rdx
mulsd  (%rdx,%r15,8),%xmm1
subsd  %xmm1,%xmm0
jmp    1310 <dludcmp+0x1d0>
xchg   %ax,%ax
neg    %rax
xorpd  %xmm0,%xmm0
mov    %r15,%rcx
jmp    1415 <dludcmp+0x2d5>
nopl   0x0(%rax)
movsd  %xmm1,(%rdx,%r15,8)
andpd  %xmm4,%xmm1
mulsd  (%r9,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
movapd %xmm0,%xmm2
cmpnlesd %xmm1,%xmm2
movapd %xmm2,%xmm3
andnpd %xmm1,%xmm3
andpd  %xmm0,%xmm2
orpd   %xmm3,%xmm2
cmovae %rcx,%r14
lea    0x1(%rcx),%rdx
movapd %xmm2,%xmm0
cmp    %r13,%rcx
mov    %rdx,%rcx
jge    14a0 <dludcmp+0x360>
mov    (%rbx,%rcx,8),%rdx
movsd  (%rdx,%r15,8),%xmm1
cmp    $0x2,%r15
jb     13d0 <dludcmp+0x290>
mov    $0x1,%esi
cmp    $0x1,%r12
je     147a <dludcmp+0x33a>
mov    $0x1,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx,%rsi,8),%xmm2
movsd  0x8(%rdx,%rsi,8),%xmm3
mov    (%rbx,%rsi,8),%rdi
mulsd  (%rdi,%r15,8),%xmm2
mov    0x8(%rbx,%rsi,8),%rdi
subsd  %xmm2,%xmm1
mulsd  (%rdi,%r15,8),%xmm3
subsd  %xmm3,%xmm1
lea    (%rax,%rsi,1),%rdi
add    $0x2,%rdi
add    $0x2,%rsi
cmp    $0x1,%rdi
jne    1440 <dludcmp+0x300>
test   $0x1,%r12b
je     13d0 <dludcmp+0x290>
movsd  (%rdx,%rsi,8),%xmm2
mov    (%rbx,%rsi,8),%rsi
mulsd  (%rsi,%r15,8),%xmm2
subsd  %xmm2,%xmm1
jmp    13d0 <dludcmp+0x290>
nopl   0x0(%rax)
cmp    %r14,%r15
je     1501 <dludcmp+0x3c1>
mov    $0x8,%ebp
nopw   0x0(%rax,%rax,1)
mov    (%rbx,%r14,8),%rdi
add    %rbp,%rdi
mov    (%rbx,%r15,8),%rsi
add    %rbp,%rsi
call   1040 <dswap@plt>
add    $0x8,%rbp
add    $0xffffffffffffffff,%r13
jne    14b0 <dludcmp+0x370>
mov    0x10(%rsp),%rax
movsd  (%rax),%xmm0
xorps  0xb33(%rip),%xmm0
movlps %xmm0,(%rax)
mov    (%rsp),%r9
movsd  (%r9,%r15,8),%xmm0
movsd  %xmm0,(%r9,%r14,8)
mov    0x8(%rsp),%r13
movapd 0xb03(%rip),%xmm4
xorpd  %xmm5,%xmm5
mov    0x18(%rsp),%rax
mov    %r14,(%rax,%r15,8)
mov    (%rbx,%r15,8),%rax
movsd  (%rax,%r15,8),%xmm0
ucomisd %xmm5,%xmm0
jne    1532 <dludcmp+0x3f2>
jp     1532 <dludcmp+0x3f2>
movabs $0x3e7ad7f29abcaf48,%rcx
d7 7a 3e
mov    %rcx,(%rax,%r15,8)
movsd  0xaf6(%rip),%xmm0
cmp    %r13,%r15
je     12e0 <dludcmp+0x1a0>
jge    12e0 <dludcmp+0x1a0>
movsd  0xad7(%rip),%xmm1
divsd  %xmm0,%xmm1
mov    %r15,%rax
mov    0x8(%rbx,%rax,8),%rcx
add    $0x1,%rax
movsd  (%rcx,%r15,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%r15,8)
cmp    %r13,%rax
jl     1550 <dludcmp+0x410>
jmp    12e0 <dludcmp+0x1a0>
mov    %r9,%rdi
mov    $0x1,%esi
mov    %r13,%rdx
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1060 <free_dvector@plt>

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

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