float determinant(float b[50][50], int m) {
  int i, j, p, sign = 1, k;
  float c[50][50];
  float sum = 0;
  if (m == 2) {
    sum = b[0][0] * b[1][1] - b[0][1] * b[1][0];
    return sum;
  }
  for (p = 0; p < m; p++) {
    int h = 0, k = 0;
    for (i = 1; i < m; i++) {
      for (j = 0; j < m; j++) {
        if (j == p)
          continue;
        c[h][k] = b[i][j];
        k++;
        if (k == m - 1) {
          h++;
          k = 0;
        }
      }
    }
    for (k = 0; k < p; k++)
      sign = sign * (-1);
    sum = sum + b[0][p] * sign * determinant(c, m - 1);
  }
  return sum;
}


x86 assembly (clang, -O0):
<determinant>:
push   %rbp
mov    %rsp,%rbp
sub    $0x2760,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movl   $0x1,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2744(%rbp)
ff
cmpl   $0x2,-0x14(%rbp)
jne    1196 <determinant+0x86>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  0xcc(%rax),%xmm2
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x10(%rbp),%rax
mulss  0xc8(%rax),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2744(%rbp)
ff
movss  -0x2744(%rbp),%xmm0
ff
movss  %xmm0,-0x4(%rbp)
jmp    1347 <determinant+0x237>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    133a <determinant+0x22a>
movl   $0x0,-0x2748(%rbp)
movl   $0x0,-0x274c(%rbp)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1295 <determinant+0x185>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1282 <determinant+0x172>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    11f4 <determinant+0xe4>
jmp    1274 <determinant+0x164>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0xc8,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x2748(%rbp),%rcx
lea    -0x2740(%rbp),%rax
imul   $0xc8,%rcx,%rcx
add    %rcx,%rax
movslq -0x274c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x274c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x274c(%rbp)
mov    -0x274c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jne    126f <determinant+0x15f>
mov    -0x2748(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2748(%rbp)
movl   $0x0,-0x274c(%rbp)
jmp    1274 <determinant+0x164>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11d7 <determinant+0xc7>
jmp    1287 <determinant+0x177>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11c4 <determinant+0xb4>
movl   $0x0,-0x274c(%rbp)
mov    -0x274c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    12c9 <determinant+0x1b9>
imul   $0xffffffff,-0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x274c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x274c(%rbp)
jmp    129f <determinant+0x18f>
movss  -0x2744(%rbp),%xmm0
ff
movss  %xmm0,-0x2754(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x2750(%rbp)
ff
lea    -0x2740(%rbp),%rdi
mov    -0x14(%rbp),%esi
sub    $0x1,%esi
call   1030 <determinant@plt>
movss  -0x2754(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movss  -0x2750(%rbp),%xmm0
ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2744(%rbp)
ff
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    119d <determinant+0x8d>
movss  -0x2744(%rbp),%xmm0
ff
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x2760,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<determinant>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x2738,%rsp
mov    %rdi,%r12
cmp    $0x2,%esi
jne    1153 <determinant+0x43>
movss  0xc8(%r12),%xmm1
mulss  0x4(%r12),%xmm1
movss  (%r12),%xmm0
mulss  0xcc(%r12),%xmm0
subss  %xmm1,%xmm0
jmp    12b3 <determinant+0x1a3>
test   %esi,%esi
jle    12ad <determinant+0x19d>
lea    -0x1(%rsi),%ebp
mov    %esi,%ebx
lea    0xc8(%r12),%rax
mov    %rax,0x10(%rsp)
mov    $0x1,%r15d
xorps  %xmm1,%xmm1
xor    %r13d,%r13d
xor    %r14d,%r14d
mov    %rsi,0x18(%rsp)
jmp    11d2 <determinant+0xc2>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %r15d,%xmm0
mulss  (%r12,%r14,4),%xmm0
movss  %xmm0,0x8(%rsp)
lea    0x20(%rsp),%rdi
mov    %ebp,%esi
call   1030 <determinant@plt>
mulss  0x8(%rsp),%xmm0
movss  0xc(%rsp),%xmm1
addss  %xmm0,%xmm1
add    $0x1,%r14
cmp    %rbx,%r14
mov    0x18(%rsp),%rsi
je     12b0 <determinant+0x1a0>
cmp    $0x1,%esi
jne    1230 <determinant+0x120>
test   %r14,%r14
movss  %xmm1,0xc(%rsp)
je     1190 <determinant+0x80>
lea    -0x1(%r14),%rax
cmp    $0x7,%rax
jb     1206 <determinant+0xf6>
mov    %r14,%rax
and    $0xfffffffffffffff8,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0xfffffffffffffff8,%rax
jne    1200 <determinant+0xf0>
test   $0x7,%r14b
je     1190 <determinant+0x80>
mov    %r14d,%eax
and    $0x7,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
neg    %r15d
add    $0xffffffffffffffff,%rax
jne    1220 <determinant+0x110>
jmp    1190 <determinant+0x80>
xchg   %ax,%ax
xor    %eax,%eax
mov    $0x1,%r8d
mov    0x10(%rsp),%rdx
xor    %esi,%esi
jmp    1264 <determinant+0x154>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
add    $0xc8,%rdx
cmp    %rbx,%r8
je     11d7 <determinant+0xc7>
xor    %edi,%edi
jmp    1279 <determinant+0x169>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %rdi,%rbx
je     1250 <determinant+0x140>
cmp    %rdi,%r14
je     1270 <determinant+0x160>
movss  (%rdx,%rdi,4),%xmm0
movslq %esi,%rsi
cltq
imul   $0xc8,%rsi,%rcx
add    %rsp,%rcx
add    $0x20,%rcx
movss  %xmm0,(%rcx,%rax,4)
add    $0x1,%eax
xor    %ecx,%ecx
cmp    %ebp,%eax
sete   %cl
cmove  %r13d,%eax
add    %ecx,%esi
jmp    1270 <determinant+0x160>
xorps  %xmm1,%xmm1
movaps %xmm1,%xmm0
add    $0x2738,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float @determinant([50 x float]* nocapture %0, i32 %1) {
%3 = alloca [50 x [50 x float]]
%4 = bitcast [50 x [50 x float]]* %3 to i8*
%5 = icmp eq i32 %1, 2
br i1 %5, label %13, label %6
6: ; preds = %2
%7 = icmp sgt i32 %1, 0
br i1 %7, label %8, label %87
8: ; preds = %6
%9 = icmp eq i32 %1, 1
%10 = add nsw i32 %1, -1
%11 = getelementptr inbounds [50 x [50 x float]], [50 x [50 x float]]* %3, i64 0, i64 0
%12 = zext i32 %1 to i64
br label %25
13: ; preds = %2
%14 = getelementptr inbounds [50 x float], [50 x float]* %0, i64 0, i64 0
%15 = load float, float* %14
%16 = getelementptr inbounds [50 x float], [50 x float]* %0, i64 1, i64 1
%17 = load float, float* %16
%18 = getelementptr inbounds [50 x float], [50 x float]* %0, i64 0, i64 1
%19 = load float, float* %18
%20 = getelementptr inbounds [50 x float], [50 x float]* %0, i64 1, i64 0
%21 = load float, float* %20
%22 = fneg float %19
%23 = fmul float %21, %22
%24 = tail call float @llvm.fmuladd.f32(float %15, float %17, float %23)
br label %87
25: ; preds = %8, %77
%26 = phi i64 [ 0, %8 ], [ %85, %77 ]
%27 = phi float [ 0.000000e+00, %8 ], [ %84, %77 ]
%28 = phi i32 [ 1, %8 ], [ %78, %77 ]
%29 = add nsw i64 %26, -1
br i1 %9, label %30, label %37
30: ; preds = %62, %25
%31 = icmp eq i64 %26, 0
br i1 %31, label %77, label %32
32: ; preds = %30
%33 = and i64 %26, 7
%34 = icmp ult i64 %29, 7
br i1 %34, label %69, label %35
35: ; preds = %32
%36 = and i64 %26, 9223372036854775800
br label %65
37: ; preds = %25, %62
%38 = phi i64 [ %63, %62 ], [ 1, %25 ]
%39 = phi i32 [ %59, %62 ], [ 0, %25 ]
%40 = phi i32 [ %58, %62 ], [ 0, %25 ]
br label %41
41: ; preds = %37, %57
%42 = phi i64 [ 0, %37 ], [ %60, %57 ]
%43 = phi i32 [ %39, %37 ], [ %59, %57 ]
%44 = phi i32 [ %40, %37 ], [ %58, %57 ]
%45 = icmp eq i64 %42, %26
br i1 %45, label %57, label %46
46: ; preds = %41
%47 = getelementptr inbounds [50 x float], [50 x float]* %0, i64 %38, i64 %42
%48 = load float, float* %47
%49 = sext i32 %44 to i64
%50 = sext i32 %43 to i64
%51 = getelementptr inbounds [50 x [50 x float]], [50 x [50 x float]]* %3, i64 0, i64 %49, i64 %50
store float %48, float* %51
%52 = add nsw i32 %43, 1
%53 = icmp eq i32 %52, %10
%54 = zext i1 %53 to i32
%55 = add nsw i32 %44, %54
%56 = select i1 %53, i32 0, i32 %52
br label %57
57: ; preds = %46, %41
%58 = phi i32 [ %44, %41 ], [ %55, %46 ]
%59 = phi i32 [ %43, %41 ], [ %56, %46 ]
%60 = add nuw nsw i64 %42, 1
%61 = icmp eq i64 %60, %12
br i1 %61, label %62, label %41
62: ; preds = %57
%63 = add nuw nsw i64 %38, 1
%64 = icmp eq i64 %63, %12
br i1 %64, label %30, label %37
65: ; preds = %65, %35
%66 = phi i64 [ 0, %35 ], [ %67, %65 ]
%67 = add i64 %66, 8
%68 = icmp eq i64 %67, %36
br i1 %68, label %69, label %65
69: ; preds = %65, %32
%70 = icmp eq i64 %33, 0
br i1 %70, label %77, label %71
71: ; preds = %69, %71
%72 = phi i32 [ %74, %71 ], [ %28, %69 ]
%73 = phi i64 [ %75, %71 ], [ 0, %69 ]
%74 = sub nsw i32 0, %72
%75 = add i64 %73, 1
%76 = icmp eq i64 %75, %33
br i1 %76, label %77, label %71
77: ; preds = %69, %71, %30
%78 = phi i32 [ %28, %30 ], [ %28, %69 ], [ %74, %71 ]
%79 = getelementptr inbounds [50 x float], [50 x float]* %0, i64 0, i64 %26
%80 = load float, float* %79
%81 = sitofp i32 %78 to float
%82 = fmul float %80, %81
%83 = call float @determinant([50 x float]* %11, i32 %10)
%84 = tail call float @llvm.fmuladd.f32(float %82, float %83, float %27)
%85 = add nuw nsw i64 %26, 1
%86 = icmp eq i64 %85, %12
br i1 %86, label %87, label %25
87: ; preds = %77, %6, %13
%88 = phi float [ %24, %13 ], [ 0.000000e+00, %6 ], [ %84, %77 ]
ret float %88
}