#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int atoi(char *);
double dgemm_gflops(int, int, int);
int omp_get_max_threads();
int /*<<< orphan*/ printf(char *, int, ...);
int main(int argc, char *argv[]) {
  int maxm = (argc > 1) ? atoi(argv[1]) : 1000;
  int incm = (argc > 2) ? atoi(argv[2]) : 100;
  int maxk = (argc > 3) ? atoi(argv[3]) : 1000;
  int inck = (argc > 4) ? atoi(argv[4]) : 100;
  printf("testing DGEMM %d threads \n", omp_get_max_threads());
  printf("maxm=%d, incm=%d, maxk=%d, inck=%d\n", maxm, incm, maxk, inck);
  for (int m = incm; m <= maxm; m += incm) {
    int n = m;
    for (int k = inck; k <= maxk; k += inck) {
      double eff_peak = dgemm_gflops(m, n, k);
      printf("m=n=%5d, k=%5d = %lf GF/s \n", m, k, eff_peak);
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
cmpl   $0x1,-0x8(%rbp)
jle    1175 <main+0x35>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rdi
call   1060 <atoi@plt>
mov    %eax,-0x3c(%rbp)
jmp    1182 <main+0x42>
mov    $0x3e8,%eax
mov    %eax,-0x3c(%rbp)
jmp    1182 <main+0x42>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x2,-0x8(%rbp)
jle    11a7 <main+0x67>
mov    -0x10(%rbp),%rax
mov    0x10(%rax),%rdi
call   1060 <atoi@plt>
mov    %eax,-0x40(%rbp)
jmp    11b4 <main+0x74>
mov    $0x64,%eax
mov    %eax,-0x40(%rbp)
jmp    11b4 <main+0x74>
mov    -0x40(%rbp),%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x3,-0x8(%rbp)
jle    11d9 <main+0x99>
mov    -0x10(%rbp),%rax
mov    0x18(%rax),%rdi
call   1060 <atoi@plt>
mov    %eax,-0x44(%rbp)
jmp    11e6 <main+0xa6>
mov    $0x3e8,%eax
mov    %eax,-0x44(%rbp)
jmp    11e6 <main+0xa6>
mov    -0x44(%rbp),%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x4,-0x8(%rbp)
jle    120b <main+0xcb>
mov    -0x10(%rbp),%rax
mov    0x20(%rax),%rdi
call   1060 <atoi@plt>
mov    %eax,-0x48(%rbp)
jmp    1218 <main+0xd8>
mov    $0x64,%eax
mov    %eax,-0x48(%rbp)
jmp    1218 <main+0xd8>
mov    -0x48(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    $0x0,%al
call   1050 <omp_get_max_threads@plt>
mov    %eax,%esi
lea    0xdd2(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x14(%rbp),%esi
mov    -0x18(%rbp),%edx
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
lea    0xdd2(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x18(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     12c7 <main+0x187>
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     12b4 <main+0x174>
mov    -0x24(%rbp),%edi
mov    -0x28(%rbp),%esi
mov    -0x2c(%rbp),%edx
call   1030 <dgemm_gflops@plt>
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%esi
mov    -0x2c(%rbp),%edx
movsd  -0x38(%rbp),%xmm0
lea    0xda0(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    -0x20(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    126e <main+0x12e>
jmp    12b9 <main+0x179>
mov    -0x18(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    1256 <main+0x116>
xor    %eax,%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    $0x64,%r14d
mov    $0x3e8,%r12d
cmp    $0x2,%edi
jl     11a1 <main+0x61>
mov    %rsi,%r15
mov    %edi,%ebx
mov    0x8(%rsi),%rdi
call   1060 <atoi@plt>
mov    %eax,%r13d
cmp    $0x2,%ebx
je     11a7 <main+0x67>
mov    0x10(%r15),%rdi
call   1060 <atoi@plt>
mov    %eax,%ebp
cmp    $0x4,%ebx
jb     11ac <main+0x6c>
mov    0x18(%r15),%rdi
call   1060 <atoi@plt>
mov    %eax,%r12d
cmp    $0x4,%ebx
je     11ac <main+0x6c>
mov    0x20(%r15),%rdi
call   1060 <atoi@plt>
mov    %eax,%r14d
jmp    11ac <main+0x6c>
mov    $0x3e8,%r13d
mov    $0x64,%ebp
xor    %eax,%eax
call   1050 <omp_get_max_threads@plt>
lea    0xe46(%rip),%rdi
mov    %eax,%esi
xor    %eax,%eax
call   1040 <printf@plt>
lea    0xe51(%rip),%rdi
mov    %r13d,%esi
mov    %ebp,%edx
mov    %r12d,%ecx
mov    %r14d,%r8d
xor    %eax,%eax
call   1040 <printf@plt>
mov    %ebp,0x4(%rsp)
cmp    %r13d,%ebp
jle    11f6 <main+0xb6>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0xe42(%rip),%rbp
mov    0x4(%rsp),%ebx
jmp    121c <main+0xdc>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    0x4(%rsp),%ebx
mov    %r15d,%r13d
cmp    %r15d,%ebx
jg     11e5 <main+0xa5>
mov    %r13d,%r15d
mov    %r14d,%r13d
cmp    %r12d,%r14d
jg     1210 <main+0xd0>
nopw   0x0(%rax,%rax,1)
mov    %ebx,%edi
mov    %ebx,%esi
mov    %r13d,%edx
call   1030 <dgemm_gflops@plt>
mov    %rbp,%rdi
mov    %ebx,%esi
mov    %r13d,%edx
mov    $0x1,%al
call   1040 <printf@plt>
add    %r14d,%r13d
cmp    %r12d,%r13d
jle    1230 <main+0xf0>
jmp    1210 <main+0xd0>

LLVM IR (-O0):
@.str = [27 x i8] c"testing DGEMM %d threads \0A\00"
@.str.1 = [36 x i8] c"maxm=%d, incm=%d, maxk=%d, inck=%d\0A\00"
@.str.2 = [28 x i8] c"m=n=%5d, k=%5d = %lf GF/s \0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
%14 = load i32, i32* %4
%15 = icmp sgt i32 %14, 1
br i1 %15, label %16, label %21
16: ; preds = %2
%17 = load i8**, i8*** %5
%18 = getelementptr inbounds i8*, i8** %17, i64 1
%19 = load i8*, i8** %18
%20 = call i32 @atoi(i8* %19)
br label %22
21: ; preds = %2
br label %22
22: ; preds = %21, %16
%23 = phi i32 [ %20, %16 ], [ 1000, %21 ]
store i32 %23, i32* %6
%24 = load i32, i32* %4
%25 = icmp sgt i32 %24, 2
br i1 %25, label %26, label %31
26: ; preds = %22
%27 = load i8**, i8*** %5
%28 = getelementptr inbounds i8*, i8** %27, i64 2
%29 = load i8*, i8** %28
%30 = call i32 @atoi(i8* %29)
br label %32
31: ; preds = %22
br label %32
32: ; preds = %31, %26
%33 = phi i32 [ %30, %26 ], [ 100, %31 ]
store i32 %33, i32* %7
%34 = load i32, i32* %4
%35 = icmp sgt i32 %34, 3
br i1 %35, label %36, label %41
36: ; preds = %32
%37 = load i8**, i8*** %5
%38 = getelementptr inbounds i8*, i8** %37, i64 3
%39 = load i8*, i8** %38
%40 = call i32 @atoi(i8* %39)
br label %42
41: ; preds = %32
br label %42
42: ; preds = %41, %36
%43 = phi i32 [ %40, %36 ], [ 1000, %41 ]
store i32 %43, i32* %8
%44 = load i32, i32* %4
%45 = icmp sgt i32 %44, 4
br i1 %45, label %46, label %51
46: ; preds = %42
%47 = load i8**, i8*** %5
%48 = getelementptr inbounds i8*, i8** %47, i64 4
%49 = load i8*, i8** %48
%50 = call i32 @atoi(i8* %49)
br label %52
51: ; preds = %42
br label %52
52: ; preds = %51, %46
%53 = phi i32 [ %50, %46 ], [ 100, %51 ]
store i32 %53, i32* %9
%54 = call i32 (...) @omp_get_max_threads()
%55 = call i32 (i8*, i32, ...) @printf(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str, i64 0, i64 0), i32 %54)
%56 = load i32, i32* %6
%57 = load i32, i32* %7
%58 = load i32, i32* %8
%59 = load i32, i32* %9
%60 = call i32 (i8*, i32, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.1, i64 0, i64 0), i32 %56, i32 %57, i32 %58, i32 %59)
%61 = load i32, i32* %7
store i32 %61, i32* %10
br label %62
62: ; preds = %87, %52
%63 = load i32, i32* %10
%64 = load i32, i32* %6
%65 = icmp sle i32 %63, %64
br i1 %65, label %66, label %91
66: ; preds = %62
%67 = load i32, i32* %10
store i32 %67, i32* %11
%68 = load i32, i32* %9
store i32 %68, i32* %12
br label %69
69: ; preds = %82, %66
%70 = load i32, i32* %12
%71 = load i32, i32* %8
%72 = icmp sle i32 %70, %71
br i1 %72, label %73, label %86
73: ; preds = %69
%74 = load i32, i32* %10
%75 = load i32, i32* %11
%76 = load i32, i32* %12
%77 = call double @dgemm_gflops(i32 %74, i32 %75, i32 %76)
store double %77, double* %13
%78 = load i32, i32* %10
%79 = load i32, i32* %12
%80 = load double, double* %13
%81 = call i32 (i8*, i32, ...) @printf(i8* getelementptr inbounds ([28 x i8], [28 x i8]* @.str.2, i64 0, i64 0), i32 %78, i32 %79, double %80)
br label %82
82: ; preds = %73
%83 = load i32, i32* %9
%84 = load i32, i32* %12
%85 = add nsw i32 %84, %83
store i32 %85, i32* %12
br label %69
86: ; preds = %69
br label %87
87: ; preds = %86
%88 = load i32, i32* %7
%89 = load i32, i32* %10
%90 = add nsw i32 %89, %88
store i32 %90, i32* %10
br label %62
91: ; preds = %62
ret i32 0
}

LLVM IR (-O2):
@.str = [27 x i8] c"testing DGEMM %d threads \0A\00"
@.str.1 = [36 x i8] c"maxm=%d, incm=%d, maxk=%d, inck=%d\0A\00"
@.str.2 = [28 x i8] c"m=n=%5d, k=%5d = %lf GF/s \0A\00"
define i32 @main(i32 %0, i8** nocapture %1) {
%3 = icmp sgt i32 %0, 1
br i1 %3, label %4, label %23
4: ; preds = %2
%5 = getelementptr inbounds i8*, i8** %1, i64 1
%6 = load i8*, i8** %5
%7 = tail call i32 @atoi(i8* %6)
%8 = icmp eq i32 %0, 2
br i1 %8, label %23, label %9
9: ; preds = %4
%10 = getelementptr inbounds i8*, i8** %1, i64 2
%11 = load i8*, i8** %10
%12 = tail call i32 @atoi(i8* %11)
%13 = icmp ugt i32 %0, 3
br i1 %13, label %14, label %23
14: ; preds = %9
%15 = getelementptr inbounds i8*, i8** %1, i64 3
%16 = load i8*, i8** %15
%17 = tail call i32 @atoi(i8* %16)
%18 = icmp eq i32 %0, 4
br i1 %18, label %23, label %19
19: ; preds = %14
%20 = getelementptr inbounds i8*, i8** %1, i64 4
%21 = load i8*, i8** %20
%22 = tail call i32 @atoi(i8* %21)
br label %23
23: ; preds = %2, %4, %9, %14, %19
%24 = phi i32 [ %17, %19 ], [ %17, %14 ], [ 1000, %9 ], [ 1000, %4 ], [ 1000, %2 ]
%25 = phi i32 [ %7, %19 ], [ %7, %14 ], [ %7, %9 ], [ %7, %4 ], [ 1000, %2 ]
%26 = phi i32 [ %12, %19 ], [ %12, %14 ], [ %12, %9 ], [ 100, %4 ], [ 100, %2 ]
%27 = phi i32 [ %22, %19 ], [ 100, %14 ], [ 100, %9 ], [ 100, %4 ], [ 100, %2 ]
%28 = tail call i32 (...) @omp_get_max_threads()
%29 = tail call i32 (i8*, i32, ...) @printf(i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str, i64 0, i64 0), i32 %28)
%30 = tail call i32 (i8*, i32, ...) @printf(i8* getelementptr inbounds ([36 x i8], [36 x i8]* @.str.1, i64 0, i64 0), i32 %25, i32 %26, i32 %24, i32 %27)
%31 = icmp sgt i32 %26, %25
br i1 %31, label %36, label %32
32: ; preds = %23
%33 = icmp sgt i32 %27, %24
br label %34
34: ; preds = %32, %37
%35 = phi i32 [ %26, %32 ], [ %38, %37 ]
br i1 %33, label %37, label %40
36: ; preds = %37, %23
ret i32 0
37: ; preds = %40, %34
%38 = add nsw i32 %35, %26
%39 = icmp sgt i32 %38, %25
br i1 %39, label %36, label %34
40: ; preds = %34, %40
%41 = phi i32 [ %44, %40 ], [ %27, %34 ]
%42 = tail call double @dgemm_gflops(i32 %35, i32 %35, i32 %41)
%43 = tail call i32 (i8*, i32, ...) @printf(i8* getelementptr inbounds ([28 x i8], [28 x i8]* @.str.2, i64 0, i64 0), i32 %35, i32 %41, double %42)
%44 = add nsw i32 %41, %27
%45 = icmp sgt i32 %44, %24
br i1 %45, label %37, label %40
}