int f01aaf(mata, matb)
double mata[4][4], matb[4][4];
{
  float **a, **b, d, y[4];
  int i, j, *indx;
  a = (float **)matrix(0, 3, 0, 3);
  b = (float **)matrix(0, 3, 0, 3);
  indx = (int *)ivector(0, 3);
  for (i = 1; i <= 3; i++) {
    for (j = 1; j <= 3; j++) {
      a[i][j] = (float)mata[i][j];
    }
  }
  for (i = 1; i <= 3; i++) {
    for (j = 1; j <= 3; j++) {
      b[i][j] = (float)0.0;
    }
    b[i][i] = (float)1.0;
  }
  if (!ludcmp(a, (int)3, indx, &d))
    return (0);
  for (i = 1; i <= 3; i++) {
    for (j = 1; j <= 3; j++) {
      y[j] = b[i][j];
    }
    lubksb(a, 3, indx, y);
    for (j = 1; j <= 3; j++) {
      b[i][j] = y[j];
    }
  }
  for (i = 1; i <= 3; i++) {
    for (j = 1; j <= 3; j++) {
      matb[i][j] = (double)b[j][i];
    }
  }
  free_matrix(a, 0, 3, 0, 3);
  free_matrix(b, 0, 3, 0, 3);
  free_ivector(indx, 0, 3);
  return (1);
}


