int lhsbnd_fill_upper(double *lhsbnd, double *G, double *w, long *kk, int n,
                      int nk) {
  int j, k;
  long i, i4, kkk;
  for (i = 0; i < n; ++i) {
    i4 = i * 4;
    kkk = kk[i] + 3 * nk;
    for (j = 0; j < 4; ++j)
      for (k = 0; k <= j; ++k)
        lhsbnd[kkk - nk * (j - k) + j] += w[i] * G[i4 + j] * G[i4 + k];
  }
  return 0;
}


x86 assembly (clang, -O0):
<lhsbnd_fill_upper>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
cmp    %rcx,%rax
jge    121e <lhsbnd_fill_upper+0x11e>
mov    -0x38(%rbp),%rax
shl    $0x2,%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
imul   $0x3,-0x28(%rbp),%ecx
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x2c(%rbp)
cmpl   $0x4,-0x2c(%rbp)
jge    1208 <lhsbnd_fill_upper+0x108>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     11f5 <lhsbnd_fill_upper+0xf5>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movslq -0x2c(%rbp),%rdx
add    %rdx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movslq -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    -0x28(%rbp),%edx
mov    -0x2c(%rbp),%esi
sub    -0x30(%rbp),%esi
imul   %esi,%edx
movslq %edx,%rdx
sub    %rdx,%rcx
movslq -0x2c(%rbp),%rdx
add    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1173 <lhsbnd_fill_upper+0x73>
jmp    11fa <lhsbnd_fill_upper+0xfa>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1162 <lhsbnd_fill_upper+0x62>
jmp    120d <lhsbnd_fill_upper+0x10d>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1124 <lhsbnd_fill_upper+0x24>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<lhsbnd_fill_upper>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %r8d,%r8d
jle    129f <lhsbnd_fill_upper+0x19f>
movslq %r8d,%r8
lea    (%r9,%r9,2),%eax
movslq %eax,%r14
movslq %r9d,%r9
lea    (%r9,%r9,1),%r10
lea    (%r9,%r9,2),%r11
neg    %r11
add    $0x18,%rsi
xor    %r13d,%r13d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rcx,%r13,8),%r15
lea    (%r15,%r14,1),%r12
movsd  -0x18(%rsi),%xmm0
movsd  (%rdx,%r13,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%r12,8),%xmm1
movsd  %xmm1,(%rdi,%r12,8)
lea    (%r15,%r14,1),%rax
add    $0x1,%rax
movsd  (%rdx,%r13,8),%xmm0
mulsd  -0x10(%rsi),%xmm0
sub    %r9,%rax
mulsd  -0x18(%rsi),%xmm0
addsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
movsd  -0x10(%rsi),%xmm0
movsd  (%rdx,%r13,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rdi,%r12,8),%xmm1
movsd  %xmm1,0x8(%rdi,%r12,8)
lea    0x2(%r15,%r14,1),%rbx
movsd  (%rdx,%r13,8),%xmm0
mulsd  -0x8(%rsi),%xmm0
mov    %rbx,%rax
mulsd  -0x18(%rsi),%xmm0
sub    %r10,%rax
addsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
movsd  (%rdx,%r13,8),%xmm0
mulsd  -0x8(%rsi),%xmm0
sub    %r9,%rbx
mulsd  -0x10(%rsi),%xmm0
addsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rdi,%rbx,8)
movsd  -0x8(%rsi),%xmm0
movsd  (%rdx,%r13,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%rdi,%r12,8),%xmm1
movsd  %xmm1,0x10(%rdi,%r12,8)
movsd  (%rdx,%r13,8),%xmm0
mulsd  (%rsi),%xmm0
lea    0x3(%r15,%r14,1),%rax
lea    0x3(%r11,%r12,1),%rbx
mulsd  -0x18(%rsi),%xmm0
addsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rdi,%rbx,8)
movsd  (%rdx,%r13,8),%xmm0
mulsd  (%rsi),%xmm0
mov    %rax,%rbx
mulsd  -0x10(%rsi),%xmm0
sub    %r10,%rbx
addsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rdi,%rbx,8)
movsd  (%rdx,%r13,8),%xmm0
mulsd  (%rsi),%xmm0
sub    %r9,%rax
mulsd  -0x8(%rsi),%xmm0
addsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
movsd  (%rsi),%xmm0
movsd  (%rdx,%r13,8),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  0x18(%rdi,%r12,8),%xmm1
movsd  %xmm1,0x18(%rdi,%r12,8)
add    $0x1,%r13
add    $0x20,%rsi
cmp    %r13,%r8
jne    1140 <lhsbnd_fill_upper+0x40>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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