#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 *crear_vector(int);
int /*<<< orphan*/ eliminate_mod(double **, double *, int, double *, double,
                                 int *);
double fabs(double);
int /*<<< orphan*/ sustitucion(double **, double *, int, double *);
void gauss_mod(double **a, double *b, int n, double *x, double tol) {
  int i, j, er;
  double *s;
  s = crear_vector(n);
  er = 0;
  for (i = 0; i < n; i++) {
    s[i] = fabs(a[i][0]);
    for (j = 1; j < n; j++) {
      if (fabs(a[i][j]) > s[i]) {
        s[i] = fabs(a[i][j]);
      }
    }
  }
  eliminate_mod(a, s, n, b, tol, &er);
  if (er != -1)
    sustitucion(a, x, n, b);
}


x86 assembly (clang, -O0):
<gauss_mod>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%edi
call   1040 <crear_vector@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1233 <gauss_mod+0x103>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movsd  (%rax),%xmm0
movaps 0xe77(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1220 <gauss_mod+0xf0>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe37(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    120d <gauss_mod+0xdd>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe04(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1212 <gauss_mod+0xe2>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11a1 <gauss_mod+0x71>
jmp    1225 <gauss_mod+0xf5>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1166 <gauss_mod+0x36>
mov    -0x8(%rbp),%rdi
mov    -0x40(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x10(%rbp),%rcx
movsd  -0x28(%rbp),%xmm0
lea    -0x34(%rbp),%r8
call   1030 <eliminate_mod@plt>
cmpl   $0xffffffff,-0x34(%rbp)
je     126e <gauss_mod+0x13e>
mov    -0x8(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x14(%rbp),%edx
mov    -0x10(%rbp),%rcx
call   1050 <sustitucion@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gauss_mod>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
movsd  %xmm0,0x8(%rsp)
mov    %rcx,%r14
mov    %edx,%ebp
mov    %rsi,%r15
mov    %rdi,%r12
mov    %edx,%edi
call   1040 <crear_vector@plt>
movl   $0x0,0x4(%rsp)
test   %ebp,%ebp
jle    1235 <gauss_mod+0x105>
mov    %ebp,%r9d
lea    -0x1(%r9),%r8
mov    %r8,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
xor    %edi,%edi
movapd 0xe81(%rip),%xmm0
jmp    119d <gauss_mod+0x6d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %r9,%rdi
je     1235 <gauss_mod+0x105>
mov    (%r12,%rdi,8),%rdx
movsd  (%rdx),%xmm1
andpd  %xmm0,%xmm1
movlpd %xmm1,(%rax,%rdi,8)
cmp    $0x1,%ebp
je     1190 <gauss_mod+0x60>
mov    $0x1,%ebx
cmp    $0x2,%ebp
jne    11e0 <gauss_mod+0xb0>
test   $0x1,%r8b
je     1190 <gauss_mod+0x60>
movsd  (%rdx,%rbx,8),%xmm2
andpd  %xmm0,%xmm2
ucomisd %xmm1,%xmm2
jbe    1190 <gauss_mod+0x60>
movsd  %xmm2,(%rax,%rdi,8)
jmp    1190 <gauss_mod+0x60>
nopl   0x0(%rax)
mov    $0x1,%ebx
jmp    1202 <gauss_mod+0xd2>
nopw   0x0(%rax,%rax,1)
lea    (%rsi,%rbx,1),%rcx
add    $0x2,%rcx
add    $0x2,%rbx
cmp    $0x1,%rcx
je     11bd <gauss_mod+0x8d>
movsd  (%rdx,%rbx,8),%xmm2
andpd  %xmm0,%xmm2
ucomisd %xmm1,%xmm2
jbe    121a <gauss_mod+0xea>
movsd  %xmm2,(%rax,%rdi,8)
movapd %xmm2,%xmm1
movsd  0x8(%rdx,%rbx,8),%xmm2
andpd  %xmm0,%xmm2
ucomisd %xmm1,%xmm2
jbe    11f0 <gauss_mod+0xc0>
movsd  %xmm2,(%rax,%rdi,8)
movapd %xmm2,%xmm1
jmp    11f0 <gauss_mod+0xc0>
lea    0x4(%rsp),%r8
mov    %r12,%rdi
mov    %rax,%rsi
mov    %ebp,%edx
mov    %r15,%rcx
movsd  0x8(%rsp),%xmm0
call   1030 <eliminate_mod@plt>
cmpl   $0xffffffff,0x4(%rsp)
je     1267 <gauss_mod+0x137>
mov    %r12,%rdi
mov    %r14,%rsi
mov    %ebp,%edx
mov    %r15,%rcx
call   1050 <sustitucion@plt>
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

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