void sustitucion(double **a, double *x, int filas, double *b) {
  int i, j, n;
  float sum;
  n = filas;
  x[n - 1] = b[n - 1] / a[n - 1][n - 1];
  for (i = n - 2; i >= 0; i--) {
    sum = b[i];
    for (j = i + 1; j <= n - 1; j++) {
      sum = sum - a[i][j] * x[j];
    }
    x[i] = sum / a[i][i];
  }
}


x86 assembly (clang, -O0):
<sustitucion>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     1239 <sustitucion+0x139>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1200 <sustitucion+0x100>
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    118e <sustitucion+0x8e>
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1165 <sustitucion+0x65>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sustitucion>:
movslq %edx,%rax
movsd  -0x8(%rcx,%rax,8),%xmm0
mov    -0x8(%rdi,%rax,8),%rdx
divsd  -0x8(%rdx,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rsi,%rax,8)
cmp    $0x2,%eax
jl     11fe <sustitucion+0xfe>
lea    -0x1(%rax),%r9
add    $0xfffffffe,%eax
mov    %eax,%r10d
xor    %r8d,%r8d
jmp    1164 <sustitucion+0x64>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cvtss2sd %xmm0,%xmm0
divsd  (%rax,%r10,8),%xmm0
movsd  %xmm0,(%rsi,%r10,8)
lea    -0x1(%r10),%rax
add    $0x1,%r8d
test   %r10,%r10
mov    %rax,%r10
jle    11fe <sustitucion+0xfe>
movsd  (%rcx,%r10,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    (%rdi,%r10,8),%rax
cmp    %r9,%r10
jge    1140 <sustitucion+0x40>
test   $0x1,%r8b
jne    11b0 <sustitucion+0xb0>
lea    0x1(%r10),%rdx
movsd  0x8(%rax,%r10,8),%xmm1
mulsd  0x8(%rsi,%r10,8),%xmm1
cvtss2sd %xmm0,%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
test   %r8d,%r8d
jne    11c0 <sustitucion+0xc0>
jmp    1140 <sustitucion+0x40>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r10,%rdx
test   %r8d,%r8d
je     1140 <sustitucion+0x40>
nopl   0x0(%rax,%rax,1)
cvtss2sd %xmm0,%xmm0
movsd  0x8(%rax,%rdx,8),%xmm1
mulsd  0x8(%rsi,%rdx,8),%xmm1
movsd  0x10(%rax,%rdx,8),%xmm2
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  0x10(%rsi,%rdx,8),%xmm2
add    $0x2,%rdx
subsd  %xmm2,%xmm0
cvtsd2ss %xmm0,%xmm0
cmp    %edx,%r9d
jne    11c0 <sustitucion+0xc0>
jmp    1140 <sustitucion+0x40>
ret

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

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