void reduction(float a[][6], int size, int pivot, int col) {
  int i, j;
  float factor;
  factor = a[pivot][col];
  for (i = 0; i < 2 * size; i++)
    a[pivot][i] /= factor;
  for (i = 0; i < size; i++)
    if (i != pivot) {
      factor = a[i][col];
      for (j = 0; j < 2 * size; j++)
        a[i][j] = a[i][j] - a[pivot][j] * factor;
    }
}


x86 assembly (clang, -O0):
<reduction>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0xc(%rbp),%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jge    117a <reduction+0x7a>
movss  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1135 <reduction+0x35>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1256 <reduction+0x156>
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
je     1243 <reduction+0x143>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0xc(%rbp),%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jge    123e <reduction+0x13e>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x20(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11bd <reduction+0xbd>
jmp    1243 <reduction+0x143>
jmp    1248 <reduction+0x148>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1181 <reduction+0x81>
pop    %rbp
ret

x86 assembly (clang, -O2):
<reduction>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %esi,%esi
jle    1331 <reduction+0x231>
movslq %edx,%rax
movslq %ecx,%r8
lea    (%rax,%rax,2),%r12
lea    (%rdi,%r12,8),%r15
movss  (%r15,%r8,4),%xmm0
lea    (%rsi,%rsi,1),%eax
cmp    $0x2,%eax
mov    $0x1,%r13d
cmovge %eax,%r13d
cmp    $0x4,%r13d
jae    1143 <reduction+0x43>
xor    %eax,%eax
jmp    11c0 <reduction+0xc0>
mov    %r13d,%eax
and    $0x7ffffffc,%eax
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x4(%rax),%rcx
mov    %rcx,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rcx,%rcx
je     133c <reduction+0x23c>
mov    %r9,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%rdi,%r12,8),%rbp
add    $0x10,%rbp
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movups -0x10(%rbp,%rcx,4),%xmm2
movups 0x0(%rbp,%rcx,4),%xmm3
divps  %xmm1,%xmm2
movups %xmm2,-0x10(%rbp,%rcx,4)
divps  %xmm1,%xmm3
movups %xmm3,0x0(%rbp,%rcx,4)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1180 <reduction+0x80>
test   $0x1,%r9b
je     11b7 <reduction+0xb7>
movups (%r15,%rcx,4),%xmm2
divps  %xmm1,%xmm2
movups %xmm2,(%r15,%rcx,4)
cmp    %r13,%rax
je     11d9 <reduction+0xd9>
nopl   0x0(%rax)
movss  (%r15,%rax,4),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,(%r15,%rax,4)
add    $0x1,%rax
cmp    %rax,%r13
jne    11c0 <reduction+0xc0>
test   %esi,%esi
jle    1331 <reduction+0x231>
mov    %edx,%r10d
mov    %esi,%r14d
lea    (%r15,%r13,4),%rax
mov    %rax,-0x8(%rsp)
mov    %r13d,%r9d
and    $0x7ffffff8,%r9d
mov    %r13,%r11
neg    %r11
lea    0x10(%rdi),%rbx
lea    (%rdi,%r12,8),%rsi
add    $0x10,%rsi
lea    0x4(%rdi),%rcx
xor    %edx,%edx
jmp    1235 <reduction+0x135>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdx
add    $0x18,%rbx
add    $0x18,%rcx
cmp    %r14,%rdx
je     1331 <reduction+0x231>
cmp    %r10,%rdx
je     1220 <reduction+0x120>
lea    (%rdx,%rdx,2),%rax
lea    (%rdi,%rax,8),%rax
movss  (%rax,%r8,4),%xmm0
cmp    $0x8,%r13d
jb     1266 <reduction+0x166>
cmp    -0x8(%rsp),%rax
jae    12e0 <reduction+0x1e0>
lea    (%rax,%r13,4),%rbp
cmp    %rbp,%r15
jae    12e0 <reduction+0x1e0>
xor    %r12d,%r12d
mov    %r12,%rbp
test   $0x1,%r13b
je     1293 <reduction+0x193>
movss  (%rax,%r12,4),%xmm1
movss  (%r15,%r12,4),%xmm2
mulss  %xmm0,%xmm2
subss  %xmm2,%xmm1
movss  %xmm1,(%rax,%r12,4)
mov    %r12,%rbp
or     $0x1,%rbp
not    %r12
cmp    %r11,%r12
je     1220 <reduction+0x120>
nopl   0x0(%rax,%rax,1)
movss  -0x10(%rsi,%rbp,4),%xmm1
mulss  %xmm0,%xmm1
movss  -0x4(%rcx,%rbp,4),%xmm2
subss  %xmm1,%xmm2
movss  (%rcx,%rbp,4),%xmm1
movss  %xmm2,-0x4(%rcx,%rbp,4)
movss  -0xc(%rsi,%rbp,4),%xmm2
mulss  %xmm0,%xmm2
subss  %xmm2,%xmm1
movss  %xmm1,(%rcx,%rbp,4)
add    $0x2,%rbp
cmp    %rbp,%r13
jne    12a0 <reduction+0x1a0>
jmp    1220 <reduction+0x120>
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
xor    %ebp,%ebp
nopl   0x0(%rax)
movups -0x10(%rbx,%rbp,4),%xmm2
movups (%rbx,%rbp,4),%xmm3
movups -0x10(%rsi,%rbp,4),%xmm4
movups (%rsi,%rbp,4),%xmm5
mulps  %xmm1,%xmm4
subps  %xmm4,%xmm2
mulps  %xmm1,%xmm5
subps  %xmm5,%xmm3
movups %xmm2,-0x10(%rbx,%rbp,4)
movups %xmm3,(%rbx,%rbp,4)
add    $0x8,%rbp
cmp    %rbp,%r9
jne    12f0 <reduction+0x1f0>
mov    %r9,%r12
cmp    %r13,%r9
je     1220 <reduction+0x120>
jmp    1269 <reduction+0x169>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ecx,%ecx
test   $0x1,%r9b
jne    11aa <reduction+0xaa>
jmp    11b7 <reduction+0xb7>

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

