#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*/ gauss(int, double *, double *, int, int);
int /*<<< orphan*/ perm(int, double *, double *, int, double, double *);
void reduction(int n, double *L, double *D, double *Z) {
  int i, j, k;
  double del;
  j = n - 2;
  k = n - 2;
  while (j >= 0) {
    if (j <= k)
      for (i = j + 1; i < n; i++)
        gauss(n, L, Z, i, j);
    del = D[j] + L[j + 1 + j * n] * L[j + 1 + j * n] * D[j + 1];
    if (del + 1E-6 < D[j + 1]) {
      perm(n, L, D, j, del, Z);
      k = j;
      j = n - 2;
    } else
      j--;
  }
}


x86 assembly (clang, -O0):
<reduction>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x28(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x28(%rbp)
jl     1267 <reduction+0x147>
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     119e <reduction+0x7e>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1199 <reduction+0x79>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%r8d
call   1030 <gauss@plt>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1168 <reduction+0x48>
jmp    119e <reduction+0x7e>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xdf8(%rip),%xmm1
addsd  -0x38(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1259 <reduction+0x139>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x28(%rbp),%ecx
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%r8
call   1040 <perm@plt>
mov    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x28(%rbp)
jmp    1262 <reduction+0x142>
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
jmp    1149 <reduction+0x29>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<reduction>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
cmp    $0x2,%edi
jl     121c <reduction+0xfc>
mov    %rcx,%r15
mov    %rdx,%rbp
mov    %rsi,%r12
mov    %edi,%ebx
lea    -0x2(%rbx),%eax
mov    %eax,%ecx
mov    %eax,0x8(%rsp)
mov    %eax,%r13d
mov    %rdx,0x10(%rsp)
jmp    117c <reduction+0x5c>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0xffffffff,%r13d
mov    %r13d,%eax
mov    0xc(%rsp),%ecx
mov    %ecx,%r13d
mov    %r13d,%ecx
mov    %eax,%r13d
test   %eax,%eax
js     121c <reduction+0xfc>
lea    0x1(%r13),%r14d
mov    %ecx,0xc(%rsp)
cmp    %ecx,%r13d
jg     11b9 <reduction+0x99>
cmp    %ebx,%r14d
jge    11b9 <reduction+0x99>
mov    %r14d,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %ebx,%edi
mov    %r12,%rsi
mov    %r15,%rdx
mov    %ebp,%ecx
mov    %r13d,%r8d
call   1030 <gauss@plt>
add    $0x1,%ebp
cmp    %ebp,%ebx
jne    11a0 <reduction+0x80>
movslq %r13d,%rax
mov    %r13d,%ecx
imul   %ebx,%ecx
add    %r14d,%ecx
movslq %ecx,%rcx
movsd  (%r12,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
movslq %r14d,%rcx
mov    0x10(%rsp),%rbp
movsd  0x0(%rbp,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
addsd  0x0(%rbp,%rax,8),%xmm0
movapd %xmm0,%xmm2
addsd  0xe0a(%rip),%xmm2
ucomisd %xmm2,%xmm1
jbe    1160 <reduction+0x40>
mov    %ebx,%edi
mov    %r12,%rsi
mov    %rbp,%rdx
mov    %r13d,%ecx
mov    %r15,%r8
call   1040 <perm@plt>
mov    0x8(%rsp),%eax
jmp    116e <reduction+0x4e>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @reduction(i32 %0, double* %1, double* %2, double* %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
store i32 %0, i32* %5
store double* %1, double** %6
store double* %2, double** %7
store double* %3, double** %8
%13 = load i32, i32* %5
%14 = sub nsw i32 %13, 2
store i32 %14, i32* %10
%15 = load i32, i32* %5
%16 = sub nsw i32 %15, 2
store i32 %16, i32* %11
br label %17
17: ; preds = %99, %4
%18 = load i32, i32* %10
%19 = icmp sge i32 %18, 0
br i1 %19, label %20, label %100
20: ; preds = %17
%21 = load i32, i32* %10
%22 = load i32, i32* %11
%23 = icmp sle i32 %21, %22
br i1 %23, label %24, label %42
24: ; preds = %20
%25 = load i32, i32* %10
%26 = add nsw i32 %25, 1
store i32 %26, i32* %9
br label %27
27: ; preds = %38, %24
%28 = load i32, i32* %9
%29 = load i32, i32* %5
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %41
31: ; preds = %27
%32 = load i32, i32* %5
%33 = load double*, double** %6
%34 = load double*, double** %8
%35 = load i32, i32* %9
%36 = load i32, i32* %10
%37 = call i32 @gauss(i32 %32, double* %33, double* %34, i32 %35, i32 %36)
br label %38
38: ; preds = %31
%39 = load i32, i32* %9
%40 = add nsw i32 %39, 1
store i32 %40, i32* %9
br label %27
41: ; preds = %27
br label %42
42: ; preds = %41, %20
%43 = load double*, double** %7
%44 = load i32, i32* %10
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %43, i64 %45
%47 = load double, double* %46
%48 = load double*, double** %6
%49 = load i32, i32* %10
%50 = add nsw i32 %49, 1
%51 = load i32, i32* %10
%52 = load i32, i32* %5
%53 = mul nsw i32 %51, %52
%54 = add nsw i32 %50, %53
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %48, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %6
%59 = load i32, i32* %10
%60 = add nsw i32 %59, 1
%61 = load i32, i32* %10
%62 = load i32, i32* %5
%63 = mul nsw i32 %61, %62
%64 = add nsw i32 %60, %63
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %58, i64 %65
%67 = load double, double* %66
%68 = fmul double %57, %67
%69 = load double*, double** %7
%70 = load i32, i32* %10
%71 = add nsw i32 %70, 1
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %69, i64 %72
%74 = load double, double* %73
%75 = call double @llvm.fmuladd.f64(double %68, double %74, double %47)
store double %75, double* %12
%76 = load double, double* %12
%77 = fadd double %76, 0x3EB0C6F7A0B5ED8D
%78 = load double*, double** %7
%79 = load i32, i32* %10
%80 = add nsw i32 %79, 1
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %78, i64 %81
%83 = load double, double* %82
%84 = fcmp olt double %77, %83
br i1 %84, label %85, label %96
85: ; preds = %42
%86 = load i32, i32* %5
%87 = load double*, double** %6
%88 = load double*, double** %7
%89 = load i32, i32* %10
%90 = load double, double* %12
%91 = load double*, double** %8
%92 = call i32 @perm(i32 %86, double* %87, double* %88, i32 %89, double %90, double* %91)
%93 = load i32, i32* %10
store i32 %93, i32* %11
%94 = load i32, i32* %5
%95 = sub nsw i32 %94, 2
store i32 %95, i32* %10
br label %99
96: ; preds = %42
%97 = load i32, i32* %10
%98 = add nsw i32 %97, -1
store i32 %98, i32* %10
br label %99
99: ; preds = %96, %85
br label %17
100: ; preds = %17
ret void
}

LLVM IR (-O2):
define void @reduction(i32 %0, double* %1, double* %2, double* %3) {
%5 = add nsw i32 %0, -2
%6 = icmp sgt i32 %0, 1
br i1 %6, label %7, label %43
7: ; preds = %4, %39
%8 = phi i32 [ %41, %39 ], [ %5, %4 ]
%9 = phi i32 [ %40, %39 ], [ %5, %4 ]
%10 = icmp sle i32 %9, %8
%11 = add nuw nsw i32 %9, 1
%12 = icmp slt i32 %11, %0
%13 = select i1 %10, i1 %12, i1 false
br i1 %13, label %14, label %19
14: ; preds = %7, %14
%15 = phi i32 [ %17, %14 ], [ %11, %7 ]
%16 = tail call i32 @gauss(i32 %0, double* %1, double* %3, i32 %15, i32 %9)
%17 = add i32 %15, 1
%18 = icmp eq i32 %17, %0
br i1 %18, label %19, label %14
19: ; preds = %14, %7
%20 = sext i32 %9 to i64
%21 = getelementptr inbounds double, double* %2, i64 %20
%22 = load double, double* %21
%23 = mul nsw i32 %9, %0
%24 = add nsw i32 %11, %23
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double, double* %1, i64 %25
%27 = load double, double* %26
%28 = fmul double %27, %27
%29 = sext i32 %11 to i64
%30 = getelementptr inbounds double, double* %2, i64 %29
%31 = load double, double* %30
%32 = tail call double @llvm.fmuladd.f64(double %28, double %31, double %22)
%33 = fadd double %32, 0x3EB0C6F7A0B5ED8D
%34 = fcmp olt double %33, %31
br i1 %34, label %35, label %37
35: ; preds = %19
%36 = tail call i32 @perm(i32 %0, double* %1, double* %2, i32 %9, double %32, double* %3)
br label %39
37: ; preds = %19
%38 = add nsw i32 %9, -1
br label %39
39: ; preds = %37, %35
%40 = phi i32 [ %5, %35 ], [ %38, %37 ]
%41 = phi i32 [ %9, %35 ], [ %8, %37 ]
%42 = icmp sgt i32 %40, -1
br i1 %42, label %7, label %43
43: ; preds = %39, %4
ret void
}