int ffu4fi4(unsigned long *input, long ntodo, double scale, double zero,
            int *output, int *status) {
  long ii;
  double dvalue;
  if (scale == 1. && zero == 2147483648. && sizeof(long) == 4) {
    for (ii = 0; ii < ntodo; ii++)
      output[ii] = (*(long *)&input[ii]) ^ 0x80000000;
  } else if (scale == 1. && zero == 0.) {
    for (ii = 0; ii < ntodo; ii++) {
      if (input[ii] > 2147483647) {
        *status = -11;
        output[ii] = 2147483647;
      } else
        output[ii] = input[ii];
    }
  } else {
    for (ii = 0; ii < ntodo; ii++) {
      dvalue = (input[ii] - zero) / scale;
      if (dvalue < -2147483648.49) {
        *status = -11;
        output[ii] = (-2147483647 - 1);
      } else if (dvalue > 2147483647.49) {
        *status = -11;
        output[ii] = 2147483647;
      } else {
        if (dvalue >= 0)
          output[ii] = (int)(dvalue + .5);
        else
          output[ii] = (int)(dvalue - .5);
      }
    }
  }
  return (*status);
}


x86 assembly (clang, -O0):
<ffu4fi4>:
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    11b6 <ffu4fi4+0xb6>
jp     11b6 <ffu4fi4+0xb6>
movsd  -0x20(%rbp),%xmm0
movsd  0xec0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    11b6 <ffu4fi4+0xb6>
jp     11b6 <ffu4fi4+0xb6>
xor    %eax,%eax
test   $0x1,%al
jne    1167 <ffu4fi4+0x67>
jmp    11b6 <ffu4fi4+0xb6>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11b1 <ffu4fi4+0xb1>
mov    -0x8(%rbp),%rcx
mov    -0x38(%rbp),%rdx
movabs $0x80000000,%rax
xor    (%rcx,%rdx,8),%rax
mov    %eax,%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    116f <ffu4fi4+0x6f>
jmp    13a0 <ffu4fi4+0x2a0>
movsd  -0x18(%rbp),%xmm0
movsd  0xe3d(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1269 <ffu4fi4+0x169>
jp     1269 <ffu4fi4+0x169>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1269 <ffu4fi4+0x169>
jp     1269 <ffu4fi4+0x169>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1264 <ffu4fi4+0x164>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
cmpq   $0x7fffffff,(%rax,%rcx,8)
7f
jbe    1235 <ffu4fi4+0x135>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x7fffffff,(%rax,%rcx,4)
jmp    124e <ffu4fi4+0x14e>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %eax,%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1253 <ffu4fi4+0x153>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11f3 <ffu4fi4+0xf3>
jmp    139b <ffu4fi4+0x29b>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1396 <ffu4fi4+0x296>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movq   (%rax,%rcx,8),%xmm0
movaps 0xd9d(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xda1(%rip),%xmm1
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xd4b(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    12ee <ffu4fi4+0x1ee>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x80000000,(%rax,%rcx,4)
jmp    1380 <ffu4fi4+0x280>
movsd  -0x40(%rbp),%xmm0
movsd  0xd1d(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1323 <ffu4fi4+0x223>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x7fffffff,(%rax,%rcx,4)
jmp    137b <ffu4fi4+0x27b>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1356 <ffu4fi4+0x256>
movsd  0xce3(%rip),%xmm0
addsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1376 <ffu4fi4+0x276>
movsd  -0x40(%rbp),%xmm0
movsd  0xcbd(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    137b <ffu4fi4+0x27b>
jmp    1380 <ffu4fi4+0x280>
jmp    1385 <ffu4fi4+0x285>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1271 <ffu4fi4+0x171>
jmp    139b <ffu4fi4+0x29b>
jmp    13a0 <ffu4fi4+0x2a0>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<ffu4fi4>:
ucomisd 0xef8(%rip),%xmm0
jne    1155 <ffu4fi4+0x55>
jp     1155 <ffu4fi4+0x55>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    1155 <ffu4fi4+0x55>
jp     1155 <ffu4fi4+0x55>
test   %rsi,%rsi
jle    121a <ffu4fi4+0x11a>
cmp    $0x1,%rsi
jne    121d <ffu4fi4+0x11d>
xor    %eax,%eax
test   $0x1,%sil
je     121a <ffu4fi4+0x11a>
mov    (%rdi,%rax,8),%rsi
cmp    $0x7fffffff,%rsi
jbe    114f <ffu4fi4+0x4f>
movl   $0xfffffff5,(%rcx)
mov    $0x7fffffff,%esi
mov    %esi,(%rdx,%rax,4)
mov    (%rcx),%eax
ret
test   %rsi,%rsi
jle    121a <ffu4fi4+0x11a>
xor    %eax,%eax
movapd 0xec8(%rip),%xmm2
movapd 0xed0(%rip),%xmm3
movsd  0xe90(%rip),%xmm4
movsd  0xe90(%rip),%xmm5
xorpd  %xmm10,%xmm10
movsd  0xe92(%rip),%xmm8
movsd  0xe81(%rip),%xmm9
jmp    11b9 <ffu4fi4+0xb9>
nopl   0x0(%rax)
movl   $0xfffffff5,(%rcx)
mov    $0x80000000,%r8d
mov    %r8d,(%rdx,%rax,4)
add    $0x1,%rax
cmp    %rax,%rsi
je     121a <ffu4fi4+0x11a>
movsd  (%rdi,%rax,8),%xmm6
unpcklps %xmm2,%xmm6
subpd  %xmm3,%xmm6
movapd %xmm6,%xmm7
unpckhpd %xmm6,%xmm7
addsd  %xmm6,%xmm7
subsd  %xmm1,%xmm7
divsd  %xmm0,%xmm7
ucomisd %xmm7,%xmm4
ja     11a0 <ffu4fi4+0xa0>
ucomisd %xmm5,%xmm7
jbe    1200 <ffu4fi4+0x100>
movl   $0xfffffff5,(%rcx)
mov    $0x7fffffff,%r8d
jmp    11ac <ffu4fi4+0xac>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
ucomisd %xmm10,%xmm7
jae    120e <ffu4fi4+0x10e>
addsd  %xmm8,%xmm7
jmp    1213 <ffu4fi4+0x113>
addsd  %xmm9,%xmm7
cvttsd2si %xmm7,%r8d
jmp    11ac <ffu4fi4+0xac>
mov    (%rcx),%eax
ret
mov    %rsi,%r8
and    $0xfffffffffffffffe,%r8
xor    %eax,%eax
jmp    1242 <ffu4fi4+0x142>
nopl   0x0(%rax,%rax,1)
mov    %r9d,0x4(%rdx,%rax,4)
add    $0x2,%rax
cmp    %rax,%r8
je     112d <ffu4fi4+0x2d>
mov    (%rdi,%rax,8),%r9
test   $0xffffffff80000000,%r9
je     125b <ffu4fi4+0x15b>
movl   $0xfffffff5,(%rcx)
mov    $0x7fffffff,%r9d
mov    %r9d,(%rdx,%rax,4)
mov    0x8(%rdi,%rax,8),%r9
cmp    $0x7fffffff,%r9
jbe    1230 <ffu4fi4+0x130>
movl   $0xfffffff5,(%rcx)
mov    $0x7fffffff,%r9d
jmp    1230 <ffu4fi4+0x130>

LLVM IR (-O0):
define i32 @ffu4fi4(i64* %0, i64 %1, double %2, double %3, i32* %4, i32* %5) {
%7 = alloca i64*
%8 = alloca i64
%9 = alloca double
%10 = alloca double
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i64
%14 = alloca double
store i64* %0, i64** %7
store i64 %1, i64* %8
store double %2, double* %9
store double %3, double* %10
store i32* %4, i32** %11
store i32* %5, i32** %12
%15 = load double, double* %9
%16 = fcmp oeq double %15, 1.000000e+00
br i1 %16, label %17, label %40
17: ; preds = %6
%18 = load double, double* %10
%19 = fcmp oeq double %18, 0x41E0000000000000
br i1 %19, label %20, label %40
20: ; preds = %17
br i1 false, label %21, label %40
21: ; preds = %20
store i64 0, i64* %13
br label %22
22: ; preds = %36, %21
%23 = load i64, i64* %13
%24 = load i64, i64* %8
%25 = icmp slt i64 %23, %24
br i1 %25, label %26, label %39
26: ; preds = %22
%27 = load i64*, i64** %7
%28 = load i64, i64* %13
%29 = getelementptr inbounds i64, i64* %27, i64 %28
%30 = load i64, i64* %29
%31 = xor i64 %30, 2147483648
%32 = trunc i64 %31 to i32
%33 = load i32*, i32** %11
%34 = load i64, i64* %13
%35 = getelementptr inbounds i32, i32* %33, i64 %34
store i32 %32, i32* %35
br label %36
36: ; preds = %26
%37 = load i64, i64* %13
%38 = add nsw i64 %37, 1
store i64 %38, i64* %13
br label %22
39: ; preds = %22
br label %131
40: ; preds = %20, %17, %6
%41 = load double, double* %9
%42 = fcmp oeq double %41, 1.000000e+00
br i1 %42, label %43, label %76
43: ; preds = %40
%44 = load double, double* %10
%45 = fcmp oeq double %44, 0.000000e+00
br i1 %45, label %46, label %76
46: ; preds = %43
store i64 0, i64* %13
br label %47
47: ; preds = %72, %46
%48 = load i64, i64* %13
%49 = load i64, i64* %8
%50 = icmp slt i64 %48, %49
br i1 %50, label %51, label %75
51: ; preds = %47
%52 = load i64*, i64** %7
%53 = load i64, i64* %13
%54 = getelementptr inbounds i64, i64* %52, i64 %53
%55 = load i64, i64* %54
%56 = icmp ugt i64 %55, 2147483647
br i1 %56, label %57, label %62
57: ; preds = %51
%58 = load i32*, i32** %12
store i32 -11, i32* %58
%59 = load i32*, i32** %11
%60 = load i64, i64* %13
%61 = getelementptr inbounds i32, i32* %59, i64 %60
store i32 2147483647, i32* %61
br label %71
62: ; preds = %51
%63 = load i64*, i64** %7
%64 = load i64, i64* %13
%65 = getelementptr inbounds i64, i64* %63, i64 %64
%66 = load i64, i64* %65
%67 = trunc i64 %66 to i32
%68 = load i32*, i32** %11
%69 = load i64, i64* %13
%70 = getelementptr inbounds i32, i32* %68, i64 %69
store i32 %67, i32* %70
br label %71
71: ; preds = %62, %57
br label %72
72: ; preds = %71
%73 = load i64, i64* %13
%74 = add nsw i64 %73, 1
store i64 %74, i64* %13
br label %47
75: ; preds = %47
br label %130
76: ; preds = %43, %40
store i64 0, i64* %13
br label %77
77: ; preds = %126, %76
%78 = load i64, i64* %13
%79 = load i64, i64* %8
%80 = icmp slt i64 %78, %79
br i1 %80, label %81, label %129
81: ; preds = %77
%82 = load i64*, i64** %7
%83 = load i64, i64* %13
%84 = getelementptr inbounds i64, i64* %82, i64 %83
%85 = load i64, i64* %84
%86 = uitofp i64 %85 to double
%87 = load double, double* %10
%88 = fsub double %86, %87
%89 = load double, double* %9
%90 = fdiv double %88, %89
store double %90, double* %14
%91 = load double, double* %14
%92 = fcmp olt double %91, 0xC1E00000000FAE14
br i1 %92, label %93, label %98
93: ; preds = %81
%94 = load i32*, i32** %12
store i32 -11, i32* %94
%95 = load i32*, i32** %11
%96 = load i64, i64* %13
%97 = getelementptr inbounds i32, i32* %95, i64 %96
store i32 -2147483648, i32* %97
br label %125
98: ; preds = %81
%99 = load double, double* %14
%100 = fcmp ogt double %99, 0x41DFFFFFFFDF5C29
br i1 %100, label %101, label %106
101: ; preds = %98
%102 = load i32*, i32** %12
store i32 -11, i32* %102
%103 = load i32*, i32** %11
%104 = load i64, i64* %13
%105 = getelementptr inbounds i32, i32* %103, i64 %104
store i32 2147483647, i32* %105
br label %124
106: ; preds = %98
%107 = load double, double* %14
%108 = fcmp oge double %107, 0.000000e+00
br i1 %108, label %109, label %116
109: ; preds = %106
%110 = load double, double* %14
%111 = fadd double %110, 5.000000e-01
%112 = fptosi double %111 to i32
%113 = load i32*, i32** %11
%114 = load i64, i64* %13
%115 = getelementptr inbounds i32, i32* %113, i64 %114
store i32 %112, i32* %115
br label %123
116: ; preds = %106
%117 = load double, double* %14
%118 = fsub double %117, 5.000000e-01
%119 = fptosi double %118 to i32
%120 = load i32*, i32** %11
%121 = load i64, i64* %13
%122 = getelementptr inbounds i32, i32* %120, i64 %121
store i32 %119, i32* %122
br label %123
123: ; preds = %116, %109
br label %124
124: ; preds = %123, %101
br label %125
125: ; preds = %124, %93
br label %126
126: ; preds = %125
%127 = load i64, i64* %13
%128 = add nsw i64 %127, 1
store i64 %128, i64* %13
br label %77
129: ; preds = %77
br label %130
130: ; preds = %129, %75
br label %131
131: ; preds = %130, %39
%132 = load i32*, i32** %12
%133 = load i32, i32* %132
ret i32 %133
}

LLVM IR (-O2):
define i32 @ffu4fi4(i64* nocapture %0, i64 %1, double %2, double %3, i32* 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 %12, label %11
11: ; preds = %6
br i1 %10, label %43, label %81
12: ; preds = %6
br i1 %10, label %13, label %81
13: ; preds = %12
%14 = and i64 %1, 1
%15 = icmp eq i64 %1, 1
br i1 %15, label %68, label %16
16: ; preds = %13
%17 = and i64 %1, -2
br label %18
18: ; preds = %37, %16
%19 = phi i64 [ 0, %16 ], [ %40, %37 ]
%20 = phi i64 [ 0, %16 ], [ %41, %37 ]
%21 = getelementptr inbounds i64, i64* %0, i64 %19
%22 = load i64, i64* %21
%23 = icmp ugt i64 %22, 2147483647
br i1 %23, label %24, label %25
24: ; preds = %18
store i32 -11, i32* %5
br label %27
25: ; preds = %18
%26 = trunc i64 %22 to i32
br label %27
27: ; preds = %24, %25
%28 = phi i32 [ 2147483647, %24 ], [ %26, %25 ]
%29 = getelementptr inbounds i32, i32* %4, i64 %19
store i32 %28, i32* %29
%30 = or i64 %19, 1
%31 = getelementptr inbounds i64, i64* %0, i64 %30
%32 = load i64, i64* %31
%33 = icmp ugt i64 %32, 2147483647
br i1 %33, label %36, label %34
34: ; preds = %27
%35 = trunc i64 %32 to i32
br label %37
36: ; preds = %27
store i32 -11, i32* %5
br label %37
37: ; preds = %36, %34
%38 = phi i32 [ 2147483647, %36 ], [ %35, %34 ]
%39 = getelementptr inbounds i32, i32* %4, i64 %30
store i32 %38, i32* %39
%40 = add nuw nsw i64 %19, 2
%41 = add i64 %20, 2
%42 = icmp eq i64 %41, %17
br i1 %42, label %68, label %18
43: ; preds = %11, %63
%44 = phi i64 [ %66, %63 ], [ 0, %11 ]
%45 = getelementptr inbounds i64, i64* %0, i64 %44
%46 = load i64, i64* %45
%47 = uitofp i64 %46 to double
%48 = fsub double %47, %3
%49 = fdiv double %48, %2
%50 = fcmp olt double %49, 0xC1E00000000FAE14
br i1 %50, label %51, label %52
51: ; preds = %43
store i32 -11, i32* %5
br label %63
52: ; preds = %43
%53 = fcmp ogt double %49, 0x41DFFFFFFFDF5C29
br i1 %53, label %54, label %55
54: ; preds = %52
store i32 -11, i32* %5
br label %63
55: ; preds = %52
%56 = fcmp ult double %49, 0.000000e+00
br i1 %56, label %60, label %57
57: ; preds = %55
%58 = fadd double %49, 5.000000e-01
%59 = fptosi double %58 to i32
br label %63
60: ; preds = %55
%61 = fadd double %49, -5.000000e-01
%62 = fptosi double %61 to i32
br label %63
63: ; preds = %51, %57, %60, %54
%64 = phi i32 [ -2147483648, %51 ], [ %59, %57 ], [ %62, %60 ], [ 2147483647, %54 ]
%65 = getelementptr inbounds i32, i32* %4, i64 %44
store i32 %64, i32* %65
%66 = add nuw nsw i64 %44, 1
%67 = icmp eq i64 %66, %1
br i1 %67, label %81, label %43
68: ; preds = %37, %13
%69 = phi i64 [ 0, %13 ], [ %40, %37 ]
%70 = icmp eq i64 %14, 0
br i1 %70, label %81, label %71
71: ; preds = %68
%72 = getelementptr inbounds i64, i64* %0, i64 %69
%73 = load i64, i64* %72
%74 = icmp ugt i64 %73, 2147483647
br i1 %74, label %77, label %75
75: ; preds = %71
%76 = trunc i64 %73 to i32
br label %78
77: ; preds = %71
store i32 -11, i32* %5
br label %78
78: ; preds = %77, %75
%79 = phi i32 [ 2147483647, %77 ], [ %76, %75 ]
%80 = getelementptr inbounds i32, i32* %4, i64 %69
store i32 %79, i32* %80
br label %81
81: ; preds = %63, %78, %68, %11, %12
%82 = load i32, i32* %5
ret i32 %82
}