void PutColumn(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++) {
    Image[Index] = (float)Column[j];
    Index += Nx;
  }
}


x86 assembly (clang, -O0):
<PutColumn>:
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 <PutColumn+0x72>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
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 <PutColumn+0x2a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<PutColumn>:
push   %rbp
push   %r14
push   %rbx
test   %r9,%r9
jle    1283 <PutColumn+0x183>
cmp    $0x4,%r9
jae    111b <PutColumn+0x1b>
xor    %r10d,%r10d
jmp    1216 <PutColumn+0x116>
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    1216 <PutColumn+0x116>
test   %r14b,%r14b
jne    1216 <PutColumn+0x116>
and    %r11b,%al
jne    1216 <PutColumn+0x116>
lea    -0x80000000(%rcx),%rdx
movabs $0xffffffff00000000,%rax
ff ff ff
cmp    %rax,%rdx
jb     1216 <PutColumn+0x116>
lea    -0x80000000(%rsi),%rdx
cmp    %rax,%rdx
jb     1216 <PutColumn+0x116>
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)
movupd (%r8,%rbx,8),%xmm0
movupd 0x10(%r8,%rbx,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
movss  %xmm0,(%rax)
shufps $0x55,%xmm0,%xmm0
movss  %xmm0,(%rcx)
movss  %xmm1,(%rax,%rsi,8)
shufps $0x55,%xmm1,%xmm1
movss  %xmm1,(%rax,%r14,1)
add    $0x4,%rbx
add    %r11,%rax
add    %r11,%rcx
cmp    %rbx,%r10
jne    11d0 <PutColumn+0xd0>
cmp    %r9,%r10
je     1283 <PutColumn+0x183>
jmp    1219 <PutColumn+0x119>
mov    %rcx,%rdx
mov    %r10,%rax
not    %rax
test   $0x1,%r9b
je     123e <PutColumn+0x13e>
movsd  (%r8,%r10,8),%xmm0
cvtsd2ss %xmm0,%xmm0
movslq %edx,%rdx
movss  %xmm0,(%rdi,%rdx,4)
add    %rsi,%rdx
or     $0x1,%r10
add    %r9,%rax
je     1283 <PutColumn+0x183>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%r8,%r10,8),%xmm0
cvtsd2ss %xmm0,%xmm0
movslq %edx,%rax
movss  %xmm0,(%rdi,%rax,4)
add    %esi,%eax
movsd  0x8(%r8,%r10,8),%xmm0
cvtsd2ss %xmm0,%xmm0
movslq %eax,%rdx
movss  %xmm0,(%rdi,%rdx,4)
add    %rsi,%rdx
add    $0x2,%r10
cmp    %r10,%r9
jne    1250 <PutColumn+0x150>
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define void @PutColumn(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 double*, double** %11
%25 = load i32, i32* %13
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = load double, double* %27
%29 = fptrunc double %28 to float
%30 = load float*, float** %7
%31 = load i32, i32* %14
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %30, i64 %32
store float %29, float* %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 @PutColumn(float* nocapture writeonly %0, i64 %1, i64 %2, i64 %3, double* nocapture %4, i64 %5) {
%7 = icmp sgt i64 %5, 0
br i1 %7, label %8, label %110
8: ; preds = %6
%9 = icmp ult i64 %5, 4
br i1 %9, label %70, 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 %70, 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 ], [ %66, %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 double, double* %4, i64 %44
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = getelementptr inbounds double, double* %50, i64 2
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = fptrunc <2 x double> %52 to <2 x float>
%57 = fptrunc <2 x double> %55 to <2 x float>
%58 = getelementptr inbounds float, float* %0, i64 %46
%59 = getelementptr inbounds float, float* %0, i64 %47
%60 = getelementptr inbounds float, float* %0, i64 %48
%61 = getelementptr inbounds float, float* %0, i64 %49
%62 = extractelement <2 x float> %56, i64 0
store float %62, float* %58
%63 = extractelement <2 x float> %56, i64 1
store float %63, float* %59
%64 = extractelement <2 x float> %57, i64 0
store float %64, float* %60
%65 = extractelement <2 x float> %57, i64 1
store float %65, float* %61
%66 = add nuw i64 %44, 4
%67 = icmp eq i64 %66, %38
br i1 %67, label %68, label %43
68: ; preds = %43
%69 = icmp eq i64 %38, %5
br i1 %69, label %110, label %70
70: ; preds = %10, %8, %68
%71 = phi i64 [ 0, %10 ], [ 0, %8 ], [ %38, %68 ]
%72 = phi i64 [ %3, %10 ], [ %3, %8 ], [ %40, %68 ]
%73 = xor i64 %71, -1
%74 = and i64 %5, 1
%75 = icmp eq i64 %74, 0
br i1 %75, label %85, label %76
76: ; preds = %70
%77 = getelementptr inbounds double, double* %4, i64 %71
%78 = load double, double* %77
%79 = fptrunc double %78 to float
%80 = shl i64 %72, 32
%81 = ashr exact i64 %80, 32
%82 = getelementptr inbounds float, float* %0, i64 %81
store float %79, float* %82
%83 = add nsw i64 %81, %1
%84 = or i64 %71, 1
br label %85
85: ; preds = %76, %70
%86 = phi i64 [ %71, %70 ], [ %84, %76 ]
%87 = phi i64 [ %72, %70 ], [ %83, %76 ]
%88 = sub i64 0, %5
%89 = icmp eq i64 %73, %88
br i1 %89, label %110, label %90
90: ; preds = %85, %90
%91 = phi i64 [ %108, %90 ], [ %86, %85 ]
%92 = phi i64 [ %107, %90 ], [ %87, %85 ]
%93 = getelementptr inbounds double, double* %4, i64 %91
%94 = load double, double* %93
%95 = fptrunc double %94 to float
%96 = shl i64 %92, 32
%97 = ashr exact i64 %96, 32
%98 = getelementptr inbounds float, float* %0, i64 %97
store float %95, float* %98
%99 = add i64 %92, %1
%100 = add nuw nsw i64 %91, 1
%101 = getelementptr inbounds double, double* %4, i64 %100
%102 = load double, double* %101
%103 = fptrunc double %102 to float
%104 = shl i64 %99, 32
%105 = ashr exact i64 %104, 32
%106 = getelementptr inbounds float, float* %0, i64 %105
store float %103, float* %106
%107 = add nsw i64 %105, %1
%108 = add nuw nsw i64 %91, 2
%109 = icmp eq i64 %108, %5
br i1 %109, label %110, label %90
110: ; preds = %85, %90, %68, %6
ret void
}