void ExtractImage(double *src, unsigned char *dst, int align_rows,
                  int align_cols, int rows, int cols) {
  int pad_top, pad_bottom, pad_left, pad_right;
  unsigned char *pd;
  double *ps;
  int i, j;
  pad_top = (align_rows - rows) >> 1;
  pad_bottom = align_rows - rows - pad_top;
  pad_left = (align_cols - cols) >> 1;
  pad_right = align_cols - cols - pad_left;
  ps = src + pad_top * align_cols + pad_left;
  pd = dst;
  for (i = 0; i < rows; i++) {
    for (j = 0; j < cols; j++)
      *pd++ = (unsigned char)*ps++;
    ps += pad_right + pad_left;
  }
}


x86 assembly (clang, -O0):
<ExtractImage>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x1c(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x1c(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11f6 <ExtractImage+0xf6>
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11d0 <ExtractImage+0xd0>
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
cvttsd2si (%rax),%eax
mov    %al,%cl
mov    -0x38(%rbp),%rax
mov    %rax,%rdx
add    $0x1,%rdx
mov    %rdx,-0x38(%rbp)
mov    %cl,(%rax)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1190 <ExtractImage+0x90>
mov    -0x30(%rbp),%ecx
add    -0x2c(%rbp),%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    117d <ExtractImage+0x7d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ExtractImage>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %r8d,%r8d
jle    1321 <ExtractImage+0x221>
sub    %r8d,%edx
sar    %edx
imul   %ecx,%edx
sub    %r9d,%ecx
movslq %edx,%rax
lea    (%rdi,%rax,8),%rdx
movslq %ecx,%rax
sar    %ecx
movslq %ecx,%rcx
lea    (%rdx,%rcx,8),%rcx
lea    -0x1(%r9),%r10d
lea    0x1(%r10),%r13
mov    %r13,%r12
and    $0xfffffffffffffffc,%r12
lea    -0x4(%r12),%r15
mov    %r15,-0x8(%rsp)
shr    $0x2,%r15
add    $0x1,%r15
mov    %r15,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x10(%rsp)
xor    %ebx,%ebx
jmp    1180 <ExtractImage+0x80>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%rcx,%rax,8),%rcx
add    $0x1,%ebx
cmp    %r8d,%ebx
je     1321 <ExtractImage+0x221>
test   %r9d,%r9d
jle    1170 <ExtractImage+0x70>
xor    %edx,%edx
cmp    $0x3,%r10d
jb     1290 <ExtractImage+0x190>
lea    (%rcx,%r13,8),%rdi
cmp    %rdi,%rsi
jae    11a7 <ExtractImage+0xa7>
lea    (%rsi,%r13,1),%rdi
cmp    %rcx,%rdi
ja     1290 <ExtractImage+0x190>
cmpq   $0x0,-0x8(%rsp)
je     1310 <ExtractImage+0x210>
mov    -0x10(%rsp),%rdx
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movupd (%rcx,%rbp,8),%xmm0
movupd 0x10(%rcx,%rbp,8),%xmm1
cvttpd2dq %xmm0,%xmm0
packuswb %xmm0,%xmm0
packuswb %xmm0,%xmm0
cvttpd2dq %xmm1,%xmm1
packuswb %xmm1,%xmm1
packuswb %xmm1,%xmm1
movd   %xmm0,%edi
mov    %di,(%rsi,%rbp,1)
movd   %xmm1,%edi
mov    %di,0x2(%rsi,%rbp,1)
movupd 0x20(%rcx,%rbp,8),%xmm0
movupd 0x30(%rcx,%rbp,8),%xmm1
cvttpd2dq %xmm0,%xmm0
packuswb %xmm0,%xmm0
packuswb %xmm0,%xmm0
cvttpd2dq %xmm1,%xmm1
packuswb %xmm1,%xmm1
packuswb %xmm1,%xmm1
movd   %xmm0,%edi
mov    %di,0x4(%rsi,%rbp,1)
movd   %xmm1,%edi
mov    %di,0x6(%rsi,%rbp,1)
add    $0x8,%rbp
add    $0xfffffffffffffffe,%rdx
jne    11c0 <ExtractImage+0xc0>
test   $0x1,%r15b
je     126e <ExtractImage+0x16e>
movupd (%rcx,%rbp,8),%xmm0
movupd 0x10(%rcx,%rbp,8),%xmm1
cvttpd2dq %xmm0,%xmm0
packuswb %xmm0,%xmm0
packuswb %xmm0,%xmm0
cvttpd2dq %xmm1,%xmm1
packuswb %xmm1,%xmm1
packuswb %xmm1,%xmm1
movd   %xmm0,%edx
mov    %dx,(%rsi,%rbp,1)
movd   %xmm1,%edx
mov    %dx,0x2(%rsi,%rbp,1)
add    %r12,%rsi
lea    (%rcx,%r12,8),%rcx
mov    %r12d,%edx
cmp    %r12,%r13
je     1170 <ExtractImage+0x70>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r9d,%edi
sub    %edx,%edi
mov    %edx,%ebp
not    %ebp
add    %r9d,%ebp
and    $0x3,%edi
je     12cc <ExtractImage+0x1cc>
neg    %edi
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
cvttsd2si (%rcx),%r11d
add    $0x8,%rcx
mov    %r11b,(%rsi)
add    $0x1,%rsi
add    $0xffffffffffffffff,%r14
cmp    %r14d,%edi
jne    12b0 <ExtractImage+0x1b0>
sub    %r14d,%edx
cmp    $0x3,%ebp
jb     1170 <ExtractImage+0x70>
mov    %r9d,%edi
sub    %edx,%edi
nopw   0x0(%rax,%rax,1)
cvttsd2si (%rcx),%edx
mov    %dl,(%rsi)
cvttsd2si 0x8(%rcx),%edx
mov    %dl,0x1(%rsi)
cvttsd2si 0x10(%rcx),%edx
mov    %dl,0x2(%rsi)
cvttsd2si 0x18(%rcx),%edx
add    $0x20,%rcx
mov    %dl,0x3(%rsi)
add    $0x4,%rsi
add    $0xfffffffc,%edi
jne    12e0 <ExtractImage+0x1e0>
jmp    1170 <ExtractImage+0x70>
xor    %ebp,%ebp
test   $0x1,%r15b
jne    123a <ExtractImage+0x13a>
jmp    126e <ExtractImage+0x16e>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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