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


x86 assembly (clang, -O0):
<ffuintfr4>:
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    119c <ffuintfr4+0x9c>
jp     119c <ffuintfr4+0x9c>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    119c <ffuintfr4+0x9c>
jp     119c <ffuintfr4+0x9c>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1197 <ffuintfr4+0x97>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cvtsi2ss %rax,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    115b <ffuintfr4+0x5b>
jmp    11f3 <ffuintfr4+0xf3>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11ee <ffuintfr4+0xee>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cvtsi2sd %rax,%xmm0
subsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11a4 <ffuintfr4+0xa4>
jmp    11f3 <ffuintfr4+0xf3>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<ffuintfr4>:
ucomisd 0xef8(%rip),%xmm0
jne    112f <ffuintfr4+0x2f>
jp     112f <ffuintfr4+0x2f>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    112f <ffuintfr4+0x2f>
jp     112f <ffuintfr4+0x2f>
test   %rsi,%rsi
jle    1457 <ffuintfr4+0x357>
cmp    $0x8,%rsi
jae    114a <ffuintfr4+0x4a>
xor    %r11d,%r11d
jmp    12d0 <ffuintfr4+0x1d0>
test   %rsi,%rsi
jle    1457 <ffuintfr4+0x357>
cmp    $0x4,%rsi
jae    12f0 <ffuintfr4+0x1f0>
xor    %r11d,%r11d
jmp    1430 <ffuintfr4+0x330>
mov    %rsi,%r11
and    $0xfffffffffffffff8,%r11
lea    -0x8(%r11),%r8
mov    %r8,%r9
shr    $0x3,%r9
add    $0x1,%r9
test   %r8,%r8
je     1469 <ffuintfr4+0x369>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r8d,%r8d
movdqa 0xea5(%rip),%xmm0
movdqa 0xead(%rip),%xmm1
movdqa 0xeb5(%rip),%xmm2
movaps 0xebe(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movdqu (%rdi,%r8,4),%xmm4
movdqu 0x10(%rdi,%r8,4),%xmm5
movdqa %xmm4,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrld  $0x10,%xmm4
por    %xmm2,%xmm4
subps  %xmm3,%xmm4
addps  %xmm6,%xmm4
movdqa %xmm5,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrld  $0x10,%xmm5
por    %xmm2,%xmm5
subps  %xmm3,%xmm5
addps  %xmm6,%xmm5
movups %xmm4,(%rdx,%r8,4)
movups %xmm5,0x10(%rdx,%r8,4)
movdqu 0x20(%rdi,%r8,4),%xmm4
movdqu 0x30(%rdi,%r8,4),%xmm5
movdqa %xmm4,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrld  $0x10,%xmm4
por    %xmm2,%xmm4
subps  %xmm3,%xmm4
addps  %xmm6,%xmm4
movdqa %xmm5,%xmm6
pand   %xmm0,%xmm6
por    %xmm1,%xmm6
psrld  $0x10,%xmm5
por    %xmm2,%xmm5
subps  %xmm3,%xmm5
addps  %xmm6,%xmm5
movups %xmm4,0x20(%rdx,%r8,4)
movups %xmm5,0x30(%rdx,%r8,4)
add    $0x10,%r8
add    $0xfffffffffffffffe,%r10
jne    11a0 <ffuintfr4+0xa0>
test   $0x1,%r9b
je     12bb <ffuintfr4+0x1bb>
movdqu (%rdi,%r8,4),%xmm0
movdqu 0x10(%rdi,%r8,4),%xmm1
movdqa 0xdb9(%rip),%xmm2
movdqa %xmm0,%xmm3
pand   %xmm2,%xmm3
movdqa 0xdb9(%rip),%xmm4
por    %xmm4,%xmm3
psrld  $0x10,%xmm0
movdqa 0xdb8(%rip),%xmm5
por    %xmm5,%xmm0
movaps 0xdbd(%rip),%xmm6
subps  %xmm6,%xmm0
addps  %xmm3,%xmm0
pand   %xmm1,%xmm2
por    %xmm4,%xmm2
psrld  $0x10,%xmm1
por    %xmm5,%xmm1
subps  %xmm6,%xmm1
addps  %xmm2,%xmm1
movups %xmm0,(%rdx,%r8,4)
movups %xmm1,0x10(%rdx,%r8,4)
cmp    %rsi,%r11
je     1457 <ffuintfr4+0x357>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    (%rdi,%r11,4),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movss  %xmm0,(%rdx,%r11,4)
add    $0x1,%r11
cmp    %r11,%rsi
jne    12d0 <ffuintfr4+0x1d0>
jmp    1457 <ffuintfr4+0x357>
mov    %rsi,%r11
and    $0xfffffffffffffffc,%r11
movapd %xmm1,%xmm3
unpcklpd %xmm1,%xmm3
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
lea    -0x4(%r11),%r8
mov    %r8,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %r8,%r8
je     145a <ffuintfr4+0x35a>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r8d,%r8d
xorpd  %xmm4,%xmm4
movapd 0xcdb(%rip),%xmm5
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rdi,%r8,4),%xmm6
movapd %xmm6,%xmm7
unpcklps %xmm4,%xmm7
orpd   %xmm5,%xmm7
subpd  %xmm5,%xmm7
unpckhps %xmm4,%xmm6
orpd   %xmm5,%xmm6
subpd  %xmm5,%xmm6
subpd  %xmm3,%xmm6
subpd  %xmm3,%xmm7
divpd  %xmm2,%xmm6
cvtpd2ps %xmm6,%xmm6
divpd  %xmm2,%xmm7
cvtpd2ps %xmm7,%xmm7
unpcklpd %xmm6,%xmm7
movupd %xmm7,(%rdx,%r8,4)
movupd 0x10(%rdi,%r8,4),%xmm6
movapd %xmm6,%xmm7
unpcklps %xmm4,%xmm7
orpd   %xmm5,%xmm7
subpd  %xmm5,%xmm7
unpckhps %xmm4,%xmm6
orpd   %xmm5,%xmm6
subpd  %xmm5,%xmm6
subpd  %xmm3,%xmm6
subpd  %xmm3,%xmm7
divpd  %xmm2,%xmm7
divpd  %xmm2,%xmm6
cvtpd2ps %xmm6,%xmm6
cvtpd2ps %xmm7,%xmm7
unpcklpd %xmm6,%xmm7
movupd %xmm7,0x10(%rdx,%r8,4)
add    $0x8,%r8
add    $0xfffffffffffffffe,%r10
jne    1340 <ffuintfr4+0x240>
test   $0x1,%r9b
je     1428 <ffuintfr4+0x328>
movupd (%rdi,%r8,4),%xmm4
xorpd  %xmm5,%xmm5
movapd %xmm4,%xmm6
unpcklps %xmm5,%xmm6
movapd 0xc1d(%rip),%xmm7
orpd   %xmm7,%xmm6
subpd  %xmm7,%xmm6
unpckhps %xmm5,%xmm4
orpd   %xmm7,%xmm4
subpd  %xmm7,%xmm4
subpd  %xmm3,%xmm4
subpd  %xmm3,%xmm6
divpd  %xmm2,%xmm6
divpd  %xmm2,%xmm4
cvtpd2ps %xmm4,%xmm2
cvtpd2ps %xmm6,%xmm3
unpcklpd %xmm2,%xmm3
movupd %xmm3,(%rdx,%r8,4)
cmp    %rsi,%r11
je     1457 <ffuintfr4+0x357>
nopl   (%rax)
mov    (%rdi,%r11,4),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %rax,%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%rdx,%r11,4)
add    $0x1,%r11
cmp    %r11,%rsi
jne    1430 <ffuintfr4+0x330>
mov    (%rcx),%eax
ret
xor    %r8d,%r8d
test   $0x1,%r9b
jne    13da <ffuintfr4+0x2da>
jmp    1428 <ffuintfr4+0x328>
xor    %r8d,%r8d
test   $0x1,%r9b
jne    1252 <ffuintfr4+0x152>
jmp    12bb <ffuintfr4+0x1bb>

