void solve_uy(int n, double *u, double *var, double *r) {
  int i, j;
  for (i = n - 1; i >= 0; i--) {
    var[i] = r[i];
    for (j = i + 1; j < n; j++)
      var[i] -= (u[j * n + i] * var[j]);
    var[i] /= u[i * n + i];
  }
}


x86 assembly (clang, -O0):
<solve_uy>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     11f1 <solve_uy+0xf1>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b7 <solve_uy+0xb7>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x24(%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    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1149 <solve_uy+0x49>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    111c <solve_uy+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<solve_uy>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
jle    11fb <solve_uy+0xfb>
mov    %edi,%r15d
lea    -0x1(%r15),%r8
lea    (%r15,%r15,1),%r9
xor    %r10d,%r10d
mov    %rsi,%r11
mov    %r15,%rbx
jmp    1156 <solve_uy+0x56>
nopl   0x0(%rax,%rax,1)
mov    %ebx,%eax
imul   %edi,%eax
cltq
add    %rbx,%rax
divsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rbx,8)
add    $0x1,%r10d
add    $0xfffffffffffffff8,%r11
cmp    $0x1,%r14
jle    11fb <solve_uy+0xfb>
mov    %rbx,%r14
add    $0xffffffffffffffff,%rbx
movsd  -0x8(%rcx,%r14,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%r14,8)
cmp    %r15,%r14
jge    1130 <solve_uy+0x30>
mov    %r14,%rax
test   $0x1,%r10b
je     1198 <solve_uy+0x98>
imul   %r15,%rax
add    %rbx,%rax
movsd  (%rsi,%rax,8),%xmm1
mulsd  (%rdx,%r14,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rbx,8)
lea    0x1(%r14),%rax
cmp    $0x1,%r10d
je     1130 <solve_uy+0x30>
mov    %r15,%rbp
imul   %rax,%rbp
lea    (%r11,%rbp,8),%r12
lea    0x1(%rax),%rbp
imul   %r15,%rbp
lea    (%r11,%rbp,8),%r13
mov    %r8,%rbp
nopl   0x0(%rax,%rax,1)
movsd  (%r12,%rbp,8),%xmm1
mulsd  (%rdx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rbx,8)
movsd  0x0(%r13,%rbp,8),%xmm1
mulsd  0x8(%rdx,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rbx,8)
add    $0x2,%rax
add    %r9,%rbp
cmp    %eax,%r15d
jne    11c0 <solve_uy+0xc0>
jmp    1130 <solve_uy+0x30>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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