void GetRow(float *Image, long Nx, long Ny, long RowNb, double *Row,
            long RowSize) {
  int i;
  int BaseIndex;
  Ny;
  BaseIndex = RowNb * Nx;
  for (i = 0L; i < RowSize; i++)
    Row[i] = (double)Image[BaseIndex + i];
}


x86 assembly (clang, -O0):
<GetRow>:
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
imul   -0x10(%rbp),%rax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
movslq -0x34(%rbp),%rax
cmp    -0x30(%rbp),%rax
jge    116e <GetRow+0x6e>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112f <GetRow+0x2f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<GetRow>:
test   %r9,%r9
jle    11c8 <GetRow+0xc8>
imul   %esi,%ecx
movslq %ecx,%r10
cmp    $0x4,%r9
jae    111c <GetRow+0x1c>
xor    %eax,%eax
jmp    11a7 <GetRow+0xa7>
mov    %r9,%rax
and    $0xfffffffffffffffc,%rax
lea    -0x4(%rax),%rcx
mov    %rcx,%r11
shr    $0x2,%r11
add    $0x1,%r11
test   %rcx,%rcx
je     11c9 <GetRow+0xc9>
mov    %r11,%rdx
and    $0xfffffffffffffffe,%rdx
lea    (%rdi,%r10,4),%rcx
add    $0x18,%rcx
xor    %esi,%esi
nopl   0x0(%rax)
cvtps2pd -0x18(%rcx,%rsi,4),%xmm0
cvtps2pd -0x10(%rcx,%rsi,4),%xmm1
movups %xmm0,(%r8,%rsi,8)
movups %xmm1,0x10(%r8,%rsi,8)
cvtps2pd -0x8(%rcx,%rsi,4),%xmm0
cvtps2pd (%rcx,%rsi,4),%xmm1
movups %xmm0,0x20(%r8,%rsi,8)
movups %xmm1,0x30(%r8,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdx
jne    1150 <GetRow+0x50>
test   $0x1,%r11b
je     11a2 <GetRow+0xa2>
lea    (%rsi,%r10,1),%rcx
cvtps2pd (%rdi,%rcx,4),%xmm0
cvtps2pd 0x8(%rdi,%rcx,4),%xmm1
movups %xmm0,(%r8,%rsi,8)
movups %xmm1,0x10(%r8,%rsi,8)
cmp    %r9,%rax
je     11c8 <GetRow+0xc8>
lea    (%rdi,%r10,4),%rcx
nopl   0x0(%rax,%rax,1)
movss  (%rcx,%rax,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%r8,%rax,8)
add    $0x1,%rax
cmp    %rax,%r9
jne    11b0 <GetRow+0xb0>
ret
xor    %esi,%esi
test   $0x1,%r11b
jne    118a <GetRow+0x8a>
jmp    11a2 <GetRow+0xa2>

LLVM IR (-O0):
define void @GetRow(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 = load i64, i64* %8
%18 = mul nsw i64 %16, %17
%19 = trunc i64 %18 to i32
store i32 %19, i32* %14
store i32 0, i32* %13
br label %20
20: ; preds = %38, %6
%21 = load i32, i32* %13
%22 = sext i32 %21 to i64
%23 = load i64, i64* %12
%24 = icmp slt i64 %22, %23
br i1 %24, label %25, label %41
25: ; preds = %20
%26 = load float*, float** %7
%27 = load i32, i32* %14
%28 = load i32, i32* %13
%29 = add nsw i32 %27, %28
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %26, i64 %30
%32 = load float, float* %31
%33 = fpext float %32 to double
%34 = load double*, double** %11
%35 = load i32, i32* %13
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
store double %33, double* %37
br label %38
38: ; preds = %25
%39 = load i32, i32* %13
%40 = add nsw i32 %39, 1
store i32 %40, i32* %13
br label %20
41: ; preds = %20
ret void
}

LLVM IR (-O2):
define void @GetRow(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 %85
8: ; preds = %6
%9 = mul nsw i64 %3, %1
%10 = shl i64 %9, 32
%11 = ashr exact i64 %10, 32
%12 = icmp ult i64 %5, 4
br i1 %12, label %74, label %13
13: ; preds = %8
%14 = and i64 %5, -4
%15 = add i64 %14, -4
%16 = lshr exact i64 %15, 2
%17 = add nuw nsw i64 %16, 1
%18 = and i64 %17, 1
%19 = icmp eq i64 %15, 0
br i1 %19, label %55, label %20
20: ; preds = %13
%21 = and i64 %17, 9223372036854775806
br label %22
22: ; preds = %22, %20
%23 = phi i64 [ 0, %20 ], [ %52, %22 ]
%24 = phi i64 [ 0, %20 ], [ %53, %22 ]
%25 = add nsw i64 %23, %11
%26 = getelementptr inbounds float, float* %0, i64 %25
%27 = bitcast float* %26 to <2 x float>*
%28 = load <2 x float>, <2 x float>* %27
%29 = getelementptr inbounds float, float* %26, i64 2
%30 = bitcast float* %29 to <2 x float>*
%31 = load <2 x float>, <2 x float>* %30
%32 = fpext <2 x float> %28 to <2 x double>
%33 = fpext <2 x float> %31 to <2 x double>
%34 = getelementptr inbounds double, double* %4, i64 %23
%35 = bitcast double* %34 to <2 x double>*
store <2 x double> %32, <2 x double>* %35
%36 = getelementptr inbounds double, double* %34, i64 2
%37 = bitcast double* %36 to <2 x double>*
store <2 x double> %33, <2 x double>* %37
%38 = or i64 %23, 4
%39 = add nsw i64 %38, %11
%40 = getelementptr inbounds float, float* %0, i64 %39
%41 = bitcast float* %40 to <2 x float>*
%42 = load <2 x float>, <2 x float>* %41
%43 = getelementptr inbounds float, float* %40, i64 2
%44 = bitcast float* %43 to <2 x float>*
%45 = load <2 x float>, <2 x float>* %44
%46 = fpext <2 x float> %42 to <2 x double>
%47 = fpext <2 x float> %45 to <2 x double>
%48 = getelementptr inbounds double, double* %4, i64 %38
%49 = bitcast double* %48 to <2 x double>*
store <2 x double> %46, <2 x double>* %49
%50 = getelementptr inbounds double, double* %48, i64 2
%51 = bitcast double* %50 to <2 x double>*
store <2 x double> %47, <2 x double>* %51
%52 = add nuw i64 %23, 8
%53 = add i64 %24, 2
%54 = icmp eq i64 %53, %21
br i1 %54, label %55, label %22
55: ; preds = %22, %13
%56 = phi i64 [ 0, %13 ], [ %52, %22 ]
%57 = icmp eq i64 %18, 0
br i1 %57, label %72, label %58
58: ; preds = %55
%59 = add nsw i64 %56, %11
%60 = getelementptr inbounds float, float* %0, i64 %59
%61 = bitcast float* %60 to <2 x float>*
%62 = load <2 x float>, <2 x float>* %61
%63 = getelementptr inbounds float, float* %60, i64 2
%64 = bitcast float* %63 to <2 x float>*
%65 = load <2 x float>, <2 x float>* %64
%66 = fpext <2 x float> %62 to <2 x double>
%67 = fpext <2 x float> %65 to <2 x double>
%68 = getelementptr inbounds double, double* %4, i64 %56
%69 = bitcast double* %68 to <2 x double>*
store <2 x double> %66, <2 x double>* %69
%70 = getelementptr inbounds double, double* %68, i64 2
%71 = bitcast double* %70 to <2 x double>*
store <2 x double> %67, <2 x double>* %71
br label %72
72: ; preds = %55, %58
%73 = icmp eq i64 %14, %5
br i1 %73, label %85, label %74
74: ; preds = %8, %72
%75 = phi i64 [ 0, %8 ], [ %14, %72 ]
br label %76
76: ; preds = %74, %76
%77 = phi i64 [ %83, %76 ], [ %75, %74 ]
%78 = add nsw i64 %77, %11
%79 = getelementptr inbounds float, float* %0, i64 %78
%80 = load float, float* %79
%81 = fpext float %80 to double
%82 = getelementptr inbounds double, double* %4, i64 %77
store double %81, double* %82
%83 = add nuw nsw i64 %77, 1
%84 = icmp eq i64 %83, %5
br i1 %84, label %85, label %76
85: ; preds = %76, %72, %6
ret void
}