void CoFactor(double **a, int n, double **b) {
  int i, j, ii, jj, i1, j1;
  double det;
  double **c;
  c = malloc((n - 1) * sizeof(double *));
  for (i = 0; i < n - 1; i++)
    c[i] = malloc((n - 1) * sizeof(double));
  for (j = 0; j < n; j++) {
    for (i = 0; i < n; i++) {
      i1 = 0;
      for (ii = 0; ii < n; ii++) {
        if (ii == i)
          continue;
        j1 = 0;
        for (jj = 0; jj < n; jj++) {
          if (jj == j)
            continue;
          c[i1][j1] = a[ii][jj];
          j1++;
        }
        i1++;
      }
      det = Determinant(c, n - 1);
      b[i][j] = pow(-1.0, i + j + 2.0) * det;
    }
  }
  for (i = 0; i < n - 1; i++)
    free(c[i]);
  free(c);
}


x86 assembly (clang, -O0):
<CoFactor>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11b0 <CoFactor+0x70>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,%rdx
mov    -0x40(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1170 <CoFactor+0x30>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12fc <CoFactor+0x1bc>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12e9 <CoFactor+0x1a9>
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1284 <CoFactor+0x144>
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jne    1201 <CoFactor+0xc1>
jmp    1276 <CoFactor+0x136>
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    126d <CoFactor+0x12d>
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    122c <CoFactor+0xec>
jmp    125f <CoFactor+0x11f>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    120f <CoFactor+0xcf>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e4 <CoFactor+0xa4>
mov    -0x40(%rbp),%rdi
mov    -0xc(%rbp),%esi
sub    $0x1,%esi
mov    $0x0,%al
call   1040 <Determinant@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xd58(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xd44(%rip),%xmm0
call   1060 <pow@plt>
mulsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11ca <CoFactor+0x8a>
jmp    12ee <CoFactor+0x1ae>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11b7 <CoFactor+0x77>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1333 <CoFactor+0x1f3>
mov    -0x40(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
call   1030 <free@plt>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1303 <CoFactor+0x1c3>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<CoFactor>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,0x18(%rsp)
mov    %rdi,%r14
lea    -0x1(%rsi),%ebx
mov    %rsi,0x8(%rsp)
movslq %esi,%rbp
lea    -0x8(,%rbp,8),%r15
ff
mov    %r15,%rdi
call   1050 <malloc@plt>
mov    %rax,%r12
mov    %ebx,0x4(%rsp)
mov    %ebx,%ebx
cmp    $0x2,%ebp
jl     11a5 <CoFactor+0x65>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r15,%rdi
call   1050 <malloc@plt>
mov    %rax,(%r12,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%rbx
jne    1190 <CoFactor+0x50>
mov    %rbx,0x10(%rsp)
mov    0x8(%rsp),%r9
test   %r9d,%r9d
jle    1321 <CoFactor+0x1e1>
mov    %r9d,%r15d
mov    %r15d,%ebp
and    $0xfffffffe,%ebp
xor    %r13d,%r13d
jmp    11dd <CoFactor+0x9d>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r13
cmp    %r15,%r13
je     1321 <CoFactor+0x1e1>
xor    %r10d,%r10d
jmp    1258 <CoFactor+0x118>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r12,%rdi
mov    0x4(%rsp),%esi
xor    %eax,%eax
mov    %r10,%rbx
call   1040 <Determinant@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x20(%rsp)
lea    (%rbx,%r13,1),%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  0xddf(%rip),%xmm1
movsd  0xddf(%rip),%xmm0
call   1060 <pow@plt>
mov    %rbx,%r10
mov    0x8(%rsp),%r9
mulsd  0x20(%rsp),%xmm0
mov    0x18(%rsp),%rax
mov    (%rax,%rbx,8),%rax
movsd  %xmm0,(%rax,%r13,8)
add    $0x1,%r10
cmp    %r15,%r10
je     11d0 <CoFactor+0x90>
xor    %eax,%eax
xor    %r8d,%r8d
jmp    126d <CoFactor+0x12d>
nop
add    $0x1,%r8d
add    $0x1,%rax
cmp    %r15,%rax
je     11f0 <CoFactor+0xb0>
cmp    %r10,%rax
je     1264 <CoFactor+0x124>
movslq %r8d,%rdx
cmp    $0x1,%r9d
jne    1290 <CoFactor+0x150>
xor    %ecx,%ecx
xor    %esi,%esi
test   $0x1,%r15b
jne    12fe <CoFactor+0x1be>
jmp    1260 <CoFactor+0x120>
nopw   0x0(%rax,%rax,1)
xor    %edi,%edi
xor    %esi,%esi
jmp    12a9 <CoFactor+0x169>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %rbp,%rdi
je     12f0 <CoFactor+0x1b0>
mov    %rdi,%rcx
cmp    %rdi,%r13
je     12c9 <CoFactor+0x189>
mov    (%r14,%rax,8),%rdi
movsd  (%rdi,%rcx,8),%xmm0
mov    (%r12,%rdx,8),%rdi
movslq %esi,%rsi
movsd  %xmm0,(%rdi,%rsi,8)
add    $0x1,%esi
lea    0x1(%rcx),%rdi
cmp    %r13,%rdi
je     12a0 <CoFactor+0x160>
mov    (%r14,%rax,8),%rbx
movsd  0x8(%rbx,%rcx,8),%xmm0
mov    (%r12,%rdx,8),%rbx
movslq %esi,%rsi
movsd  %xmm0,(%rbx,%rsi,8)
add    $0x1,%esi
jmp    12a0 <CoFactor+0x160>
nopl   (%rax)
add    $0x2,%rcx
test   $0x1,%r15b
je     1260 <CoFactor+0x120>
cmp    %r13,%rcx
je     1260 <CoFactor+0x120>
mov    (%r14,%rax,8),%rdi
movsd  (%rdi,%rcx,8),%xmm0
mov    (%r12,%rdx,8),%rcx
movslq %esi,%rdx
movsd  %xmm0,(%rcx,%rdx,8)
jmp    1260 <CoFactor+0x120>
cmp    $0x2,%r9d
mov    0x10(%rsp),%rbp
jl     1342 <CoFactor+0x202>
xor    %ebx,%ebx
xchg   %ax,%ax
mov    (%r12,%rbx,8),%rdi
call   1030 <free@plt>
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    1330 <CoFactor+0x1f0>
mov    %r12,%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @CoFactor(double** %0, i32 %1, double** %2) {
%4 = alloca double**
%5 = alloca i32
%6 = alloca double**
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double**
store double** %0, double*** %4
store i32 %1, i32* %5
store double** %2, double*** %6
%15 = load i32, i32* %5
%16 = sub nsw i32 %15, 1
%17 = sext i32 %16 to i64
%18 = mul i64 %17, 8
%19 = call i8* @malloc(i64 %18)
%20 = bitcast i8* %19 to double**
store double** %20, double*** %14
store i32 0, i32* %7
br label %21
21: ; preds = %37, %3
%22 = load i32, i32* %7
%23 = load i32, i32* %5
%24 = sub nsw i32 %23, 1
%25 = icmp slt i32 %22, %24
br i1 %25, label %26, label %40
26: ; preds = %21
%27 = load i32, i32* %5
%28 = sub nsw i32 %27, 1
%29 = sext i32 %28 to i64
%30 = mul i64 %29, 8
%31 = call i8* @malloc(i64 %30)
%32 = bitcast i8* %31 to double*
%33 = load double**, double*** %14
%34 = load i32, i32* %7
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double*, double** %33, i64 %35
store double* %32, double** %36
br label %37
37: ; preds = %26
%38 = load i32, i32* %7
%39 = add nsw i32 %38, 1
store i32 %39, i32* %7
br label %21
40: ; preds = %21
store i32 0, i32* %8
br label %41
41: ; preds = %125, %40
%42 = load i32, i32* %8
%43 = load i32, i32* %5
%44 = icmp slt i32 %42, %43
br i1 %44, label %45, label %128
45: ; preds = %41
store i32 0, i32* %7
br label %46
46: ; preds = %121, %45
%47 = load i32, i32* %7
%48 = load i32, i32* %5
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %124
50: ; preds = %46
store i32 0, i32* %11
store i32 0, i32* %9
br label %51
51: ; preds = %96, %50
%52 = load i32, i32* %9
%53 = load i32, i32* %5
%54 = icmp slt i32 %52, %53
br i1 %54, label %55, label %99
55: ; preds = %51
%56 = load i32, i32* %9
%57 = load i32, i32* %7
%58 = icmp eq i32 %56, %57
br i1 %58, label %59, label %60
59: ; preds = %55
br label %96
60: ; preds = %55
store i32 0, i32* %12
store i32 0, i32* %10
br label %61
61: ; preds = %90, %60
%62 = load i32, i32* %10
%63 = load i32, i32* %5
%64 = icmp slt i32 %62, %63
br i1 %64, label %65, label %93
65: ; preds = %61
%66 = load i32, i32* %10
%67 = load i32, i32* %8
%68 = icmp eq i32 %66, %67
br i1 %68, label %69, label %70
69: ; preds = %65
br label %90
70: ; preds = %65
%71 = load double**, double*** %4
%72 = load i32, i32* %9
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double*, double** %71, i64 %73
%75 = load double*, double** %74
%76 = load i32, i32* %10
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = load double**, double*** %14
%81 = load i32, i32* %11
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds double*, double** %80, i64 %82
%84 = load double*, double** %83
%85 = load i32, i32* %12
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
store double %79, double* %87
%88 = load i32, i32* %12
%89 = add nsw i32 %88, 1
store i32 %89, i32* %12
br label %90
90: ; preds = %70, %69
%91 = load i32, i32* %10
%92 = add nsw i32 %91, 1
store i32 %92, i32* %10
br label %61
93: ; preds = %61
%94 = load i32, i32* %11
%95 = add nsw i32 %94, 1
store i32 %95, i32* %11
br label %96
96: ; preds = %93, %59
%97 = load i32, i32* %9
%98 = add nsw i32 %97, 1
store i32 %98, i32* %9
br label %51
99: ; preds = %51
%100 = load double**, double*** %14
%101 = load i32, i32* %5
%102 = sub nsw i32 %101, 1
%103 = call i32 (double**, i32, ...) bitcast (i32 (...)* @Determinant to i32 (double**, i32, ...)*)(double** %100, i32 %102)
%104 = sitofp i32 %103 to double
store double %104, double* %13
%105 = load i32, i32* %7
%106 = load i32, i32* %8
%107 = add nsw i32 %105, %106
%108 = sitofp i32 %107 to double
%109 = fadd double %108, 2.000000e+00
%110 = call double @pow(double -1.000000e+00, double %109)
%111 = load double, double* %13
%112 = fmul double %110, %111
%113 = load double**, double*** %6
%114 = load i32, i32* %7
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds double*, double** %113, i64 %115
%117 = load double*, double** %116
%118 = load i32, i32* %8
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %117, i64 %119
store double %112, double* %120
br label %121
121: ; preds = %99
%122 = load i32, i32* %7
%123 = add nsw i32 %122, 1
store i32 %123, i32* %7
br label %46
124: ; preds = %46
br label %125
125: ; preds = %124
%126 = load i32, i32* %8
%127 = add nsw i32 %126, 1
store i32 %127, i32* %8
br label %41
128: ; preds = %41
store i32 0, i32* %7
br label %129
129: ; preds = %141, %128
%130 = load i32, i32* %7
%131 = load i32, i32* %5
%132 = sub nsw i32 %131, 1
%133 = icmp slt i32 %130, %132
br i1 %133, label %134, label %144
134: ; preds = %129
%135 = load double**, double*** %14
%136 = load i32, i32* %7
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double*, double** %135, i64 %137
%139 = load double*, double** %138
%140 = bitcast double* %139 to i8*
br label %141
141: ; preds = %134
%142 = load i32, i32* %7
%143 = add nsw i32 %142, 1
store i32 %143, i32* %7
br label %129
144: ; preds = %129
%145 = load double**, double*** %14
%146 = bitcast double** %145 to i8*
ret void
}

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