void fmt_convert_from_double(double *dest[7], const void *vsrc, int nch,
                             int n) {
  int i, j, ch;
  const double *src = vsrc;
  for (ch = 0; ch < nch; ch++) {
    double *dest_ch = dest[ch];
    const double *src_ch = src + ch;
    for (i = 0, j = 0; i < n; i++, j += nch)
      dest_ch[i] = (double)src_ch[j];
  }
}


x86 assembly (clang, -O0):
<fmt_convert_from_double>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ae <fmt_convert_from_double+0xae>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    119b <fmt_convert_from_double+0x9b>
mov    -0x40(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
jmp    115e <fmt_convert_from_double+0x5e>
jmp    11a0 <fmt_convert_from_double+0xa0>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1121 <fmt_convert_from_double+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fmt_convert_from_double>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %ecx,-0x3c(%rsp)
mov    %rsi,-0x38(%rsp)
mov    %rdi,-0x10(%rsp)
test   %edx,%edx
jle    12e5 <fmt_convert_from_double+0x1e5>
mov    %edx,%eax
mov    -0x3c(%rsp),%r13d
lea    -0x1(%r13),%rcx
imul   %rax,%rcx
mov    %rcx,-0x8(%rsp)
mov    %r13d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,-0x20(%rsp)
imul   %rax,%rcx
mov    %rcx,-0x28(%rsp)
mov    %r13d,%ecx
and    $0x3,%ecx
mov    %rcx,-0x18(%rsp)
lea    0x0(,%r13,8),%rcx
and    $0xffffffffffffffe0,%rcx
mov    %rcx,-0x30(%rsp)
mov    %rax,%rbx
shl    $0x5,%rbx
lea    0x0(,%rax,8),%rcx
lea    (%rcx,%rcx,2),%r15
mov    %rax,%r12
shl    $0x4,%r12
xor    %edx,%edx
mov    -0x38(%rsp),%rbp
jmp    11a1 <fmt_convert_from_double+0xa1>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdx
add    $0x8,%rbp
cmp    %rax,%rdx
je     12e5 <fmt_convert_from_double+0x1e5>
mov    -0x38(%rsp),%rdi
lea    (%rdi,%rdx,8),%rsi
mov    -0x8(%rsp),%rcx
add    %rdx,%rcx
lea    (%rdi,%rcx,8),%rdi
cmp    %rdi,%rsi
mov    %rdi,%rcx
cmovb  %rsi,%rcx
cmova  %rsi,%rdi
cmpl   $0x0,-0x3c(%rsp)
jle    1190 <fmt_convert_from_double+0x90>
mov    -0x10(%rsp),%rsi
mov    (%rsi,%rdx,8),%r10
cmpl   $0x4,-0x3c(%rsp)
jb     11f5 <fmt_convert_from_double+0xf5>
add    $0x8,%rdi
cmp    %r10,%rdi
jbe    1293 <fmt_convert_from_double+0x193>
lea    (%r10,%r13,8),%rsi
cmp    %rsi,%rcx
jae    1293 <fmt_convert_from_double+0x193>
xor    %r9d,%r9d
xor    %r11d,%r11d
mov    %r11,%rdi
not    %rdi
add    %r13,%rdi
mov    -0x18(%rsp),%rcx
test   %rcx,%rcx
je     122a <fmt_convert_from_double+0x12a>
xchg   %ax,%ax
movsd  0x0(%rbp,%r9,8),%xmm0
movsd  %xmm0,(%r10,%r11,8)
add    $0x1,%r11
add    %rax,%r9
add    $0xffffffffffffffff,%rcx
jne    1210 <fmt_convert_from_double+0x110>
cmp    $0x3,%rdi
jb     1190 <fmt_convert_from_double+0x90>
lea    0x0(,%r9,8),%rdi
lea    (%r15,%r9,8),%r8
lea    (%r12,%r9,8),%r14
add    %rax,%r9
mov    %rbp,%rcx
nopw   0x0(%rax,%rax,1)
movsd  (%rcx,%rdi,1),%xmm0
movsd  %xmm0,(%r10,%r11,8)
movsd  (%rcx,%r9,8),%xmm0
movsd  %xmm0,0x8(%r10,%r11,8)
movsd  (%rcx,%r14,1),%xmm0
movsd  %xmm0,0x10(%r10,%r11,8)
movsd  (%rcx,%r8,1),%xmm0
movsd  %xmm0,0x18(%r10,%r11,8)
add    $0x4,%r11
add    %rbx,%rcx
cmp    %r11,%r13
jne    1250 <fmt_convert_from_double+0x150>
jmp    1190 <fmt_convert_from_double+0x90>
mov    %rbp,%rdi
xor    %ecx,%ecx
mov    -0x30(%rsp),%rsi
nopl   (%rax)
movsd  (%rdi),%xmm0
movhps (%rdi,%rax,8),%xmm0
movsd  (%rdi,%r12,1),%xmm1
movhps (%rdi,%r15,1),%xmm1
movups %xmm0,(%r10,%rcx,1)
movups %xmm1,0x10(%r10,%rcx,1)
add    $0x20,%rcx
add    %rbx,%rdi
cmp    %rcx,%rsi
jne    12a0 <fmt_convert_from_double+0x1a0>
mov    -0x28(%rsp),%r9
mov    -0x20(%rsp),%rcx
mov    %rcx,%r11
cmp    %r13,%rcx
je     1190 <fmt_convert_from_double+0x90>
jmp    11fb <fmt_convert_from_double+0xfb>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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