void finishSolvingmakeLowerDiag(double aug_matrix[][10 * 10 + 1], int n) {
  int i, j;
  double ans, knownterms;
  for (i = n - 1; i >= 0; i--) {
    knownterms = 0;
    for (j = i + 1; j < n; j++)
      knownterms += (aug_matrix[i][j] * aug_matrix[j][n]);
    ans = (aug_matrix[i][n] - knownterms) / aug_matrix[i][i];
    aug_matrix[i][n] = ans;
    for (j = 0; j < n; j++) {
      if (i != j)
        aug_matrix[i][j] = 0;
      else
        aug_matrix[i][j] = 1;
    }
  }
}


x86 assembly (clang, -O0):
<finishSolvingmakeLowerDiag>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x10(%rbp)
cmpl   $0x0,-0x10(%rbp)
jl     127c <finishSolvingmakeLowerDiag+0x17c>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1191 <finishSolvingmakeLowerDiag+0x91>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    112f <finishSolvingmakeLowerDiag+0x2f>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0x10(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0xc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1269 <finishSolvingmakeLowerDiag+0x169>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
je     1233 <finishSolvingmakeLowerDiag+0x133>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1256 <finishSolvingmakeLowerDiag+0x156>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
imul   $0x328,%rcx,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
movsd  0xdaf(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    125b <finishSolvingmakeLowerDiag+0x15b>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11f8 <finishSolvingmakeLowerDiag+0xf8>
jmp    126e <finishSolvingmakeLowerDiag+0x16e>
mov    -0x10(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x10(%rbp)
jmp    1114 <finishSolvingmakeLowerDiag+0x14>
pop    %rbp
ret

x86 assembly (clang, -O2):
<finishSolvingmakeLowerDiag>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %esi,%esi
jle    12f2 <finishSolvingmakeLowerDiag+0x1f2>
mov    %esi,%r15d
lea    -0x1(%r15),%r8
mov    %r15d,%r13d
and    $0x3,%r13d
mov    %r15d,%r12d
and    $0xfffffffc,%r12d
imul   $0x328,%r8,%rcx
add    %rdi,%rcx
lea    (%rdi,%r15,8),%r9
xor    %r10d,%r10d
xorpd  %xmm0,%xmm0
movsd  0xebd(%rip),%xmm1
mov    %r8,%r11
mov    %r15,%rbx
jmp    1169 <finishSolvingmakeLowerDiag+0x69>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10d
add    $0xfffffffffffffcd8,%rcx
add    $0xffffffffffffffff,%r11
cmp    $0x2,%r14
jl     12f2 <finishSolvingmakeLowerDiag+0x1f2>
mov    %rbx,%r14
add    $0xffffffffffffffff,%rbx
xorpd  %xmm2,%xmm2
cmp    %r15,%r14
jge    1203 <finishSolvingmakeLowerDiag+0x103>
xorpd  %xmm3,%xmm3
mov    %r14,%rax
test   $0x1,%r10b
je     11b6 <finishSolvingmakeLowerDiag+0xb6>
imul   $0x328,%rbx,%rax
add    %rdi,%rax
movsd  (%rax,%r14,8),%xmm2
imul   $0x328,%r14,%rax
add    %rdi,%rax
mulsd  (%rax,%r15,8),%xmm2
addsd  %xmm0,%xmm2
lea    0x1(%r14),%rax
movapd %xmm2,%xmm3
cmp    $0x1,%r10d
je     1203 <finishSolvingmakeLowerDiag+0x103>
imul   $0x328,%rax,%rsi
add    %r9,%rsi
cs nopw 0x0(%rax,%rax,1)
movsd  (%rcx,%rax,8),%xmm4
mulsd  (%rsi),%xmm4
movsd  0x8(%rcx,%rax,8),%xmm2
addsd  %xmm3,%xmm4
mulsd  0x328(%rsi),%xmm2
addsd  %xmm4,%xmm2
add    $0x2,%rax
add    $0x650,%rsi
movapd %xmm2,%xmm3
cmp    %eax,%r15d
jne    11d0 <finishSolvingmakeLowerDiag+0xd0>
imul   $0x328,%rbx,%rax
add    %rdi,%rax
movsd  (%rax,%r15,8),%xmm3
subsd  %xmm2,%xmm3
divsd  (%rax,%rbx,8),%xmm3
movsd  %xmm3,(%rax,%r15,8)
cmp    $0x3,%r8
jae    1280 <finishSolvingmakeLowerDiag+0x180>
xor    %eax,%eax
test   %r13,%r13
je     1150 <finishSolvingmakeLowerDiag+0x50>
lea    (%rcx,%rax,8),%rsi
mov    %r11,%rdx
sub    %rax,%rdx
xor    %eax,%eax
jmp    1262 <finishSolvingmakeLowerDiag+0x162>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,(%rsi,%rax,8)
add    $0x1,%rax
cmp    %rax,%r13
je     1150 <finishSolvingmakeLowerDiag+0x50>
movapd %xmm1,%xmm2
cmp    %rax,%rdx
je     1250 <finishSolvingmakeLowerDiag+0x150>
xorpd  %xmm2,%xmm2
jmp    1250 <finishSolvingmakeLowerDiag+0x150>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    -0x1(%r14),%rsi
xor    %eax,%eax
jmp    129f <finishSolvingmakeLowerDiag+0x19f>
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,0x18(%rcx,%rax,8)
add    $0x4,%rax
cmp    %r12,%rax
je     122a <finishSolvingmakeLowerDiag+0x12a>
movapd %xmm1,%xmm2
cmp    %rax,%rsi
je     12ac <finishSolvingmakeLowerDiag+0x1ac>
xorpd  %xmm2,%xmm2
movsd  %xmm2,(%rcx,%rax,8)
lea    0x1(%rax),%rdx
movapd %xmm1,%xmm2
cmp    %rdx,%rbx
je     12c2 <finishSolvingmakeLowerDiag+0x1c2>
xorpd  %xmm2,%xmm2
movsd  %xmm2,0x8(%rcx,%rax,8)
lea    0x2(%rax),%rdx
movapd %xmm1,%xmm2
cmp    %rdx,%rbx
je     12d9 <finishSolvingmakeLowerDiag+0x1d9>
xorpd  %xmm2,%xmm2
movsd  %xmm2,0x10(%rcx,%rax,8)
lea    0x3(%rax),%rdx
movapd %xmm1,%xmm2
cmp    %rdx,%rbx
je     1290 <finishSolvingmakeLowerDiag+0x190>
xorpd  %xmm2,%xmm2
jmp    1290 <finishSolvingmakeLowerDiag+0x190>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @finishSolvingmakeLowerDiag([101 x double]* nocapture %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %91
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = add nsw i64 %5, -1
%7 = and i64 %5, 3
%8 = icmp ult i64 %6, 3
%9 = and i64 %5, 4294967292
%10 = icmp eq i64 %7, 0
br label %25
11: ; preds = %70, %62
%12 = phi i64 [ 0, %62 ], [ %88, %70 ]
br i1 %10, label %22, label %13
13: ; preds = %11, %13
%14 = phi i64 [ %19, %13 ], [ %12, %11 ]
%15 = phi i64 [ %20, %13 ], [ 0, %11 ]
%16 = icmp eq i64 %28, %14
%17 = select i1 %16, double 1.000000e+00, double 0.000000e+00
%18 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %14
store double %17, double* %18
%19 = add nuw nsw i64 %14, 1
%20 = add i64 %15, 1
%21 = icmp eq i64 %20, %7
br i1 %21, label %22, label %13
22: ; preds = %13, %11
%23 = icmp sgt i64 %27, 1
%24 = add i32 %26, 1
br i1 %23, label %25, label %91
25: ; preds = %4, %22
%26 = phi i32 [ 0, %4 ], [ %24, %22 ]
%27 = phi i64 [ %5, %4 ], [ %28, %22 ]
%28 = add nsw i64 %27, -1
%29 = icmp slt i64 %27, %5
br i1 %29, label %30, label %62
30: ; preds = %25
%31 = and i32 %26, 1
%32 = icmp eq i32 %31, 0
br i1 %32, label %40, label %33
33: ; preds = %30
%34 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %27
%35 = load double, double* %34
%36 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %27, i64 %5
%37 = load double, double* %36
%38 = tail call double @llvm.fmuladd.f64(double %35, double %37, double 0.000000e+00)
%39 = add nsw i64 %27, 1
br label %40
40: ; preds = %33, %30
%41 = phi double [ undef, %30 ], [ %38, %33 ]
%42 = phi i64 [ %27, %30 ], [ %39, %33 ]
%43 = phi double [ 0.000000e+00, %30 ], [ %38, %33 ]
%44 = icmp eq i32 %26, 1
br i1 %44, label %62, label %45
45: ; preds = %40, %45
%46 = phi i64 [ %59, %45 ], [ %42, %40 ]
%47 = phi double [ %58, %45 ], [ %43, %40 ]
%48 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %46
%49 = load double, double* %48
%50 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %46, i64 %5
%51 = load double, double* %50
%52 = tail call double @llvm.fmuladd.f64(double %49, double %51, double %47)
%53 = add nsw i64 %46, 1
%54 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %53
%55 = load double, double* %54
%56 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %53, i64 %5
%57 = load double, double* %56
%58 = tail call double @llvm.fmuladd.f64(double %55, double %57, double %52)
%59 = add nsw i64 %46, 2
%60 = trunc i64 %59 to i32
%61 = icmp eq i32 %60, %1
br i1 %61, label %62, label %45
62: ; preds = %40, %45, %25
%63 = phi double [ 0.000000e+00, %25 ], [ %41, %40 ], [ %58, %45 ]
%64 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %5
%65 = load double, double* %64
%66 = fsub double %65, %63
%67 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %28
%68 = load double, double* %67
%69 = fdiv double %66, %68
store double %69, double* %64
br i1 %8, label %11, label %70
70: ; preds = %62, %70
%71 = phi i64 [ %88, %70 ], [ 0, %62 ]
%72 = phi i64 [ %89, %70 ], [ 0, %62 ]
%73 = icmp eq i64 %28, %71
%74 = select i1 %73, double 1.000000e+00, double 0.000000e+00
%75 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %71
store double %74, double* %75
%76 = or i64 %71, 1
%77 = icmp eq i64 %28, %76
%78 = select i1 %77, double 1.000000e+00, double 0.000000e+00
%79 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %76
store double %78, double* %79
%80 = or i64 %71, 2
%81 = icmp eq i64 %28, %80
%82 = select i1 %81, double 1.000000e+00, double 0.000000e+00
%83 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %80
store double %82, double* %83
%84 = or i64 %71, 3
%85 = icmp eq i64 %28, %84
%86 = select i1 %85, double 1.000000e+00, double 0.000000e+00
%87 = getelementptr inbounds [101 x double], [101 x double]* %0, i64 %28, i64 %84
store double %86, double* %87
%88 = add nuw nsw i64 %71, 4
%89 = add i64 %72, 4
%90 = icmp eq i64 %89, %9
br i1 %90, label %11, label %70
91: ; preds = %22, %2
ret void
}