void d_rep_mat(int reps, int row, int col, double *A, int lda, double *B,
               int ldb) {
  int i, j, l;
  for (l = 0; l < reps; l++) {
    for (j = 0; j < col; j++) {
      for (i = 0; i < row; i++) {
        B[i + j * ldb + l * row * col] = A[i + j * lda];
      }
    }
  }
}


x86 assembly (clang, -O0):
<d_rep_mat>:
push   %rbp
mov    %rsp,%rbp
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,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11c6 <d_rep_mat+0xc6>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b3 <d_rep_mat+0xb3>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11a0 <d_rep_mat+0xa0>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   -0x8(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1149 <d_rep_mat+0x49>
jmp    11a5 <d_rep_mat+0xa5>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1136 <d_rep_mat+0x36>
jmp    11b8 <d_rep_mat+0xb8>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1123 <d_rep_mat+0x23>
pop    %rbp
ret

x86 assembly (clang, -O2):
<d_rep_mat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %rcx,-0x70(%rsp)
mov    %edx,-0x7c(%rsp)
mov    %esi,-0x80(%rsp)
mov    %edi,-0x50(%rsp)
test   %edi,%edi
jle    146d <d_rep_mat+0x36d>
mov    -0x7c(%rsp),%eax
mov    %eax,%edx
mov    -0x80(%rsp),%ecx
imul   %ecx,%edx
mov    %rdx,-0x48(%rsp)
movslq %r8d,%rdx
mov    %eax,%eax
mov    %rax,0x18(%rsp)
mov    %ecx,%eax
lea    -0x1(%rax),%rcx
mov    %rcx,(%rsp)
lea    (%r9,%rax,8),%rcx
mov    %rcx,-0x8(%rsp)
mov    %eax,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,-0x10(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,-0x18(%rsp)
mov    %rcx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %rsi,%rcx
mov    %rsi,-0x20(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x28(%rsp)
mov    %rax,-0x58(%rsp)
neg    %rax
mov    %rax,-0x38(%rsp)
mov    -0x70(%rsp),%rax
add    $0x30,%rax
mov    %rax,-0x30(%rsp)
mov    %rdx,0x8(%rsp)
lea    0x0(,%rdx,8),%rax
mov    %rax,0x10(%rsp)
movl   $0x4,-0x74(%rsp)
xor    %eax,%eax
mov    %rax,-0x60(%rsp)
movl   $0x1,-0x78(%rsp)
xor    %eax,%eax
mov    %rax,-0x68(%rsp)
jmp    1211 <d_rep_mat+0x111>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x68(%rsp),%rdx
add    $0x1,%rdx
mov    -0x48(%rsp),%rax
add    %eax,-0x74(%rsp)
mov    -0x60(%rsp),%rcx
add    %eax,%ecx
mov    %rcx,-0x60(%rsp)
add    %eax,-0x78(%rsp)
mov    %rdx,-0x68(%rsp)
cmp    -0x50(%rsp),%edx
je     146d <d_rep_mat+0x36d>
cmpl   $0x0,-0x7c(%rsp)
jle    11e0 <d_rep_mat+0xe0>
mov    -0x68(%rsp),%rcx
mov    %rcx,%rdx
mov    -0x48(%rsp),%rax
imul   %rax,%rdx
mov    %rdx,-0x40(%rsp)
imul   %ecx,%eax
mov    %eax,-0x4c(%rsp)
mov    -0x78(%rsp),%r11d
mov    -0x60(%rsp),%rbp
mov    -0x70(%rsp),%rcx
mov    -0x30(%rsp),%r14
mov    -0x74(%rsp),%r8d
xor    %r15d,%r15d
jmp    1286 <d_rep_mat+0x186>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r15
mov    0x58(%rsp),%eax
add    %eax,%r8d
mov    0x10(%rsp),%rdx
add    %rdx,%r14
add    %rdx,%rcx
add    %eax,%ebp
add    %eax,%r11d
cmp    0x18(%rsp),%r15
je     11e0 <d_rep_mat+0xe0>
cmpl   $0x0,-0x80(%rsp)
jle    1260 <d_rep_mat+0x160>
mov    %r15,%r13
imul   0x8(%rsp),%r13
mov    %r15d,%r10d
imul   0x58(%rsp),%r10d
add    -0x4c(%rsp),%r10d
cmpl   $0x4,-0x80(%rsp)
jb     1320 <d_rep_mat+0x220>
mov    %r15d,%eax
mov    0x58(%rsp),%edx
imul   %edx,%eax
add    -0x40(%rsp),%eax
mov    (%rsp),%rdx
lea    (%rax,%rdx,1),%esi
cmp    %eax,%esi
jl     1320 <d_rep_mat+0x220>
mov    %rdx,%rax
shr    $0x20,%rax
mov    $0x0,%eax
jne    1322 <d_rep_mat+0x222>
mov    %r15d,%eax
mov    0x58(%rsp),%edx
imul   %edx,%eax
add    -0x40(%rsp),%eax
mov    -0x70(%rsp),%rdx
lea    (%rdx,%r13,8),%rsi
cltq
mov    -0x8(%rsp),%rdi
lea    (%rdi,%rax,8),%rdi
cmp    %rdi,%rsi
jae    13bb <d_rep_mat+0x2bb>
mov    -0x58(%rsp),%rsi
add    %r13,%rsi
lea    (%rdx,%rsi,8),%rsi
lea    (%r9,%rax,8),%rax
cmp    %rsi,%rax
jae    13bb <d_rep_mat+0x2bb>
cs nopw 0x0(%rax,%rax,1)
xor    %eax,%eax
mov    %rax,%r12
testb  $0x1,-0x58(%rsp)
je     134d <d_rep_mat+0x24d>
add    %rax,%r13
mov    -0x70(%rsp),%rdx
movsd  (%rdx,%r13,8),%xmm0
add    %eax,%r10d
movslq %r10d,%rsi
movsd  %xmm0,(%r9,%rsi,8)
mov    %rax,%r12
or     $0x1,%r12
not    %rax
cmp    -0x38(%rsp),%rax
je     1260 <d_rep_mat+0x160>
mov    -0x38(%rsp),%rax
lea    (%rax,%r12,1),%rbx
lea    (%rcx,%r12,8),%rax
lea    (%r12,%rbp,1),%edi
add    %r11d,%r12d
mov    $0x1,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  -0x8(%rax,%rsi,8),%xmm0
movslq %edi,%rdi
movsd  %xmm0,(%r9,%rdi,8)
movsd  (%rax,%rsi,8),%xmm0
movslq %r12d,%r12
movsd  %xmm0,(%r9,%r12,8)
lea    (%rbx,%rsi,1),%rdx
add    $0x2,%rdx
add    $0x2,%rsi
add    $0x2,%edi
add    $0x2,%r12d
cmp    $0x1,%rdx
jne    1380 <d_rep_mat+0x280>
jmp    1260 <d_rep_mat+0x160>
cmpq   $0x0,-0x18(%rsp)
je     1462 <d_rep_mat+0x362>
mov    -0x28(%rsp),%rax
mov    %r8d,%edi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups -0x30(%r14,%rbx,8),%xmm0
movups -0x20(%r14,%rbx,8),%xmm1
lea    -0x4(%rdi),%esi
movslq %esi,%rsi
movups %xmm0,(%r9,%rsi,8)
movups %xmm1,0x10(%r9,%rsi,8)
movups -0x10(%r14,%rbx,8),%xmm0
movups (%r14,%rbx,8),%xmm1
movslq %edi,%rdi
movups %xmm0,(%r9,%rdi,8)
movups %xmm1,0x10(%r9,%rdi,8)
add    $0x8,%rbx
add    $0x8,%edi
add    $0xfffffffffffffffe,%rax
jne    13e0 <d_rep_mat+0x2e0>
testb  $0x1,-0x20(%rsp)
je     144a <d_rep_mat+0x34a>
mov    %rbx,%rax
add    %r13,%rax
movups (%rdx,%rax,8),%xmm0
movups 0x10(%rdx,%rax,8),%xmm1
lea    (%r10,%rbx,1),%eax
cltq
movups %xmm0,(%r9,%rax,8)
movups %xmm1,0x10(%r9,%rax,8)
mov    -0x10(%rsp),%rdx
mov    %rdx,%rax
cmp    -0x58(%rsp),%rdx
je     1260 <d_rep_mat+0x160>
jmp    1322 <d_rep_mat+0x222>
xor    %ebx,%ebx
testb  $0x1,-0x20(%rsp)
jne    142a <d_rep_mat+0x32a>
jmp    144a <d_rep_mat+0x34a>
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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