#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__ fabs(double);
int /*<<< orphan*/ printf(char *);
int /*<<< orphan*/ swap_rows(double *, size_t, size_t, size_t const);
void gaussian_elimination(double *a, const size_t rows, const size_t cols) {
  size_t row = 0;
  for (size_t col = 0; col < cols - 1; ++col) {
    size_t pivot = row;
    for (size_t i = row + 1; i < rows; ++i) {
      if (fabs(a[i * cols + col]) > fabs(a[pivot * cols + col])) {
        pivot = i;
      }
    }
    if (a[pivot * cols + col] == 0) {
      printf("The matrix is singular.\n");
      continue;
    }
    if (col != pivot) {
      swap_rows(a, col, pivot, cols);
    }
    for (size_t i = row + 1; i < rows; ++i) {
      double scale = a[i * cols + col] / a[row * cols + col];
      for (size_t j = col + 1; j < cols; ++j) {
        a[i * cols + j] -= a[row * cols + j] * scale;
      }
      a[i * cols + col] = 0;
    }
    row++;
  }
}


x86 assembly (clang, -O0):
<gaussian_elimination>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x18(%rbp),%rcx
sub    $0x1,%rcx
cmp    %rcx,%rax
jae    1362 <gaussian_elimination+0x232>
mov    -0x20(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    11f3 <gaussian_elimination+0xc3>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,-0x58(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,%rcx
mov    -0x58(%rbp),%rax
cmp    %rcx,%rax
jle    11dd <gaussian_elimination+0xad>
mov    -0x38(%rbp),%rax
mov    %rax,-0x30(%rbp)
jmp    11e2 <gaussian_elimination+0xb2>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    117d <gaussian_elimination+0x4d>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    122d <gaussian_elimination+0xfd>
jp     122d <gaussian_elimination+0xfd>
lea    0xddd(%rip),%rdi
call   1030 <printf@plt>
jmp    1351 <gaussian_elimination+0x221>
mov    -0x28(%rbp),%rax
cmp    -0x30(%rbp),%rax
je     1250 <gaussian_elimination+0x120>
mov    -0x8(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    -0x30(%rbp),%rdx
mov    -0x18(%rbp),%rcx
call   1050 <swap_rows@plt>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    1345 <gaussian_elimination+0x215>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x28(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x18(%rbp),%rax
jae    131b <gaussian_elimination+0x1eb>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x50(%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    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    12a7 <gaussian_elimination+0x177>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
imul   -0x18(%rbp),%rcx
add    -0x28(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    125c <gaussian_elimination+0x12c>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1154 <gaussian_elimination+0x24>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gaussian_elimination>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %rsi,0x8(%rsp)
mov    %rdx,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,0x20(%rsp)
je     1555 <gaussian_elimination+0x425>
mov    %rdi,%r12
mov    %rdx,%rdi
mov    0x8(%rsp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,0x58(%rsp)
mov    %rdx,%rax
shl    $0x4,%rax
mov    %rax,0x80(%rsp)
lea    0x18(%r12),%rax
mov    %rax,0x48(%rsp)
lea    0x0(,%rdx,8),%rax
mov    %rax,0x78(%rsp)
xor    %r11d,%r11d
mov    0x20(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    %r12,0x50(%rsp)
xor    %r14d,%r14d
mov    %rdx,(%rsp)
jmp    1207 <gaussian_elimination+0xd7>
nopl   (%rax)
mov    %r11,%rax
imul   (%rsp),%rax
add    %r14,%rax
movsd  (%r12,%rax,8),%xmm0
ucomisd 0xe37(%rip),%xmm0
jne    1360 <gaussian_elimination+0x230>
jp     1360 <gaussian_elimination+0x230>
lea    0xe2c(%rip),%rdi
call   1030 <puts@plt>
mov    0x10(%rsp),%r11
add    $0x1,%r14
addq   $0x8,0x50(%rsp)
addq   $0xffffffffffffffff,0x40(%rsp)
addq   $0x8,0x48(%rsp)
cmp    0x20(%rsp),%r14
je     1555 <gaussian_elimination+0x425>
lea    0x1(%r11),%rax
mov    %r11,0x10(%rsp)
mov    %rax,0x28(%rsp)
cmp    0x8(%rsp),%rax
jae    11b0 <gaussian_elimination+0x80>
mov    0x10(%rsp),%r13
mov    %r13d,%eax
not    %eax
mov    0x8(%rsp),%rcx
add    %ecx,%eax
mov    0x28(%rsp),%r15
mov    %r15,%rbp
mov    %r13,%rbx
test   $0x1,%al
je     127e <gaussian_elimination+0x14e>
mov    %r15,%rax
mov    (%rsp),%rbx
imul   %rbx,%rax
add    %r14,%rax
movsd  (%r12,%rax,8),%xmm0
call   1040 <fabs@plt>
mov    %rax,%rbp
mov    %r13,%rax
imul   %rbx,%rax
add    %r14,%rax
movsd  (%r12,%rax,8),%xmm0
call   1040 <fabs@plt>
cmp    %rax,%rbp
mov    %r13,%r11
cmovg  %r15,%r11
lea    0x2(%r13),%rbp
mov    %r11,%rbx
cmp    %r13,0x58(%rsp)
je     11b0 <gaussian_elimination+0x80>
lea    0x1(%rbp),%rax
mov    (%rsp),%rcx
imul   %rcx,%rax
mov    0x50(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
mov    %rax,0x38(%rsp)
mov    %rcx,%rax
imul   %rbp,%rax
lea    (%rdx,%rax,8),%rax
mov    %rax,0x30(%rsp)
xor    %r13d,%r13d
jmp    12dd <gaussian_elimination+0x1ad>
nopl   0x0(%rax,%rax,1)
mov    %rbx,%r11
add    0x80(%rsp),%r13
add    $0x2,%rbp
mov    %r11,%rbx
cmp    0x8(%rsp),%rbp
je     11b0 <gaussian_elimination+0x80>
mov    0x30(%rsp),%rax
movsd  (%rax,%r13,1),%xmm0
call   1040 <fabs@plt>
mov    %rax,0x18(%rsp)
mov    %rbx,%rax
mov    (%rsp),%r15
imul   %r15,%rax
add    %r14,%rax
movsd  (%r12,%rax,8),%xmm0
call   1040 <fabs@plt>
mov    %rbp,%rcx
cmp    %rax,0x18(%rsp)
jg     1318 <gaussian_elimination+0x1e8>
mov    %rbx,%rcx
mov    0x38(%rsp),%rax
movsd  (%rax,%r13,1),%xmm0
mov    %rcx,%rbx
call   1040 <fabs@plt>
mov    %rax,0x18(%rsp)
mov    %rbx,%rax
imul   %r15,%rax
add    %r14,%rax
movsd  (%r12,%rax,8),%xmm0
call   1040 <fabs@plt>
cmp    %rax,0x18(%rsp)
jle    12c0 <gaussian_elimination+0x190>
lea    0x1(%rbp),%r11
jmp    12c3 <gaussian_elimination+0x193>
nopl   0x0(%rax)
cmp    %r11,%r14
je     1377 <gaussian_elimination+0x247>
mov    %r12,%rdi
mov    %r14,%rsi
mov    %r11,%rdx
mov    (%rsp),%rcx
call   1050 <swap_rows@plt>
mov    0x28(%rsp),%r11
cmp    0x8(%rsp),%r11
mov    0x10(%rsp),%rbp
jae    11e6 <gaussian_elimination+0xb6>
mov    0x20(%rsp),%r15
sub    %r14,%r15
mov    0x40(%rsp),%rcx
and    $0xfffffffffffffffc,%rcx
mov    %r14,%rax
add    $0x1,%r14
lea    (%r12,%r15,8),%r8
lea    0x10(,%rbp,8),%rsi
mov    (%rsp),%rdi
imul   %rdi,%rbp
lea    (%rax,%rbp,1),%rbx
mov    %rbx,0x18(%rsp)
mov    %r11,%rbx
imul   %rdi,%rbx
lea    (%r14,%rbx,1),%rdx
mov    %rdx,0x38(%rsp)
imul   %rdi,%rsi
mov    %rsi,0x10(%rsp)
lea    (%r14,%rbp,1),%rsi
lea    (%r12,%rsi,8),%rdx
mov    %rdx,0x70(%rsp)
lea    (%r8,%rsi,8),%rdx
mov    %rdx,0x30(%rsp)
mov    %r15,%rsi
and    $0xfffffffffffffffc,%rsi
mov    %rbp,%rdx
mov    %rax,%r13
mov    %rsi,0x68(%rsp)
add    %rsi,%rax
add    $0x1,%rax
mov    %rax,0x60(%rsp)
mov    0x48(%rsp),%rax
lea    (%rax,%rbp,8),%rsi
lea    (%rax,%rbx,8),%rbp
lea    (%r12,%rdx,8),%r10
lea    (%r12,%rbx,8),%r8
xor    %r9d,%r9d
jmp    1456 <gaussian_elimination+0x326>
nopl   0x0(%rax,%rax,1)
movq   $0x0,(%r12,%rbx,8)
add    $0x1,%r11
add    $0x1,%r9
mov    0x78(%rsp),%rax
add    %rax,%rbp
add    %rax,%r8
cmp    0x8(%rsp),%r11
je     1548 <gaussian_elimination+0x418>
mov    %r11,%rbx
imul   %rdi,%rbx
add    %r13,%rbx
cmp    %rdi,%r14
jae    1430 <gaussian_elimination+0x300>
movsd  (%r12,%rbx,8),%xmm0
mov    0x18(%rsp),%rax
divsd  (%r12,%rax,8),%xmm0
mov    %r14,%rdx
cmp    $0x4,%r15
jb     1520 <gaussian_elimination+0x3f0>
mov    %r9,%rdx
imul   %rdi,%rdx
mov    %r12,%rdi
mov    0x38(%rsp),%rax
lea    (%rax,%rdx,1),%r12
lea    (%rdi,%r12,8),%rax
mov    %rdi,%r12
mov    (%rsp),%rdi
cmp    0x30(%rsp),%rax
jae    14bb <gaussian_elimination+0x38b>
lea    (%r12,%rdx,8),%rax
add    0x10(%rsp),%rax
mov    %r14,%rdx
cmp    0x70(%rsp),%rax
ja     1520 <gaussian_elimination+0x3f0>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movupd -0x10(%rsi,%rdx,8),%xmm2
movupd (%rsi,%rdx,8),%xmm3
movupd -0x10(%rbp,%rdx,8),%xmm4
movupd 0x0(%rbp,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
subpd  %xmm2,%xmm4
mulpd  %xmm1,%xmm3
subpd  %xmm3,%xmm5
movupd %xmm4,-0x10(%rbp,%rdx,8)
movupd %xmm5,0x0(%rbp,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    14d0 <gaussian_elimination+0x3a0>
mov    0x60(%rsp),%rdx
cmp    0x68(%rsp),%r15
je     1430 <gaussian_elimination+0x300>
nopl   0x0(%rax)
movsd  (%r8,%rdx,8),%xmm1
movsd  (%r10,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,(%r8,%rdx,8)
add    $0x1,%rdx
cmp    %rdi,%rdx
jb     1520 <gaussian_elimination+0x3f0>
jmp    1430 <gaussian_elimination+0x300>
mov    0x28(%rsp),%r11
mov    %r13,%r14
jmp    11e6 <gaussian_elimination+0xb6>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [25 x i8] c"The matrix is singular.\0A\00"
define void @gaussian_elimination(double* %0, i64 %1, i64 %2) {
%4 = alloca double*
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca double
%13 = alloca i64
store double* %0, double** %4
store i64 %1, i64* %5
store i64 %2, i64* %6
store i64 0, i64* %7
store i64 0, i64* %8
br label %14
14: ; preds = %143, %3
%15 = load i64, i64* %8
%16 = load i64, i64* %6
%17 = sub i64 %16, 1
%18 = icmp ult i64 %15, %17
br i1 %18, label %19, label %146
19: ; preds = %14
%20 = load i64, i64* %7
store i64 %20, i64* %9
%21 = load i64, i64* %7
%22 = add i64 %21, 1
store i64 %22, i64* %10
br label %23
23: ; preds = %50, %19
%24 = load i64, i64* %10
%25 = load i64, i64* %5
%26 = icmp ult i64 %24, %25
br i1 %26, label %27, label %53
27: ; preds = %23
%28 = load double*, double** %4
%29 = load i64, i64* %10
%30 = load i64, i64* %6
%31 = mul i64 %29, %30
%32 = load i64, i64* %8
%33 = add i64 %31, %32
%34 = getelementptr inbounds double, double* %28, i64 %33
%35 = load double, double* %34
%36 = call i64 @fabs(double %35)
%37 = load double*, double** %4
%38 = load i64, i64* %9
%39 = load i64, i64* %6
%40 = mul i64 %38, %39
%41 = load i64, i64* %8
%42 = add i64 %40, %41
%43 = getelementptr inbounds double, double* %37, i64 %42
%44 = load double, double* %43
%45 = call i64 @fabs(double %44)
%46 = icmp sgt i64 %36, %45
br i1 %46, label %47, label %49
47: ; preds = %27
%48 = load i64, i64* %10
store i64 %48, i64* %9
br label %49
49: ; preds = %47, %27
br label %50
50: ; preds = %49
%51 = load i64, i64* %10
%52 = add i64 %51, 1
store i64 %52, i64* %10
br label %23
53: ; preds = %23
%54 = load double*, double** %4
%55 = load i64, i64* %9
%56 = load i64, i64* %6
%57 = mul i64 %55, %56
%58 = load i64, i64* %8
%59 = add i64 %57, %58
%60 = getelementptr inbounds double, double* %54, i64 %59
%61 = load double, double* %60
%62 = fcmp oeq double %61, 0.000000e+00
br i1 %62, label %63, label %65
63: ; preds = %53
%64 = call i32 @printf(i8* getelementptr inbounds ([25 x i8], [25 x i8]* @.str, i64 0, i64 0))
br label %143
65: ; preds = %53
%66 = load i64, i64* %8
%67 = load i64, i64* %9
%68 = icmp ne i64 %66, %67
br i1 %68, label %69, label %75
69: ; preds = %65
%70 = load double*, double** %4
%71 = load i64, i64* %8
%72 = load i64, i64* %9
%73 = load i64, i64* %6
%74 = call i32 @swap_rows(double* %70, i64 %71, i64 %72, i64 %73)
br label %75
75: ; preds = %69, %65
%76 = load i64, i64* %7
%77 = add i64 %76, 1
store i64 %77, i64* %11
br label %78
78: ; preds = %137, %75
%79 = load i64, i64* %11
%80 = load i64, i64* %5
%81 = icmp ult i64 %79, %80
br i1 %81, label %82, label %140
82: ; preds = %78
%83 = load double*, double** %4
%84 = load i64, i64* %11
%85 = load i64, i64* %6
%86 = mul i64 %84, %85
%87 = load i64, i64* %8
%88 = add i64 %86, %87
%89 = getelementptr inbounds double, double* %83, i64 %88
%90 = load double, double* %89
%91 = load double*, double** %4
%92 = load i64, i64* %7
%93 = load i64, i64* %6
%94 = mul i64 %92, %93
%95 = load i64, i64* %8
%96 = add i64 %94, %95
%97 = getelementptr inbounds double, double* %91, i64 %96
%98 = load double, double* %97
%99 = fdiv double %90, %98
store double %99, double* %12
%100 = load i64, i64* %8
%101 = add i64 %100, 1
store i64 %101, i64* %13
br label %102
102: ; preds = %126, %82
%103 = load i64, i64* %13
%104 = load i64, i64* %6
%105 = icmp ult i64 %103, %104
br i1 %105, label %106, label %129
106: ; preds = %102
%107 = load double*, double** %4
%108 = load i64, i64* %7
%109 = load i64, i64* %6
%110 = mul i64 %108, %109
%111 = load i64, i64* %13
%112 = add i64 %110, %111
%113 = getelementptr inbounds double, double* %107, i64 %112
%114 = load double, double* %113
%115 = load double, double* %12
%116 = load double*, double** %4
%117 = load i64, i64* %11
%118 = load i64, i64* %6
%119 = mul i64 %117, %118
%120 = load i64, i64* %13
%121 = add i64 %119, %120
%122 = getelementptr inbounds double, double* %116, i64 %121
%123 = load double, double* %122
%124 = fneg double %114
%125 = call double @llvm.fmuladd.f64(double %124, double %115, double %123)
store double %125, double* %122
br label %126
126: ; preds = %106
%127 = load i64, i64* %13
%128 = add i64 %127, 1
store i64 %128, i64* %13
br label %102
129: ; preds = %102
%130 = load double*, double** %4
%131 = load i64, i64* %11
%132 = load i64, i64* %6
%133 = mul i64 %131, %132
%134 = load i64, i64* %8
%135 = add i64 %133, %134
%136 = getelementptr inbounds double, double* %130, i64 %135
store double 0.000000e+00, double* %136
br label %137
137: ; preds = %129
%138 = load i64, i64* %11
%139 = add i64 %138, 1
store i64 %139, i64* %11
br label %78
140: ; preds = %78
%141 = load i64, i64* %7
%142 = add i64 %141, 1
store i64 %142, i64* %7
br label %143
143: ; preds = %140, %63
%144 = load i64, i64* %8
%145 = add i64 %144, 1
store i64 %145, i64* %8
br label %14
146: ; preds = %14
ret void
}

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