LLVM IR (-O0):
define i32 @ffuintfr4(i32* %0, i64 %1, double %2, double %3, float* %4, i32* %5) {
%7 = alloca i32*
%8 = alloca i64
%9 = alloca double
%10 = alloca double
%11 = alloca float*
%12 = alloca i32*
%13 = alloca i64
store i32* %0, i32** %7
store i64 %1, i64* %8
store double %2, double* %9
store double %3, double* %10
store float* %4, float** %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 i32*, i32** %7
%26 = load i64, i64* %13
%27 = getelementptr inbounds i32, i32* %25, i64 %26
%28 = load i32, i32* %27
%29 = uitofp i32 %28 to float
%30 = load float*, float** %11
%31 = load i64, i64* %13
%32 = getelementptr inbounds float, float* %30, i64 %31
store float %29, float* %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 %60
37: ; preds = %16, %6
store i64 0, i64* %13
br label %38
38: ; preds = %56, %37
%39 = load i64, i64* %13
%40 = load i64, i64* %8
%41 = icmp slt i64 %39, %40
br i1 %41, label %42, label %59
42: ; preds = %38
%43 = load i32*, i32** %7
%44 = load i64, i64* %13
%45 = getelementptr inbounds i32, i32* %43, i64 %44
%46 = load i32, i32* %45
%47 = uitofp i32 %46 to double
%48 = load double, double* %10
%49 = fsub double %47, %48
%50 = load double, double* %9
%51 = fdiv double %49, %50
%52 = fptrunc double %51 to float
%53 = load float*, float** %11
%54 = load i64, i64* %13
%55 = getelementptr inbounds float, float* %53, i64 %54
store float %52, float* %55
br label %56
56: ; preds = %42
%57 = load i64, i64* %13
%58 = add nsw i64 %57, 1
store i64 %58, i64* %13
br label %38
59: ; preds = %38
br label %60
60: ; preds = %59, %36
%61 = load i32*, i32** %12
%62 = load i32, i32* %61
ret i32 %62
}

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