int reduceRhs(double *ma, double *vec, int i0, int dim, int *profile) {
  int i, k, counter;
  double matrix_i_0;
  counter = 0;
  for (i = 0; i < dim; i++) {
    matrix_i_0 = ma[counter];
    counter++;
    for (k = 1; k < profile[i]; k++) {
      vec[i0 * dim + i] -= ma[counter] * vec[i0 * dim + i - k];
      counter++;
    }
    if (matrix_i_0 > 0.0) {
      vec[i0 * dim + i] /= matrix_i_0;
    } else {
      vec[i0 * dim + i] = 0.0;
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<reduceRhs>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1244 <reduceRhs+0x144>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%rcx
movslq -0x24(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jge    11dd <reduceRhs+0xdd>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x24(%rbp),%ecx
sub    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%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    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1152 <reduceRhs+0x52>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1218 <reduceRhs+0x118>
movsd  -0x38(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1231 <reduceRhs+0x131>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1236 <reduceRhs+0x136>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1124 <reduceRhs+0x24>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<reduceRhs>:
push   %r15
push   %r14
push   %r12
push   %rbx
test   %ecx,%ecx
jle    1216 <reduceRhs+0x116>
imul   %ecx,%edx
movslq %edx,%r10
mov    %ecx,%r9d
lea    (%rsi,%r10,8),%r15
xor    %r11d,%r11d
xor    %eax,%eax
jmp    1146 <reduceRhs+0x46>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x1,%r11
add    $0x8,%r15
cmp    %r9,%r11
je     1216 <reduceRhs+0x116>
cltq
movsd  (%rdi,%rax,8),%xmm0
mov    (%r8,%r11,4),%r14d
add    $0x1,%eax
cmp    $0x2,%r14d
jl     11f6 <reduceRhs+0xf6>
mov    %r14d,%edx
lea    (%r11,%r10,1),%rcx
cltq
movsd  (%rsi,%rcx,8),%xmm1
add    $0xffffffffffffffff,%r14
mov    $0x1,%ebx
cmp    $0x2,%edx
je     11d3 <reduceRhs+0xd3>
mov    %r14,%r12
and    $0xfffffffffffffffe,%r12
mov    $0xffffffffffffffff,%rbx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm2
mulsd  (%r15,%rbx,8),%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
movsd  0x8(%rdi,%rax,8),%xmm2
mulsd  -0x8(%r15,%rbx,8),%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x2,%rax
lea    (%r12,%rbx,1),%rdx
add    $0xfffffffffffffffe,%rdx
add    $0xfffffffffffffffe,%rbx
cmp    $0xffffffffffffffff,%rdx
jne    1190 <reduceRhs+0x90>
neg    %rbx
test   $0x1,%r14b
je     11f6 <reduceRhs+0xf6>
movsd  (%rdi,%rax,8),%xmm2
mov    %rcx,%rdx
sub    %rbx,%rdx
mulsd  (%rsi,%rdx,8),%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x1,%rax
lea    (%r11,%r10,1),%rcx
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1130 <reduceRhs+0x30>
movsd  (%rsi,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
jmp    1130 <reduceRhs+0x30>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

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

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