#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 */
double im2col_get_pixel(double const *, int const, int const, int const,
                        int const, int const, int const);
void im2col(const double *data_im, const int channels, const int height,
            const int width, const int ksize, const int stride, const int pad,
            double *data_col) {
  const int height_col = (height + 2 * pad - ksize) / stride + 1;
  const int width_col = (width + 2 * pad - ksize) / stride + 1;
  const int channels_col = channels * ksize * ksize;
  for (int c = 0; c < channels_col; ++c) {
    const int w_offset = c % ksize;
    const int h_offset = (c / ksize) % ksize;
    const int c_im = c / ksize / ksize;
    for (int h = 0; h < height_col; ++h) {
      for (int w = 0; w < width_col; ++w) {
        const int im_row = h_offset + h * stride;
        const int im_col = w_offset + w * stride;
        const int col_index = (c * height_col + h) * width_col + w;
        data_col[col_index] =
            im2col_get_pixel(data_im, height, width, im_row, im_col, c_im, pad);
      }
    }
  }
}


x86 assembly (clang, -O0):
<im2col>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
mov    -0x10(%rbp),%eax
mov    0x10(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
sub    -0x18(%rbp),%eax
cltd
idivl  -0x1c(%rbp)
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    0x10(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
sub    -0x18(%rbp),%eax
cltd
idivl  -0x1c(%rbp)
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0xc(%rbp),%eax
imul   -0x18(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1264 <im2col+0x154>
mov    -0x2c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
mov    %edx,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
cltd
idivl  -0x18(%rbp)
mov    %edx,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
cltd
idivl  -0x18(%rbp)
cltd
idivl  -0x18(%rbp)
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1251 <im2col+0x141>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    123e <im2col+0x12e>
mov    -0x34(%rbp),%eax
mov    -0x3c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x44(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x40(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
imul   -0x20(%rbp),%eax
add    -0x3c(%rbp),%eax
imul   -0x24(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%esi
mov    -0x14(%rbp),%edx
mov    -0x44(%rbp),%ecx
mov    -0x48(%rbp),%r8d
mov    -0x38(%rbp),%r9d
mov    0x10(%rbp),%eax
mov    %eax,(%rsp)
call   1030 <im2col_get_pixel@plt>
mov    0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11c5 <im2col+0xb5>
jmp    1243 <im2col+0x133>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11b2 <im2col+0xa2>
jmp    1256 <im2col+0x146>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1179 <im2col+0x69>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<im2col>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %edx,%r14d
mov    0xb0(%rsp),%r11d
mov    %edx,%eax
sub    %r8d,%eax
lea    (%rax,%r11,2),%eax
cltd
idiv   %r9d
mov    %rax,0x28(%rsp)
mov    %ecx,%eax
sub    %r8d,%eax
lea    (%rax,%r11,2),%eax
cltd
idiv   %r9d
mov    %rax,0x38(%rsp)
imul   %r8d,%esi
mov    %r8d,0x14(%rsp)
imul   %r8d,%esi
test   %esi,%esi
jle    12a8 <im2col+0x198>
mov    %ecx,%ebp
mov    %rdi,%r13
mov    0xb8(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    0x28(%rsp),%rax
add    $0x1,%eax
mov    0x38(%rsp),%rcx
add    $0x1,%ecx
movslq %ecx,%rcx
cltq
mov    %esi,%edx
mov    %rdx,0x40(%rsp)
mov    %eax,%edx
mov    %rdx,0x58(%rsp)
mov    %ecx,%esi
imul   %rcx,%rax
shl    $0x3,%rax
mov    %rax,0x48(%rsp)
shl    $0x3,%rcx
mov    %rcx,0x60(%rsp)
xor    %ecx,%ecx
mov    %r9d,0x20(%rsp)
mov    %rsi,0x68(%rsp)
jmp    11de <im2col+0xce>
nopl   (%rax)
mov    0x50(%rsp),%rcx
add    $0x1,%rcx
mov    0x48(%rsp),%rax
add    %rax,0x30(%rsp)
cmp    0x40(%rsp),%rcx
je     12a8 <im2col+0x198>
mov    %rcx,0x50(%rsp)
mov    %ecx,%eax
cltd
mov    0x14(%rsp),%ecx
idiv   %ecx
mov    %edx,0x1c(%rsp)
cltd
idiv   %ecx
mov    %edx,0x18(%rsp)
cmpl   $0x0,0x28(%rsp)
js     11c0 <im2col+0xb0>
mov    %eax,%r12d
mov    0x30(%rsp),%rcx
xor    %edx,%edx
mov    %rdx,0x8(%rsp)
jmp    122d <im2col+0x11d>
nop
mov    0x8(%rsp),%rbx
add    $0x1,%rbx
add    0x60(%rsp),%rcx
mov    %rbx,%rdx
mov    %rbx,0x8(%rsp)
cmp    0x58(%rsp),%rbx
je     11c0 <im2col+0xb0>
cmpl   $0x0,0x38(%rsp)
js     1210 <im2col+0x100>
mov    0x8(%rsp),%rdx
imul   %r9d,%edx
add    0x18(%rsp),%edx
mov    %edx,0x24(%rsp)
mov    0x1c(%rsp),%r15d
xor    %ebx,%ebx
mov    %rcx,0x70(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r11d,(%rsp)
mov    %r13,%rdi
mov    %r14d,%esi
mov    %ebp,%edx
mov    0x24(%rsp),%ecx
mov    %r15d,%r8d
mov    %r12d,%r9d
call   1030 <im2col_get_pixel@plt>
mov    0x70(%rsp),%rcx
mov    0x68(%rsp),%rsi
mov    0xb0(%rsp),%r11d
mov    0x20(%rsp),%r9d
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x1,%rbx
add    %r9d,%r15d
cmp    %rbx,%rsi
jne    1260 <im2col+0x150>
jmp    1210 <im2col+0x100>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @im2col(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, double* %7) {
%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 i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
store double* %0, double** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store double* %7, double** %16
%29 = load i32, i32* %11
%30 = load i32, i32* %15
%31 = mul nsw i32 2, %30
%32 = add nsw i32 %29, %31
%33 = load i32, i32* %13
%34 = sub nsw i32 %32, %33
%35 = load i32, i32* %14
%36 = sdiv i32 %34, %35
%37 = add nsw i32 %36, 1
store i32 %37, i32* %17
%38 = load i32, i32* %12
%39 = load i32, i32* %15
%40 = mul nsw i32 2, %39
%41 = add nsw i32 %38, %40
%42 = load i32, i32* %13
%43 = sub nsw i32 %41, %42
%44 = load i32, i32* %14
%45 = sdiv i32 %43, %44
%46 = add nsw i32 %45, 1
store i32 %46, i32* %18
%47 = load i32, i32* %10
%48 = load i32, i32* %13
%49 = mul nsw i32 %47, %48
%50 = load i32, i32* %13
%51 = mul nsw i32 %49, %50
store i32 %51, i32* %19
store i32 0, i32* %20
br label %52
52: ; preds = %119, %8
%53 = load i32, i32* %20
%54 = load i32, i32* %19
%55 = icmp slt i32 %53, %54
br i1 %55, label %56, label %122
56: ; preds = %52
%57 = load i32, i32* %20
%58 = load i32, i32* %13
%59 = srem i32 %57, %58
store i32 %59, i32* %21
%60 = load i32, i32* %20
%61 = load i32, i32* %13
%62 = sdiv i32 %60, %61
%63 = load i32, i32* %13
%64 = srem i32 %62, %63
store i32 %64, i32* %22
%65 = load i32, i32* %20
%66 = load i32, i32* %13
%67 = sdiv i32 %65, %66
%68 = load i32, i32* %13
%69 = sdiv i32 %67, %68
store i32 %69, i32* %23
store i32 0, i32* %24
br label %70
70: ; preds = %115, %56
%71 = load i32, i32* %24
%72 = load i32, i32* %17
%73 = icmp slt i32 %71, %72
br i1 %73, label %74, label %118
74: ; preds = %70
store i32 0, i32* %25
br label %75
75: ; preds = %111, %74
%76 = load i32, i32* %25
%77 = load i32, i32* %18
%78 = icmp slt i32 %76, %77
br i1 %78, label %79, label %114
79: ; preds = %75
%80 = load i32, i32* %22
%81 = load i32, i32* %24
%82 = load i32, i32* %14
%83 = mul nsw i32 %81, %82
%84 = add nsw i32 %80, %83
store i32 %84, i32* %26
%85 = load i32, i32* %21
%86 = load i32, i32* %25
%87 = load i32, i32* %14
%88 = mul nsw i32 %86, %87
%89 = add nsw i32 %85, %88
store i32 %89, i32* %27
%90 = load i32, i32* %20
%91 = load i32, i32* %17
%92 = mul nsw i32 %90, %91
%93 = load i32, i32* %24
%94 = add nsw i32 %92, %93
%95 = load i32, i32* %18
%96 = mul nsw i32 %94, %95
%97 = load i32, i32* %25
%98 = add nsw i32 %96, %97
store i32 %98, i32* %28
%99 = load double*, double** %9
%100 = load i32, i32* %11
%101 = load i32, i32* %12
%102 = load i32, i32* %26
%103 = load i32, i32* %27
%104 = load i32, i32* %23
%105 = load i32, i32* %15
%106 = call double @im2col_get_pixel(double* %99, i32 %100, i32 %101, i32 %102, i32 %103, i32 %104, i32 %105)
%107 = load double*, double** %16
%108 = load i32, i32* %28
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %107, i64 %109
store double %106, double* %110
br label %111
111: ; preds = %79
%112 = load i32, i32* %25
%113 = add nsw i32 %112, 1
store i32 %113, i32* %25
br label %75
114: ; preds = %75
br label %115
115: ; preds = %114
%116 = load i32, i32* %24
%117 = add nsw i32 %116, 1
store i32 %117, i32* %24
br label %70
118: ; preds = %70
br label %119
119: ; preds = %118
%120 = load i32, i32* %20
%121 = add nsw i32 %120, 1
store i32 %121, i32* %20
br label %52
122: ; preds = %52
ret void
}

LLVM IR (-O2):
define void @im2col(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, double* nocapture writeonly %7) {
%9 = shl nsw i32 %6, 1
%10 = sub i32 %2, %4
%11 = add i32 %10, %9
%12 = sdiv i32 %11, %5
%13 = add i32 %12, 1
%14 = sub i32 %3, %4
%15 = add i32 %14, %9
%16 = sdiv i32 %15, %5
%17 = add i32 %16, 1
%18 = mul nsw i32 %4, %1
%19 = mul i32 %18, %4
%20 = icmp sgt i32 %19, 0
br i1 %20, label %21, label %29
21: ; preds = %8
%22 = icmp slt i32 %12, 0
%23 = icmp slt i32 %16, 0
%24 = sext i32 %17 to i64
%25 = sext i32 %13 to i64
%26 = zext i32 %19 to i64
%27 = zext i32 %13 to i64
%28 = zext i32 %17 to i64
br label %30
29: ; preds = %47, %8
ret void
30: ; preds = %21, %47
%31 = phi i64 [ 0, %21 ], [ %48, %47 ]
%32 = trunc i64 %31 to i32
%33 = sdiv i32 %32, %4
%34 = srem i32 %32, %4
%35 = sdiv i32 %33, %4
%36 = srem i32 %33, %4
br i1 %22, label %47, label %37
37: ; preds = %30
%38 = mul nsw i64 %31, %25
br label %39
39: ; preds = %37, %50
%40 = phi i64 [ 0, %37 ], [ %51, %50 ]
br i1 %23, label %50, label %41
41: ; preds = %39
%42 = add nsw i64 %40, %38
%43 = mul nsw i64 %42, %24
%44 = trunc i64 %40 to i32
%45 = mul i32 %44, %5
%46 = add i32 %45, %36
br label %53
47: ; preds = %50, %30
%48 = add nuw nsw i64 %31, 1
%49 = icmp eq i64 %48, %26
br i1 %49, label %29, label %30
50: ; preds = %53, %39
%51 = add nuw nsw i64 %40, 1
%52 = icmp eq i64 %51, %27
br i1 %52, label %47, label %39
53: ; preds = %41, %53
%54 = phi i64 [ 0, %41 ], [ %61, %53 ]
%55 = add nsw i64 %54, %43
%56 = trunc i64 %54 to i32
%57 = mul i32 %56, %5
%58 = add i32 %57, %34
%59 = tail call double @im2col_get_pixel(double* %0, i32 %2, i32 %3, i32 %46, i32 %58, i32 %35, i32 %6)
%60 = getelementptr inbounds double, double* %7, i64 %55
store double %59, double* %60
%61 = add nuw nsw i64 %54, 1
%62 = icmp eq i64 %61, %28
br i1 %62, label %50, label %53
}