void LeastSquaresAddTerms(double **matrix, double **vectors,
                          const double *terms, const double *results,
                          const unsigned long rank, const unsigned long nvecs) {
  register unsigned long i, j;
  for (j = 0; j < rank; j++) {
    for (i = 0; i < rank; i++)
      matrix[i][j] += terms[i] * terms[j];
    for (i = 0; i < nvecs; i++)
      vectors[i][j] += results[i] * terms[j];
  }
  return;
}


x86 assembly (clang, -O0):
<LeastSquaresAddTerms>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    120e <LeastSquaresAddTerms+0x10e>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    1195 <LeastSquaresAddTerms+0x95>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    113a <LeastSquaresAddTerms+0x3a>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x30(%rbp),%rax
jae    11f8 <LeastSquaresAddTerms+0xf8>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    119d <LeastSquaresAddTerms+0x9d>
jmp    11fd <LeastSquaresAddTerms+0xfd>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1124 <LeastSquaresAddTerms+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<LeastSquaresAddTerms>:
push   %r14
push   %rbx
test   %r8,%r8
je     122b <LeastSquaresAddTerms+0x12b>
mov    %r8,%r11
and    $0xfffffffffffffffe,%r11
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r14d,%r14d
jmp    112d <LeastSquaresAddTerms+0x2d>
nop
add    $0x1,%r14
cmp    %r8,%r14
je     122b <LeastSquaresAddTerms+0x12b>
xor    %ebx,%ebx
cmp    $0x1,%r8
je     1181 <LeastSquaresAddTerms+0x81>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx,%rbx,8),%xmm0
mov    (%rdi,%rbx,8),%rax
mulsd  (%rdx,%r14,8),%xmm0
addsd  (%rax,%r14,8),%xmm0
movsd  %xmm0,(%rax,%r14,8)
movsd  0x8(%rdx,%rbx,8),%xmm0
mov    0x8(%rdi,%rbx,8),%rax
mulsd  (%rdx,%r14,8),%xmm0
addsd  (%rax,%r14,8),%xmm0
movsd  %xmm0,(%rax,%r14,8)
add    $0x2,%rbx
cmp    %rbx,%r11
jne    1140 <LeastSquaresAddTerms+0x40>
test   $0x1,%r8b
je     11a2 <LeastSquaresAddTerms+0xa2>
movsd  (%rdx,%rbx,8),%xmm0
mulsd  (%rdx,%r14,8),%xmm0
mov    (%rdi,%rbx,8),%rbx
addsd  (%rbx,%r14,8),%xmm0
movsd  %xmm0,(%rbx,%r14,8)
test   %r9,%r9
je     1120 <LeastSquaresAddTerms+0x20>
xor    %ebx,%ebx
cmp    $0x1,%r9
je     1201 <LeastSquaresAddTerms+0x101>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rcx,%rbx,8),%xmm0
mov    (%rsi,%rbx,8),%rax
mulsd  (%rdx,%r14,8),%xmm0
addsd  (%rax,%r14,8),%xmm0
movsd  %xmm0,(%rax,%r14,8)
movsd  0x8(%rcx,%rbx,8),%xmm0
mov    0x8(%rsi,%rbx,8),%rax
mulsd  (%rdx,%r14,8),%xmm0
addsd  (%rax,%r14,8),%xmm0
movsd  %xmm0,(%rax,%r14,8)
add    $0x2,%rbx
cmp    %rbx,%r10
jne    11c0 <LeastSquaresAddTerms+0xc0>
test   $0x1,%r9b
je     1120 <LeastSquaresAddTerms+0x20>
movsd  (%rcx,%rbx,8),%xmm0
mulsd  (%rdx,%r14,8),%xmm0
mov    (%rsi,%rbx,8),%rax
addsd  (%rax,%r14,8),%xmm0
movsd  %xmm0,(%rax,%r14,8)
jmp    1120 <LeastSquaresAddTerms+0x20>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @LeastSquaresAddTerms(double** %0, double** %1, double* %2, double* %3, i64 %4, i64 %5) {
%7 = alloca double**
%8 = alloca double**
%9 = alloca double*
%10 = alloca double*
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
store double** %0, double*** %7
store double** %1, double*** %8
store double* %2, double** %9
store double* %3, double** %10
store i64 %4, i64* %11
store i64 %5, i64* %12
store i64 0, i64* %14
br label %15
15: ; preds = %70, %6
%16 = load i64, i64* %14
%17 = load i64, i64* %11
%18 = icmp ult i64 %16, %17
br i1 %18, label %19, label %73
19: ; preds = %15
store i64 0, i64* %13
br label %20
20: ; preds = %41, %19
%21 = load i64, i64* %13
%22 = load i64, i64* %11
%23 = icmp ult i64 %21, %22
br i1 %23, label %24, label %44
24: ; preds = %20
%25 = load double*, double** %9
%26 = load i64, i64* %13
%27 = getelementptr inbounds double, double* %25, i64 %26
%28 = load double, double* %27
%29 = load double*, double** %9
%30 = load i64, i64* %14
%31 = getelementptr inbounds double, double* %29, i64 %30
%32 = load double, double* %31
%33 = load double**, double*** %7
%34 = load i64, i64* %13
%35 = getelementptr inbounds double*, double** %33, i64 %34
%36 = load double*, double** %35
%37 = load i64, i64* %14
%38 = getelementptr inbounds double, double* %36, i64 %37
%39 = load double, double* %38
%40 = call double @llvm.fmuladd.f64(double %28, double %32, double %39)
store double %40, double* %38
br label %41
41: ; preds = %24
%42 = load i64, i64* %13
%43 = add i64 %42, 1
store i64 %43, i64* %13
br label %20
44: ; preds = %20
store i64 0, i64* %13
br label %45
45: ; preds = %66, %44
%46 = load i64, i64* %13
%47 = load i64, i64* %12
%48 = icmp ult i64 %46, %47
br i1 %48, label %49, label %69
49: ; preds = %45
%50 = load double*, double** %10
%51 = load i64, i64* %13
%52 = getelementptr inbounds double, double* %50, i64 %51
%53 = load double, double* %52
%54 = load double*, double** %9
%55 = load i64, i64* %14
%56 = getelementptr inbounds double, double* %54, i64 %55
%57 = load double, double* %56
%58 = load double**, double*** %8
%59 = load i64, i64* %13
%60 = getelementptr inbounds double*, double** %58, i64 %59
%61 = load double*, double** %60
%62 = load i64, i64* %14
%63 = getelementptr inbounds double, double* %61, i64 %62
%64 = load double, double* %63
%65 = call double @llvm.fmuladd.f64(double %53, double %57, double %64)
store double %65, double* %63
br label %66
66: ; preds = %49
%67 = load i64, i64* %13
%68 = add i64 %67, 1
store i64 %68, i64* %13
br label %45
69: ; preds = %45
br label %70
70: ; preds = %69
%71 = load i64, i64* %14
%72 = add i64 %71, 1
store i64 %72, i64* %14
br label %15
73: ; preds = %15
ret void
}

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