LLVM IR (-O2):
define void @reduction([6 x float]* nocapture %0, i32 %1, i32 %2, i32 %3) {
%5 = sext i32 %2 to i64
%6 = sext i32 %3 to i64
%7 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %6
%8 = load float, float* %7
%9 = shl i32 %1, 1
%10 = icmp sgt i32 %1, 0
br i1 %10, label %11, label %151
11: ; preds = %4
%12 = call i32 @llvm.smax.i32(i32 %9, i32 1)
%13 = zext i32 %12 to i64
%14 = icmp ult i32 %12, 4
br i1 %14, label %54, label %15
15: ; preds = %11
%16 = and i64 %13, 2147483644
%17 = insertelement <4 x float> poison, float %8, i64 0
%18 = shufflevector <4 x float> %17, <4 x float> poison, <4 x i32> zeroinitializer
%19 = add nsw i64 %16, -4
%20 = lshr exact i64 %19, 2
%21 = add nuw nsw i64 %20, 1
%22 = and i64 %21, 1
%23 = icmp eq i64 %19, 0
br i1 %23, label %43, label %24
24: ; preds = %15
%25 = and i64 %21, 9223372036854775806
br label %26
26: ; preds = %26, %24
%27 = phi i64 [ 0, %24 ], [ %40, %26 ]
%28 = phi i64 [ 0, %24 ], [ %41, %26 ]
%29 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %27
%30 = bitcast float* %29 to <4 x float>*
%31 = load <4 x float>, <4 x float>* %30
%32 = fdiv <4 x float> %31, %18
%33 = bitcast float* %29 to <4 x float>*
store <4 x float> %32, <4 x float>* %33
%34 = or i64 %27, 4
%35 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %34
%36 = bitcast float* %35 to <4 x float>*
%37 = load <4 x float>, <4 x float>* %36
%38 = fdiv <4 x float> %37, %18
%39 = bitcast float* %35 to <4 x float>*
store <4 x float> %38, <4 x float>* %39
%40 = add nuw i64 %27, 8
%41 = add i64 %28, 2
%42 = icmp eq i64 %41, %25
br i1 %42, label %43, label %26
43: ; preds = %26, %15
%44 = phi i64 [ 0, %15 ], [ %40, %26 ]
%45 = icmp eq i64 %22, 0
br i1 %45, label %52, label %46
46: ; preds = %43
%47 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %44
%48 = bitcast float* %47 to <4 x float>*
%49 = load <4 x float>, <4 x float>* %48
%50 = fdiv <4 x float> %49, %18
%51 = bitcast float* %47 to <4 x float>*
store <4 x float> %50, <4 x float>* %51
br label %52
52: ; preds = %43, %46
%53 = icmp eq i64 %16, %13
br i1 %53, label %56, label %54
54: ; preds = %11, %52
%55 = phi i64 [ 0, %11 ], [ %16, %52 ]
br label %70
56: ; preds = %70, %52
br i1 %10, label %57, label %151
57: ; preds = %56
%58 = call i32 @llvm.smax.i32(i32 %9, i32 1)
%59 = zext i32 %2 to i64
%60 = zext i32 %1 to i64
%61 = zext i32 %58 to i64
%62 = getelementptr [6 x float], [6 x float]* %0, i64 %5, i64 0
%63 = getelementptr [6 x float], [6 x float]* %0, i64 %5, i64 %13
%64 = icmp ult i32 %12, 8
%65 = and i64 %13, 2147483640
%66 = icmp eq i64 %65, %13
%67 = and i64 %13, 1
%68 = icmp eq i64 %67, 0
%69 = sub nsw i64 0, %13
br label %77
70: ; preds = %54, %70
%71 = phi i64 [ %75, %70 ], [ %55, %54 ]
%72 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %71
%73 = load float, float* %72
%74 = fdiv float %73, %8
store float %74, float* %72
%75 = add nuw nsw i64 %71, 1
%76 = icmp eq i64 %75, %13
br i1 %76, label %56, label %70
77: ; preds = %57, %148
%78 = phi i64 [ 0, %57 ], [ %149, %148 ]
%79 = getelementptr [6 x float], [6 x float]* %0, i64 %78, i64 0
%80 = getelementptr [6 x float], [6 x float]* %0, i64 %78, i64 %13
%81 = icmp eq i64 %78, %59
br i1 %81, label %148, label %82
82: ; preds = %77
%83 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %78, i64 %6
%84 = load float, float* %83
br i1 %64, label %117, label %85
85: ; preds = %82
%86 = icmp ult float* %79, %63
%87 = icmp ult float* %62, %80
%88 = and i1 %86, %87
br i1 %88, label %117, label %89
89: ; preds = %85
%90 = insertelement <4 x float> poison, float %84, i64 0
%91 = shufflevector <4 x float> %90, <4 x float> poison, <4 x i32> zeroinitializer
%92 = insertelement <4 x float> poison, float %84, i64 0
%93 = shufflevector <4 x float> %92, <4 x float> poison, <4 x i32> zeroinitializer
br label %94
94: ; preds = %94, %89
%95 = phi i64 [ 0, %89 ], [ %114, %94 ]
%96 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %78, i64 %95
%97 = bitcast float* %96 to <4 x float>*
%98 = load <4 x float>, <4 x float>* %97
%99 = getelementptr inbounds float, float* %96, i64 4
%100 = bitcast float* %99 to <4 x float>*
%101 = load <4 x float>, <4 x float>* %100
%102 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %95
%103 = bitcast float* %102 to <4 x float>*
%104 = load <4 x float>, <4 x float>* %103
%105 = getelementptr inbounds float, float* %102, i64 4
%106 = bitcast float* %105 to <4 x float>*
%107 = load <4 x float>, <4 x float>* %106
%108 = fneg <4 x float> %104
%109 = fneg <4 x float> %107
%110 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %108, <4 x float> %91, <4 x float> %98)
%111 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %109, <4 x float> %93, <4 x float> %101)
%112 = bitcast float* %96 to <4 x float>*
store <4 x float> %110, <4 x float>* %112
%113 = bitcast float* %99 to <4 x float>*
store <4 x float> %111, <4 x float>* %113
%114 = add nuw i64 %95, 8
%115 = icmp eq i64 %114, %65
br i1 %115, label %116, label %94
116: ; preds = %94
br i1 %66, label %148, label %117
117: ; preds = %85, %82, %116
%118 = phi i64 [ 0, %85 ], [ 0, %82 ], [ %65, %116 ]
%119 = xor i64 %118, -1
br i1 %68, label %128, label %120
120: ; preds = %117
%121 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %78, i64 %118
%122 = load float, float* %121
%123 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %118
%124 = load float, float* %123
%125 = fneg float %124
%126 = tail call float @llvm.fmuladd.f32(float %125, float %84, float %122)
store float %126, float* %121
%127 = or i64 %118, 1
br label %128
128: ; preds = %120, %117
%129 = phi i64 [ %118, %117 ], [ %127, %120 ]
%130 = icmp eq i64 %119, %69
br i1 %130, label %148, label %131
131: ; preds = %128, %131
%132 = phi i64 [ %146, %131 ], [ %129, %128 ]
%133 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %78, i64 %132
%134 = load float, float* %133
%135 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %132
%136 = load float, float* %135
%137 = fneg float %136
%138 = tail call float @llvm.fmuladd.f32(float %137, float %84, float %134)
store float %138, float* %133
%139 = add nuw nsw i64 %132, 1
%140 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %78, i64 %139
%141 = load float, float* %140
%142 = getelementptr inbounds [6 x float], [6 x float]* %0, i64 %5, i64 %139
%143 = load float, float* %142
%144 = fneg float %143
%145 = tail call float @llvm.fmuladd.f32(float %144, float %84, float %141)
store float %145, float* %140
%146 = add nuw nsw i64 %132, 2
%147 = icmp eq i64 %146, %61
br i1 %147, label %148, label %131
148: ; preds = %128, %131, %116, %77
%149 = add nuw nsw i64 %78, 1
%150 = icmp eq i64 %149, %60
br i1 %150, label %151, label %77
151: ; preds = %148, %4, %56
ret void
}