void mmadd(double alpha, double *A, int ALda, int Ai, int Aj, double beta,
           double *B, int BLda, int Bi, int Bj, double *C, int CLda, int Ci,
           int Cj, int m, int n) {
  int i, j;
  int offsetA = Aj * ALda + Ai;
  int offsetB = Bj * BLda + Bi;
  int offsetC = Cj * CLda + Ci;
  for (j = 0; j < n; ++j) {
    for (i = 0; i < m; ++i) {
      C[offsetC + i] = alpha * A[offsetA + i] + beta * B[offsetB + i];
    }
    offsetA += ALda;
    offsetB += BLda;
    offsetC += CLda;
  }
}


x86 assembly (clang, -O0):
<mmadd>:
push   %rbp
mov    %rsp,%rbp
mov    0x48(%rbp),%eax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movsd  %xmm1,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
mov    -0x1c(%rbp),%eax
imul   -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    0x18(%rbp),%eax
imul   -0x34(%rbp),%eax
add    0x10(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    0x38(%rbp),%eax
imul   0x28(%rbp),%eax
add    0x30(%rbp),%eax
mov    %eax,-0x48(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x48(%rbp),%eax
jge    1208 <mmadd+0x108>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x40(%rbp),%eax
jge    11df <mmadd+0xdf>
movsd  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mov    -0x30(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x20(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    117d <mmadd+0x7d>
mov    -0x14(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    0x28(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    116a <mmadd+0x6a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mmadd>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %r8,-0x70(%rsp)
mov    %rdi,-0x78(%rsp)
cmpl   $0x0,0x80(%rsp)
jle    147d <mmadd+0x37d>
mov    0x78(%rsp),%ebx
mov    0x60(%rsp),%eax
movslq %eax,%r14
mov    %r14,-0x48(%rsp)
imul   0x70(%rsp),%eax
add    0x68(%rsp),%eax
mov    0x58(%rsp),%r8
movslq %r9d,%r11
imul   0x50(%rsp),%r9d
add    0x48(%rsp),%r9d
imul   %esi,%ecx
add    %edx,%ecx
cltq
movslq %r9d,%rbp
movslq %ecx,%rcx
movslq %esi,%rdi
mov    %rax,%rsi
mov    %ebx,%ebx
add    %rbx,%rax
mov    %rax,-0x20(%rsp)
lea    (%rcx,%rbx,1),%rax
mov    %rax,-0x28(%rsp)
lea    (%rbx,%rbp,1),%rax
mov    %rax,-0x30(%rsp)
mov    %ebx,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x38(%rsp)
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
mov    %rbx,%rax
neg    %rax
mov    %rax,-0x18(%rsp)
lea    0x0(,%rbx,8),%rax
and    $0xffffffffffffffe0,%rax
mov    %rax,-0x40(%rsp)
mov    -0x78(%rsp),%rdx
lea    (%rdx,%rcx,8),%r9
add    $0x10,%r9
mov    %rdi,-0x58(%rsp)
lea    0x0(,%rdi,8),%rax
mov    %rax,0x8(%rsp)
mov    -0x70(%rsp),%rax
lea    0x10(%rax,%rbp,8),%r10
mov    %r11,-0x50(%rsp)
lea    0x0(,%r11,8),%rdi
mov    %rdi,(%rsp)
lea    0x10(%r8,%rsi,8),%r11
lea    0x0(,%r14,8),%rdi
mov    %rdi,-0x8(%rsp)
mov    %rsi,%r15
lea    0x8(%r8,%rsi,8),%r8
lea    (%rax,%rbp,8),%rdi
mov    %rcx,-0x80(%rsp)
lea    (%rdx,%rcx,8),%r12
mov    %rbp,%rdx
xor    %eax,%eax
mov    %rax,-0x68(%rsp)
jmp    1284 <mmadd+0x184>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    -0x80(%rsp),%rax
add    -0x58(%rsp),%rax
mov    %rax,-0x80(%rsp)
add    -0x50(%rsp),%rdx
add    -0x48(%rsp),%r15
mov    -0x68(%rsp),%rsi
add    $0x1,%rsi
mov    0x8(%rsp),%rax
add    %rax,%r9
mov    (%rsp),%rbp
add    %rbp,%r10
mov    -0x8(%rsp),%rcx
add    %rcx,%r11
add    %rcx,%r8
add    %rbp,%rdi
add    %rax,%r12
mov    %rsi,-0x68(%rsp)
cmp    0x80(%rsp),%esi
je     147d <mmadd+0x37d>
cmpl   $0x0,0x78(%rsp)
jle    1230 <mmadd+0x130>
cmpl   $0x4,0x78(%rsp)
jae    12a0 <mmadd+0x1a0>
xor    %r13d,%r13d
jmp    13c0 <mmadd+0x2c0>
nopw   0x0(%rax,%rax,1)
mov    -0x68(%rsp),%rax
mov    %rax,%rsi
imul   -0x48(%rsp),%rsi
add    -0x20(%rsp),%rsi
mov    %rax,%rbp
imul   -0x58(%rsp),%rbp
add    -0x28(%rsp),%rbp
mov    %rdx,-0x10(%rsp)
mov    %rax,%rdx
imul   -0x50(%rsp),%rdx
add    -0x30(%rsp),%rdx
mov    0x58(%rsp),%r14
mov    %r15,-0x60(%rsp)
lea    (%r14,%r15,8),%rax
lea    (%r14,%rsi,8),%rsi
mov    -0x78(%rsp),%rcx
lea    (%rcx,%rbp,8),%rbp
cmp    %rbp,%rax
mov    -0x80(%rsp),%rbp
lea    (%rcx,%rbp,8),%rbp
setb   %r13b
cmp    %rsi,%rbp
mov    -0x70(%rsp),%rbp
lea    0x0(%rbp,%rdx,8),%rdx
setb   %r14b
cmp    %rdx,%rax
mov    -0x10(%rsp),%rdx
lea    0x0(%rbp,%rdx,8),%rax
setb   %bpl
cmp    %rsi,%rax
setb   %al
test   %r14b,%r13b
jne    13ab <mmadd+0x2ab>
mov    $0x0,%r13d
and    %al,%bpl
mov    -0x60(%rsp),%r15
jne    13c0 <mmadd+0x2c0>
xor    %ebp,%ebp
mov    -0x40(%rsp),%rax
nopw   0x0(%rax,%rax,1)
movupd -0x10(%r9,%rbp,1),%xmm4
movupd (%r9,%rbp,1),%xmm5
movupd -0x10(%r10,%rbp,1),%xmm6
movupd (%r10,%rbp,1),%xmm7
mulpd  %xmm2,%xmm6
mulpd  %xmm2,%xmm7
mulpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
mulpd  %xmm3,%xmm5
addpd  %xmm7,%xmm5
movupd %xmm4,-0x10(%r11,%rbp,1)
movupd %xmm5,(%r11,%rbp,1)
add    $0x20,%rbp
cmp    %rbp,%rax
jne    1350 <mmadd+0x250>
mov    -0x38(%rsp),%rbp
mov    %rbp,%r13
cmp    %rbx,%rbp
je     1230 <mmadd+0x130>
jmp    13c0 <mmadd+0x2c0>
xor    %r13d,%r13d
mov    -0x60(%rsp),%r15
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r13,%r14
test   $0x1,%bl
je     1414 <mmadd+0x314>
mov    -0x80(%rsp),%rax
add    %r13,%rax
mov    %rdx,%rcx
lea    0x0(%r13,%rdx,1),%rdx
mov    -0x70(%rsp),%rsi
movsd  (%rsi,%rdx,8),%xmm4
mov    %r15,%rbp
mulsd  %xmm1,%xmm4
mov    -0x78(%rsp),%rdx
movsd  (%rdx,%rax,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
lea    0x0(%r13,%r15,1),%rax
mov    0x58(%rsp),%rdx
movsd  %xmm5,(%rdx,%rax,8)
mov    %rcx,%rdx
mov    %r13,%r14
or     $0x1,%r14
not    %r13
cmp    -0x18(%rsp),%r13
je     1230 <mmadd+0x130>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdi,%r14,8),%xmm4
mulsd  %xmm1,%xmm4
movsd  (%r12,%r14,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,-0x8(%r8,%r14,8)
movsd  0x8(%rdi,%r14,8),%xmm4
mulsd  %xmm1,%xmm4
movsd  0x8(%r12,%r14,8),%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,(%r8,%r14,8)
add    $0x2,%r14
cmp    %r14,%rbx
jne    1430 <mmadd+0x330>
jmp    1230 <mmadd+0x130>
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @mmadd(double %0, double* nocapture %1, i32 %2, i32 %3, i32 %4, double %5, double* nocapture %6, i32 %7, i32 %8, i32 %9, double* nocapture writeonly %10, i32 %11, i32 %12, i32 %13, i32 %14, i32 %15) {
%17 = icmp sgt i32 %15, 0
br i1 %17, label %18, label %155
18: ; preds = %16
%19 = mul nsw i32 %13, %11
%20 = add i32 %19, %12
%21 = mul nsw i32 %9, %7
%22 = add i32 %21, %8
%23 = mul nsw i32 %4, %2
%24 = add i32 %23, %3
%25 = icmp sgt i32 %14, 0
%26 = sext i32 %20 to i64
%27 = sext i32 %11 to i64
%28 = sext i32 %22 to i64
%29 = sext i32 %7 to i64
%30 = sext i32 %24 to i64
%31 = sext i32 %2 to i64
%32 = zext i32 %14 to i64
%33 = add nsw i64 %26, %32
%34 = add nsw i64 %30, %32
%35 = add nsw i64 %28, %32
%36 = icmp ult i32 %14, 4
%37 = and i64 %32, 4294967292
%38 = insertelement <2 x double> poison, double %5, i64 0
%39 = shufflevector <2 x double> %38, <2 x double> poison, <2 x i32> zeroinitializer
%40 = insertelement <2 x double> poison, double %5, i64 0
%41 = shufflevector <2 x double> %40, <2 x double> poison, <2 x i32> zeroinitializer
%42 = insertelement <2 x double> poison, double %0, i64 0
%43 = shufflevector <2 x double> %42, <2 x double> poison, <2 x i32> zeroinitializer
%44 = insertelement <2 x double> poison, double %0, i64 0
%45 = shufflevector <2 x double> %44, <2 x double> poison, <2 x i32> zeroinitializer
%46 = icmp eq i64 %37, %32
%47 = and i64 %32, 1
%48 = icmp eq i64 %47, 0
%49 = sub nsw i64 0, %32
br label %50
50: ; preds = %18, %148
%51 = phi i64 [ 0, %18 ], [ %154, %148 ]
%52 = phi i64 [ %30, %18 ], [ %149, %148 ]
%53 = phi i64 [ %28, %18 ], [ %150, %148 ]
%54 = phi i64 [ %26, %18 ], [ %151, %148 ]
%55 = phi i32 [ 0, %18 ], [ %152, %148 ]
%56 = getelementptr double, double* %10, i64 %54
%57 = mul i64 %51, %27
%58 = add i64 %33, %57
%59 = getelementptr double, double* %10, i64 %58
%60 = getelementptr double, double* %1, i64 %52
%61 = mul i64 %51, %31
%62 = add i64 %34, %61
%63 = getelementptr double, double* %1, i64 %62
%64 = getelementptr double, double* %6, i64 %53
%65 = mul i64 %51, %29
%66 = add i64 %35, %65
%67 = getelementptr double, double* %6, i64 %66
br i1 %25, label %68, label %148
68: ; preds = %50
br i1 %36, label %105, label %69
69: ; preds = %68
%70 = icmp ult double* %56, %63
%71 = icmp ult double* %60, %59
%72 = and i1 %70, %71
%73 = icmp ult double* %56, %67
%74 = icmp ult double* %64, %59
%75 = and i1 %73, %74
%76 = or i1 %72, %75
br i1 %76, label %105, label %77
77: ; preds = %69, %77
%78 = phi i64 [ %102, %77 ], [ 0, %69 ]
%79 = add nsw i64 %78, %52
%80 = getelementptr inbounds double, double* %1, i64 %79
%81 = bitcast double* %80 to <2 x double>*
%82 = load <2 x double>, <2 x double>* %81
%83 = getelementptr inbounds double, double* %80, i64 2
%84 = bitcast double* %83 to <2 x double>*
%85 = load <2 x double>, <2 x double>* %84
%86 = add nsw i64 %78, %53
%87 = getelementptr inbounds double, double* %6, i64 %86
%88 = bitcast double* %87 to <2 x double>*
%89 = load <2 x double>, <2 x double>* %88
%90 = getelementptr inbounds double, double* %87, i64 2
%91 = bitcast double* %90 to <2 x double>*
%92 = load <2 x double>, <2 x double>* %91
%93 = fmul <2 x double> %89, %39
%94 = fmul <2 x double> %92, %41
%95 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %43, <2 x double> %82, <2 x double> %93)
%96 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %45, <2 x double> %85, <2 x double> %94)
%97 = add nsw i64 %78, %54
%98 = getelementptr inbounds double, double* %10, i64 %97
%99 = bitcast double* %98 to <2 x double>*
store <2 x double> %95, <2 x double>* %99
%100 = getelementptr inbounds double, double* %98, i64 2
%101 = bitcast double* %100 to <2 x double>*
store <2 x double> %96, <2 x double>* %101
%102 = add nuw i64 %78, 4
%103 = icmp eq i64 %102, %37
br i1 %103, label %104, label %77
104: ; preds = %77
br i1 %46, label %148, label %105
105: ; preds = %69, %68, %104
%106 = phi i64 [ 0, %69 ], [ 0, %68 ], [ %37, %104 ]
%107 = xor i64 %106, -1
br i1 %48, label %120, label %108
108: ; preds = %105
%109 = add nsw i64 %106, %52
%110 = getelementptr inbounds double, double* %1, i64 %109
%111 = load double, double* %110
%112 = add nsw i64 %106, %53
%113 = getelementptr inbounds double, double* %6, i64 %112
%114 = load double, double* %113
%115 = fmul double %114, %5
%116 = tail call double @llvm.fmuladd.f64(double %0, double %111, double %115)
%117 = add nsw i64 %106, %54
%118 = getelementptr inbounds double, double* %10, i64 %117
store double %116, double* %118
%119 = or i64 %106, 1
br label %120
120: ; preds = %108, %105
%121 = phi i64 [ %106, %105 ], [ %119, %108 ]
%122 = icmp eq i64 %107, %49
br i1 %122, label %148, label %123
123: ; preds = %120, %123
%124 = phi i64 [ %146, %123 ], [ %121, %120 ]
%125 = add nsw i64 %124, %52
%126 = getelementptr inbounds double, double* %1, i64 %125
%127 = load double, double* %126
%128 = add nsw i64 %124, %53
%129 = getelementptr inbounds double, double* %6, i64 %128
%130 = load double, double* %129
%131 = fmul double %130, %5
%132 = tail call double @llvm.fmuladd.f64(double %0, double %127, double %131)
%133 = add nsw i64 %124, %54
%134 = getelementptr inbounds double, double* %10, i64 %133
store double %132, double* %134
%135 = add nuw nsw i64 %124, 1
%136 = add nsw i64 %135, %52
%137 = getelementptr inbounds double, double* %1, i64 %136
%138 = load double, double* %137
%139 = add nsw i64 %135, %53
%140 = getelementptr inbounds double, double* %6, i64 %139
%141 = load double, double* %140
%142 = fmul double %141, %5
%143 = tail call double @llvm.fmuladd.f64(double %0, double %138, double %142)
%144 = add nsw i64 %135, %54
%145 = getelementptr inbounds double, double* %10, i64 %144
store double %143, double* %145
%146 = add nuw nsw i64 %124, 2
%147 = icmp eq i64 %146, %32
br i1 %147, label %148, label %123
148: ; preds = %120, %123, %104, %50
%149 = add i64 %52, %31
%150 = add i64 %53, %29
%151 = add i64 %54, %27
%152 = add nuw nsw i32 %55, 1
%153 = icmp eq i32 %152, %15
%154 = add i64 %51, 1
br i1 %153, label %155, label %50
155: ; preds = %148, %16
ret void
}