void multAB(double *C, double *A, double *B, const int rA, const int cA,
            const int cB) {
  int i, j, k;
  double *a, *c, tmp;
  for (i = 0; i < cB; i++) {
    for (k = 0; k < cA; k++) {
      tmp = B[i * cA + k];
      a = A + k * rA;
      c = C + i * rA;
      for (j = 0; j < rA; j++) {
        c[j] += tmp * a[j];
      }
    }
  }
}


x86 assembly (clang, -O0):
<multAB>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1201 <multAB+0x101>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11ee <multAB+0xee>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11db <multAB+0xdb>
movsd  -0x48(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1195 <multAB+0x95>
jmp    11e0 <multAB+0xe0>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1135 <multAB+0x35>
jmp    11f3 <multAB+0xf3>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1122 <multAB+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<multAB>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,-0x74(%rsp)
mov    %rdi,-0x68(%rsp)
test   %r9d,%r9d
jle    13e4 <multAB+0x2e4>
movslq %ecx,%r11
movslq -0x74(%rsp),%rax
mov    %r9d,%edi
mov    %rdi,-0x50(%rsp)
mov    %rax,-0x58(%rsp)
mov    %eax,%eax
mov    %rax,-0x8(%rsp)
mov    %r11d,%ebx
mov    %ebx,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x30(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x38(%rsp)
mov    %rax,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %rdi,%rax
mov    %rdi,-0x40(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,-0x48(%rsp)
mov    %rbx,%rax
neg    %rax
mov    %rax,-0x18(%rsp)
lea    0x30(%rsi),%rax
mov    %rax,-0x60(%rsp)
lea    0x0(,%r11,8),%r10
mov    -0x68(%rsp),%rax
lea    0x30(%rax),%r9
lea    0x8(%rax),%r13
xor    %eax,%eax
mov    %rax,-0x70(%rsp)
mov    %r11,-0x28(%rsp)
jmp    11d2 <multAB+0xd2>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x70(%rsp),%rdi
add    $0x1,%rdi
add    %r10,%r9
add    %r10,%r13
mov    %rdi,%rax
mov    %rdi,-0x70(%rsp)
cmp    -0x50(%rsp),%rdi
je     13e4 <multAB+0x2e4>
cmpl   $0x0,-0x74(%rsp)
jle    11b0 <multAB+0xb0>
mov    -0x70(%rsp),%rbp
mov    %rbp,%rax
imul   %r11,%rax
mov    -0x68(%rsp),%rdi
lea    (%rdi,%rax,8),%r14
add    %rbx,%rax
lea    (%rdi,%rax,8),%rax
mov    %rax,-0x20(%rsp)
imul   -0x58(%rsp),%rbp
mov    %rbp,-0x10(%rsp)
mov    %rsi,%r8
mov    -0x60(%rsp),%r15
xor    %r12d,%r12d
jmp    1235 <multAB+0x135>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r12
add    %r10,%r15
add    %r10,%r8
cmp    -0x8(%rsp),%r12
je     11b0 <multAB+0xb0>
test   %ecx,%ecx
jle    1220 <multAB+0x120>
mov    %r12,%rax
imul   %r11,%rax
mov    -0x10(%rsp),%rdi
add    %r12,%rdi
lea    (%rsi,%rax,8),%rbp
movsd  (%rdx,%rdi,8),%xmm0
cmp    $0x4,%ecx
jb     126d <multAB+0x16d>
add    %rbx,%rax
lea    (%rsi,%rax,8),%rax
cmp    %rax,%r14
jae    12df <multAB+0x1df>
cmp    -0x20(%rsp),%rbp
jae    12df <multAB+0x1df>
xor    %edi,%edi
mov    %rdi,%rax
test   $0x1,%bl
je     1294 <multAB+0x194>
movsd  0x0(%rbp,%rdi,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%r14,%rdi,8),%xmm1
movsd  %xmm1,(%r14,%rdi,8)
mov    %rdi,%rax
or     $0x1,%rax
not    %rdi
cmp    -0x18(%rsp),%rdi
je     1220 <multAB+0x120>
xchg   %ax,%ax
movsd  (%r8,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  -0x8(%r13,%rax,8),%xmm1
movsd  %xmm1,-0x8(%r13,%rax,8)
movsd  0x8(%r8,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x0(%r13,%rax,8),%xmm1
movsd  %xmm1,0x0(%r13,%rax,8)
add    $0x2,%rax
cmp    %rax,%rbx
jne    12a0 <multAB+0x1a0>
jmp    1220 <multAB+0x120>
mov    %rdx,%rax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,-0x38(%rsp)
je     13d8 <multAB+0x2d8>
mov    -0x48(%rsp),%rdx
xor    %r11d,%r11d
xchg   %ax,%ax
movupd -0x30(%r15,%r11,8),%xmm2
movupd -0x20(%r15,%r11,8),%xmm3
mulpd  %xmm1,%xmm2
movupd -0x30(%r9,%r11,8),%xmm4
addpd  %xmm2,%xmm4
movupd -0x20(%r9,%r11,8),%xmm2
movupd -0x10(%r9,%r11,8),%xmm5
movupd (%r9,%r11,8),%xmm6
mulpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm4,-0x30(%r9,%r11,8)
movupd %xmm3,-0x20(%r9,%r11,8)
movupd -0x10(%r15,%r11,8),%xmm2
movupd (%r15,%r11,8),%xmm3
mulpd  %xmm1,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm6,%xmm3
movupd %xmm2,-0x10(%r9,%r11,8)
movupd %xmm3,(%r9,%r11,8)
add    $0x8,%r11
add    $0xfffffffffffffffe,%rdx
jne    1300 <multAB+0x200>
testb  $0x1,-0x40(%rsp)
je     13ba <multAB+0x2ba>
movupd 0x0(%rbp,%r11,8),%xmm2
movupd 0x10(%rbp,%r11,8),%xmm3
movupd (%r14,%r11,8),%xmm4
movupd 0x10(%r14,%r11,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,(%r14,%r11,8)
movupd %xmm3,0x10(%r14,%r11,8)
mov    -0x30(%rsp),%rdx
mov    %rdx,%rdi
cmp    %rbx,%rdx
mov    %rax,%rdx
mov    -0x28(%rsp),%r11
je     1220 <multAB+0x120>
jmp    126f <multAB+0x16f>
xor    %r11d,%r11d
testb  $0x1,-0x40(%rsp)
jne    1382 <multAB+0x282>
jmp    13ba <multAB+0x2ba>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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