void ch_fwd_sub(const double *_rr, double *_x, int _k, int _n) {
  int i;
  for (i = 0; i < _k; i++) {
    int j;
    _x[i] /= _rr[(((((_n) * ((_n) + 1) >> 1)) -
                   ((((_n) - (i)) * (((_n) - (i)) + 1) >> 1))) +
                  (i) - (i))];
    for (j = i + 1; j < _k; j++)
      _x[j] -= _rr[(((((_n) * ((_n) + 1) >> 1)) -
                     ((((_n) - (i)) * (((_n) - (i)) + 1) >> 1))) +
                    (j) - (i))] *
               _x[i];
  }
}


x86 assembly (clang, -O0):
<ch_fwd_sub>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    121f <ch_fwd_sub+0x11f>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x18(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
sar    $0x1,%ecx
mov    -0x18(%rbp),%edx
sub    -0x1c(%rbp),%edx
mov    -0x18(%rbp),%esi
sub    -0x1c(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
sar    $0x1,%edx
sub    %edx,%ecx
add    -0x1c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    120c <ch_fwd_sub+0x10c>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    -0x18(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
sar    $0x1,%ecx
mov    -0x18(%rbp),%edx
sub    -0x1c(%rbp),%edx
mov    -0x18(%rbp),%esi
sub    -0x1c(%rbp),%esi
add    $0x1,%esi
imul   %esi,%edx
sar    $0x1,%edx
sub    %edx,%ecx
add    -0x20(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x20(%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    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    117c <ch_fwd_sub+0x7c>
jmp    1211 <ch_fwd_sub+0x111>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <ch_fwd_sub+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ch_fwd_sub>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
test   %edx,%edx
jle    11ed <ch_fwd_sub+0xed>
lea    0x1(%rcx),%r9d
imul   %ecx,%r9d
sar    %r9d
mov    %edx,%r10d
movslq %r9d,%r8
mov    $0x1,%r11d
xor    %r14d,%r14d
mov    %rsi,%rdx
mov    %r10,%r15
jmp    1155 <ch_fwd_sub+0x55>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r11
add    $0xffffffffffffffff,%r15
add    $0x8,%rdx
cmp    %r10,%r14
je     11ed <ch_fwd_sub+0xed>
mov    %r14,%r12
mov    %ecx,%eax
sub    %r12d,%eax
lea    0x1(%rax),%ebp
imul   %eax,%ebp
sar    %ebp
mov    %r9d,%eax
sub    %ebp,%eax
movslq %eax,%r14
movsd  (%rsi,%r12,8),%xmm0
divsd  (%rdi,%r14,8),%xmm0
movsd  %xmm0,(%rsi,%r12,8)
lea    0x1(%r12),%r14
cmp    %r10,%r14
jae    1140 <ch_fwd_sub+0x40>
sub    %r12d,%eax
add    %r11d,%eax
cltq
movsd  (%rsi,%r11,8),%xmm1
mulsd  (%rdi,%rax,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%r11,8)
lea    0x1(%r11),%rax
cmp    %r10,%rax
je     1140 <ch_fwd_sub+0x40>
neg    %ebp
movslq %ebp,%rbp
add    %r8,%rbp
mov    $0x2,%eax
nopl   0x0(%rax)
lea    (%rax,%rbp,1),%ebx
movslq %ebx,%rbx
movsd  (%rdi,%rbx,8),%xmm0
movsd  (%rdx,%rax,8),%xmm1
mulsd  (%rsi,%r12,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rax,8)
add    $0x1,%rax
cmp    %rax,%r15
jne    11c0 <ch_fwd_sub+0xc0>
jmp    1140 <ch_fwd_sub+0x40>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

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