void MyGemm(int m, int n, int k, double *A, int ldA, double *B, int ldB,
            double *C, int ldC) {
  for (int j = 0; j < n; j += 4) {
    int jb = ((n - j) < (4) ? (n - j) : (4));
    for (int i = 0; i < m; i += 4) {
      int ib = ((m - i) < (4) ? (m - i) : (4));
      Gemm_P_Ger(ib, jb, k, &A[(0) * ldA + (i)], ldA, &B[(j)*ldB + (0)], ldB,
                 &C[(j)*ldC + (i)], ldC);
    }
  }
}


x86 assembly (clang, -O0):
<MyGemm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1248 <MyGemm+0x138>
mov    -0x8(%rbp),%eax
sub    -0x2c(%rbp),%eax
cmp    $0x4,%eax
jge    1167 <MyGemm+0x57>
mov    -0x8(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
jmp    1174 <MyGemm+0x64>
mov    $0x4,%eax
mov    %eax,-0x3c(%rbp)
jmp    1174 <MyGemm+0x64>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1235 <MyGemm+0x125>
mov    -0x4(%rbp),%eax
sub    -0x34(%rbp),%eax
cmp    $0x4,%eax
jge    11aa <MyGemm+0x9a>
mov    -0x4(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x40(%rbp)
jmp    11b7 <MyGemm+0xa7>
mov    $0x4,%eax
mov    %eax,-0x40(%rbp)
jmp    11b7 <MyGemm+0xa7>
mov    -0x40(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%edi
mov    -0x30(%rbp),%esi
mov    -0xc(%rbp),%edx
mov    -0x18(%rbp),%rcx
imul   $0x0,-0x1c(%rbp),%eax
add    -0x34(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rcx
mov    -0x1c(%rbp),%r8d
mov    -0x28(%rbp),%r9
mov    -0x2c(%rbp),%eax
imul   0x10(%rbp),%eax
add    $0x0,%eax
cltq
shl    $0x3,%rax
add    %rax,%r9
mov    0x10(%rbp),%r11d
mov    0x18(%rbp),%r10
mov    -0x2c(%rbp),%eax
imul   0x20(%rbp),%eax
add    -0x34(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%r10
mov    0x20(%rbp),%eax
mov    %r11d,(%rsp)
mov    %r10,0x8(%rsp)
mov    %eax,0x10(%rsp)
mov    $0x0,%al
call   1030 <Gemm_P_Ger@plt>
mov    -0x34(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x34(%rbp)
jmp    1181 <MyGemm+0x71>
jmp    123a <MyGemm+0x12a>
mov    -0x2c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x2c(%rbp)
jmp    113e <MyGemm+0x2e>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<MyGemm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r9,0x28(%rsp)
mov    %rcx,0x20(%rsp)
mov    %edi,0x4(%rsp)
test   %esi,%esi
jle    1252 <MyGemm+0x142>
mov    %r8d,%r14d
mov    %edx,%ebp
mov    0xa0(%rsp),%r10d
mov    0x98(%rsp),%rax
mov    %rax,0x8(%rsp)
mov    0x90(%rsp),%r11d
movslq 0x4(%rsp),%rdi
mov    %esi,%eax
movslq %r11d,%rcx
mov    %rcx,0x18(%rsp)
movslq %r10d,%rcx
shl    $0x5,%rcx
mov    %rcx,0x30(%rsp)
xor    %ecx,%ecx
mov    %rdi,0x40(%rsp)
mov    %rax,0x38(%rsp)
jmp    11a1 <MyGemm+0x91>
nop
mov    0x10(%rsp),%rcx
add    $0x4,%rcx
mov    0x30(%rsp),%rax
add    %rax,0x8(%rsp)
mov    0x38(%rsp),%rax
cmp    %rax,%rcx
jge    1252 <MyGemm+0x142>
mov    %rax,%rsi
mov    %rcx,0x10(%rsp)
sub    %rcx,%rsi
cmp    $0x4,%rsi
mov    $0x4,%eax
cmovge %eax,%esi
cmpl   $0x0,0x4(%rsp)
jle    1180 <MyGemm+0x70>
mov    0x10(%rsp),%rax
imul   0x18(%rsp),%rax
mov    0x28(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,0x48(%rsp)
mov    0x20(%rsp),%r15
mov    %rdi,%r12
mov    0x8(%rsp),%rbx
xor    %r13d,%r13d
mov    %rsi,0x50(%rsp)
nopl   (%rax)
cmp    $0x4,%r12
mov    $0x4,%edi
cmovl  %r12d,%edi
sub    $0x8,%rsp
mov    0x58(%rsp),%rsi
mov    %ebp,%edx
mov    %r15,%rcx
mov    %r14d,%r8d
mov    0x50(%rsp),%r9
xor    %eax,%eax
push   %r10
push   %rbx
push   %r11
call   1030 <Gemm_P_Ger@plt>
mov    0x60(%rsp),%rdi
mov    0xb0(%rsp),%r11d
mov    0xc0(%rsp),%r10d
add    $0x20,%rsp
add    $0x4,%r13
add    $0x20,%rbx
add    $0xfffffffffffffffc,%r12
add    $0x20,%r15
cmp    %rdi,%r13
jl     11f0 <MyGemm+0xe0>
jmp    1180 <MyGemm+0x70>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @MyGemm(i32 %0, i32 %1, i32 %2, double* %3, i32 %4, double* %5, i32 %6, double* %7, i32 %8) {
%10 = icmp sgt i32 %1, 0
br i1 %10, label %11, label %17
11: ; preds = %9
%12 = icmp sgt i32 %0, 0
%13 = sext i32 %0 to i64
%14 = zext i32 %1 to i64
%15 = sext i32 %6 to i64
%16 = sext i32 %8 to i64
br label %18
17: ; preds = %28, %9
ret void
18: ; preds = %11, %28
%19 = phi i64 [ 0, %11 ], [ %29, %28 ]
%20 = sub nsw i64 %14, %19
%21 = icmp slt i64 %20, 4
%22 = select i1 %21, i64 %20, i64 4
%23 = trunc i64 %22 to i32
br i1 %12, label %24, label %28
24: ; preds = %18
%25 = mul nsw i64 %19, %15
%26 = getelementptr inbounds double, double* %5, i64 %25
%27 = mul nsw i64 %19, %16
br label %31
28: ; preds = %31, %18
%29 = add nuw i64 %19, 4
%30 = icmp slt i64 %29, %14
br i1 %30, label %18, label %17
31: ; preds = %24, %31
%32 = phi i64 [ 0, %24 ], [ %41, %31 ]
%33 = sub nsw i64 %13, %32
%34 = icmp slt i64 %33, 4
%35 = select i1 %34, i64 %33, i64 4
%36 = trunc i64 %35 to i32
%37 = getelementptr inbounds double, double* %3, i64 %32
%38 = add nsw i64 %32, %27
%39 = getelementptr inbounds double, double* %7, i64 %38
%40 = tail call i32 (i32, i32, i32, double*, i32, double*, i32, double*, i32, ...) bitcast (i32 (...)* @Gemm_P_Ger to i32 (i32, i32, i32, double*, i32, double*, i32, double*, i32, ...)*)(i32 %36, i32 %23, i32 %2, double* %37, i32 %4, double* %26, i32 %6, double* %39, i32 %8)
%41 = add nuw i64 %32, 4
%42 = icmp slt i64 %41, %13
br i1 %42, label %31, label %28
}