int cimg2dimg(unsigned char *src, double *dst, int w, int h, int dp,
              double norm) {
  int i, j, k;
  int lw;
  lw = w * dp;
  for (j = 0; j < h; j++) {
    for (i = 0; i < w; i++) {
      for (k = 0; k < dp; k++)
        dst[j * lw + i * dp + k] = norm * src[j * lw + i * dp + k];
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<cimg2dimg>:
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)
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11dc <cimg2dimg+0xdc>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c9 <cimg2dimg+0xc9>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11b6 <cimg2dimg+0xb6>
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x38(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movzbl (%rax,%rcx,1),%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x38(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1152 <cimg2dimg+0x52>
jmp    11bb <cimg2dimg+0xbb>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    113f <cimg2dimg+0x3f>
jmp    11ce <cimg2dimg+0xce>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    112c <cimg2dimg+0x2c>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<cimg2dimg>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %ecx,-0x58(%rsp)
test   %ecx,%ecx
jle    1451 <cimg2dimg+0x351>
mov    %r8d,%r9d
lea    -0x1(%r9),%r10
mov    %r8d,%eax
imul   %edx,%eax
mov    %eax,-0x64(%rsp)
lea    (%rsi,%r9,8),%rax
mov    %rax,-0x18(%rsp)
lea    (%rdi,%r9,1),%r12
mov    %r9d,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x38(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x40(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %rcx,%rax
mov    %rcx,-0x48(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x50(%rsp)
mov    %r9,%rax
neg    %rax
mov    %rax,-0x10(%rsp)
movl   $0x4,-0x68(%rsp)
movl   $0x1,-0x6c(%rsp)
xor    %eax,%eax
mov    %rax,-0x60(%rsp)
pxor   %xmm2,%xmm2
movl   $0x0,-0x70(%rsp)
mov    %r9,-0x20(%rsp)
mov    %r10,-0x28(%rsp)
mov    %r12,-0x30(%rsp)
jmp    11dd <cimg2dimg+0xdd>
nopl   0x0(%rax,%rax,1)
mov    -0x70(%rsp),%ebp
add    $0x1,%ebp
mov    -0x64(%rsp),%eax
add    %eax,-0x68(%rsp)
add    %eax,-0x6c(%rsp)
mov    -0x60(%rsp),%rcx
add    %eax,%ecx
mov    %rcx,-0x60(%rsp)
mov    %ebp,-0x70(%rsp)
cmp    -0x58(%rsp),%ebp
je     1451 <cimg2dimg+0x351>
test   %edx,%edx
jle    11b0 <cimg2dimg+0xb0>
mov    -0x64(%rsp),%ecx
mov    -0x70(%rsp),%eax
imul   %eax,%ecx
mov    %ecx,-0x54(%rsp)
imul   %edx,%eax
mov    %rax,-0x8(%rsp)
xor    %r11d,%r11d
mov    -0x60(%rsp),%r15
mov    -0x6c(%rsp),%eax
mov    %eax,%r14d
mov    -0x68(%rsp),%r13d
jmp    1222 <cimg2dimg+0x122>
xchg   %ax,%ax
add    $0x1,%r11d
add    %r8d,%r13d
add    %r8d,%r14d
add    %r8d,%r15d
cmp    %edx,%r11d
je     11b0 <cimg2dimg+0xb0>
test   %r8d,%r8d
jle    1210 <cimg2dimg+0x110>
mov    -0x8(%rsp),%rax
lea    (%r11,%rax,1),%ebx
imul   %r8d,%ebx
cmp    $0x4,%r8d
jb     1290 <cimg2dimg+0x190>
mov    %r11d,%ecx
imul   %r8d,%ecx
add    -0x54(%rsp),%ecx
lea    (%rcx,%r10,1),%eax
cmp    %ecx,%eax
jl     1290 <cimg2dimg+0x190>
mov    %r10,%rax
shr    $0x20,%rax
mov    $0x0,%eax
jne    1292 <cimg2dimg+0x192>
movslq %ecx,%rax
lea    (%r12,%rax,1),%rcx
lea    (%rsi,%rax,8),%rbp
cmp    %rbp,%rcx
jbe    131f <cimg2dimg+0x21f>
lea    (%rdi,%rax,1),%rcx
mov    -0x18(%rsp),%rbp
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
cmp    %rax,%rcx
jae    131f <cimg2dimg+0x21f>
nopl   0x0(%rax)
xor    %eax,%eax
mov    %rax,%rcx
test   $0x1,%r9b
je     12bb <cimg2dimg+0x1bb>
add    %eax,%ebx
movslq %ebx,%rcx
movzbl (%rdi,%rcx,1),%ebp
xorps  %xmm3,%xmm3
cvtsi2sd %ebp,%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
mov    %rax,%rcx
or     $0x1,%rcx
not    %rax
cmp    -0x10(%rsp),%rax
je     1210 <cimg2dimg+0x110>
lea    (%r14,%rcx,1),%ebp
mov    %r9,%rax
sub    %rcx,%rax
add    %r15d,%ecx
cs nopw 0x0(%rax,%rax,1)
movslq %ecx,%rcx
movzbl (%rdi,%rcx,1),%ebx
xorps  %xmm3,%xmm3
cvtsi2sd %ebx,%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
movslq %ebp,%rbp
movzbl (%rdi,%rbp,1),%ebx
xorps  %xmm3,%xmm3
cvtsi2sd %ebx,%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,(%rsi,%rbp,8)
add    $0x2,%ebp
add    $0x2,%ecx
add    $0xfffffffffffffffe,%rax
jne    12e0 <cimg2dimg+0x1e0>
jmp    1210 <cimg2dimg+0x110>
cmpq   $0x0,-0x40(%rsp)
je     1445 <cimg2dimg+0x345>
mov    -0x50(%rsp),%r10
mov    %r13d,%r9d
xor    %r12d,%r12d
cs nopw 0x0(%rax,%rax,1)
lea    -0x4(%r9),%eax
cltq
movzwl (%rdi,%rax,1),%ecx
movd   %ecx,%xmm3
movzwl 0x2(%rdi,%rax,1),%ecx
movd   %ecx,%xmm4
punpcklbw %xmm2,%xmm3
punpcklwd %xmm2,%xmm3
cvtdq2pd %xmm3,%xmm3
punpcklbw %xmm2,%xmm4
punpcklwd %xmm2,%xmm4
cvtdq2pd %xmm4,%xmm4
mulpd  %xmm1,%xmm3
mulpd  %xmm1,%xmm4
movupd %xmm3,(%rsi,%rax,8)
movupd %xmm4,0x10(%rsi,%rax,8)
movslq %r9d,%r9
movzwl (%rdi,%r9,1),%eax
movd   %eax,%xmm3
movzwl 0x2(%rdi,%r9,1),%eax
movd   %eax,%xmm4
punpcklbw %xmm2,%xmm3
punpcklwd %xmm2,%xmm3
cvtdq2pd %xmm3,%xmm3
punpcklbw %xmm2,%xmm4
punpcklwd %xmm2,%xmm4
cvtdq2pd %xmm4,%xmm4
mulpd  %xmm1,%xmm3
mulpd  %xmm1,%xmm4
movupd %xmm3,(%rsi,%r9,8)
movupd %xmm4,0x10(%rsi,%r9,8)
add    $0x8,%r12
add    $0x8,%r9d
add    $0xfffffffffffffffe,%r10
jne    1340 <cimg2dimg+0x240>
testb  $0x1,-0x48(%rsp)
je     1420 <cimg2dimg+0x320>
add    %ebx,%r12d
movslq %r12d,%rax
movzwl (%rdi,%rax,1),%ecx
movd   %ecx,%xmm3
movzwl 0x2(%rdi,%rax,1),%ecx
movd   %ecx,%xmm4
punpcklbw %xmm2,%xmm3
punpcklwd %xmm2,%xmm3
cvtdq2pd %xmm3,%xmm3
punpcklbw %xmm2,%xmm4
punpcklwd %xmm2,%xmm4
cvtdq2pd %xmm4,%xmm4
mulpd  %xmm1,%xmm3
mulpd  %xmm1,%xmm4
movupd %xmm3,(%rsi,%rax,8)
movupd %xmm4,0x10(%rsi,%rax,8)
mov    -0x38(%rsp),%rcx
mov    %rcx,%rax
mov    -0x20(%rsp),%r9
cmp    %r9,%rcx
mov    -0x28(%rsp),%r10
mov    -0x30(%rsp),%r12
je     1210 <cimg2dimg+0x110>
jmp    1292 <cimg2dimg+0x192>
xor    %r12d,%r12d
testb  $0x1,-0x48(%rsp)
jne    13de <cimg2dimg+0x2de>
jmp    1420 <cimg2dimg+0x320>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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