int ffr8fi4(double *input, long ntodo, double scale, double zero, int *output,
            int *status) {
  long ii;
  double dvalue;
  if (scale == 1. && zero == 0.) {
    for (ii = 0; ii < ntodo; ii++) {
      if (input[ii] < -2147483648.49) {
        *status = -11;
        output[ii] = (-2147483647 - 1);
      } else if (input[ii] > 2147483647.49) {
        *status = -11;
        output[ii] = 2147483647;
      } else
        output[ii] = (int)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):
<ffr8fi4>:
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    1217 <ffr8fi4+0x117>
jp     1217 <ffr8fi4+0x117>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1217 <ffr8fi4+0x117>
jp     1217 <ffr8fi4+0x117>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1212 <ffr8fi4+0x112>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  0xe8f(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    11a2 <ffr8fi4+0xa2>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x80000000,(%rax,%rcx,4)
jmp    11fc <ffr8fi4+0xfc>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe59(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11df <ffr8fi4+0xdf>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x7fffffff,(%rax,%rcx,4)
jmp    11f7 <ffr8fi4+0xf7>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    11fc <ffr8fi4+0xfc>
jmp    1201 <ffr8fi4+0x101>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    115b <ffr8fi4+0x5b>
jmp    1327 <ffr8fi4+0x227>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1322 <ffr8fi4+0x222>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xdb7(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    127a <ffr8fi4+0x17a>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x80000000,(%rax,%rcx,4)
jmp    130c <ffr8fi4+0x20c>
movsd  -0x40(%rbp),%xmm0
movsd  0xd89(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12af <ffr8fi4+0x1af>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movl   $0x7fffffff,(%rax,%rcx,4)
jmp    1307 <ffr8fi4+0x207>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     12e2 <ffr8fi4+0x1e2>
movsd  0xd4f(%rip),%xmm0
addsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1302 <ffr8fi4+0x202>
movsd  -0x40(%rbp),%xmm0
movsd  0xd29(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    1307 <ffr8fi4+0x207>
jmp    130c <ffr8fi4+0x20c>
jmp    1311 <ffr8fi4+0x211>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    121f <ffr8fi4+0x11f>
jmp    1327 <ffr8fi4+0x227>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<ffr8fi4>:
ucomisd 0xef8(%rip),%xmm0
jne    1187 <ffr8fi4+0x87>
jp     1187 <ffr8fi4+0x87>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    1187 <ffr8fi4+0x87>
jp     1187 <ffr8fi4+0x87>
test   %rsi,%rsi
jle    1217 <ffr8fi4+0x117>
xor    %eax,%eax
movsd  0xedd(%rip),%xmm0
movsd  0xedd(%rip),%xmm1
jmp    115d <ffr8fi4+0x5d>
cs nopw 0x0(%rax,%rax,1)
nop
movl   $0xfffffff5,(%rcx)
mov    $0x80000000,%r8d
mov    %r8d,(%rdx,%rax,4)
add    $0x1,%rax
cmp    %rax,%rsi
je     1217 <ffr8fi4+0x117>
movsd  (%rdi,%rax,8),%xmm2
ucomisd %xmm2,%xmm0
ja     1140 <ffr8fi4+0x40>
ucomisd %xmm1,%xmm2
jbe    1180 <ffr8fi4+0x80>
movl   $0xfffffff5,(%rcx)
mov    $0x7fffffff,%r8d
jmp    114c <ffr8fi4+0x4c>
nopl   0x0(%rax)
cvttsd2si %xmm2,%r8d
jmp    114c <ffr8fi4+0x4c>
test   %rsi,%rsi
jle    1217 <ffr8fi4+0x117>
xor    %eax,%eax
movsd  0xe6e(%rip),%xmm2
movsd  0xe6e(%rip),%xmm3
xorpd  %xmm4,%xmm4
movsd  0xe72(%rip),%xmm5
movsd  0xe62(%rip),%xmm6
jmp    11d9 <ffr8fi4+0xd9>
nopl   0x0(%rax,%rax,1)
movl   $0xfffffff5,(%rcx)
mov    $0x80000000,%r8d
mov    %r8d,(%rdx,%rax,4)
add    $0x1,%rax
cmp    %rax,%rsi
je     1217 <ffr8fi4+0x117>
movsd  (%rdi,%rax,8),%xmm7
subsd  %xmm1,%xmm7
divsd  %xmm0,%xmm7
ucomisd %xmm7,%xmm2
ja     11c0 <ffr8fi4+0xc0>
ucomisd %xmm3,%xmm7
jbe    1200 <ffr8fi4+0x100>
movl   $0xfffffff5,(%rcx)
mov    $0x7fffffff,%r8d
jmp    11cc <ffr8fi4+0xcc>
ucomisd %xmm4,%xmm7
jae    120c <ffr8fi4+0x10c>
addsd  %xmm5,%xmm7
jmp    1210 <ffr8fi4+0x110>
addsd  %xmm6,%xmm7
cvttsd2si %xmm7,%r8d
jmp    11cc <ffr8fi4+0xcc>
mov    (%rcx),%eax
ret

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

LLVM IR (-O2):
define i32 @ffr8fi4(double* 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 %29, label %53
12: ; preds = %6
br i1 %10, label %13, label %53
13: ; preds = %12, %24
%14 = phi i64 [ %27, %24 ], [ 0, %12 ]
%15 = getelementptr inbounds double, double* %0, i64 %14
%16 = load double, double* %15
%17 = fcmp olt double %16, 0xC1E00000000FAE14
br i1 %17, label %18, label %19
18: ; preds = %13
store i32 -11, i32* %5
br label %24
19: ; preds = %13
%20 = fcmp ogt double %16, 0x41DFFFFFFFDF5C29
br i1 %20, label %21, label %22
21: ; preds = %19
store i32 -11, i32* %5
br label %24
22: ; preds = %19
%23 = fptosi double %16 to i32
br label %24
24: ; preds = %18, %22, %21
%25 = phi i32 [ -2147483648, %18 ], [ %23, %22 ], [ 2147483647, %21 ]
%26 = getelementptr inbounds i32, i32* %4, i64 %14
store i32 %25, i32* %26
%27 = add nuw nsw i64 %14, 1
%28 = icmp eq i64 %27, %1
br i1 %28, label %53, label %13
29: ; preds = %11, %48
%30 = phi i64 [ %51, %48 ], [ 0, %11 ]
%31 = getelementptr inbounds double, double* %0, i64 %30
%32 = load double, double* %31
%33 = fsub double %32, %3
%34 = fdiv double %33, %2
%35 = fcmp olt double %34, 0xC1E00000000FAE14
br i1 %35, label %36, label %37
36: ; preds = %29
store i32 -11, i32* %5
br label %48
37: ; preds = %29
%38 = fcmp ogt double %34, 0x41DFFFFFFFDF5C29
br i1 %38, label %39, label %40
39: ; preds = %37
store i32 -11, i32* %5
br label %48
40: ; preds = %37
%41 = fcmp ult double %34, 0.000000e+00
br i1 %41, label %45, label %42
42: ; preds = %40
%43 = fadd double %34, 5.000000e-01
%44 = fptosi double %43 to i32
br label %48
45: ; preds = %40
%46 = fadd double %34, -5.000000e-01
%47 = fptosi double %46 to i32
br label %48
48: ; preds = %36, %42, %45, %39
%49 = phi i32 [ -2147483648, %36 ], [ %44, %42 ], [ %47, %45 ], [ 2147483647, %39 ]
%50 = getelementptr inbounds i32, i32* %4, i64 %30
store i32 %49, i32* %50
%51 = add nuw nsw i64 %30, 1
%52 = icmp eq i64 %51, %1
br i1 %52, label %53, label %29
53: ; preds = %48, %24, %11, %12
%54 = load i32, i32* %5
ret i32 %54
}