void copy_strided(double *out, const double *in, unsigned g, unsigned s,
                  unsigned n) {
  if (g == 1)
    for (; n; --n, in += s)
      *out++ = *in;
  else {
    s *= g;
    for (; n; --n, in += s)
      memcpy(out, in, g * sizeof(double)), out += g;
  }
}


x86 assembly (clang, -O0):
<copy_strided>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
cmpl   $0x1,-0x14(%rbp)
jne    1185 <copy_strided+0x75>
jmp    1139 <copy_strided+0x29>
cmpl   $0x0,-0x1c(%rbp)
je     1180 <copy_strided+0x70>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
movsd  %xmm0,(%rax)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,%ecx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
jmp    1139 <copy_strided+0x29>
jmp    11ea <copy_strided+0xda>
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x0,-0x1c(%rbp)
je     11e5 <copy_strided+0xd5>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%eax
mov    %eax,%edx
shl    $0x3,%rdx
call   1030 <memcpy@plt>
mov    -0x14(%rbp),%ecx
mov    -0x8(%rbp),%rax
mov    %ecx,%ecx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,%ecx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
jmp    118f <copy_strided+0x7f>
jmp    11ea <copy_strided+0xda>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<copy_strided>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r8d,%r15d
mov    %rsi,%rbx
mov    %rdi,%r14
cmp    $0x1,%edx
jne    1213 <copy_strided+0x103>
test   %r15d,%r15d
je     12f9 <copy_strided+0x1e9>
mov    %ecx,%eax
lea    -0x1(%r15),%ecx
mov    %r15d,%edx
and    $0x7,%edx
je     117a <copy_strided+0x6a>
lea    0x0(,%rax,8),%rdi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rbx),%xmm0
movsd  %xmm0,(%r14)
add    $0x8,%r14
add    $0x1,%esi
add    %rdi,%rbx
cmp    %esi,%edx
jne    1160 <copy_strided+0x50>
sub    %esi,%r15d
cmp    $0x7,%ecx
jb     12f9 <copy_strided+0x1e9>
mov    %r15d,%ecx
shl    $0x3,%rax
xor    %edx,%edx
nopl   0x0(%rax)
movsd  (%rbx),%xmm0
movsd  %xmm0,(%r14,%rdx,8)
movsd  (%rbx,%rax,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x8(%r14,%rdx,8)
movsd  (%rax,%rbx,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x10(%r14,%rdx,8)
movsd  (%rax,%rbx,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x18(%r14,%rdx,8)
movsd  (%rax,%rbx,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x20(%r14,%rdx,8)
movsd  (%rax,%rbx,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x28(%r14,%rdx,8)
movsd  (%rax,%rbx,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x30(%r14,%rdx,8)
movsd  (%rax,%rbx,1),%xmm0
add    %rax,%rbx
movsd  %xmm0,0x38(%r14,%rdx,8)
add    $0x8,%rdx
add    %rax,%rbx
cmp    %edx,%ecx
jne    1190 <copy_strided+0x80>
jmp    12f9 <copy_strided+0x1e9>
test   %r15d,%r15d
je     12f9 <copy_strided+0x1e9>
imul   %edx,%ecx
mov    %edx,%r13d
shl    $0x3,%r13
lea    -0x1(%r15),%eax
mov    %r15d,%r12d
and    $0x3,%r12d
mov    %rcx,0x10(%rsp)
je     1274 <copy_strided+0x164>
mov    %eax,(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,0x8(%rsp)
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
mov    %r14,%rdi
mov    %rbx,%rsi
mov    %r13,%rdx
call   1030 <memcpy@plt>
add    $0x1,%ebp
add    %r13,%r14
add    0x8(%rsp),%rbx
cmp    %ebp,%r12d
jne    1250 <copy_strided+0x140>
sub    %ebp,%r15d
mov    (%rsp),%eax
cmp    $0x3,%eax
mov    0x10(%rsp),%rbp
jb     12f9 <copy_strided+0x1e9>
mov    %rbp,%r12
shl    $0x5,%r12
lea    0x0(,%rbp,8),%rax
mov    %rax,0x8(%rsp)
lea    (%rax,%rax,2),%rax
mov    %rax,(%rsp)
shl    $0x4,%rbp
xchg   %ax,%ax
mov    %r14,%rdi
mov    %rbx,%rsi
mov    %r13,%rdx
call   1030 <memcpy@plt>
mov    0x8(%rsp),%rax
lea    (%rbx,%rax,1),%rsi
add    %r13,%r14
mov    %r14,%rdi
mov    %r13,%rdx
call   1030 <memcpy@plt>
lea    (%rbx,%rbp,1),%rsi
add    %r13,%r14
mov    %r14,%rdi
mov    %r13,%rdx
call   1030 <memcpy@plt>
mov    (%rsp),%rax
lea    (%rbx,%rax,1),%rsi
add    %r13,%r14
mov    %r14,%rdi
mov    %r13,%rdx
call   1030 <memcpy@plt>
add    %r12,%rbx
add    %r13,%r14
add    $0xfffffffc,%r15d
jne    12a0 <copy_strided+0x190>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @copy_strided(double* nocapture writeonly %0, double* nocapture %1, i32 %2, i32 %3, i32 %4) {
%6 = icmp eq i32 %2, 1
%7 = icmp eq i32 %4, 0
br i1 %6, label %8, label %60
8: ; preds = %5
br i1 %7, label %108, label %9
9: ; preds = %8
%10 = zext i32 %3 to i64
%11 = add i32 %4, -1
%12 = and i32 %4, 7
%13 = icmp eq i32 %12, 0
br i1 %13, label %25, label %14
14: ; preds = %9, %14
%15 = phi i32 [ %21, %14 ], [ %4, %9 ]
%16 = phi double* [ %20, %14 ], [ %0, %9 ]
%17 = phi double* [ %22, %14 ], [ %1, %9 ]
%18 = phi i32 [ %23, %14 ], [ 0, %9 ]
%19 = load double, double* %17
%20 = getelementptr inbounds double, double* %16, i64 1
store double %19, double* %16
%21 = add i32 %15, -1
%22 = getelementptr inbounds double, double* %17, i64 %10
%23 = add i32 %18, 1
%24 = icmp eq i32 %23, %12
br i1 %24, label %25, label %14
25: ; preds = %14, %9
%26 = phi i32 [ %4, %9 ], [ %21, %14 ]
%27 = phi double* [ %0, %9 ], [ %20, %14 ]
%28 = phi double* [ %1, %9 ], [ %22, %14 ]
%29 = icmp ult i32 %11, 7
br i1 %29, label %108, label %30
30: ; preds = %25, %30
%31 = phi i32 [ %57, %30 ], [ %26, %25 ]
%32 = phi double* [ %56, %30 ], [ %27, %25 ]
%33 = phi double* [ %58, %30 ], [ %28, %25 ]
%34 = load double, double* %33
%35 = getelementptr inbounds double, double* %32, i64 1
store double %34, double* %32
%36 = getelementptr inbounds double, double* %33, i64 %10
%37 = load double, double* %36
%38 = getelementptr inbounds double, double* %32, i64 2
store double %37, double* %35
%39 = getelementptr inbounds double, double* %36, i64 %10
%40 = load double, double* %39
%41 = getelementptr inbounds double, double* %32, i64 3
store double %40, double* %38
%42 = getelementptr inbounds double, double* %39, i64 %10
%43 = load double, double* %42
%44 = getelementptr inbounds double, double* %32, i64 4
store double %43, double* %41
%45 = getelementptr inbounds double, double* %42, i64 %10
%46 = load double, double* %45
%47 = getelementptr inbounds double, double* %32, i64 5
store double %46, double* %44
%48 = getelementptr inbounds double, double* %45, i64 %10
%49 = load double, double* %48
%50 = getelementptr inbounds double, double* %32, i64 6
store double %49, double* %47
%51 = getelementptr inbounds double, double* %48, i64 %10
%52 = load double, double* %51
%53 = getelementptr inbounds double, double* %32, i64 7
store double %52, double* %50
%54 = getelementptr inbounds double, double* %51, i64 %10
%55 = load double, double* %54
%56 = getelementptr inbounds double, double* %32, i64 8
store double %55, double* %53
%57 = add i32 %31, -8
%58 = getelementptr inbounds double, double* %54, i64 %10
%59 = icmp eq i32 %57, 0
br i1 %59, label %108, label %30
60: ; preds = %5
br i1 %7, label %108, label %61
61: ; preds = %60
%62 = mul i32 %3, %2
%63 = zext i32 %2 to i64
%64 = shl nuw nsw i64 %63, 3
%65 = zext i32 %62 to i64
%66 = add i32 %4, -1
%67 = and i32 %4, 3
%68 = icmp eq i32 %67, 0
br i1 %68, label %81, label %69
69: ; preds = %61, %69
%70 = phi i32 [ %77, %69 ], [ %4, %61 ]
%71 = phi double* [ %76, %69 ], [ %0, %61 ]
%72 = phi double* [ %78, %69 ], [ %1, %61 ]
%73 = phi i32 [ %79, %69 ], [ 0, %61 ]
%74 = bitcast double* %71 to i8*
%75 = bitcast double* %72 to i8*
tail call void @llvm.memcpy.p0i8.p0i8.i64(i8* align 8 %74, i8* align 8 %75, i64 %64, i1 false)
%76 = getelementptr inbounds double, double* %71, i64 %63
%77 = add i32 %70, -1
%78 = getelementptr inbounds double, double* %72, i64 %65
%79 = add i32 %73, 1
%80 = icmp eq i32 %79, %67
br i1 %80, label %81, label %69
81: ; preds = %69, %61
%82 = phi i32 [ %4, %61 ], [ %77, %69 ]
%83 = phi double* [ %0, %61 ], [ %76, %69 ]
%84 = phi double* [ %1, %61 ], [ %78, %69 ]
%85 = icmp ult i32 %66, 3
br i1 %85, label %108, label %86
86: ; preds = %81, %86
%87 = phi i32 [ %105, %86 ], [ %82, %81 ]
%88 = phi double* [ %104, %86 ], [ %83, %81 ]
%89 = phi double* [ %106, %86 ], [ %84, %81 ]
%90 = bitcast double* %88 to i8*
%91 = bitcast double* %89 to i8*
tail call void @llvm.memcpy.p0i8.p0i8.i64(i8* align 8 %90, i8* align 8 %91, i64 %64, i1 false)
%92 = getelementptr inbounds double, double* %88, i64 %63
%93 = getelementptr inbounds double, double* %89, i64 %65
%94 = bitcast double* %92 to i8*
%95 = bitcast double* %93 to i8*
tail call void @llvm.memcpy.p0i8.p0i8.i64(i8* align 8 %94, i8* align 8 %95, i64 %64, i1 false)
%96 = getelementptr inbounds double, double* %92, i64 %63
%97 = getelementptr inbounds double, double* %93, i64 %65
%98 = bitcast double* %96 to i8*
%99 = bitcast double* %97 to i8*
tail call void @llvm.memcpy.p0i8.p0i8.i64(i8* align 8 %98, i8* align 8 %99, i64 %64, i1 false)
%100 = getelementptr inbounds double, double* %96, i64 %63
%101 = getelementptr inbounds double, double* %97, i64 %65
%102 = bitcast double* %100 to i8*
%103 = bitcast double* %101 to i8*
tail call void @llvm.memcpy.p0i8.p0i8.i64(i8* align 8 %102, i8* align 8 %103, i64 %64, i1 false)
%104 = getelementptr inbounds double, double* %100, i64 %63
%105 = add i32 %87, -4
%106 = getelementptr inbounds double, double* %101, i64 %65
%107 = icmp eq i32 %105, 0
br i1 %107, label %108, label %86
108: ; preds = %81, %86, %25, %30, %60, %8
ret void
}