#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 */
int /*<<< orphan*/ ludcmp(double **, int, double *);
void det(double **a, int n, double *d) {
  int j;
  ludcmp(a, n, d);
  for (j = 0; j < n; j++)
    *d *= a[j][j];
}


x86 assembly (clang, -O0):
<det>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x18(%rbp),%rdx
call   1030 <ludcmp@plt>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1175 <det+0x65>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    113a <det+0x2a>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<det>:
push   %rbp
push   %r14
push   %rbx
mov    %rdx,%r14
mov    %esi,%ebp
mov    %rdi,%rbx
call   1030 <ludcmp@plt>
test   %ebp,%ebp
jle    11b8 <det+0xa8>
mov    %ebp,%edx
movsd  (%r14),%xmm0
lea    -0x1(%rdx),%rcx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1143 <det+0x33>
xor    %ecx,%ecx
jmp    1197 <det+0x87>
and    $0xfffffffc,%edx
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
mov    (%rbx,%rcx,8),%rsi
mulsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%r14)
mov    0x8(%rbx,%rcx,8),%rsi
mulsd  0x8(%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%r14)
mov    0x10(%rbx,%rcx,8),%rsi
mulsd  0x10(%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%r14)
mov    0x18(%rbx,%rcx,8),%rsi
mulsd  0x18(%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%r14)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1150 <det+0x40>
test   %rax,%rax
je     11b8 <det+0xa8>
nopl   0x0(%rax)
mov    (%rbx,%rcx,8),%rdx
mulsd  (%rdx,%rcx,8),%xmm0
movsd  %xmm0,(%r14)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rax
jne    11a0 <det+0x90>
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define void @det(double** %0, i32 %1, double* %2) {
%4 = alloca double**
%5 = alloca i32
%6 = alloca double*
%7 = alloca i32
store double** %0, double*** %4
store i32 %1, i32* %5
store double* %2, double** %6
%8 = load double**, double*** %4
%9 = load i32, i32* %5
%10 = load double*, double** %6
%11 = call i32 @ludcmp(double** %8, i32 %9, double* %10)
store i32 0, i32* %7
br label %12
12: ; preds = %29, %3
%13 = load i32, i32* %7
%14 = load i32, i32* %5
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %32
16: ; preds = %12
%17 = load double**, double*** %4
%18 = load i32, i32* %7
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double*, double** %17, i64 %19
%21 = load double*, double** %20
%22 = load i32, i32* %7
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %21, i64 %23
%25 = load double, double* %24
%26 = load double*, double** %6
%27 = load double, double* %26
%28 = fmul double %27, %25
store double %28, double* %26
br label %29
29: ; preds = %16
%30 = load i32, i32* %7
%31 = add nsw i32 %30, 1
store i32 %31, i32* %7
br label %12
32: ; preds = %12
ret void
}

LLVM IR (-O2):
define void @det(double** %0, i32 %1, double* %2) {
%4 = tail call i32 @ludcmp(double** %0, i32 %1, double* %2)
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %60
6: ; preds = %3
%7 = zext i32 %1 to i64
%8 = load double, double* %2
%9 = add nsw i64 %7, -1
%10 = and i64 %7, 3
%11 = icmp ult i64 %9, 3
br i1 %11, label %44, label %12
12: ; preds = %6
%13 = and i64 %7, 4294967292
br label %14
14: ; preds = %14, %12
%15 = phi double [ %8, %12 ], [ %40, %14 ]
%16 = phi i64 [ 0, %12 ], [ %41, %14 ]
%17 = phi i64 [ 0, %12 ], [ %42, %14 ]
%18 = getelementptr inbounds double*, double** %0, i64 %16
%19 = load double*, double** %18
%20 = getelementptr inbounds double, double* %19, i64 %16
%21 = load double, double* %20
%22 = fmul double %21, %15
store double %22, double* %2
%23 = or i64 %16, 1
%24 = getelementptr inbounds double*, double** %0, i64 %23
%25 = load double*, double** %24
%26 = getelementptr inbounds double, double* %25, i64 %23
%27 = load double, double* %26
%28 = fmul double %27, %22
store double %28, double* %2
%29 = or i64 %16, 2
%30 = getelementptr inbounds double*, double** %0, i64 %29
%31 = load double*, double** %30
%32 = getelementptr inbounds double, double* %31, i64 %29
%33 = load double, double* %32
%34 = fmul double %33, %28
store double %34, double* %2
%35 = or i64 %16, 3
%36 = getelementptr inbounds double*, double** %0, i64 %35
%37 = load double*, double** %36
%38 = getelementptr inbounds double, double* %37, i64 %35
%39 = load double, double* %38
%40 = fmul double %39, %34
store double %40, double* %2
%41 = add nuw nsw i64 %16, 4
%42 = add i64 %17, 4
%43 = icmp eq i64 %42, %13
br i1 %43, label %44, label %14
44: ; preds = %14, %6
%45 = phi double [ %8, %6 ], [ %40, %14 ]
%46 = phi i64 [ 0, %6 ], [ %41, %14 ]
%47 = icmp eq i64 %10, 0
br i1 %47, label %60, label %48
48: ; preds = %44, %48
%49 = phi double [ %56, %48 ], [ %45, %44 ]
%50 = phi i64 [ %57, %48 ], [ %46, %44 ]
%51 = phi i64 [ %58, %48 ], [ 0, %44 ]
%52 = getelementptr inbounds double*, double** %0, i64 %50
%53 = load double*, double** %52
%54 = getelementptr inbounds double, double* %53, i64 %50
%55 = load double, double* %54
%56 = fmul double %55, %49
store double %56, double* %2
%57 = add nuw nsw i64 %50, 1
%58 = add i64 %51, 1
%59 = icmp eq i64 %58, %10
br i1 %59, label %60, label %48
60: ; preds = %44, %48, %3
ret void
}