#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
unsigned long SHRT_MAX;
short SHRT_MIN;
int ffu4fi2(unsigned long *input, long ntodo, double scale, double zero,
            short *output, int *status) {
  long ii;
  double dvalue;
  if (scale == 1. && zero == 0.) {
    for (ii = 0; ii < ntodo; ii++) {
      if (input[ii] > SHRT_MAX) {
        *status = -11;
        output[ii] = SHRT_MAX;
      } else
        output[ii] = input[ii];
    }
  } else {
    for (ii = 0; ii < ntodo; ii++) {
      dvalue = (input[ii] - zero) / scale;
      if (dvalue < -32768.49) {
        *status = -11;
        output[ii] = SHRT_MIN;
      } else if (dvalue > 32767.49) {
        *status = -11;
        output[ii] = SHRT_MAX;
      } else {
        if (dvalue >= 0)
          output[ii] = (short)(dvalue + .5);
        else
          output[ii] = (short)(dvalue - .5);
      }
    }
  }
  return (*status);
}


x86 assembly (clang, -O0):
<ffu4fi2>:
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    11e3 <ffu4fi2+0xe3>
jp     11e3 <ffu4fi2+0xe3>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11e3 <ffu4fi2+0xe3>
jp     11e3 <ffu4fi2+0xe3>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11de <ffu4fi2+0xde>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x2e5c(%rip),%rcx
cmp    (%rcx),%rax
jbe    11ad <ffu4fi2+0xad>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    0x2e42(%rip),%rax
mov    (%rax),%rax
mov    %ax,%dx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    11c8 <ffu4fi2+0xc8>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %ax,%dx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    11cd <ffu4fi2+0xcd>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    115b <ffu4fi2+0x5b>
jmp    1332 <ffu4fi2+0x232>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    132d <ffu4fi2+0x22d>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movq   (%rax,%rcx,8),%xmm0
movaps 0xe23(%rip),%xmm1
punpckldq %xmm1,%xmm0
movapd 0xe27(%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  0xdc9(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    126f <ffu4fi2+0x16f>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    0x2d85(%rip),%rax
mov    (%rax),%dx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    1317 <ffu4fi2+0x217>
movsd  -0x40(%rbp),%xmm0
movsd  0xd94(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12ae <ffu4fi2+0x1ae>
mov    -0x30(%rbp),%rax
movl   $0xfffffff5,(%rax)
mov    0x2d41(%rip),%rax
mov    (%rax),%rax
mov    %ax,%dx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    1312 <ffu4fi2+0x212>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     12e9 <ffu4fi2+0x1e9>
movsd  -0x40(%rbp),%xmm0
movsd  0xd53(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    130d <ffu4fi2+0x20d>
movsd  -0x40(%rbp),%xmm0
movsd  0xd22(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    1312 <ffu4fi2+0x212>
jmp    1317 <ffu4fi2+0x217>
jmp    131c <ffu4fi2+0x21c>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11eb <ffu4fi2+0xeb>
jmp    1332 <ffu4fi2+0x232>
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<ffu4fi2>:
ucomisd 0xef8(%rip),%xmm0
jne    115a <ffu4fi2+0x5a>
jp     115a <ffu4fi2+0x5a>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    115a <ffu4fi2+0x5a>
jp     115a <ffu4fi2+0x5a>
test   %rsi,%rsi
jle    121a <ffu4fi2+0x11a>
mov    0x2eb0(%rip),%rax
mov    (%rax),%r8
cmp    $0x1,%rsi
jne    121d <ffu4fi2+0x11d>
xor    %eax,%eax
test   $0x1,%sil
je     121a <ffu4fi2+0x11a>
mov    (%rdi,%rax,8),%rsi
cmp    %r8,%rsi
jbe    1153 <ffu4fi2+0x53>
movl   $0xfffffff5,(%rcx)
mov    %r8,%rsi
mov    %si,(%rdx,%rax,2)
mov    (%rcx),%eax
ret
test   %rsi,%rsi
jle    121a <ffu4fi2+0x11a>
mov    0x2e6e(%rip),%rax
movzwl (%rax),%r8d
xor    %eax,%eax
movapd 0xeb8(%rip),%xmm2
movapd 0xec0(%rip),%xmm3
movsd  0xe80(%rip),%xmm4
mov    0x2e51(%rip),%r9
movsd  0xe79(%rip),%xmm5
xorpd  %xmm10,%xmm10
movsd  0xe7b(%rip),%xmm8
movsd  0xe6a(%rip),%xmm9
jmp    11c8 <ffu4fi2+0xc8>
movl   $0xfffffff5,(%rcx)
movzwl (%r9),%r10d
mov    %r10w,(%rdx,%rax,2)
add    $0x1,%rax
cmp    %rax,%rsi
je     121a <ffu4fi2+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     11b0 <ffu4fi2+0xb0>
ucomisd %xmm5,%xmm7
jbe    1200 <ffu4fi2+0x100>
movl   $0xfffffff5,(%rcx)
mov    %r8d,%r10d
jmp    11ba <ffu4fi2+0xba>
nop
ucomisd %xmm10,%xmm7
jae    120e <ffu4fi2+0x10e>
addsd  %xmm8,%xmm7
jmp    1213 <ffu4fi2+0x113>
addsd  %xmm9,%xmm7
cvttsd2si %xmm7,%r10d
jmp    11ba <ffu4fi2+0xba>
mov    (%rcx),%eax
ret
mov    %rsi,%r9
and    $0xfffffffffffffffe,%r9
xor    %eax,%eax
jmp    1243 <ffu4fi2+0x143>
nopl   0x0(%rax,%rax,1)
mov    %r10w,0x2(%rdx,%rax,2)
add    $0x2,%rax
cmp    %rax,%r9
je     1137 <ffu4fi2+0x37>
mov    (%rdi,%rax,8),%r10
cmp    %r8,%r10
jbe    1255 <ffu4fi2+0x155>
movl   $0xfffffff5,(%rcx)
mov    %r8,%r10
mov    %r10w,(%rdx,%rax,2)
mov    0x8(%rdi,%rax,8),%r10
cmp    %r8,%r10
jbe    1230 <ffu4fi2+0x130>
movl   $0xfffffff5,(%rcx)
mov    %r8,%r10
jmp    1230 <ffu4fi2+0x130>

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

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