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[i * n + j] * 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    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%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   %r14
push   %rbx
test   %edi,%edi
jle    1208 <solve_uy+0x108>
mov    %edi,%r8d
lea    -0x1(%r8),%rax
imul   %r8,%rax
lea    (%rsi,%rax,8),%r11
add    $0x8,%r11
lea    0x0(,%r8,8),%r9
neg    %r9
xor    %r10d,%r10d
mov    %r8,%rax
jmp    1166 <solve_uy+0x66>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %eax,%ebx
imul   %edi,%ebx
movslq %ebx,%rbx
add    %rax,%rbx
divsd  (%rsi,%rbx,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x1,%r10d
add    %r9,%r11
cmp    $0x1,%r14
jle    1208 <solve_uy+0x108>
mov    %rax,%r14
add    $0xffffffffffffffff,%rax
movsd  -0x8(%rcx,%r14,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%r14,8)
cmp    %r8,%r14
jge    1140 <solve_uy+0x40>
test   $0x1,%r10b
jne    11a0 <solve_uy+0xa0>
mov    %r14,%rbx
cmp    $0x1,%r10d
jne    11d0 <solve_uy+0xd0>
jmp    1140 <solve_uy+0x40>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rax,%rbx
imul   %r8,%rbx
add    %r14,%rbx
movsd  (%rsi,%rbx,8),%xmm1
mulsd  (%rdx,%r14,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rax,8)
lea    0x1(%r14),%rbx
cmp    $0x1,%r10d
je     1140 <solve_uy+0x40>
nopl   0x0(%rax)
movsd  -0x8(%r11,%rbx,8),%xmm1
mulsd  (%rdx,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rax,8)
movsd  (%r11,%rbx,8),%xmm1
mulsd  0x8(%rdx,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x2,%rbx
cmp    %ebx,%r8d
jne    11d0 <solve_uy+0xd0>
jmp    1140 <solve_uy+0x40>
pop    %rbx
pop    %r14
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* %9
%35 = load i32, i32* %5
%36 = mul nsw i32 %34, %35
%37 = load i32, i32* %10
%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 %66
6: ; preds = %4
%7 = zext i32 %0 to i64
br label %8
8: ; preds = %6, %56
%9 = phi i32 [ 0, %6 ], [ %65, %56 ]
%10 = phi i64 [ %7, %6 ], [ %11, %56 ]
%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 %56
17: ; preds = %8
%18 = mul nsw i64 %11, %7
%19 = and i32 %9, 1
%20 = icmp eq i32 %19, 0
br i1 %20, label %30, label %21
21: ; preds = %17
%22 = add nsw i64 %10, %18
%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 = %21, %17
%31 = phi double [ undef, %17 ], [ %28, %21 ]
%32 = phi double [ %14, %17 ], [ %28, %21 ]
%33 = phi i64 [ %10, %17 ], [ %29, %21 ]
%34 = icmp eq i32 %9, 1
br i1 %34, label %56, label %35
35: ; preds = %30, %35
%36 = phi double [ %52, %35 ], [ %32, %30 ]
%37 = phi i64 [ %53, %35 ], [ %33, %30 ]
%38 = add nsw i64 %37, %18
%39 = getelementptr inbounds double, double* %1, i64 %38
%40 = load double, double* %39
%41 = getelementptr inbounds double, double* %2, i64 %37
%42 = load double, double* %41
%43 = fneg double %40
%44 = tail call double @llvm.fmuladd.f64(double %43, double %42, double %36)
store double %44, double* %15
%45 = add nsw i64 %37, 1
%46 = add nsw i64 %45, %18
%47 = getelementptr inbounds double, double* %1, i64 %46
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %2, i64 %45
%50 = load double, double* %49
%51 = fneg double %48
%52 = tail call double @llvm.fmuladd.f64(double %51, double %50, double %44)
store double %52, double* %15
%53 = add nsw i64 %37, 2
%54 = trunc i64 %53 to i32
%55 = icmp eq i32 %54, %0
br i1 %55, label %56, label %35
56: ; preds = %30, %35, %8
%57 = phi double [ %14, %8 ], [ %31, %30 ], [ %52, %35 ]
%58 = mul nsw i32 %12, %0
%59 = sext i32 %58 to i64
%60 = add nsw i64 %11, %59
%61 = getelementptr inbounds double, double* %1, i64 %60
%62 = load double, double* %61
%63 = fdiv double %57, %62
store double %63, double* %15
%64 = icmp sgt i64 %10, 1
%65 = add i32 %9, 1
br i1 %64, label %8, label %66
66: ; preds = %56, %4
ret void
}