void pfuslv(int neqns, double **penv, double *diag, double *rhs) {
  int i, k;
  double s, *ptr;
  for (i = neqns - 1; i >= 0; i--) {
    if (rhs[i] == 0.0e0)
      continue;
    s = rhs[i];
    k = i - 1;
    for (ptr = penv[i + 1] - 1; ptr >= penv[i]; ptr--, k--) {
      rhs[k] -= (*ptr * s);
    }
  }
  return;
}


x86 assembly (clang, -O0):
<pfuslv>:
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     1200 <pfuslv+0x100>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    114b <pfuslv+0x4b>
jp     114b <pfuslv+0x4b>
jmp    11f2 <pfuslv+0xf2>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
cmp    (%rcx,%rdx,8),%rax
jb     11ed <pfuslv+0xed>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x28(%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    -0x38(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
jmp    117f <pfuslv+0x7f>
jmp    11f2 <pfuslv+0xf2>
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    111c <pfuslv+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<pfuslv>:
test   %edi,%edi
jle    1183 <pfuslv+0x83>
movslq %edi,%rax
mov    %edi,%r8d
lea    (%rcx,%rax,8),%r10
add    $0xfffffffffffffff0,%r10
xorpd  %xmm0,%xmm0
jmp    112a <pfuslv+0x2a>
nopl   0x0(%rax,%rax,1)
add    $0xfffffffffffffff8,%r10
cmp    $0x1,%r9
jle    1183 <pfuslv+0x83>
mov    %r8,%r9
add    $0xffffffffffffffff,%r8
mov    %r8d,%eax
movsd  (%rcx,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1141 <pfuslv+0x41>
jnp    1120 <pfuslv+0x20>
mov    (%rsi,%r9,8),%rdx
mov    (%rsi,%rax,8),%rdi
add    $0xfffffffffffffff8,%rdx
cmp    %rdi,%rdx
jb     1120 <pfuslv+0x20>
mov    %r10,%rax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rax),%xmm2
movsd  (%rdx),%xmm3
mulsd  %xmm1,%xmm3
subsd  %xmm3,%xmm2
movsd  %xmm2,(%rax)
add    $0xfffffffffffffff8,%rdx
add    $0xfffffffffffffff8,%rax
cmp    %rdi,%rdx
jae    1160 <pfuslv+0x60>
jmp    1120 <pfuslv+0x20>
ret

LLVM IR (-O0):
define void @pfuslv(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
%11 = alloca double
%12 = alloca double*
store i32 %0, i32* %5
store double** %1, double*** %6
store double* %2, double** %7
store double* %3, double** %8
%13 = load i32, i32* %5
%14 = sub nsw i32 %13, 1
store i32 %14, i32* %9
br label %15
15: ; preds = %66, %4
%16 = load i32, i32* %9
%17 = icmp sge i32 %16, 0
br i1 %17, label %18, label %69
18: ; preds = %15
%19 = load double*, double** %8
%20 = load i32, i32* %9
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %19, i64 %21
%23 = load double, double* %22
%24 = fcmp oeq double %23, 0.000000e+00
br i1 %24, label %25, label %26
25: ; preds = %18
br label %66
26: ; preds = %18
%27 = load double*, double** %8
%28 = load i32, i32* %9
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
store double %31, double* %11
%32 = load i32, i32* %9
%33 = sub nsw i32 %32, 1
store i32 %33, i32* %10
%34 = load double**, double*** %6
%35 = load i32, i32* %9
%36 = add nsw i32 %35, 1
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double*, double** %34, i64 %37
%39 = load double*, double** %38
%40 = getelementptr inbounds double, double* %39, i64 -1
store double* %40, double** %12
br label %41
41: ; preds = %60, %26
%42 = load double*, double** %12
%43 = load double**, double*** %6
%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 = icmp uge double* %42, %47
br i1 %48, label %49, label %65
49: ; preds = %41
%50 = load double*, double** %12
%51 = load double, double* %50
%52 = load double, double* %11
%53 = load double*, double** %8
%54 = load i32, i32* %10
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = fneg double %51
%59 = call double @llvm.fmuladd.f64(double %58, double %52, double %57)
store double %59, double* %56
br label %60
60: ; preds = %49
%61 = load double*, double** %12
%62 = getelementptr inbounds double, double* %61, i32 -1
store double* %62, double** %12
%63 = load i32, i32* %10
%64 = add nsw i32 %63, -1
store i32 %64, i32* %10
br label %41
65: ; preds = %41
br label %66
66: ; preds = %65, %25
%67 = load i32, i32* %9
%68 = add nsw i32 %67, -1
store i32 %68, i32* %9
br label %15
69: ; preds = %15
ret void
}

LLVM IR (-O2):
define void @pfuslv(i32 %0, double** nocapture %1, double* nocapture readnone %2, double* nocapture %3) {
%5 = icmp sgt i32 %0, 0
br i1 %5, label %6, label %39
6: ; preds = %4
%7 = add nsw i32 %0, -2
%8 = sext i32 %7 to i64
%9 = zext i32 %0 to i64
br label %10
10: ; preds = %6, %36
%11 = phi i64 [ %9, %6 ], [ %13, %36 ]
%12 = phi i64 [ %8, %6 ], [ %38, %36 ]
%13 = add nsw i64 %11, -1
%14 = and i64 %13, 4294967295
%15 = getelementptr inbounds double, double* %3, i64 %14
%16 = load double, double* %15
%17 = fcmp oeq double %16, 0.000000e+00
br i1 %17, label %36, label %18
18: ; preds = %10
%19 = getelementptr inbounds double*, double** %1, i64 %11
%20 = load double*, double** %19
%21 = getelementptr inbounds double*, double** %1, i64 %14
%22 = load double*, double** %21
%23 = getelementptr inbounds double, double* %20, i64 -1
%24 = icmp ult double* %23, %22
br i1 %24, label %36, label %25
25: ; preds = %18, %25
%26 = phi i64 [ %33, %25 ], [ %12, %18 ]
%27 = phi double* [ %34, %25 ], [ %23, %18 ]
%28 = load double, double* %27
%29 = getelementptr inbounds double, double* %3, i64 %26
%30 = load double, double* %29
%31 = fneg double %28
%32 = tail call double @llvm.fmuladd.f64(double %31, double %16, double %30)
store double %32, double* %29
%33 = add i64 %26, -1
%34 = getelementptr inbounds double, double* %27, i64 -1
%35 = icmp ult double* %34, %22
br i1 %35, label %36, label %25
36: ; preds = %25, %18, %10
%37 = icmp sgt i64 %11, 1
%38 = add nsw i64 %12, -1
br i1 %37, label %10, label %39
39: ; preds = %36, %4
ret void
}