#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 */
double Det3x3(double, double, double, double, double, double, double, double,
              double);
double Determinant44(double a[4][4]) {
  double a1, a2, a3, a4, b1, b2, b3, b4, c1, c2, c3, c4, d1, d2, d3, d4;
  a1 = a[0][0];
  b1 = a[0][1];
  c1 = a[0][2];
  d1 = a[0][3];
  a2 = a[1][0];
  b2 = a[1][1];
  c2 = a[1][2];
  d2 = a[1][3];
  a3 = a[2][0];
  b3 = a[2][1];
  c3 = a[2][2];
  d3 = a[2][3];
  a4 = a[3][0];
  b4 = a[3][1];
  c4 = a[3][2];
  d4 = a[3][3];
  return (a1 * Det3x3(b2, b3, b4, c2, c3, c4, d2, d3, d4) -
          b1 * Det3x3(a2, a3, a4, c2, c3, c4, d2, d3, d4) +
          c1 * Det3x3(a2, a3, a4, b2, b3, b4, d2, d3, d4) -
          d1 * Det3x3(a2, a3, a4, b2, b3, b4, c2, c3, c4));
}


x86 assembly (clang, -O0):
<Determinant44>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
mov    %rdi,-0x8(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x48(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x50(%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x58(%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x60(%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x68(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x70(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x78(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm3
movsd  -0x60(%rbp),%xmm4
movsd  -0x68(%rbp),%xmm5
movsd  -0x78(%rbp),%xmm6
movsd  -0x80(%rbp),%xmm7
movsd  -0x88(%rbp),%xmm8
ff ff
movsd  %xmm8,(%rsp)
call   1030 <Det3x3@plt>
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm3
movsd  -0x60(%rbp),%xmm4
movsd  -0x68(%rbp),%xmm5
movsd  -0x78(%rbp),%xmm6
movsd  -0x80(%rbp),%xmm7
movsd  -0x88(%rbp),%xmm8
ff ff
movsd  %xmm8,(%rsp)
call   1030 <Det3x3@plt>
movsd  -0xc0(%rbp),%xmm1
ff
movsd  -0xb8(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0xb0(%rbp),%xmm0
ff
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm4
movsd  -0x48(%rbp),%xmm5
movsd  -0x78(%rbp),%xmm6
movsd  -0x80(%rbp),%xmm7
movsd  -0x88(%rbp),%xmm8
ff ff
movsd  %xmm8,(%rsp)
call   1030 <Det3x3@plt>
movsd  -0xa8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xa0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm4
movsd  -0x48(%rbp),%xmm5
movsd  -0x58(%rbp),%xmm6
movsd  -0x60(%rbp),%xmm7
movsd  -0x68(%rbp),%xmm8
movsd  %xmm8,(%rsp)
call   1030 <Det3x3@plt>
movsd  -0x98(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x90(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Determinant44>:
sub    $0x98,%rsp
movsd  (%rdi),%xmm0
movsd  %xmm0,0x80(%rsp)
movsd  0x8(%rdi),%xmm0
movsd  %xmm0,0x78(%rsp)
movsd  0x10(%rdi),%xmm0
movsd  %xmm0,0x88(%rsp)
movsd  0x18(%rdi),%xmm0
movsd  %xmm0,0x90(%rsp)
movsd  0x20(%rdi),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0x28(%rdi),%xmm0
movsd  %xmm0,0x68(%rsp)
movsd  0x30(%rdi),%xmm3
movsd  %xmm3,0x70(%rsp)
movsd  0x38(%rdi),%xmm6
movsd  %xmm6,0x40(%rsp)
movsd  0x40(%rdi),%xmm1
movsd  %xmm1,0x20(%rsp)
movsd  0x48(%rdi),%xmm1
movsd  %xmm1,0x58(%rsp)
movsd  0x50(%rdi),%xmm4
movsd  %xmm4,0x60(%rsp)
movsd  0x58(%rdi),%xmm7
movsd  %xmm7,0x38(%rsp)
movsd  0x60(%rdi),%xmm2
movsd  %xmm2,0x18(%rsp)
movsd  0x68(%rdi),%xmm2
movsd  %xmm2,0x50(%rsp)
movsd  0x70(%rdi),%xmm5
movsd  %xmm5,0x48(%rsp)
movsd  0x78(%rdi),%xmm8
movsd  %xmm8,0x30(%rsp)
movsd  %xmm8,(%rsp)
call   1030 <Det3x3@plt>
movsd  %xmm0,0x10(%rsp)
movsd  0x30(%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x28(%rsp),%xmm0
movsd  0x20(%rsp),%xmm1
movsd  0x18(%rsp),%xmm2
movsd  0x70(%rsp),%xmm3
movsd  0x60(%rsp),%xmm4
movsd  0x48(%rsp),%xmm5
movsd  0x40(%rsp),%xmm6
movsd  0x38(%rsp),%xmm7
call   1030 <Det3x3@plt>
mulsd  0x78(%rsp),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  0x80(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  0x30(%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x28(%rsp),%xmm0
movsd  0x20(%rsp),%xmm1
movsd  0x18(%rsp),%xmm2
movsd  0x68(%rsp),%xmm3
movsd  0x58(%rsp),%xmm4
movsd  0x50(%rsp),%xmm5
movsd  0x40(%rsp),%xmm6
movsd  0x38(%rsp),%xmm7
call   1030 <Det3x3@plt>
mulsd  0x88(%rsp),%xmm0
addsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x48(%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x28(%rsp),%xmm0
movsd  0x20(%rsp),%xmm1
movsd  0x18(%rsp),%xmm2
movsd  0x68(%rsp),%xmm3
movsd  0x58(%rsp),%xmm4
movsd  0x50(%rsp),%xmm5
movsd  0x70(%rsp),%xmm6
movsd  0x60(%rsp),%xmm7
call   1030 <Det3x3@plt>
mulsd  0x90(%rsp),%xmm0
movsd  0x10(%rsp),%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
add    $0x98,%rsp
ret

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

LLVM IR (-O2):
define double @Determinant44([4 x double]* nocapture %0) {
%2 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 0
%3 = load double, double* %2
%4 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 1
%5 = load double, double* %4
%6 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 2
%7 = load double, double* %6
%8 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 3
%9 = load double, double* %8
%10 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 0
%11 = load double, double* %10
%12 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 1
%13 = load double, double* %12
%14 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 2
%15 = load double, double* %14
%16 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 3
%17 = load double, double* %16
%18 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 0
%19 = load double, double* %18
%20 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 1
%21 = load double, double* %20
%22 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 2
%23 = load double, double* %22
%24 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 3
%25 = load double, double* %24
%26 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 0
%27 = load double, double* %26
%28 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 1
%29 = load double, double* %28
%30 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 2
%31 = load double, double* %30
%32 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 3, i64 3
%33 = load double, double* %32
%34 = tail call double @Det3x3(double %13, double %21, double %29, double %15, double %23, double %31, double %17, double %25, double %33)
%35 = tail call double @Det3x3(double %11, double %19, double %27, double %15, double %23, double %31, double %17, double %25, double %33)
%36 = fneg double %5
%37 = fmul double %35, %36
%38 = tail call double @llvm.fmuladd.f64(double %3, double %34, double %37)
%39 = tail call double @Det3x3(double %11, double %19, double %27, double %13, double %21, double %29, double %17, double %25, double %33)
%40 = tail call double @llvm.fmuladd.f64(double %7, double %39, double %38)
%41 = tail call double @Det3x3(double %11, double %19, double %27, double %13, double %21, double %29, double %15, double %23, double %31)
%42 = fneg double %9
%43 = tail call double @llvm.fmuladd.f64(double %42, double %41, double %40)
ret double %43
}