x86 assembly (clang, -O0):
<f01aaf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
xor    %edx,%edx
mov    $0x3,%ecx
mov    %edx,%edi
mov    %ecx,%esi
mov    $0x0,%al
call   1060 <matrix@plt>
cltq
mov    %rax,-0x20(%rbp)
xor    %edx,%edx
mov    $0x3,%ecx
mov    %edx,%edi
mov    %ecx,%esi
mov    $0x0,%al
call   1060 <matrix@plt>
cltq
mov    %rax,-0x28(%rbp)
xor    %edi,%edi
mov    $0x3,%esi
mov    $0x0,%al
call   1080 <ivector@plt>
cltq
mov    %rax,-0x50(%rbp)
movl   $0x1,-0x44(%rbp)
cmpl   $0x3,-0x44(%rbp)
jg     1228 <f01aaf+0xc8>
movl   $0x1,-0x48(%rbp)
cmpl   $0x3,-0x48(%rbp)
jg     1215 <f01aaf+0xb5>
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    11cc <f01aaf+0x6c>
jmp    121a <f01aaf+0xba>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11bb <f01aaf+0x5b>
movl   $0x1,-0x44(%rbp)
cmpl   $0x3,-0x44(%rbp)
jg     129b <f01aaf+0x13b>
movl   $0x1,-0x48(%rbp)
cmpl   $0x3,-0x48(%rbp)
jg     1270 <f01aaf+0x110>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1240 <f01aaf+0xe0>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
movss  0xd78(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    122f <f01aaf+0xcf>
mov    -0x20(%rbp),%rdi
mov    -0x50(%rbp),%rdx
mov    $0x3,%esi
lea    -0x2c(%rbp),%rcx
mov    $0x0,%al
call   1030 <ludcmp@plt>
cmp    $0x0,%eax
jne    12c8 <f01aaf+0x168>
movl   $0x0,-0x4(%rbp)
jmp    143d <f01aaf+0x2dd>
movl   $0x1,-0x44(%rbp)
cmpl   $0x3,-0x44(%rbp)
jg     1380 <f01aaf+0x220>
movl   $0x1,-0x48(%rbp)
cmpl   $0x3,-0x48(%rbp)
jg     1317 <f01aaf+0x1b7>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x48(%rbp),%rax
movss  %xmm0,-0x40(%rbp,%rax,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    12e0 <f01aaf+0x180>
mov    -0x20(%rbp),%rdi
mov    -0x50(%rbp),%rdx
lea    -0x40(%rbp),%rcx
mov    $0x3,%esi
mov    $0x0,%al
call   1050 <lubksb@plt>
movl   $0x1,-0x48(%rbp)
cmpl   $0x3,-0x48(%rbp)
jg     136d <f01aaf+0x20d>
movslq -0x48(%rbp),%rax
movss  -0x40(%rbp,%rax,4),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1336 <f01aaf+0x1d6>
jmp    1372 <f01aaf+0x212>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12cf <f01aaf+0x16f>
movl   $0x1,-0x44(%rbp)
cmpl   $0x3,-0x44(%rbp)
jg     13f4 <f01aaf+0x294>
movl   $0x1,-0x48(%rbp)
cmpl   $0x3,-0x48(%rbp)
jg     13e1 <f01aaf+0x281>
mov    -0x28(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1398 <f01aaf+0x238>
jmp    13e6 <f01aaf+0x286>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1387 <f01aaf+0x227>
mov    -0x20(%rbp),%rdi
xor    %ecx,%ecx
mov    $0x3,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
mov    $0x0,%al
call   1040 <free_matrix@plt>
mov    -0x28(%rbp),%rdi
xor    %ecx,%ecx
mov    $0x3,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
mov    $0x0,%al
call   1040 <free_matrix@plt>
mov    -0x50(%rbp),%rdi
xor    %esi,%esi
mov    $0x3,%edx
mov    $0x0,%al
call   1070 <free_ivector@plt>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<f01aaf>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%r14
mov    %rdi,%rbp
xor    %r13d,%r13d
xor    %edi,%edi
mov    $0x3,%esi
xor    %edx,%edx
mov    $0x3,%ecx
xor    %eax,%eax
call   1060 <matrix@plt>
movslq %eax,%r15
xor    %edi,%edi
mov    $0x3,%esi
xor    %edx,%edx
mov    $0x3,%ecx
xor    %eax,%eax
call   1060 <matrix@plt>
movslq %eax,%rbx
xor    %edi,%edi
mov    $0x3,%esi
xor    %eax,%eax
call   1080 <ivector@plt>
mov    0x8(%r15),%rcx
movsd  0x28(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rcx)
movsd  0x30(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x8(%rcx)
movsd  0x38(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0xc(%rcx)
mov    0x10(%r15),%rcx
movsd  0x48(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rcx)
movsd  0x50(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x8(%rcx)
movsd  0x58(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0xc(%rcx)
mov    0x18(%r15),%rcx
movsd  0x68(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rcx)
movsd  0x70(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x8(%rcx)
movsd  0x78(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movslq %eax,%r12
movss  %xmm0,0xc(%rcx)
mov    0x8(%rbx),%rax
movq   $0x0,0x8(%rax)
movl   $0x3f800000,0x4(%rax)
mov    0x10(%rbx),%rax
movq   $0x0,0x4(%rax)
movl   $0x0,0xc(%rax)
movl   $0x3f800000,0x8(%rax)
mov    0x18(%rbx),%rax
movq   $0x0,0x4(%rax)
movl   $0x3f800000,0xc(%rax)
lea    0x14(%rsp),%rcx
mov    %r15,%rdi
mov    $0x3,%esi
mov    %r12,%rdx
xor    %eax,%eax
call   1030 <ludcmp@plt>
test   %eax,%eax
je     1475 <f01aaf+0x315>
mov    0x8(%rbx),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,0x4(%rsp)
movss  0x8(%rax),%xmm0
movss  %xmm0,0x8(%rsp)
movss  0xc(%rax),%xmm0
movss  %xmm0,0xc(%rsp)
mov    %rsp,%rbp
mov    %r15,%rdi
mov    $0x3,%esi
mov    %r12,%rdx
mov    %rbp,%rcx
xor    %eax,%eax
call   1050 <lubksb@plt>
mov    0x8(%rbx),%rax
movss  0x4(%rsp),%xmm0
movss  %xmm0,0x4(%rax)
movss  0x8(%rsp),%xmm0
movss  %xmm0,0x8(%rax)
movss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rax)
mov    0x10(%rbx),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,0x4(%rsp)
movss  0x8(%rax),%xmm0
movss  %xmm0,0x8(%rsp)
movss  0xc(%rax),%xmm0
movss  %xmm0,0xc(%rsp)
mov    %r15,%rdi
mov    $0x3,%esi
mov    %r12,%rdx
mov    %rbp,%rcx
xor    %eax,%eax
call   1050 <lubksb@plt>
mov    0x10(%rbx),%rax
movss  0x4(%rsp),%xmm0
movss  %xmm0,0x4(%rax)
movss  0x8(%rsp),%xmm0
movss  %xmm0,0x8(%rax)
movss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rax)
mov    0x18(%rbx),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,0x4(%rsp)
movss  0x8(%rax),%xmm0
movss  %xmm0,0x8(%rsp)
movss  0xc(%rax),%xmm0
movss  %xmm0,0xc(%rsp)
mov    %r15,%rdi
mov    $0x3,%esi
mov    %r12,%rdx
mov    %rbp,%rcx
xor    %eax,%eax
call   1050 <lubksb@plt>
mov    0x18(%rbx),%rax
movss  0x4(%rsp),%xmm0
movss  %xmm0,0x4(%rax)
movss  0x8(%rsp),%xmm1
movss  %xmm1,0x8(%rax)
movss  0xc(%rsp),%xmm2
movss  %xmm2,0xc(%rax)
mov    0x8(%rbx),%rax
mov    0x10(%rbx),%rcx
movss  0x4(%rcx),%xmm3
movss  0x4(%rax),%xmm4
unpcklps %xmm3,%xmm4
cvtps2pd %xmm4,%xmm3
movups %xmm3,0x28(%r14)
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x38(%r14)
movss  0x8(%rcx),%xmm0
movss  0x8(%rax),%xmm3
unpcklps %xmm0,%xmm3
cvtps2pd %xmm3,%xmm0
movups %xmm0,0x48(%r14)
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
movsd  %xmm0,0x58(%r14)
movss  0xc(%rcx),%xmm0
movss  0xc(%rax),%xmm1
unpcklps %xmm0,%xmm1
cvtps2pd %xmm1,%xmm0
movups %xmm0,0x68(%r14)
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
movsd  %xmm0,0x78(%r14)
mov    %r15,%rdi
xor    %esi,%esi
mov    $0x3,%edx
xor    %ecx,%ecx
mov    $0x3,%r8d
xor    %eax,%eax
call   1040 <free_matrix@plt>
mov    %rbx,%rdi
xor    %esi,%esi
mov    $0x3,%edx
xor    %ecx,%ecx
mov    $0x3,%r8d
xor    %eax,%eax
call   1040 <free_matrix@plt>
mov    %r12,%rdi
xor    %esi,%esi
mov    $0x3,%edx
xor    %eax,%eax
call   1070 <free_ivector@plt>
mov    $0x1,%r13d
mov    %r13d,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @f01aaf([4 x double]* %0, [4 x double]* %1) {
%3 = alloca i32
%4 = alloca [4 x double]*
%5 = alloca [4 x double]*
%6 = alloca float**
%7 = alloca float**
%8 = alloca float
%9 = alloca [4 x float]
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32*
store [4 x double]* %0, [4 x double]** %4
store [4 x double]* %1, [4 x double]** %5
%13 = call i32 (i32, i32, i32, i32, ...) bitcast (i32 (...)* @matrix to i32 (i32, i32, i32, i32, ...)*)(i32 0, i32 3, i32 0, i32 3)
%14 = sext i32 %13 to i64
%15 = inttoptr i64 %14 to float**
store float** %15, float*** %6
%16 = call i32 (i32, i32, i32, i32, ...) bitcast (i32 (...)* @matrix to i32 (i32, i32, i32, i32, ...)*)(i32 0, i32 3, i32 0, i32 3)
%17 = sext i32 %16 to i64
%18 = inttoptr i64 %17 to float**
store float** %18, float*** %7
%19 = call i32 (i32, i32, ...) bitcast (i32 (...)* @ivector to i32 (i32, i32, ...)*)(i32 0, i32 3)
%20 = sext i32 %19 to i64
%21 = inttoptr i64 %20 to i32*
store i32* %21, i32** %12
store i32 1, i32* %10
br label %22
22: ; preds = %51, %2
%23 = load i32, i32* %10
%24 = icmp sle i32 %23, 3
br i1 %24, label %25, label %54
25: ; preds = %22
store i32 1, i32* %11
br label %26
26: ; preds = %47, %25
%27 = load i32, i32* %11
%28 = icmp sle i32 %27, 3
br i1 %28, label %29, label %50
29: ; preds = %26
%30 = load [4 x double]*, [4 x double]** %4
%31 = load i32, i32* %10
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds [4 x double], [4 x double]* %30, i64 %32
%34 = load i32, i32* %11
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds [4 x double], [4 x double]* %33, i64 0, i64 %35
%37 = load double, double* %36
%38 = fptrunc double %37 to float
%39 = load float**, float*** %6
%40 = load i32, i32* %10
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float*, float** %39, i64 %41
%43 = load float*, float** %42
%44 = load i32, i32* %11
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds float, float* %43, i64 %45
store float %38, float* %46
br label %47
47: ; preds = %29
%48 = load i32, i32* %11
%49 = add nsw i32 %48, 1
store i32 %49, i32* %11
br label %26
50: ; preds = %26
br label %51
51: ; preds = %50
%52 = load i32, i32* %10
%53 = add nsw i32 %52, 1
store i32 %53, i32* %10
br label %22
54: ; preds = %22
store i32 1, i32* %10
br label %55
55: ; preds = %83, %54
%56 = load i32, i32* %10
%57 = icmp sle i32 %56, 3
br i1 %57, label %58, label %86
58: ; preds = %55
store i32 1, i32* %11
br label %59
59: ; preds = %71, %58
%60 = load i32, i32* %11
%61 = icmp sle i32 %60, 3
br i1 %61, label %62, label %74
62: ; preds = %59
%63 = load float**, float*** %7
%64 = load i32, i32* %10
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds float*, float** %63, i64 %65
%67 = load float*, float** %66
%68 = load i32, i32* %11
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds float, float* %67, i64 %69
store float 0.000000e+00, float* %70
br label %71
71: ; preds = %62
%72 = load i32, i32* %11
%73 = add nsw i32 %72, 1
store i32 %73, i32* %11
br label %59
74: ; preds = %59
%75 = load float**, float*** %7
%76 = load i32, i32* %10
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float*, float** %75, i64 %77
%79 = load float*, float** %78
%80 = load i32, i32* %10
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %79, i64 %81
store float 1.000000e+00, float* %82
br label %83
83: ; preds = %74
%84 = load i32, i32* %10
%85 = add nsw i32 %84, 1
store i32 %85, i32* %10
br label %55
86: ; preds = %55
%87 = load float**, float*** %6
%88 = load i32*, i32** %12
%89 = call i32 (float**, i32, i32*, float*, ...) bitcast (i32 (...)* @ludcmp to i32 (float**, i32, i32*, float*, ...)*)(float** %87, i32 3, i32* %88, float* %8)
%90 = icmp ne i32 %89, 0
br i1 %90, label %92, label %91
91: ; preds = %86
store i32 0, i32* %3
br label %184
92: ; preds = %86
store i32 1, i32* %10
br label %93
93: ; preds = %141, %92
%94 = load i32, i32* %10
%95 = icmp sle i32 %94, 3
br i1 %95, label %96, label %144
96: ; preds = %93
store i32 1, i32* %11
br label %97
97: ; preds = %113, %96
%98 = load i32, i32* %11
%99 = icmp sle i32 %98, 3
br i1 %99, label %100, label %116
100: ; preds = %97
%101 = load float**, float*** %7
%102 = load i32, i32* %10
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds float*, float** %101, i64 %103
%105 = load float*, float** %104
%106 = load i32, i32* %11
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds float, float* %105, i64 %107
%109 = load float, float* %108
%110 = load i32, i32* %11
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds [4 x float], [4 x float]* %9, i64 0, i64 %111
store float %109, float* %112
br label %113
113: ; preds = %100
%114 = load i32, i32* %11
%115 = add nsw i32 %114, 1
store i32 %115, i32* %11
br label %97
116: ; preds = %97
%117 = load float**, float*** %6
%118 = load i32*, i32** %12
%119 = getelementptr inbounds [4 x float], [4 x float]* %9, i64 0, i64 0
%120 = call i32 (float**, i32, i32*, float*, ...) bitcast (i32 (...)* @lubksb to i32 (float**, i32, i32*, float*, ...)*)(float** %117, i32 3, i32* %118, float* %119)
store i32 1, i32* %11
br label %121
121: ; preds = %137, %116
%122 = load i32, i32* %11
%123 = icmp sle i32 %122, 3
br i1 %123, label %124, label %140
124: ; preds = %121
%125 = load i32, i32* %11
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds [4 x float], [4 x float]* %9, i64 0, i64 %126
%128 = load float, float* %127
%129 = load float**, float*** %7
%130 = load i32, i32* %10
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds float*, float** %129, i64 %131
%133 = load float*, float** %132
%134 = load i32, i32* %11
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds float, float* %133, i64 %135
store float %128, float* %136
br label %137
137: ; preds = %124
%138 = load i32, i32* %11
%139 = add nsw i32 %138, 1
store i32 %139, i32* %11
br label %121
140: ; preds = %121
br label %141
141: ; preds = %140
%142 = load i32, i32* %10
%143 = add nsw i32 %142, 1
store i32 %143, i32* %10
br label %93
144: ; preds = %93
store i32 1, i32* %10
br label %145
145: ; preds = %174, %144
%146 = load i32, i32* %10
%147 = icmp sle i32 %146, 3
br i1 %147, label %148, label %177
148: ; preds = %145
store i32 1, i32* %11
br label %149
149: ; preds = %170, %148
%150 = load i32, i32* %11
%151 = icmp sle i32 %150, 3
br i1 %151, label %152, label %173
152: ; preds = %149
%153 = load float**, float*** %7
%154 = load i32, i32* %11
%155 = sext i32 %154 to i64
%156 = getelementptr inbounds float*, float** %153, i64 %155
%157 = load float*, float** %156
%158 = load i32, i32* %10
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds float, float* %157, i64 %159
%161 = load float, float* %160
%162 = fpext float %161 to double
%163 = load [4 x double]*, [4 x double]** %5
%164 = load i32, i32* %10
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds [4 x double], [4 x double]* %163, i64 %165
%167 = load i32, i32* %11
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds [4 x double], [4 x double]* %166, i64 0, i64 %168
store double %162, double* %169
br label %170
170: ; preds = %152
%171 = load i32, i32* %11
%172 = add nsw i32 %171, 1
store i32 %172, i32* %11
br label %149
173: ; preds = %149
br label %174
174: ; preds = %173
%175 = load i32, i32* %10
%176 = add nsw i32 %175, 1
store i32 %176, i32* %10
br label %145
177: ; preds = %145
%178 = load float**, float*** %6
%179 = call i32 (float**, i32, i32, i32, i32, ...) bitcast (i32 (...)* @free_matrix to i32 (float**, i32, i32, i32, i32, ...)*)(float** %178, i32 0, i32 3, i32 0, i32 3)
%180 = load float**, float*** %7
%181 = call i32 (float**, i32, i32, i32, i32, ...) bitcast (i32 (...)* @free_matrix to i32 (float**, i32, i32, i32, i32, ...)*)(float** %180, i32 0, i32 3, i32 0, i32 3)
%182 = load i32*, i32** %12
%183 = call i32 (i32*, i32, i32, ...) bitcast (i32 (...)* @free_ivector to i32 (i32*, i32, i32, ...)*)(i32* %182, i32 0, i32 3)
store i32 1, i32* %3
br label %184
184: ; preds = %177, %91
%185 = load i32, i32* %3
ret i32 %185
}

LLVM IR (-O2):
define i32 @f01aaf([4 x double]* nocapture %0, [4 x double]* nocapture writeonly %1) {
%3 = alloca float
%4 = alloca [4 x float]
%5 = bitcast float* %3 to i8*
%6 = bitcast [4 x float]* %4 to i8*
%7 = tail call i32 (i32, i32, i32, i32, ...) bitcast (i32 (...)* @matrix to i32 (i32, i32, i32, i32, ...)*)(i32 0, i32 3, i32 0, i32 3)
%8 = sext i32 %7 to i64
%9 = inttoptr i64 %8 to float**
%10 = tail call i32 (i32, i32, i32, i32, ...) bitcast (i32 (...)* @matrix to i32 (i32, i32, i32, i32, ...)*)(i32 0, i32 3, i32 0, i32 3)
%11 = sext i32 %10 to i64
%12 = tail call i32 (i32, i32, ...) bitcast (i32 (...)* @ivector to i32 (i32, i32, ...)*)(i32 0, i32 3)
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds float*, float** %9, i64 1
%15 = load float*, float** %14
%16 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 1
%17 = load double, double* %16
%18 = fptrunc double %17 to float
%19 = getelementptr inbounds float, float* %15, i64 1
store float %18, float* %19
%20 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 2
%21 = load double, double* %20
%22 = fptrunc double %21 to float
%23 = getelementptr inbounds float, float* %15, i64 2
store float %22, float* %23
%24 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 3
%25 = load double, double* %24
%26 = fptrunc double %25 to float
%27 = getelementptr inbounds float, float* %15, i64 3
store float %26, float* %27
%28 = getelementptr inbounds float*, float** %9, i64 2
%29 = load float*, float** %28
%30 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 1
%31 = load double, double* %30
%32 = fptrunc double %31 to float
%33 = getelementptr inbounds float, float* %29, i64 1
store float %32, float* %33
%34 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 2
%35 = load double, double* %34
%36 = fptrunc double %35 to float
%37 = getelementptr inbounds float, float* %29, i64 2
store float %36, float* %37
%38 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 3
%39 = load double, double* %38
%40 = fptrunc double %39 to float
%41 = getelementptr inbounds float, float* %29, i64 3
store float %40, float* %41
%42 = getelementptr inbounds float*, float** %9, i64 3
%43 = load float*, float** %42
%44 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 1
%45 = load double, double* %44
%46 = fptrunc double %45 to float
%47 = getelementptr inbounds float, float* %43, i64 1
store float %46, float* %47
%48 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 2
%49 = load double, double* %48
%50 = fptrunc double %49 to float
%51 = getelementptr inbounds float, float* %43, i64 2
store float %50, float* %51
%52 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 3
%53 = load double, double* %52
%54 = fptrunc double %53 to float
%55 = getelementptr inbounds float, float* %43, i64 3
store float %54, float* %55
%56 = inttoptr i64 %11 to float**
%57 = getelementptr inbounds float*, float** %56, i64 1
%58 = load float*, float** %57
%59 = getelementptr float, float* %58, i64 1
%60 = getelementptr float, float* %58, i64 2
%61 = bitcast float* %60 to i64*
store i64 0, i64* %61
store float 1.000000e+00, float* %59
%62 = getelementptr inbounds float*, float** %56, i64 2
%63 = load float*, float** %62
%64 = getelementptr float, float* %63, i64 1
%65 = bitcast float* %64 to i8*
%66 = getelementptr inbounds float, float* %63, i64 2
store float 1.000000e+00, float* %66
%67 = getelementptr inbounds float*, float** %56, i64 3
%68 = load float*, float** %67
%69 = getelementptr float, float* %68, i64 1
%70 = bitcast float* %69 to i64*
store i64 0, i64* %70
%71 = getelementptr inbounds float, float* %68, i64 3
store float 1.000000e+00, float* %71
%72 = inttoptr i64 %13 to i32*
%73 = call i32 (float**, i32, i32*, float*, ...) bitcast (i32 (...)* @ludcmp to i32 (float**, i32, i32*, float*, ...)*)(float** %9, i32 3, i32* %72, float* %3)
%74 = icmp eq i32 %73, 0
br i1 %74, label %163, label %75
75: ; preds = %2
%76 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 0
%77 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 1
%78 = load float*, float** %57
%79 = getelementptr inbounds float, float* %78, i64 1
%80 = load float, float* %79
store float %80, float* %77
%81 = getelementptr inbounds float, float* %78, i64 2
%82 = load float, float* %81
%83 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 2
store float %82, float* %83
%84 = getelementptr inbounds float, float* %78, i64 3
%85 = load float, float* %84
%86 = getelementptr inbounds [4 x float], [4 x float]* %4, i64 0, i64 3
store float %85, float* %86
%87 = call i32 (float**, i32, i32*, float*, ...) bitcast (i32 (...)* @lubksb to i32 (float**, i32, i32*, float*, ...)*)(float** %9, i32 3, i32* %72, float* %76)
%88 = load float*, float** %57
%89 = load float, float* %77
%90 = getelementptr inbounds float, float* %88, i64 1
store float %89, float* %90
%91 = load float, float* %83
%92 = getelementptr inbounds float, float* %88, i64 2
store float %91, float* %92
%93 = load float, float* %86
%94 = getelementptr inbounds float, float* %88, i64 3
store float %93, float* %94
%95 = load float*, float** %62
%96 = getelementptr inbounds float, float* %95, i64 1
%97 = load float, float* %96
store float %97, float* %77
%98 = getelementptr inbounds float, float* %95, i64 2
%99 = load float, float* %98
store float %99, float* %83
%100 = getelementptr inbounds float, float* %95, i64 3
%101 = load float, float* %100
store float %101, float* %86
%102 = call i32 (float**, i32, i32*, float*, ...) bitcast (i32 (...)* @lubksb to i32 (float**, i32, i32*, float*, ...)*)(float** %9, i32 3, i32* %72, float* %76)
%103 = load float*, float** %62
%104 = load float, float* %77
%105 = getelementptr inbounds float, float* %103, i64 1
store float %104, float* %105
%106 = load float, float* %83
%107 = getelementptr inbounds float, float* %103, i64 2
store float %106, float* %107
%108 = load float, float* %86
%109 = getelementptr inbounds float, float* %103, i64 3
store float %108, float* %109
%110 = load float*, float** %67
%111 = getelementptr inbounds float, float* %110, i64 1
%112 = load float, float* %111
store float %112, float* %77
%113 = getelementptr inbounds float, float* %110, i64 2
%114 = load float, float* %113
store float %114, float* %83
%115 = getelementptr inbounds float, float* %110, i64 3
%116 = load float, float* %115
store float %116, float* %86
%117 = call i32 (float**, i32, i32*, float*, ...) bitcast (i32 (...)* @lubksb to i32 (float**, i32, i32*, float*, ...)*)(float** %9, i32 3, i32* %72, float* %76)
%118 = load float*, float** %67
%119 = load float, float* %77
%120 = getelementptr inbounds float, float* %118, i64 1
store float %119, float* %120
%121 = load float, float* %83
%122 = getelementptr inbounds float, float* %118, i64 2
store float %121, float* %122
%123 = load float, float* %86
%124 = getelementptr inbounds float, float* %118, i64 3
store float %123, float* %124
%125 = load float*, float** %57
%126 = getelementptr inbounds float, float* %125, i64 1
%127 = load float, float* %126
%128 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 1, i64 1
%129 = load float*, float** %62
%130 = getelementptr inbounds float, float* %129, i64 1
%131 = load float, float* %130
%132 = insertelement <2 x float> poison, float %127, i64 0
%133 = insertelement <2 x float> %132, float %131, i64 1
%134 = fpext <2 x float> %133 to <2 x double>
%135 = bitcast double* %128 to <2 x double>*
store <2 x double> %134, <2 x double>* %135
%136 = fpext float %119 to double
%137 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 1, i64 3
store double %136, double* %137
%138 = getelementptr inbounds float, float* %125, i64 2
%139 = load float, float* %138
%140 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 2, i64 1
%141 = getelementptr inbounds float, float* %129, i64 2
%142 = load float, float* %141
%143 = insertelement <2 x float> poison, float %139, i64 0
%144 = insertelement <2 x float> %143, float %142, i64 1
%145 = fpext <2 x float> %144 to <2 x double>
%146 = bitcast double* %140 to <2 x double>*
store <2 x double> %145, <2 x double>* %146
%147 = fpext float %121 to double
%148 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 2, i64 3
store double %147, double* %148
%149 = getelementptr inbounds float, float* %125, i64 3
%150 = load float, float* %149
%151 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 3, i64 1
%152 = getelementptr inbounds float, float* %129, i64 3
%153 = load float, float* %152
%154 = insertelement <2 x float> poison, float %150, i64 0
%155 = insertelement <2 x float> %154, float %153, i64 1
%156 = fpext <2 x float> %155 to <2 x double>
%157 = bitcast double* %151 to <2 x double>*
store <2 x double> %156, <2 x double>* %157
%158 = fpext float %123 to double
%159 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 3, i64 3
store double %158, double* %159
%160 = call i32 (float**, i32, i32, i32, i32, ...) bitcast (i32 (...)* @free_matrix to i32 (float**, i32, i32, i32, i32, ...)*)(float** %9, i32 0, i32 3, i32 0, i32 3)
%161 = call i32 (float**, i32, i32, i32, i32, ...) bitcast (i32 (...)* @free_matrix to i32 (float**, i32, i32, i32, i32, ...)*)(float** %56, i32 0, i32 3, i32 0, i32 3)
%162 = call i32 (i32*, i32, i32, ...) bitcast (i32 (...)* @free_ivector to i32 (i32*, i32, i32, ...)*)(i32* %72, i32 0, i32 3)
br label %163
163: ; preds = %2, %75
%164 = phi i32 [ 1, %75 ], [ 0, %2 ]
ret i32 %164
}