#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 */
scalar_t__ malloc(int);
int /*<<< orphan*/ printf(char *, ...);
double *gauss(int rank, double augmentedmatrix[rank][rank + 1],
              double accuracy) {
  int i, j, k = 1;
  double temp;
  do {
    i = k + 1;
    do {
      temp = augmentedmatrix[i - 1][k - 1] / augmentedmatrix[k - 1][k - 1];
      j = k;
      do
        augmentedmatrix[i - 1][j - 1] -= augmentedmatrix[k - 1][j - 1] * temp;
      while (++j <= rank + 1);
    } while (++i <= rank);
  } while (++k <= rank - 1);
  double *solutionvector = (double *)malloc(rank * sizeof(double));
  solutionvector[rank - 1] =
      augmentedmatrix[rank - 1][rank] / augmentedmatrix[rank - 1][rank - 1];
  i = rank - 1;
  do {
    j = i + 1;
    temp = augmentedmatrix[i - 1][rank];
    do
      temp -= augmentedmatrix[i - 1][j - 1] * solutionvector[j - 1];
    while (++j <= rank);
    solutionvector[i - 1] = temp / augmentedmatrix[i - 1][i - 1];
  } while (--i >= 1);
  i = 1;
  printf("GAUSS METHOD:\n\n");
  do {
    printf("%+09g\t", *(solutionvector + i - 1));
  } while (++i <= rank);
  return solutionvector;
}


