void inverse(double **a, int dim) {
  int i, j, *indx;
  double **y, d, *col;
  y = dmatrix(dim, dim);
  indx = ivector(dim);
  col = dvector(dim);
  ludcmp(a, dim, indx, &d);
  for (j = 0; j < dim; j++) {
    memset(col, 0, sizeof(double) * dim);
    col[j] = 1.0;
    lubksb(a, dim, indx, col);
    for (i = 0; i < dim; i++)
      y[i][j] = col[i];
  }
  for (i = 0; i < dim; i++)
    for (j = 0; j < dim; j++)
      a[i][j] = y[i][j];
  free_dmatrix(y, dim);
  free(col);
  free(indx);
}


x86 assembly (clang, -O0):
<inverse>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%edi
mov    -0xc(%rbp),%esi
mov    $0x0,%al
call   1060 <dmatrix@plt>
cltq
mov    %rax,-0x28(%rbp)
mov    -0xc(%rbp),%edi
mov    $0x0,%al
call   1090 <ivector@plt>
cltq
mov    %rax,-0x20(%rbp)
mov    -0xc(%rbp),%edi
mov    $0x0,%al
call   1080 <dvector@plt>
cltq
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x20(%rbp),%rdx
lea    -0x30(%rbp),%rcx
mov    $0x0,%al
call   1040 <ludcmp@plt>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    127f <inverse+0xff>
mov    -0x38(%rbp),%rdi
movslq -0xc(%rbp),%rdx
shl    $0x3,%rdx
xor    %esi,%esi
call   1050 <memset@plt>
mov    -0x38(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  0xdf2(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x20(%rbp),%rdx
mov    -0x38(%rbp),%rcx
mov    $0x0,%al
call   1070 <lubksb@plt>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    126c <inverse+0xec>
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1230 <inverse+0xb0>
jmp    1271 <inverse+0xf1>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11df <inverse+0x5f>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12f0 <inverse+0x170>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12dd <inverse+0x15d>
mov    -0x28(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1299 <inverse+0x119>
jmp    12e2 <inverse+0x162>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1286 <inverse+0x106>
mov    -0x28(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    $0x0,%al
call   10a0 <free_dmatrix@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<inverse>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %esi,%r14d
mov    %rdi,%r15
mov    %esi,%edi
xor    %eax,%eax
call   1060 <dmatrix@plt>
movslq %eax,%rbx
mov    %r14d,%edi
xor    %eax,%eax
call   1090 <ivector@plt>
movslq %eax,%r12
mov    %r14d,%edi
xor    %eax,%eax
call   1080 <dvector@plt>
movslq %eax,%rbp
lea    0x30(%rsp),%rcx
mov    %r15,0x18(%rsp)
mov    %r15,%rdi
mov    %r14d,%esi
mov    %r12,0x10(%rsp)
mov    %r12,%rdx
xor    %eax,%eax
call   1040 <ludcmp@plt>
mov    %r14d,0xc(%rsp)
test   %r14d,%r14d
jle    14bc <inverse+0x33c>
mov    0xc(%rsp),%r14d
lea    0x0(,%r14,8),%rax
mov    %rax,0x28(%rsp)
lea    -0x1(%r14),%rax
mov    %rax,0x20(%rsp)
mov    %r14d,%r15d
and    $0x3,%r15d
mov    %r14d,%r12d
and    $0xfffffffc,%r12d
xor    %r13d,%r13d
jmp    122d <inverse+0xad>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r13
cmp    %r14,%r13
je     130d <inverse+0x18d>
mov    %rbp,%rdi
xor    %esi,%esi
mov    0x28(%rsp),%rdx
call   1050 <memset@plt>
movabs $0x3ff0000000000000,%rax
mov    %rax,0x0(%rbp,%r13,8)
mov    0x18(%rsp),%rdi
mov    0xc(%rsp),%esi
mov    0x10(%rsp),%rdx
mov    %rbp,%rcx
xor    %eax,%eax
call   1070 <lubksb@plt>
cmpq   $0x3,0x20(%rsp)
jae    1270 <inverse+0xf0>
xor    %eax,%eax
jmp    12cc <inverse+0x14c>
nop
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%rax,8),%xmm0
mov    (%rbx,%rax,8),%rcx
movsd  %xmm0,(%rcx,%r13,8)
movsd  0x8(%rbp,%rax,8),%xmm0
mov    0x8(%rbx,%rax,8),%rcx
movsd  %xmm0,(%rcx,%r13,8)
movsd  0x10(%rbp,%rax,8),%xmm0
mov    0x10(%rbx,%rax,8),%rcx
movsd  %xmm0,(%rcx,%r13,8)
movsd  0x18(%rbp,%rax,8),%xmm0
mov    0x18(%rbx,%rax,8),%rcx
movsd  %xmm0,(%rcx,%r13,8)
add    $0x4,%rax
cmp    %rax,%r12
jne    1280 <inverse+0x100>
test   %r15,%r15
je     1220 <inverse+0xa0>
lea    (%rbx,%rax,8),%rcx
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm0
mov    (%rcx,%rdx,8),%rsi
movsd  %xmm0,(%rsi,%r13,8)
add    $0x1,%rdx
cmp    %rdx,%r15
jne    12f0 <inverse+0x170>
jmp    1220 <inverse+0xa0>
cmpl   $0x0,0xc(%rsp)
jle    14bc <inverse+0x33c>
mov    %r14d,%r10d
and    $0xfffffffc,%r10d
add    $0xfffffffffffffffc,%r10
mov    %r10,%r8
shr    $0x2,%r8
add    $0x1,%r8
mov    %r14d,%r11d
and    $0xfffffffc,%r11d
mov    %r8d,%r9d
and    $0x3,%r9d
and    $0xfffffffffffffffc,%r8
mov    %r14d,%r15d
and    $0x3,%r15d
mov    %r9,%r12
shl    $0x5,%r12
xor    %r13d,%r13d
jmp    136d <inverse+0x1ed>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r13
cmp    %r14,%r13
je     14bc <inverse+0x33c>
mov    (%rbx,%r13,8),%rax
mov    0x18(%rsp),%rcx
mov    (%rcx,%r13,8),%rdx
cmpl   $0x4,0xc(%rsp)
jb     139b <inverse+0x21b>
lea    (%rax,%r14,8),%rcx
cmp    %rcx,%rdx
jae    1410 <inverse+0x290>
lea    (%rdx,%r14,8),%rcx
cmp    %rcx,%rax
jae    1410 <inverse+0x290>
xor    %edi,%edi
mov    %rdi,%rsi
not    %rsi
add    %r14,%rsi
test   %r15,%r15
je     13c4 <inverse+0x244>
mov    %r15,%rcx
xchg   %ax,%ax
movsd  (%rax,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,8)
add    $0x1,%rdi
add    $0xffffffffffffffff,%rcx
jne    13b0 <inverse+0x230>
cmp    $0x3,%rsi
jb     1360 <inverse+0x1e0>
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,8)
movsd  0x8(%rax,%rdi,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rdi,8)
movsd  0x10(%rax,%rdi,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rdi,8)
movsd  0x18(%rax,%rdi,8),%xmm0
movsd  %xmm0,0x18(%rdx,%rdi,8)
add    $0x4,%rdi
cmp    %rdi,%r14
jne    13d0 <inverse+0x250>
jmp    1360 <inverse+0x1e0>
nopl   0x0(%rax)
cmp    $0xc,%r10
jae    141a <inverse+0x29a>
xor    %edi,%edi
jmp    1478 <inverse+0x2f8>
mov    %r8,%rsi
xor    %edi,%edi
nop
movups (%rax,%rdi,8),%xmm0
movups 0x10(%rax,%rdi,8),%xmm1
movups %xmm0,(%rdx,%rdi,8)
movups %xmm1,0x10(%rdx,%rdi,8)
movups 0x20(%rax,%rdi,8),%xmm0
movups 0x30(%rax,%rdi,8),%xmm1
movups %xmm0,0x20(%rdx,%rdi,8)
movups %xmm1,0x30(%rdx,%rdi,8)
movups 0x40(%rax,%rdi,8),%xmm0
movups 0x50(%rax,%rdi,8),%xmm1
movups %xmm0,0x40(%rdx,%rdi,8)
movups %xmm1,0x50(%rdx,%rdi,8)
movups 0x60(%rax,%rdi,8),%xmm0
movups 0x70(%rax,%rdi,8),%xmm1
movups %xmm0,0x60(%rdx,%rdi,8)
movups %xmm1,0x70(%rdx,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1420 <inverse+0x2a0>
test   %r9,%r9
je     14ab <inverse+0x32b>
lea    (%rdx,%rdi,8),%rsi
add    $0x10,%rsi
lea    (%rax,%rdi,8),%rcx
add    $0x10,%rcx
xor    %edi,%edi
nop
movups -0x10(%rcx,%rdi,1),%xmm0
movups (%rcx,%rdi,1),%xmm1
movups %xmm0,-0x10(%rsi,%rdi,1)
movups %xmm1,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%r12
jne    1490 <inverse+0x310>
mov    %r11,%rdi
cmp    %r14,%r11
je     1360 <inverse+0x1e0>
jmp    139d <inverse+0x21d>
mov    %rbx,%rdi
mov    0xc(%rsp),%esi
xor    %eax,%eax
call   10a0 <free_dmatrix@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    0x10(%rsp),%rdi
call   1030 <free@plt>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @inverse(double** %0, i32 %1) {
%3 = alloca double**
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32*
%8 = alloca double**
%9 = alloca double
%10 = alloca double*
store double** %0, double*** %3
store i32 %1, i32* %4
%11 = load i32, i32* %4
%12 = load i32, i32* %4
%13 = call i32 (i32, i32, ...) bitcast (i32 (...)* @dmatrix to i32 (i32, i32, ...)*)(i32 %11, i32 %12)
%14 = sext i32 %13 to i64
%15 = inttoptr i64 %14 to double**
store double** %15, double*** %8
%16 = load i32, i32* %4
%17 = call i32 (i32, ...) bitcast (i32 (...)* @ivector to i32 (i32, ...)*)(i32 %16)
%18 = sext i32 %17 to i64
%19 = inttoptr i64 %18 to i32*
store i32* %19, i32** %7
%20 = load i32, i32* %4
%21 = call i32 (i32, ...) bitcast (i32 (...)* @dvector to i32 (i32, ...)*)(i32 %20)
%22 = sext i32 %21 to i64
%23 = inttoptr i64 %22 to double*
store double* %23, double** %10
%24 = load double**, double*** %3
%25 = load i32, i32* %4
%26 = load i32*, i32** %7
%27 = call i32 (double**, i32, i32*, double*, ...) bitcast (i32 (...)* @ludcmp to i32 (double**, i32, i32*, double*, ...)*)(double** %24, i32 %25, i32* %26, double* %9)
store i32 0, i32* %6
br label %28
28: ; preds = %69, %2
%29 = load i32, i32* %6
%30 = load i32, i32* %4
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %72
32: ; preds = %28
%33 = load double*, double** %10
%34 = bitcast double* %33 to i8*
%35 = load i32, i32* %4
%36 = sext i32 %35 to i64
%37 = mul i64 8, %36
%38 = load double*, double** %10
%39 = load i32, i32* %6
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %38, i64 %40
store double 1.000000e+00, double* %41
%42 = load double**, double*** %3
%43 = load i32, i32* %4
%44 = load i32*, i32** %7
%45 = load double*, double** %10
%46 = call i32 (double**, i32, i32*, double*, ...) bitcast (i32 (...)* @lubksb to i32 (double**, i32, i32*, double*, ...)*)(double** %42, i32 %43, i32* %44, double* %45)
store i32 0, i32* %5
br label %47
47: ; preds = %65, %32
%48 = load i32, i32* %5
%49 = load i32, i32* %4
%50 = icmp slt i32 %48, %49
br i1 %50, label %51, label %68
51: ; preds = %47
%52 = load double*, double** %10
%53 = load i32, i32* %5
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %52, i64 %54
%56 = load double, double* %55
%57 = load double**, double*** %8
%58 = load i32, i32* %5
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double*, double** %57, i64 %59
%61 = load double*, double** %60
%62 = load i32, i32* %6
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
store double %56, double* %64
br label %65
65: ; preds = %51
%66 = load i32, i32* %5
%67 = add nsw i32 %66, 1
store i32 %67, i32* %5
br label %47
68: ; preds = %47
br label %69
69: ; preds = %68
%70 = load i32, i32* %6
%71 = add nsw i32 %70, 1
store i32 %71, i32* %6
br label %28
72: ; preds = %28
store i32 0, i32* %5
br label %73
73: ; preds = %104, %72
%74 = load i32, i32* %5
%75 = load i32, i32* %4
%76 = icmp slt i32 %74, %75
br i1 %76, label %77, label %107
77: ; preds = %73
store i32 0, i32* %6
br label %78
78: ; preds = %100, %77
%79 = load i32, i32* %6
%80 = load i32, i32* %4
%81 = icmp slt i32 %79, %80
br i1 %81, label %82, label %103
82: ; preds = %78
%83 = load double**, double*** %8
%84 = load i32, i32* %5
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double*, double** %83, i64 %85
%87 = load double*, double** %86
%88 = load i32, i32* %6
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = load double**, double*** %3
%93 = load i32, i32* %5
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double*, double** %92, i64 %94
%96 = load double*, double** %95
%97 = load i32, i32* %6
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %96, i64 %98
store double %91, double* %99
br label %100
100: ; preds = %82
%101 = load i32, i32* %6
%102 = add nsw i32 %101, 1
store i32 %102, i32* %6
br label %78
103: ; preds = %78
br label %104
104: ; preds = %103
%105 = load i32, i32* %5
%106 = add nsw i32 %105, 1
store i32 %106, i32* %5
br label %73
107: ; preds = %73
%108 = load double**, double*** %8
%109 = load i32, i32* %4
%110 = call i32 (double**, i32, ...) bitcast (i32 (...)* @free_dmatrix to i32 (double**, i32, ...)*)(double** %108, i32 %109)
%111 = load double*, double** %10
%112 = bitcast double* %111 to i8*
%113 = load i32*, i32** %7
%114 = bitcast i32* %113 to i8*
ret void
}

LLVM IR (-O2):
define void @inverse(double** %0, i32 %1) {
%3 = alloca double
%4 = bitcast double* %3 to i8*
%5 = tail call i32 (i32, i32, ...) bitcast (i32 (...)* @dmatrix to i32 (i32, i32, ...)*)(i32 %1, i32 %1)
%6 = sext i32 %5 to i64
%7 = inttoptr i64 %6 to double**
%8 = tail call i32 (i32, ...) bitcast (i32 (...)* @ivector to i32 (i32, ...)*)(i32 %1)
%9 = sext i32 %8 to i64
%10 = inttoptr i64 %9 to i32*
%11 = tail call i32 (i32, ...) bitcast (i32 (...)* @dvector to i32 (i32, ...)*)(i32 %1)
%12 = sext i32 %11 to i64
%13 = inttoptr i64 %12 to double*
%14 = call i32 (double**, i32, i32*, double*, ...) bitcast (i32 (...)* @ludcmp to i32 (double**, i32, i32*, double*, ...)*)(double** %0, i32 %1, i32* %10, double* %3)
%15 = icmp sgt i32 %1, 0
br i1 %15, label %16, label %210
16: ; preds = %2
%17 = inttoptr i64 %12 to i8*
%18 = zext i32 %1 to i64
%19 = shl nuw nsw i64 %18, 3
%20 = zext i32 %1 to i64
%21 = add nsw i64 %20, -1
%22 = and i64 %20, 3
%23 = icmp ult i64 %21, 3
%24 = and i64 %20, 4294967292
%25 = icmp eq i64 %22, 0
br label %42
26: ; preds = %88
br i1 %15, label %27, label %210
27: ; preds = %26
%28 = zext i32 %1 to i64
%29 = and i64 %28, 4294967292
%30 = add nsw i64 %29, -4
%31 = lshr exact i64 %30, 2
%32 = add nuw nsw i64 %31, 1
%33 = icmp ult i32 %1, 4
%34 = and i64 %20, 4294967292
%35 = and i64 %32, 3
%36 = icmp ult i64 %30, 12
%37 = and i64 %32, 9223372036854775804
%38 = icmp eq i64 %35, 0
%39 = icmp eq i64 %34, %20
%40 = and i64 %28, 3
%41 = icmp eq i64 %40, 0
br label %91
42: ; preds = %88, %16
%43 = phi i64 [ 0, %16 ], [ %89, %88 ]
%44 = getelementptr inbounds double, double* %13, i64 %43
store double 1.000000e+00, double* %44
%45 = call i32 (double**, i32, i32*, double*, ...) bitcast (i32 (...)* @lubksb to i32 (double**, i32, i32*, double*, ...)*)(double** %0, i32 %1, i32* %10, double* %13)
br i1 %23, label %75, label %46
46: ; preds = %42, %46
%47 = phi i64 [ %72, %46 ], [ 0, %42 ]
%48 = phi i64 [ %73, %46 ], [ 0, %42 ]
%49 = getelementptr inbounds double, double* %13, i64 %47
%50 = load double, double* %49
%51 = getelementptr inbounds double*, double** %7, i64 %47
%52 = load double*, double** %51
%53 = getelementptr inbounds double, double* %52, i64 %43
store double %50, double* %53
%54 = or i64 %47, 1
%55 = getelementptr inbounds double, double* %13, i64 %54
%56 = load double, double* %55
%57 = getelementptr inbounds double*, double** %7, i64 %54
%58 = load double*, double** %57
%59 = getelementptr inbounds double, double* %58, i64 %43
store double %56, double* %59
%60 = or i64 %47, 2
%61 = getelementptr inbounds double, double* %13, i64 %60
%62 = load double, double* %61
%63 = getelementptr inbounds double*, double** %7, i64 %60
%64 = load double*, double** %63
%65 = getelementptr inbounds double, double* %64, i64 %43
store double %62, double* %65
%66 = or i64 %47, 3
%67 = getelementptr inbounds double, double* %13, i64 %66
%68 = load double, double* %67
%69 = getelementptr inbounds double*, double** %7, i64 %66
%70 = load double*, double** %69
%71 = getelementptr inbounds double, double* %70, i64 %43
store double %68, double* %71
%72 = add nuw nsw i64 %47, 4
%73 = add i64 %48, 4
%74 = icmp eq i64 %73, %24
br i1 %74, label %75, label %46
75: ; preds = %46, %42
%76 = phi i64 [ 0, %42 ], [ %72, %46 ]
br i1 %25, label %88, label %77
77: ; preds = %75, %77
%78 = phi i64 [ %85, %77 ], [ %76, %75 ]
%79 = phi i64 [ %86, %77 ], [ 0, %75 ]
%80 = getelementptr inbounds double, double* %13, i64 %78
%81 = load double, double* %80
%82 = getelementptr inbounds double*, double** %7, i64 %78
%83 = load double*, double** %82
%84 = getelementptr inbounds double, double* %83, i64 %43
store double %81, double* %84
%85 = add nuw nsw i64 %78, 1
%86 = add i64 %79, 1
%87 = icmp eq i64 %86, %22
br i1 %87, label %88, label %77
88: ; preds = %77, %75
%89 = add nuw nsw i64 %43, 1
%90 = icmp eq i64 %89, %20
br i1 %90, label %26, label %42
91: ; preds = %207, %27
%92 = phi i64 [ 0, %27 ], [ %208, %207 ]
%93 = getelementptr inbounds double*, double** %7, i64 %92
%94 = load double*, double** %93
%95 = getelementptr inbounds double*, double** %0, i64 %92
%96 = load double*, double** %95
br i1 %33, label %172, label %97
97: ; preds = %91
%98 = getelementptr double, double* %96, i64 %20
%99 = getelementptr double, double* %94, i64 %20
%100 = icmp ult double* %96, %99
%101 = icmp ult double* %94, %98
%102 = and i1 %100, %101
br i1 %102, label %172, label %103
103: ; preds = %97
br i1 %36, label %153, label %104
104: ; preds = %103, %104
%105 = phi i64 [ %150, %104 ], [ 0, %103 ]
%106 = phi i64 [ %151, %104 ], [ 0, %103 ]
%107 = getelementptr inbounds double, double* %94, i64 %105
%108 = bitcast double* %107 to <2 x double>*
%109 = load <2 x double>, <2 x double>* %108
%110 = getelementptr inbounds double, double* %107, i64 2
%111 = bitcast double* %110 to <2 x double>*
%112 = load <2 x double>, <2 x double>* %111
%113 = getelementptr inbounds double, double* %96, i64 %105
%114 = bitcast double* %113 to <2 x double>*
store <2 x double> %109, <2 x double>* %114
%115 = getelementptr inbounds double, double* %113, i64 2
%116 = bitcast double* %115 to <2 x double>*
store <2 x double> %112, <2 x double>* %116
%117 = or i64 %105, 4
%118 = getelementptr inbounds double, double* %94, i64 %117
%119 = bitcast double* %118 to <2 x double>*
%120 = load <2 x double>, <2 x double>* %119
%121 = getelementptr inbounds double, double* %118, i64 2
%122 = bitcast double* %121 to <2 x double>*
%123 = load <2 x double>, <2 x double>* %122
%124 = getelementptr inbounds double, double* %96, i64 %117
%125 = bitcast double* %124 to <2 x double>*
store <2 x double> %120, <2 x double>* %125
%126 = getelementptr inbounds double, double* %124, i64 2
%127 = bitcast double* %126 to <2 x double>*
store <2 x double> %123, <2 x double>* %127
%128 = or i64 %105, 8
%129 = getelementptr inbounds double, double* %94, i64 %128
%130 = bitcast double* %129 to <2 x double>*
%131 = load <2 x double>, <2 x double>* %130
%132 = getelementptr inbounds double, double* %129, i64 2
%133 = bitcast double* %132 to <2 x double>*
%134 = load <2 x double>, <2 x double>* %133
%135 = getelementptr inbounds double, double* %96, i64 %128
%136 = bitcast double* %135 to <2 x double>*
store <2 x double> %131, <2 x double>* %136
%137 = getelementptr inbounds double, double* %135, i64 2
%138 = bitcast double* %137 to <2 x double>*
store <2 x double> %134, <2 x double>* %138
%139 = or i64 %105, 12
%140 = getelementptr inbounds double, double* %94, i64 %139
%141 = bitcast double* %140 to <2 x double>*
%142 = load <2 x double>, <2 x double>* %141
%143 = getelementptr inbounds double, double* %140, i64 2
%144 = bitcast double* %143 to <2 x double>*
%145 = load <2 x double>, <2 x double>* %144
%146 = getelementptr inbounds double, double* %96, i64 %139
%147 = bitcast double* %146 to <2 x double>*
store <2 x double> %142, <2 x double>* %147
%148 = getelementptr inbounds double, double* %146, i64 2
%149 = bitcast double* %148 to <2 x double>*
store <2 x double> %145, <2 x double>* %149
%150 = add nuw i64 %105, 16
%151 = add nuw i64 %106, 4
%152 = icmp eq i64 %151, %37
br i1 %152, label %153, label %104
153: ; preds = %104, %103
%154 = phi i64 [ 0, %103 ], [ %150, %104 ]
br i1 %38, label %171, label %155
155: ; preds = %153, %155
%156 = phi i64 [ %168, %155 ], [ %154, %153 ]
%157 = phi i64 [ %169, %155 ], [ 0, %153 ]
%158 = getelementptr inbounds double, double* %94, i64 %156
%159 = bitcast double* %158 to <2 x double>*
%160 = load <2 x double>, <2 x double>* %159
%161 = getelementptr inbounds double, double* %158, i64 2
%162 = bitcast double* %161 to <2 x double>*
%163 = load <2 x double>, <2 x double>* %162
%164 = getelementptr inbounds double, double* %96, i64 %156
%165 = bitcast double* %164 to <2 x double>*
store <2 x double> %160, <2 x double>* %165
%166 = getelementptr inbounds double, double* %164, i64 2
%167 = bitcast double* %166 to <2 x double>*
store <2 x double> %163, <2 x double>* %167
%168 = add nuw i64 %156, 4
%169 = add i64 %157, 1
%170 = icmp eq i64 %169, %35
br i1 %170, label %171, label %155
171: ; preds = %155, %153
br i1 %39, label %207, label %172
172: ; preds = %97, %91, %171
%173 = phi i64 [ 0, %97 ], [ 0, %91 ], [ %34, %171 ]
%174 = xor i64 %173, -1
%175 = add nsw i64 %174, %28
br i1 %41, label %185, label %176
176: ; preds = %172, %176
%177 = phi i64 [ %182, %176 ], [ %173, %172 ]
%178 = phi i64 [ %183, %176 ], [ 0, %172 ]
%179 = getelementptr inbounds double, double* %94, i64 %177
%180 = load double, double* %179
%181 = getelementptr inbounds double, double* %96, i64 %177
store double %180, double* %181
%182 = add nuw nsw i64 %177, 1
%183 = add i64 %178, 1
%184 = icmp eq i64 %183, %40
br i1 %184, label %185, label %176
185: ; preds = %176, %172
%186 = phi i64 [ %173, %172 ], [ %182, %176 ]
%187 = icmp ult i64 %175, 3
br i1 %187, label %207, label %188
188: ; preds = %185, %188
%189 = phi i64 [ %205, %188 ], [ %186, %185 ]
%190 = getelementptr inbounds double, double* %94, i64 %189
%191 = load double, double* %190
%192 = getelementptr inbounds double, double* %96, i64 %189
store double %191, double* %192
%193 = add nuw nsw i64 %189, 1
%194 = getelementptr inbounds double, double* %94, i64 %193
%195 = load double, double* %194
%196 = getelementptr inbounds double, double* %96, i64 %193
store double %195, double* %196
%197 = add nuw nsw i64 %189, 2
%198 = getelementptr inbounds double, double* %94, i64 %197
%199 = load double, double* %198
%200 = getelementptr inbounds double, double* %96, i64 %197
store double %199, double* %200
%201 = add nuw nsw i64 %189, 3
%202 = getelementptr inbounds double, double* %94, i64 %201
%203 = load double, double* %202
%204 = getelementptr inbounds double, double* %96, i64 %201
store double %203, double* %204
%205 = add nuw nsw i64 %189, 4
%206 = icmp eq i64 %205, %28
br i1 %206, label %207, label %188
207: ; preds = %185, %188, %171
%208 = add nuw nsw i64 %92, 1
%209 = icmp eq i64 %208, %28
br i1 %209, label %210, label %91
210: ; preds = %207, %2, %26
%211 = call i32 (double**, i32, ...) bitcast (i32 (...)* @free_dmatrix to i32 (double**, i32, ...)*)(double** %7, i32 %1)
%212 = inttoptr i64 %12 to i8*
%213 = inttoptr i64 %9 to i8*
ret void
}