void GetColumn(float *Image, long Nx, long Ny, long ColumnNb, double *Column,
               long ColumnSize) {
  int j;
  int Index;
  Ny;
  Index = ColumnNb;
  for (j = 0L; j < ColumnSize; j++) {
    Column[j] = (double)Image[Index];
    Index += Nx;
  }
}


x86 assembly (clang, -O0):
<GetColumn>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
movslq -0x34(%rbp),%rax
cmp    -0x30(%rbp),%rax
jge    1172 <GetColumn+0x72>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rcx
movslq -0x38(%rbp),%rax
add    %rcx,%rax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112a <GetColumn+0x2a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<GetColumn>:
push   %rbp
push   %r14
push   %rbx
test   %r9,%r9
jle    1273 <GetColumn+0x173>
cmp    $0x4,%r9
jae    111b <GetColumn+0x1b>
xor    %r10d,%r10d
jmp    1210 <GetColumn+0x110>
mov    %esi,%eax
neg    %eax
setb   %r11b
test   %esi,%esi
cmovns %esi,%eax
lea    -0x1(%r9),%rbx
mul    %ebx
seto   %r14b
mov    %ecx,%edx
sub    %eax,%edx
add    %ecx,%eax
xor    %ebp,%ebp
cmp    %ecx,%eax
setl   %bpl
xor    %eax,%eax
cmp    %ecx,%edx
setg   %al
test   %esi,%esi
cmovs  %eax,%ebp
shr    $0x20,%rbx
setne  %al
xor    %r10d,%r10d
test   %bpl,%bpl
jne    1210 <GetColumn+0x110>
test   %r14b,%r14b
jne    1210 <GetColumn+0x110>
and    %r11b,%al
jne    1210 <GetColumn+0x110>
lea    -0x80000000(%rcx),%rdx
movabs $0xffffffff00000000,%rax
ff ff ff
cmp    %rax,%rdx
jb     1210 <GetColumn+0x110>
lea    -0x80000000(%rsi),%rdx
cmp    %rax,%rdx
jb     1210 <GetColumn+0x110>
mov    %r9,%r10
and    $0xfffffffffffffffc,%r10
mov    %r10,%rdx
imul   %rsi,%rdx
add    %rcx,%rdx
lea    (%rdi,%rcx,4),%rax
mov    %rsi,%r11
shl    $0x4,%r11
lea    0x0(,%rsi,4),%rbx
lea    (%rbx,%rbx,2),%r14
add    %rsi,%rcx
lea    (%rdi,%rcx,4),%rcx
xor    %ebx,%ebx
nopl   0x0(%rax)
movss  (%rcx),%xmm0
movss  (%rax),%xmm1
unpcklps %xmm0,%xmm1
movss  (%rax,%r14,1),%xmm0
movss  (%rax,%rsi,8),%xmm2
unpcklps %xmm0,%xmm2
cvtps2pd %xmm1,%xmm0
cvtps2pd %xmm2,%xmm1
movups %xmm0,(%r8,%rbx,8)
movups %xmm1,0x10(%r8,%rbx,8)
add    $0x4,%rbx
add    %r11,%rax
add    %r11,%rcx
cmp    %rbx,%r10
jne    11d0 <GetColumn+0xd0>
cmp    %r9,%r10
je     1273 <GetColumn+0x173>
jmp    1213 <GetColumn+0x113>
mov    %rcx,%rdx
mov    %r10,%rax
not    %rax
test   $0x1,%r9b
je     1238 <GetColumn+0x138>
movslq %edx,%rdx
movss  (%rdi,%rdx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%r8,%r10,8)
add    %rsi,%rdx
or     $0x1,%r10
add    %r9,%rax
je     1273 <GetColumn+0x173>
nopl   (%rax)
movslq %edx,%rax
movss  (%rdi,%rax,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%r8,%r10,8)
add    %esi,%eax
movslq %eax,%rdx
movss  (%rdi,%rdx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x8(%r8,%r10,8)
add    %rsi,%rdx
add    $0x2,%r10
cmp    %r10,%r9
jne    1240 <GetColumn+0x140>
pop    %rbx
pop    %r14
pop    %rbp
ret

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

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