void dither_double(const double *in, double *out, int frames, int channels) {
  int ch, k;
  for (ch = 0; ch < channels; ch++)
    for (k = ch; k < channels * frames; k += channels)
      out[k] = in[k];
}


x86 assembly (clang, -O0):
<dither_double>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1178 <dither_double+0x78>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    1165 <dither_double+0x65>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
jmp    112b <dither_double+0x2b>
jmp    116a <dither_double+0x6a>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <dither_double+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dither_double>:
test   %ecx,%ecx
jle    1134 <dither_double+0x34>
imul   %ecx,%edx
mov    %ecx,%r8d
movslq %edx,%rax
cmp    $0x1,%ecx
jne    1135 <dither_double+0x35>
xor    %r10d,%r10d
test   $0x1,%r8b
je     1134 <dither_double+0x34>
cmp    %r10d,%edx
jle    1134 <dither_double+0x34>
movsd  (%rdi,%r10,8),%xmm0
movsd  %xmm0,(%rsi,%r10,8)
add    %r8,%r10
cmp    %rax,%r10
jl     1120 <dither_double+0x20>
ret
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xor    %r10d,%r10d
jmp    1159 <dither_double+0x59>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x2,%r10
cmp    %r9,%r10
je     1115 <dither_double+0x15>
cmp    %r10d,%edx
jle    1182 <dither_double+0x82>
mov    %r10,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rsi,%rcx,8)
add    %r8,%rcx
cmp    %rax,%rcx
jl     1170 <dither_double+0x70>
mov    %r10,%rcx
or     $0x1,%rcx
cmp    %ecx,%edx
jle    1150 <dither_double+0x50>
nopl   (%rax)
movsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rsi,%rcx,8)
add    %r8,%rcx
cmp    %rax,%rcx
jl     1190 <dither_double+0x90>
jmp    1150 <dither_double+0x50>

LLVM IR (-O0):
define void @dither_double(double* %0, double* %1, i32 %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
store i32 0, i32* %9
br label %11
11: ; preds = %38, %4
%12 = load i32, i32* %9
%13 = load i32, i32* %8
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %41
15: ; preds = %11
%16 = load i32, i32* %9
store i32 %16, i32* %10
br label %17
17: ; preds = %33, %15
%18 = load i32, i32* %10
%19 = load i32, i32* %8
%20 = load i32, i32* %7
%21 = mul nsw i32 %19, %20
%22 = icmp slt i32 %18, %21
br i1 %22, label %23, label %37
23: ; preds = %17
%24 = load double*, double** %5
%25 = load i32, i32* %10
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = load double, double* %27
%29 = load double*, double** %6
%30 = load i32, i32* %10
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %29, i64 %31
store double %28, double* %32
br label %33
33: ; preds = %23
%34 = load i32, i32* %8
%35 = load i32, i32* %10
%36 = add nsw i32 %35, %34
store i32 %36, i32* %10
br label %17
37: ; preds = %17
br label %38
38: ; preds = %37
%39 = load i32, i32* %9
%40 = add nsw i32 %39, 1
store i32 %40, i32* %9
br label %11
41: ; preds = %11
ret void
}

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