#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 GaussPivot(int, int, double *, double *);
int /*<<< orphan*/ assert(int, char *);
double fabs(double);
void GaussJordan(int n, double *A, double *x, double *b) {
  int i, j, k;
  double pivot, w, _eps = 1E-16;
  for (j = 0; j < n - 1; j++) {
    pivot = GaussPivot(j, n, A, b);
    assert(fabs(pivot) > _eps, "Matriz singular.");
    for (i = j + 1; i < n; i++) {
      w = A[i * n + j] / pivot;
      for (k = j + 1; k < n; k++) {
        A[i * n + k] -= w * A[j * n + k];
      }
      b[i] -= w * b[j];
    }
  }
  for (i = n - 1; i >= 0; i--) {
    assert(fabs(pivot) > _eps, "Matriz singular.");
    x[i] = b[i];
    for (j = i + 1; j < n; j++) {
      x[i] -= A[i * n + j] * x[j];
    }
    x[i] /= A[i * n + i];
  }
  return;
}


x86 assembly (clang, -O0):
<GaussJordan>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  0xec1(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12ae <GaussJordan+0x18e>
mov    -0x28(%rbp),%edi
mov    -0x4(%rbp),%esi
mov    -0x10(%rbp),%rdx
mov    -0x20(%rbp),%rcx
call   1030 <GaussPivot@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movaps 0xe90(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x48(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%edi
lea    0xe88(%rip),%rsi
call   1040 <assert@plt>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    129b <GaussJordan+0x17b>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    124a <GaussJordan+0x12a>
movsd  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11db <GaussJordan+0xbb>
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a6 <GaussJordan+0x86>
jmp    12a0 <GaussJordan+0x180>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    114b <GaussJordan+0x2b>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     13b5 <GaussJordan+0x295>
movsd  -0x38(%rbp),%xmm0
movaps 0xd43(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x48(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%edi
lea    0xd3b(%rip),%rsi
call   1040 <assert@plt>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    137b <GaussJordan+0x25b>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    130d <GaussJordan+0x1ed>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    12b7 <GaussJordan+0x197>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<GaussJordan>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xc8,%rsp
mov    %rcx,0x18(%rsp)
mov    %rdx,%r13
mov    %rsi,%r15
mov    %edi,%eax
mov    %edi,%r12d
cmp    $0x2,%edi
jl     14a0 <GaussJordan+0x380>
mov    %r13,0x40(%rsp)
lea    -0x1(%rax),%ecx
mov    %rcx,0x58(%rsp)
lea    0x1(%r12),%rcx
mov    %rcx,0x98(%rsp)
lea    0x0(,%r12,8),%rcx
mov    %rcx,0x8(%rsp)
lea    0x8(,%r12,8),%rcx
mov    %rcx,0x38(%rsp)
lea    (%r12,%r12,1),%rcx
mov    %rcx,0x50(%rsp)
lea    0x8(%r15),%rcx
mov    %rcx,0x48(%rsp)
mov    %r12,%rcx
neg    %rcx
mov    %rcx,0x10(%rsp)
lea    -0x1(%r12),%rcx
mov    %rcx,0x30(%rsp)
lea    0x18(%r15),%r14
lea    0x18(%r15,%r12,8),%rcx
mov    %rcx,0x28(%rsp)
lea    0x8(%r15,%r12,8),%rcx
mov    %rcx,0x20(%rsp)
mov    $0x1,%ecx
mov    %rcx,(%rsp)
xor    %ebp,%ebp
mov    %r15,%rbx
mov    %rax,0x60(%rsp)
jmp    121a <GaussJordan+0xfa>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addq   $0x1,(%rsp)
addq   $0xffffffffffffffff,0x30(%rsp)
mov    0x38(%rsp),%rax
add    %rax,%r14
add    %rax,0x28(%rsp)
mov    0x8(%rsp),%rax
add    %rax,0x20(%rsp)
add    %rax,%rbx
mov    0x68(%rsp),%rbp
cmp    0x58(%rsp),%rbp
mov    0x60(%rsp),%rax
je     149b <GaussJordan+0x37b>
mov    %ebp,%edi
mov    %eax,%esi
mov    %r15,%rdx
mov    0x18(%rsp),%rcx
call   1030 <GaussPivot@plt>
movapd %xmm0,0xb0(%rsp)
movapd %xmm0,%xmm1
andpd  0xdc0(%rip),%xmm1
xor    %edi,%edi
ucomisd 0xdd6(%rip),%xmm1
seta   %dil
lea    0xdd3(%rip),%rsi
call   1040 <assert@plt>
movapd 0xb0(%rsp),%xmm7
movapd 0xda5(%rip),%xmm6
lea    0x1(%rbp),%rax
mov    %rax,0x68(%rsp)
cmp    %r12,%rax
jae    11e0 <GaussJordan+0xc0>
mov    0x30(%rsp),%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rbp,%rdx
not    %rdx
add    %r12,%rdx
mov    0x38(%rsp),%rsi
imul   %rbp,%rsi
mov    %rbp,%rdi
imul   %r12,%rdi
mov    0x50(%rsp),%rax
add    %rdi,%rax
mov    %rax,0x88(%rsp)
mov    0x48(%rsp),%rax
mov    %rsi,0xa8(%rsp)
add    %rsi,%rax
mov    %rax,0x80(%rsp)
mov    %rdi,0x90(%rsp)
lea    (%rdi,%r12,1),%rax
lea    (%r15,%rax,8),%rax
mov    %rax,0xa0(%rsp)
mov    %rdx,%rax
and    $0xfffffffffffffffc,%rax
mov    (%rsp),%r10
mov    %rax,0x78(%rsp)
add    %r10,%rax
mov    %rax,0x70(%rsp)
mov    0x20(%rsp),%r8
mov    0x28(%rsp),%r11
xor    %r9d,%r9d
jmp    1342 <GaussJordan+0x222>
cs nopw 0x0(%rax,%rax,1)
mov    0x18(%rsp),%rax
mulsd  (%rax,%rbp,8),%xmm0
addsd  (%rax,%r10,8),%xmm0
movsd  %xmm0,(%rax,%r10,8)
add    $0x1,%r10
add    $0x1,%r9
mov    0x8(%rsp),%rax
add    %rax,%r11
add    %rax,%r8
cmp    %r12,%r10
je     11e0 <GaussJordan+0xc0>
mov    %r10,%rax
imul   %r12,%rax
lea    (%rax,%rbp,1),%rsi
movsd  (%r15,%rsi,8),%xmm0
xorpd  %xmm6,%xmm0
divsd  %xmm7,%xmm0
mov    (%rsp),%rsi
cmp    $0x4,%rdx
jb     1410 <GaussJordan+0x2f0>
mov    %r9,%rsi
imul   %r12,%rsi
mov    0x98(%rsp),%rdi
add    %rsi,%rdi
lea    (%r15,%rdi,8),%rdi
add    0xa8(%rsp),%rdi
cmp    0xa0(%rsp),%rdi
jae    13ab <GaussJordan+0x28b>
add    0x88(%rsp),%rsi
lea    (%r15,%rsi,8),%rdi
mov    (%rsp),%rsi
cmp    %rdi,0x80(%rsp)
jb     1410 <GaussJordan+0x2f0>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movupd -0x10(%r14,%rsi,8),%xmm2
movupd (%r14,%rsi,8),%xmm3
movupd -0x10(%r11,%rsi,8),%xmm4
movupd (%r11,%rsi,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,-0x10(%r11,%rsi,8)
movupd %xmm3,(%r11,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%rcx
jne    13c0 <GaussJordan+0x2a0>
mov    0x70(%rsp),%rsi
cmp    0x78(%rsp),%rdx
je     1310 <GaussJordan+0x1f0>
mov    %r12d,%r13d
sub    %esi,%r13d
mov    %rsi,%rdi
test   $0x1,%r13b
je     1447 <GaussJordan+0x327>
mov    0x90(%rsp),%rdi
add    %rsi,%rdi
add    %rsi,%rax
movsd  (%r15,%rdi,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%r15,%rax,8),%xmm1
movsd  %xmm1,(%r15,%rax,8)
lea    0x1(%rsi),%rdi
not    %rsi
cmp    0x10(%rsp),%rsi
je     1310 <GaussJordan+0x1f0>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx,%rdi,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  -0x8(%r8,%rdi,8),%xmm1
movsd  %xmm1,-0x8(%r8,%rdi,8)
movsd  0x8(%rbx,%rdi,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%r8,%rdi,8),%xmm1
movsd  %xmm1,(%r8,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%r12
jne    1460 <GaussJordan+0x340>
jmp    1310 <GaussJordan+0x1f0>
mov    0x40(%rsp),%r13
test   %eax,%eax
jle    15cd <GaussJordan+0x4ad>
andpd  0xb50(%rip),%xmm7
xor    %eax,%eax
ucomisd 0xb66(%rip),%xmm7
seta   %al
mov    %eax,0x8(%rsp)
lea    -0x1(%r12),%rax
imul   %r12,%rax
lea    (%r15,%rax,8),%rbx
add    $0x8,%rbx
lea    0x0(,%r12,8),%rax
neg    %rax
mov    %rax,0x10(%rsp)
xor    %esi,%esi
mov    %r12,%r14
jmp    1515 <GaussJordan+0x3f5>
nopl   0x0(%rax)
mov    (%rsp),%esi
add    %r14,%rax
divsd  (%r15,%rax,8),%xmm0
movsd  %xmm0,0x0(%r13,%r14,8)
add    $0x1,%esi
add    0x10(%rsp),%rbx
cmp    $0x1,%rbp
jle    15cd <GaussJordan+0x4ad>
mov    %esi,(%rsp)
mov    %r14,%rbp
add    $0xffffffffffffffff,%r14
mov    0x8(%rsp),%edi
lea    0xafe(%rip),%rsi
call   1040 <assert@plt>
mov    0x18(%rsp),%rax
movsd  -0x8(%rax,%rbp,8),%xmm0
movsd  %xmm0,-0x8(%r13,%rbp,8)
mov    %r14,%rax
imul   %r12,%rax
cmp    %r12,%rbp
jge    14f0 <GaussJordan+0x3d0>
mov    (%rsp),%esi
test   $0x1,%sil
jne    1560 <GaussJordan+0x440>
mov    %rbp,%rcx
cmp    $0x1,%esi
jne    1590 <GaussJordan+0x470>
jmp    14f3 <GaussJordan+0x3d3>
lea    (%rax,%rbp,1),%rcx
movsd  (%r15,%rcx,8),%xmm1
mulsd  0x0(%r13,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%r13,%r14,8)
lea    0x1(%rbp),%rcx
cmp    $0x1,%esi
je     14f3 <GaussJordan+0x3d3>
nopl   0x0(%rax)
movsd  -0x8(%rbx,%rcx,8),%xmm1
mulsd  0x0(%r13,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%r13,%r14,8)
movsd  (%rbx,%rcx,8),%xmm1
mulsd  0x8(%r13,%rcx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%r13,%r14,8)
add    $0x2,%rcx
cmp    %ecx,%r12d
jne    1590 <GaussJordan+0x470>
jmp    14f3 <GaussJordan+0x3d3>
add    $0xc8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [17 x i8] c"Matriz singular.\00"
define void @GaussJordan(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
%13 = alloca double
%14 = alloca double
store i32 %0, i32* %5
store double* %1, double** %6
store double* %2, double** %7
store double* %3, double** %8
store double 0x3C9CD2B297D889BC, double* %14
store i32 0, i32* %10
br label %15
15: ; preds = %99, %4
%16 = load i32, i32* %10
%17 = load i32, i32* %5
%18 = sub nsw i32 %17, 1
%19 = icmp slt i32 %16, %18
br i1 %19, label %20, label %102
20: ; preds = %15
%21 = load i32, i32* %10
%22 = load i32, i32* %5
%23 = load double*, double** %6
%24 = load double*, double** %8
%25 = call double @GaussPivot(i32 %21, i32 %22, double* %23, double* %24)
store double %25, double* %12
%26 = load double, double* %12
%27 = call double @llvm.fabs.f64(double %26)
%28 = load double, double* %14
%29 = fcmp ogt double %27, %28
%30 = zext i1 %29 to i32
%31 = call i32 @assert(i32 %30, i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str, i64 0, i64 0))
%32 = load i32, i32* %10
%33 = add nsw i32 %32, 1
store i32 %33, i32* %9
br label %34
34: ; preds = %95, %20
%35 = load i32, i32* %9
%36 = load i32, i32* %5
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %98
38: ; preds = %34
%39 = load double*, double** %6
%40 = load i32, i32* %9
%41 = load i32, i32* %5
%42 = mul nsw i32 %40, %41
%43 = load i32, i32* %10
%44 = add nsw i32 %42, %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %39, i64 %45
%47 = load double, double* %46
%48 = load double, double* %12
%49 = fdiv double %47, %48
store double %49, double* %13
%50 = load i32, i32* %10
%51 = add nsw i32 %50, 1
store i32 %51, i32* %11
br label %52
52: ; preds = %78, %38
%53 = load i32, i32* %11
%54 = load i32, i32* %5
%55 = icmp slt i32 %53, %54
br i1 %55, label %56, label %81
56: ; preds = %52
%57 = load double, double* %13
%58 = load double*, double** %6
%59 = load i32, i32* %10
%60 = load i32, i32* %5
%61 = mul nsw i32 %59, %60
%62 = load i32, i32* %11
%63 = add nsw i32 %61, %62
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %58, i64 %64
%66 = load double, double* %65
%67 = load double*, double** %6
%68 = load i32, i32* %9
%69 = load i32, i32* %5
%70 = mul nsw i32 %68, %69
%71 = load i32, i32* %11
%72 = add nsw i32 %70, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %67, i64 %73
%75 = load double, double* %74
%76 = fneg double %57
%77 = call double @llvm.fmuladd.f64(double %76, double %66, double %75)
store double %77, double* %74
br label %78
78: ; preds = %56
%79 = load i32, i32* %11
%80 = add nsw i32 %79, 1
store i32 %80, i32* %11
br label %52
81: ; preds = %52
%82 = load double, double* %13
%83 = load double*, double** %8
%84 = load i32, i32* %10
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = load double*, double** %8
%89 = load i32, i32* %9
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %88, i64 %90
%92 = load double, double* %91
%93 = fneg double %82
%94 = call double @llvm.fmuladd.f64(double %93, double %87, double %92)
store double %94, double* %91
br label %95
95: ; preds = %81
%96 = load i32, i32* %9
%97 = add nsw i32 %96, 1
store i32 %97, i32* %9
br label %34
98: ; preds = %34
br label %99
99: ; preds = %98
%100 = load i32, i32* %10
%101 = add nsw i32 %100, 1
store i32 %101, i32* %10
br label %15
102: ; preds = %15
%103 = load i32, i32* %5
%104 = sub nsw i32 %103, 1
store i32 %104, i32* %9
br label %105
105: ; preds = %171, %102
%106 = load i32, i32* %9
%107 = icmp sge i32 %106, 0
br i1 %107, label %108, label %174
108: ; preds = %105
%109 = load double, double* %12
%110 = call double @llvm.fabs.f64(double %109)
%111 = load double, double* %14
%112 = fcmp ogt double %110, %111
%113 = zext i1 %112 to i32
%114 = call i32 @assert(i32 %113, i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str, i64 0, i64 0))
%115 = load double*, double** %8
%116 = load i32, i32* %9
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double, double* %115, i64 %117
%119 = load double, double* %118
%120 = load double*, double** %7
%121 = load i32, i32* %9
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %120, i64 %122
store double %119, double* %123
%124 = load i32, i32* %9
%125 = add nsw i32 %124, 1
store i32 %125, i32* %10
br label %126
126: ; preds = %152, %108
%127 = load i32, i32* %10
%128 = load i32, i32* %5
%129 = icmp slt i32 %127, %128
br i1 %129, label %130, label %155
130: ; preds = %126
%131 = load double*, double** %6
%132 = load i32, i32* %9
%133 = load i32, i32* %5
%134 = mul nsw i32 %132, %133
%135 = load i32, i32* %10
%136 = add nsw i32 %134, %135
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %131, i64 %137
%139 = load double, double* %138
%140 = load double*, double** %7
%141 = load i32, i32* %10
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %140, i64 %142
%144 = load double, double* %143
%145 = load double*, double** %7
%146 = load i32, i32* %9
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = fneg double %139
%151 = call double @llvm.fmuladd.f64(double %150, double %144, double %149)
store double %151, double* %148
br label %152
152: ; preds = %130
%153 = load i32, i32* %10
%154 = add nsw i32 %153, 1
store i32 %154, i32* %10
br label %126
155: ; preds = %126
%156 = load double*, double** %6
%157 = load i32, i32* %9
%158 = load i32, i32* %5
%159 = mul nsw i32 %157, %158
%160 = load i32, i32* %9
%161 = add nsw i32 %159, %160
%162 = sext i32 %161 to i64
%163 = getelementptr inbounds double, double* %156, i64 %162
%164 = load double, double* %163
%165 = load double*, double** %7
%166 = load i32, i32* %9
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds double, double* %165, i64 %167
%169 = load double, double* %168
%170 = fdiv double %169, %164
store double %170, double* %168
br label %171
171: ; preds = %155
%172 = load i32, i32* %9
%173 = add nsw i32 %172, -1
store i32 %173, i32* %9
br label %105
174: ; preds = %105
ret void
}

LLVM IR (-O2):
@.str = [17 x i8] c"Matriz singular.\00"
define void @GaussJordan(i32 %0, double* %1, double* nocapture %2, double* %3) {
%5 = icmp sgt i32 %0, 1
br i1 %5, label %6, label %21
6: ; preds = %4
%7 = add nsw i32 %0, -1
%8 = zext i32 %0 to i64
%9 = zext i32 %7 to i64
%10 = zext i32 %0 to i64
%11 = add nuw nsw i64 %10, 1
%12 = shl nuw nsw i64 %10, 3
%13 = add nuw nsw i64 %12, 8
%14 = shl nuw nsw i64 %10, 1
%15 = getelementptr double, double* %1, i64 1
%16 = bitcast double* %15 to i8*
%17 = sub nsw i64 0, %8
br label %29
18: ; preds = %143, %29
%19 = add nuw nsw i64 %31, 1
%20 = icmp eq i64 %47, %9
br i1 %20, label %21, label %29
21: ; preds = %18, %4
%22 = phi double [ undef, %4 ], [ %42, %18 ]
%23 = icmp sgt i32 %0, 0
br i1 %23, label %24, label %209
24: ; preds = %21
%25 = tail call double @llvm.fabs.f64(double %22)
%26 = fcmp ogt double %25, 0x3C9CD2B297D889BC
%27 = zext i1 %26 to i32
%28 = zext i32 %0 to i64
br label %151
29: ; preds = %6, %18
%30 = phi i64 [ 0, %6 ], [ %47, %18 ]
%31 = phi i64 [ 1, %6 ], [ %19, %18 ]
%32 = xor i64 %30, -1
%33 = add nsw i64 %32, %10
%34 = mul i64 %13, %30
%35 = mul i64 %30, %10
%36 = add i64 %14, %35
%37 = getelementptr i8, i8* %16, i64 %34
%38 = add i64 %35, %10
%39 = getelementptr double, double* %1, i64 %38
%40 = bitcast double* %39 to i8*
%41 = trunc i64 %30 to i32
%42 = tail call double @GaussPivot(i32 %41, i32 %0, double* %1, double* %3)
%43 = tail call double @llvm.fabs.f64(double %42)
%44 = fcmp ogt double %43, 0x3C9CD2B297D889BC
%45 = zext i1 %44 to i32
%46 = tail call i32 @assert(i32 %45, i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str, i64 0, i64 0))
%47 = add nuw nsw i64 %30, 1
%48 = icmp ult i64 %47, %8
br i1 %48, label %49, label %18
49: ; preds = %29
%50 = mul nsw i64 %30, %8
%51 = getelementptr inbounds double, double* %3, i64 %30
%52 = icmp ult i64 %33, 4
%53 = and i64 %33, -4
%54 = add i64 %31, %53
%55 = icmp eq i64 %33, %53
br label %56
56: ; preds = %143, %49
%57 = phi i64 [ %150, %143 ], [ 0, %49 ]
%58 = phi i64 [ %148, %143 ], [ %31, %49 ]
%59 = mul i64 %57, %10
%60 = mul nsw i64 %58, %8
%61 = add nuw nsw i64 %60, %30
%62 = getelementptr inbounds double, double* %1, i64 %61
%63 = load double, double* %62
%64 = fneg double %63
%65 = fdiv double %64, %42
br i1 %52, label %106, label %66
66: ; preds = %56
%67 = add i64 %36, %59
%68 = getelementptr double, double* %1, i64 %67
%69 = bitcast double* %68 to i8*
%70 = add i64 %11, %59
%71 = getelementptr double, double* %1, i64 %70
%72 = bitcast double* %71 to i8*
%73 = getelementptr i8, i8* %72, i64 %34
%74 = icmp ult i8* %73, %40
%75 = icmp ult i8* %37, %69
%76 = and i1 %74, %75
br i1 %76, label %106, label %77
77: ; preds = %66
%78 = insertelement <2 x double> poison, double %65, i64 0
%79 = shufflevector <2 x double> %78, <2 x double> poison, <2 x i32> zeroinitializer
%80 = insertelement <2 x double> poison, double %65, i64 0
%81 = shufflevector <2 x double> %80, <2 x double> poison, <2 x i32> zeroinitializer
br label %82
82: ; preds = %82, %77
%83 = phi i64 [ 0, %77 ], [ %103, %82 ]
%84 = add i64 %31, %83
%85 = add nuw nsw i64 %84, %50
%86 = getelementptr inbounds double, double* %1, i64 %85
%87 = bitcast double* %86 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = getelementptr inbounds double, double* %86, i64 2
%90 = bitcast double* %89 to <2 x double>*
%91 = load <2 x double>, <2 x double>* %90
%92 = add nuw nsw i64 %84, %60
%93 = getelementptr inbounds double, double* %1, i64 %92
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = getelementptr inbounds double, double* %93, i64 2
%97 = bitcast double* %96 to <2 x double>*
%98 = load <2 x double>, <2 x double>* %97
%99 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %79, <2 x double> %88, <2 x double> %95)
%100 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %81, <2 x double> %91, <2 x double> %98)
%101 = bitcast double* %93 to <2 x double>*
store <2 x double> %99, <2 x double>* %101
%102 = bitcast double* %96 to <2 x double>*
store <2 x double> %100, <2 x double>* %102
%103 = add nuw i64 %83, 4
%104 = icmp eq i64 %103, %53
br i1 %104, label %105, label %82
105: ; preds = %82
br i1 %55, label %143, label %106
106: ; preds = %66, %56, %105
%107 = phi i64 [ %31, %66 ], [ %31, %56 ], [ %54, %105 ]
%108 = sub i64 %8, %107
%109 = xor i64 %107, -1
%110 = and i64 %108, 1
%111 = icmp eq i64 %110, 0
br i1 %111, label %121, label %112
112: ; preds = %106
%113 = add nuw nsw i64 %107, %50
%114 = getelementptr inbounds double, double* %1, i64 %113
%115 = load double, double* %114
%116 = add nuw nsw i64 %107, %60
%117 = getelementptr inbounds double, double* %1, i64 %116
%118 = load double, double* %117
%119 = tail call double @llvm.fmuladd.f64(double %65, double %115, double %118)
store double %119, double* %117
%120 = add nuw nsw i64 %107, 1
br label %121
121: ; preds = %112, %106
%122 = phi i64 [ %107, %106 ], [ %120, %112 ]
%123 = icmp eq i64 %109, %17
br i1 %123, label %143, label %124
124: ; preds = %121, %124
%125 = phi i64 [ %141, %124 ], [ %122, %121 ]
%126 = add nuw nsw i64 %125, %50
%127 = getelementptr inbounds double, double* %1, i64 %126
%128 = load double, double* %127
%129 = add nuw nsw i64 %125, %60
%130 = getelementptr inbounds double, double* %1, i64 %129
%131 = load double, double* %130
%132 = tail call double @llvm.fmuladd.f64(double %65, double %128, double %131)
store double %132, double* %130
%133 = add nuw nsw i64 %125, 1
%134 = add nuw nsw i64 %133, %50
%135 = getelementptr inbounds double, double* %1, i64 %134
%136 = load double, double* %135
%137 = add nuw nsw i64 %133, %60
%138 = getelementptr inbounds double, double* %1, i64 %137
%139 = load double, double* %138
%140 = tail call double @llvm.fmuladd.f64(double %65, double %136, double %139)
store double %140, double* %138
%141 = add nuw nsw i64 %125, 2
%142 = icmp eq i64 %141, %10
br i1 %142, label %143, label %124
143: ; preds = %121, %124, %105
%144 = load double, double* %51
%145 = getelementptr inbounds double, double* %3, i64 %58
%146 = load double, double* %145
%147 = tail call double @llvm.fmuladd.f64(double %65, double %144, double %146)
store double %147, double* %145
%148 = add nuw nsw i64 %58, 1
%149 = icmp eq i64 %148, %10
%150 = add i64 %57, 1
br i1 %149, label %18, label %56
151: ; preds = %24, %200
%152 = phi i32 [ 0, %24 ], [ %208, %200 ]
%153 = phi i64 [ %28, %24 ], [ %155, %200 ]
%154 = phi i64 [ %28, %24 ], [ %207, %200 ]
%155 = add nsw i64 %153, -1
%156 = tail call i32 @assert(i32 %27, i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str, i64 0, i64 0))
%157 = getelementptr inbounds double, double* %3, i64 %155
%158 = load double, double* %157
%159 = getelementptr inbounds double, double* %2, i64 %155
store double %158, double* %159
%160 = icmp slt i64 %153, %28
%161 = mul nsw i64 %155, %28
br i1 %160, label %162, label %200
162: ; preds = %151
%163 = and i32 %152, 1
%164 = icmp eq i32 %163, 0
br i1 %164, label %174, label %165
165: ; preds = %162
%166 = add nsw i64 %154, %161
%167 = getelementptr inbounds double, double* %1, i64 %166
%168 = load double, double* %167
%169 = getelementptr inbounds double, double* %2, i64 %154
%170 = load double, double* %169
%171 = fneg double %168
%172 = tail call double @llvm.fmuladd.f64(double %171, double %170, double %158)
store double %172, double* %159
%173 = add nsw i64 %154, 1
br label %174
174: ; preds = %165, %162
%175 = phi double [ undef, %162 ], [ %172, %165 ]
%176 = phi i64 [ %154, %162 ], [ %173, %165 ]
%177 = phi double [ %158, %162 ], [ %172, %165 ]
%178 = icmp eq i32 %152, 1
br i1 %178, label %200, label %179
179: ; preds = %174, %179
%180 = phi i64 [ %197, %179 ], [ %176, %174 ]
%181 = phi double [ %196, %179 ], [ %177, %174 ]
%182 = add nsw i64 %180, %161
%183 = getelementptr inbounds double, double* %1, i64 %182
%184 = load double, double* %183
%185 = getelementptr inbounds double, double* %2, i64 %180
%186 = load double, double* %185
%187 = fneg double %184
%188 = tail call double @llvm.fmuladd.f64(double %187, double %186, double %181)
store double %188, double* %159
%189 = add nsw i64 %180, 1
%190 = add nsw i64 %189, %161
%191 = getelementptr inbounds double, double* %1, i64 %190
%192 = load double, double* %191
%193 = getelementptr inbounds double, double* %2, i64 %189
%194 = load double, double* %193
%195 = fneg double %192
%196 = tail call double @llvm.fmuladd.f64(double %195, double %194, double %188)
store double %196, double* %159
%197 = add nsw i64 %180, 2
%198 = trunc i64 %197 to i32
%199 = icmp eq i32 %198, %0
br i1 %199, label %200, label %179
200: ; preds = %174, %179, %151
%201 = phi double [ %158, %151 ], [ %175, %174 ], [ %196, %179 ]
%202 = add nsw i64 %161, %155
%203 = getelementptr inbounds double, double* %1, i64 %202
%204 = load double, double* %203
%205 = fdiv double %201, %204
store double %205, double* %159
%206 = icmp sgt i64 %153, 1
%207 = add nsw i64 %154, -1
%208 = add i32 %152, 1
br i1 %206, label %151, label %209
209: ; preds = %200, %21
ret void
}