x86 assembly (clang, -O0):
<gauss>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    -0x4(%rbp),%eax
add    $0x1,%eax
mov    %eax,%eax
mov    %rax,-0x40(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%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    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jle    11ab <gauss+0x8b>
jmp    1240 <gauss+0x120>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmp    -0x4(%rbp),%eax
jle    1150 <gauss+0x30>
jmp    1257 <gauss+0x137>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    1147 <gauss+0x27>
movslq -0x4(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    -0x40(%rbp),%rdx
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x40(%rbp),%rdx
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
cmp    -0x4(%rbp),%eax
jle    131a <gauss+0x1fa>
mov    -0x40(%rbp),%rdx
movsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
cmp    $0x1,%eax
jge    12e7 <gauss+0x1c7>
movl   $0x1,-0x1c(%rbp)
lea    0xc0e(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    -0x38(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  -0x8(%rax,%rcx,8),%xmm0
lea    0xc02(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmp    -0x4(%rbp),%eax
jle    13f9 <gauss+0x2d9>
mov    -0x38(%rbp),%rax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gauss>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %rsi,%rbp
lea    0x1(%rdi),%r15d
lea    -0x1(%rdi),%eax
mov    %rdi,0x18(%rsp)
movslq %edi,%rbx
cmp    $0x2,%eax
mov    $0x1,%ecx
mov    %eax,(%rsp)
cmovge %eax,%ecx
add    $0x1,%ecx
mov    %rcx,0x30(%rsp)
lea    0x1(%rbx),%rax
mov    %rax,0x28(%rsp)
lea    0x8(%rsi),%rax
mov    %rax,0x20(%rsp)
lea    0x0(,%r15,8),%r13
lea    0x8(,%r15,8),%r9
lea    0x18(%rsi),%rcx
mov    $0x1,%r8d
mov    $0x2,%edx
xor    %eax,%eax
mov    %rcx,%rdi
mov    %rbx,%r14
mov    %r15,0x68(%rsp)
mov    %rsi,0x60(%rsp)
mov    %r9,0x8(%rsp)
jmp    11fc <gauss+0xcc>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r8
mov    0x4(%rsp),%edx
add    $0x1,%edx
mov    0x40(%rsp),%rax
add    $0x1,%rax
mov    0x38(%rsp),%r14
add    $0xffffffffffffffff,%r14
mov    0x8(%rsp),%r9
add    %r9,%rdi
mov    0x48(%rsp),%rcx
add    $0x8,%rcx
add    %r13,%rsi
cmp    0x30(%rsp),%r8
je     1428 <gauss+0x2f8>
mov    %edx,0x4(%rsp)
movslq %edx,%r10
lea    -0x1(%r10),%rdx
imul   %r15,%rdx
lea    0x0(,%rdx,8),%r9
add    %rbp,%r9
mov    %rcx,0x48(%rsp)
lea    (%rcx,%rdx,8),%r11
mov    %r14,0x38(%rsp)
and    $0xfffffffffffffffc,%r14
mov    %rbx,%r12
sub    %rax,%r12
lea    (%r8,%rdx,1),%rcx
mov    %rcx,0x80(%rsp)
mov    0x28(%rsp),%rcx
add    %rcx,%rdx
mov    %rdx,0x88(%rsp)
mov    0x8(%rsp),%rdx
imul   %rax,%rdx
add    0x20(%rsp),%rdx
mov    %rdx,0x78(%rsp)
mov    %rax,0x40(%rsp)
imul   %r15,%rax
add    %rcx,%rax
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
mov    %rax,0x70(%rsp)
lea    -0x1(%r8),%rax
imul   %r15,%rax
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
mov    %rax,0x90(%rsp)
mov    %r12,0x10(%rsp)
and    $0xfffffffffffffffc,%r12
mov    %r12,0x58(%rsp)
lea    (%r8,%r12,1),%rax
mov    %rax,0x50(%rsp)
xor    %ecx,%ecx
jmp    12ca <gauss+0x19a>
nop
lea    0x1(%r10),%rax
add    $0x1,%rcx
add    %r13,%r11
add    %r13,%r9
cmp    %rbx,%r10
mov    %rax,%r10
jge    11c0 <gauss+0x90>
lea    -0x1(%r10),%rax
imul   %r15,%rax
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
movsd  -0x8(%rax,%r8,8),%xmm1
mov    0x90(%rsp),%rdx
movsd  -0x8(%rdx,%r8,8),%xmm2
movapd %xmm1,%xmm0
divsd  %xmm2,%xmm0
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,-0x8(%rax,%r8,8)
cmp    %rbx,%r8
jg     12b0 <gauss+0x180>
mov    %r8,%rdx
cmpq   $0x4,0x10(%rsp)
jae    1350 <gauss+0x220>
add    $0xffffffffffffffff,%rdx
xchg   %ax,%ax
movsd  0x8(%r9,%rdx,8),%xmm1
movsd  0x8(%rsi,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,0x8(%r9,%rdx,8)
add    $0x1,%rdx
cmp    %rbx,%rdx
jl     1320 <gauss+0x1f0>
jmp    12b0 <gauss+0x180>
nopw   0x0(%rax,%rax,1)
mov    %rcx,%rax
imul   %r15,%rax
mov    0x80(%rsp),%rdx
add    %rax,%rdx
add    0x88(%rsp),%rax
lea    0x0(,%rdx,8),%r12
add    %rbp,%r12
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
cmp    0x70(%rsp),%r12
setb   %bpl
cmp    %rax,0x78(%rsp)
mov    %r13,%r15
setb   %r13b
mov    %r8,%rdx
test   %r13b,%bpl
mov    %r15,%r13
mov    0x68(%rsp),%r15
mov    0x60(%rsp),%rbp
jne    131a <gauss+0x1ea>
mov    %r8,%rdx
cmp    %rax,%r12
jb     131a <gauss+0x1ea>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd -0x10(%rdi,%rdx,8),%xmm2
movupd (%rdi,%rdx,8),%xmm3
movupd -0x10(%r11,%rdx,8),%xmm4
movupd (%r11,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
subpd  %xmm2,%xmm4
mulpd  %xmm1,%xmm3
subpd  %xmm3,%xmm5
movupd %xmm4,-0x10(%r11,%rdx,8)
movupd %xmm5,(%r11,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%r14
jne    13d0 <gauss+0x2a0>
mov    0x50(%rsp),%rdx
mov    0x58(%rsp),%rax
cmp    %rax,0x10(%rsp)
je     12b0 <gauss+0x180>
jmp    131a <gauss+0x1ea>
lea    0x0(,%rbx,8),%edi
call   1050 <malloc@plt>
mov    %r13,%r9
mov    %rax,%r13
movslq (%rsp),%rsi
mov    %rsi,%rax
imul   %r15,%rax
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
movsd  (%rax,%rbx,8),%xmm0
divsd  (%rax,%rsi,8),%xmm0
movsd  %xmm0,0x0(%r13,%rsi,8)
lea    -0x1(%rsi),%rax
imul   %r15,%rax
add    %rsi,%rax
lea    0x0(,%rax,8),%rcx
add    %rbp,%rcx
xor    $0xfffffffffffffff8,%r9
nopl   0x0(%rax,%rax,1)
lea    -0x1(%rsi),%rdx
mov    %rdx,%rax
imul   %r15,%rax
lea    0x0(,%rax,8),%r8
add    %rbp,%r8
movsd  (%r8,%rbx,8),%xmm0
mov    %rcx,%rdi
mov    %rsi,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdi),%xmm1
mulsd  0x0(%r13,%rax,8),%xmm1
add    $0x1,%rax
subsd  %xmm1,%xmm0
add    $0x8,%rdi
cmp    %rbx,%rax
jl     14b0 <gauss+0x380>
divsd  (%r8,%rdx,8),%xmm0
movsd  %xmm0,0x0(%r13,%rdx,8)
add    %r9,%rcx
cmp    $0x1,%rsi
mov    %rdx,%rsi
jg     1480 <gauss+0x350>
lea    0xb1b(%rip),%rdi
call   1030 <puts@plt>
mov    0x18(%rsp),%rax
cmp    $0x2,%eax
mov    $0x1,%ebx
cmovge %eax,%ebx
xor    %ebp,%ebp
lea    0xaf6(%rip),%r14
nopw   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rbp,8),%xmm0
mov    %r14,%rdi
mov    $0x1,%al
call   1040 <printf@plt>
add    $0x1,%rbp
cmp    %rbp,%rbx
jne    1510 <gauss+0x3e0>
mov    %r13,%rax
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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