int ffu4fr8(unsigned long *input, long ntodo, double scale, double zero,
            double *output, int *status) {
  long ii;
  if (scale == 1. && zero == 0.) {
    for (ii = 0; ii < ntodo; ii++)
      output[ii] = (double)input[ii];
  } else {
    for (ii = 0; ii < ntodo; ii++)
      output[ii] = (input[ii] - zero) / scale;
  }
  return (*status);
}


x86 assembly (clang, -O0):
<ffu4fr8>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xed5(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    11bb <ffu4fr8+0xbb>
jp     11bb <ffu4fr8+0xbb>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11bb <ffu4fr8+0xbb>
jp     11bb <ffu4fr8+0xbb>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11b6 <ffu4fr8+0xb6>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movq   (%rax,%rcx,8),%xmm0
movaps 0xe93(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xe97(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    115b <ffu4fr8+0x5b>
jmp    122d <ffu4fr8+0x12d>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1228 <ffu4fr8+0x128>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movq   (%rax,%rcx,8),%xmm0
movaps 0xe2b(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xe2f(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11c3 <ffu4fr8+0xc3>
jmp    122d <ffu4fr8+0x12d>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<ffu4fr8>:
ucomisd 0xef8(%rip),%xmm0
jne    112e <ffu4fr8+0x2e>
jp     112e <ffu4fr8+0x2e>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    112e <ffu4fr8+0x2e>
jp     112e <ffu4fr8+0x2e>
test   %rsi,%rsi
jle    147e <ffu4fr8+0x37e>
cmp    $0x4,%rsi
jae    1148 <ffu4fr8+0x48>
xor    %eax,%eax
jmp    12d7 <ffu4fr8+0x1d7>
test   %rsi,%rsi
jle    147e <ffu4fr8+0x37e>
cmp    $0x1,%rsi
jne    131b <ffu4fr8+0x21b>
xor    %eax,%eax
jmp    143e <ffu4fr8+0x33e>
mov    %rsi,%rax
and    $0xfffffffffffffffc,%rax
lea    -0x4(%rax),%r8
mov    %r8,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %r8,%r8
je     1490 <ffu4fr8+0x390>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r8d,%r8d
movdqa 0xe97(%rip),%xmm0
movdqa 0xe9f(%rip),%xmm1
movdqa 0xea7(%rip),%xmm2
movapd 0xeaf(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movdqu (%rdi,%r8,8),%xmm4
movdqu 0x10(%rdi,%r8,8),%xmm5
movdqa %xmm4,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrlq  $0x20,%xmm4
por    %xmm2,%xmm4
subpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
movdqa %xmm5,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrlq  $0x20,%xmm5
por    %xmm2,%xmm5
subpd  %xmm3,%xmm5
addpd  %xmm6,%xmm5
movupd %xmm4,(%rdx,%r8,8)
movupd %xmm5,0x10(%rdx,%r8,8)
movdqu 0x20(%rdi,%r8,8),%xmm4
movdqu 0x30(%rdi,%r8,8),%xmm5
movdqa %xmm4,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrlq  $0x20,%xmm4
por    %xmm2,%xmm4
subpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
movdqa %xmm5,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrlq  $0x20,%xmm5
por    %xmm2,%xmm5
subpd  %xmm3,%xmm5
addpd  %xmm6,%xmm5
movupd %xmm4,0x20(%rdx,%r8,8)
movupd %xmm5,0x30(%rdx,%r8,8)
add    $0x8,%r8
add    $0xfffffffffffffffe,%r10
jne    11a0 <ffu4fr8+0xa0>
test   $0x1,%r9b
je     12ce <ffu4fr8+0x1ce>
movdqu (%rdi,%r8,8),%xmm0
movdqu 0x10(%rdi,%r8,8),%xmm1
movdqa 0xd9d(%rip),%xmm2
movdqa %xmm0,%xmm3
pand   %xmm2,%xmm3
movdqa 0xd9d(%rip),%xmm4
por    %xmm4,%xmm3
psrlq  $0x20,%xmm0
movdqa 0xd9c(%rip),%xmm5
por    %xmm5,%xmm0
movapd 0xda0(%rip),%xmm6
subpd  %xmm6,%xmm0
addpd  %xmm3,%xmm0
pand   %xmm1,%xmm2
por    %xmm4,%xmm2
psrlq  $0x20,%xmm1
por    %xmm5,%xmm1
subpd  %xmm6,%xmm1
addpd  %xmm2,%xmm1
movupd %xmm0,(%rdx,%r8,8)
movupd %xmm1,0x10(%rdx,%r8,8)
cmp    %rsi,%rax
je     147e <ffu4fr8+0x37e>
movapd 0xd71(%rip),%xmm0
movapd 0xd79(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm2
unpcklps %xmm0,%xmm2
subpd  %xmm1,%xmm2
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdx,%rax,8)
add    $0x1,%rax
cmp    %rax,%rsi
jne    12f0 <ffu4fr8+0x1f0>
jmp    147e <ffu4fr8+0x37e>
mov    %rsi,%rax
and    $0xfffffffffffffffe,%rax
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
lea    -0x2(%rax),%r8
mov    %r8,%r9
shr    %r9
add    $0x1,%r9
test   %r8,%r8
je     1481 <ffu4fr8+0x381>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r8d,%r8d
movdqa 0xcb4(%rip),%xmm8
movdqa 0xcbb(%rip),%xmm9
movdqa 0xcc3(%rip),%xmm6
movapd 0xccb(%rip),%xmm7
cs nopw 0x0(%rax,%rax,1)
nop
movdqu (%rdi,%r8,8),%xmm4
movdqa %xmm4,%xmm5
pand   %xmm8,%xmm5
por    %xmm9,%xmm5
psrlq  $0x20,%xmm4
por    %xmm6,%xmm4
subpd  %xmm7,%xmm4
addpd  %xmm5,%xmm4
subpd  %xmm2,%xmm4
divpd  %xmm3,%xmm4
movupd %xmm4,(%rdx,%r8,8)
movdqu 0x10(%rdi,%r8,8),%xmm4
movdqa %xmm4,%xmm5
pand   %xmm8,%xmm5
por    %xmm9,%xmm5
psrlq  $0x20,%xmm4
por    %xmm6,%xmm4
subpd  %xmm7,%xmm4
addpd  %xmm5,%xmm4
subpd  %xmm2,%xmm4
divpd  %xmm3,%xmm4
movupd %xmm4,0x10(%rdx,%r8,8)
add    $0x4,%r8
add    $0xfffffffffffffffe,%r10
jne    1380 <ffu4fr8+0x280>
test   $0x1,%r9b
je     1439 <ffu4fr8+0x339>
movdqu (%rdi,%r8,8),%xmm4
movdqa 0xc0a(%rip),%xmm5
pand   %xmm4,%xmm5
por    0xc0e(%rip),%xmm5
psrlq  $0x20,%xmm4
por    0xc11(%rip),%xmm4
subpd  0xc19(%rip),%xmm4
addpd  %xmm5,%xmm4
subpd  %xmm2,%xmm4
divpd  %xmm3,%xmm4
movupd %xmm4,(%rdx,%r8,8)
cmp    %rsi,%rax
je     147e <ffu4fr8+0x37e>
movapd 0xc0a(%rip),%xmm2
movapd 0xc12(%rip),%xmm3
xchg   %ax,%ax
movsd  (%rdi,%rax,8),%xmm4
unpcklps %xmm2,%xmm4
subpd  %xmm3,%xmm4
movapd %xmm4,%xmm5
unpckhpd %xmm4,%xmm5
addsd  %xmm4,%xmm5
subsd  %xmm1,%xmm5
divsd  %xmm0,%xmm5
movsd  %xmm5,(%rdx,%rax,8)
add    $0x1,%rax
cmp    %rax,%rsi
jne    1450 <ffu4fr8+0x350>
mov    (%rcx),%eax
ret
xor    %r8d,%r8d
test   $0x1,%r9b
jne    13f8 <ffu4fr8+0x2f8>
jmp    1439 <ffu4fr8+0x339>
xor    %r8d,%r8d
test   $0x1,%r9b
jne    125e <ffu4fr8+0x15e>
jmp    12ce <ffu4fr8+0x1ce>

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

LLVM IR (-O2):
define i32 @ffu4fr8(i64* nocapture %0, i64 %1, double %2, double %3, double* nocapture writeonly %4, i32* nocapture %5) {
%7 = fcmp oeq double %2, 1.000000e+00
%8 = fcmp oeq double %3, 0.000000e+00
%9 = and i1 %7, %8
%10 = icmp sgt i64 %1, 0
br i1 %9, label %66, label %11
11: ; preds = %6
br i1 %10, label %12, label %147
12: ; preds = %11
%13 = icmp eq i64 %1, 1
br i1 %13, label %64, label %14
14: ; preds = %12
%15 = and i64 %1, -2
%16 = insertelement <2 x double> poison, double %3, i64 0
%17 = shufflevector <2 x double> %16, <2 x double> poison, <2 x i32> zeroinitializer
%18 = insertelement <2 x double> poison, double %2, i64 0
%19 = shufflevector <2 x double> %18, <2 x double> poison, <2 x i32> zeroinitializer
%20 = add i64 %15, -2
%21 = lshr exact i64 %20, 1
%22 = add nuw i64 %21, 1
%23 = and i64 %22, 1
%24 = icmp eq i64 %20, 0
br i1 %24, label %50, label %25
25: ; preds = %14
%26 = and i64 %22, -2
br label %27
27: ; preds = %27, %25
%28 = phi i64 [ 0, %25 ], [ %47, %27 ]
%29 = phi i64 [ 0, %25 ], [ %48, %27 ]
%30 = getelementptr inbounds i64, i64* %0, i64 %28
%31 = bitcast i64* %30 to <2 x i64>*
%32 = load <2 x i64>, <2 x i64>* %31
%33 = uitofp <2 x i64> %32 to <2 x double>
%34 = fsub <2 x double> %33, %17
%35 = fdiv <2 x double> %34, %19
%36 = getelementptr inbounds double, double* %4, i64 %28
%37 = bitcast double* %36 to <2 x double>*
store <2 x double> %35, <2 x double>* %37
%38 = or i64 %28, 2
%39 = getelementptr inbounds i64, i64* %0, i64 %38
%40 = bitcast i64* %39 to <2 x i64>*
%41 = load <2 x i64>, <2 x i64>* %40
%42 = uitofp <2 x i64> %41 to <2 x double>
%43 = fsub <2 x double> %42, %17
%44 = fdiv <2 x double> %43, %19
%45 = getelementptr inbounds double, double* %4, i64 %38
%46 = bitcast double* %45 to <2 x double>*
store <2 x double> %44, <2 x double>* %46
%47 = add nuw i64 %28, 4
%48 = add nuw i64 %29, 2
%49 = icmp eq i64 %48, %26
br i1 %49, label %50, label %27
50: ; preds = %27, %14
%51 = phi i64 [ 0, %14 ], [ %47, %27 ]
%52 = icmp eq i64 %23, 0
br i1 %52, label %62, label %53
53: ; preds = %50
%54 = getelementptr inbounds i64, i64* %0, i64 %51
%55 = bitcast i64* %54 to <2 x i64>*
%56 = load <2 x i64>, <2 x i64>* %55
%57 = uitofp <2 x i64> %56 to <2 x double>
%58 = fsub <2 x double> %57, %17
%59 = fdiv <2 x double> %58, %19
%60 = getelementptr inbounds double, double* %4, i64 %51
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %59, <2 x double>* %61
br label %62
62: ; preds = %50, %53
%63 = icmp eq i64 %15, %1
br i1 %63, label %147, label %64
64: ; preds = %12, %62
%65 = phi i64 [ 0, %12 ], [ %15, %62 ]
br label %137
66: ; preds = %6
br i1 %10, label %67, label %147
67: ; preds = %66
%68 = icmp ult i64 %1, 4
br i1 %68, label %127, label %69
69: ; preds = %67
%70 = and i64 %1, -4
%71 = add i64 %70, -4
%72 = lshr exact i64 %71, 2
%73 = add nuw nsw i64 %72, 1
%74 = and i64 %73, 1
%75 = icmp eq i64 %71, 0
br i1 %75, label %109, label %76
76: ; preds = %69
%77 = and i64 %73, 9223372036854775806
br label %78
78: ; preds = %78, %76
%79 = phi i64 [ 0, %76 ], [ %106, %78 ]
%80 = phi i64 [ 0, %76 ], [ %107, %78 ]
%81 = getelementptr inbounds i64, i64* %0, i64 %79
%82 = bitcast i64* %81 to <2 x i64>*
%83 = load <2 x i64>, <2 x i64>* %82
%84 = getelementptr inbounds i64, i64* %81, i64 2
%85 = bitcast i64* %84 to <2 x i64>*
%86 = load <2 x i64>, <2 x i64>* %85
%87 = uitofp <2 x i64> %83 to <2 x double>
%88 = uitofp <2 x i64> %86 to <2 x double>
%89 = getelementptr inbounds double, double* %4, i64 %79
%90 = bitcast double* %89 to <2 x double>*
store <2 x double> %87, <2 x double>* %90
%91 = getelementptr inbounds double, double* %89, i64 2
%92 = bitcast double* %91 to <2 x double>*
store <2 x double> %88, <2 x double>* %92
%93 = or i64 %79, 4
%94 = getelementptr inbounds i64, i64* %0, i64 %93
%95 = bitcast i64* %94 to <2 x i64>*
%96 = load <2 x i64>, <2 x i64>* %95
%97 = getelementptr inbounds i64, i64* %94, i64 2
%98 = bitcast i64* %97 to <2 x i64>*
%99 = load <2 x i64>, <2 x i64>* %98
%100 = uitofp <2 x i64> %96 to <2 x double>
%101 = uitofp <2 x i64> %99 to <2 x double>
%102 = getelementptr inbounds double, double* %4, i64 %93
%103 = bitcast double* %102 to <2 x double>*
store <2 x double> %100, <2 x double>* %103
%104 = getelementptr inbounds double, double* %102, i64 2
%105 = bitcast double* %104 to <2 x double>*
store <2 x double> %101, <2 x double>* %105
%106 = add nuw i64 %79, 8
%107 = add i64 %80, 2
%108 = icmp eq i64 %107, %77
br i1 %108, label %109, label %78
109: ; preds = %78, %69
%110 = phi i64 [ 0, %69 ], [ %106, %78 ]
%111 = icmp eq i64 %74, 0
br i1 %111, label %125, label %112
112: ; preds = %109
%113 = getelementptr inbounds i64, i64* %0, i64 %110
%114 = bitcast i64* %113 to <2 x i64>*
%115 = load <2 x i64>, <2 x i64>* %114
%116 = getelementptr inbounds i64, i64* %113, i64 2
%117 = bitcast i64* %116 to <2 x i64>*
%118 = load <2 x i64>, <2 x i64>* %117
%119 = uitofp <2 x i64> %115 to <2 x double>
%120 = uitofp <2 x i64> %118 to <2 x double>
%121 = getelementptr inbounds double, double* %4, i64 %110
%122 = bitcast double* %121 to <2 x double>*
store <2 x double> %119, <2 x double>* %122
%123 = getelementptr inbounds double, double* %121, i64 2
%124 = bitcast double* %123 to <2 x double>*
store <2 x double> %120, <2 x double>* %124
br label %125
125: ; preds = %109, %112
%126 = icmp eq i64 %70, %1
br i1 %126, label %147, label %127
127: ; preds = %67, %125
%128 = phi i64 [ 0, %67 ], [ %70, %125 ]
br label %129
129: ; preds = %127, %129
%130 = phi i64 [ %135, %129 ], [ %128, %127 ]
%131 = getelementptr inbounds i64, i64* %0, i64 %130
%132 = load i64, i64* %131
%133 = uitofp i64 %132 to double
%134 = getelementptr inbounds double, double* %4, i64 %130
store double %133, double* %134
%135 = add nuw nsw i64 %130, 1
%136 = icmp eq i64 %135, %1
br i1 %136, label %147, label %129
137: ; preds = %64, %137
%138 = phi i64 [ %145, %137 ], [ %65, %64 ]
%139 = getelementptr inbounds i64, i64* %0, i64 %138
%140 = load i64, i64* %139
%141 = uitofp i64 %140 to double
%142 = fsub double %141, %3
%143 = fdiv double %142, %2
%144 = getelementptr inbounds double, double* %4, i64 %138
store double %143, double* %144
%145 = add nuw nsw i64 %138, 1
%146 = icmp eq i64 %145, %1
br i1 %146, label %147, label %137
147: ; preds = %137, %129, %62, %125, %11, %66
%148 = load i32, i32* %5
ret i32 %148
}