long fl2sh(n, x, iy, half_lsb, mask)
long n;
float *x;
short *iy;
double half_lsb;
short mask;
{
  register long iOvrFlw, k;
  register double y;
  iOvrFlw = 0;
  if (half_lsb == 0.0) {
    for (k = 0; k < n; k++) {
      y = x[k] * 32768;
      if (y > 32767.0) {
        y = 32767.0;
        iOvrFlw += 1;
      } else if (y < -32768.0) {
        y = -32768.0;
        iOvrFlw += 1;
      }
      iy[k] = (short)(y);
      iy[k] &= mask;
    }
  } else {
    for (k = 0; k < n; k++) {
      y = x[k] * 32768;
      if (y >= 0.0)
        y = y + half_lsb;
      else
        y = y - half_lsb;
      if (y > 32767.0) {
        y = 32767.0;
        iOvrFlw += 1;
      } else if (y < -32768.0) {
        y = -32768.0;
        iOvrFlw += 1;
      }
      if (y >= 0.0) {
        iy[k] = (short)(y);
        iy[k] &= mask;
      } else {
        iy[k] = (short)(-y);
        iy[k] &= mask;
        iy[k] = -iy[k];
      }
    }
  }
  return iOvrFlw;
}


x86 assembly (clang, -O0):
<fl2sh>:
push   %rbp
mov    %rsp,%rbp
mov    %cx,%ax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %ax,-0x22(%rbp)
movq   $0x0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    121a <fl2sh+0x11a>
jp     121a <fl2sh+0x11a>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    1215 <fl2sh+0x115>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0xe9e(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  0xe8b(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11a5 <fl2sh+0xa5>
movsd  0xe79(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11d6 <fl2sh+0xd6>
movsd  0xe63(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    11d1 <fl2sh+0xd1>
movsd  0xe50(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11d6 <fl2sh+0xd6>
movsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
movswl -0x22(%rbp),%esi
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movswl (%rax,%rcx,2),%edx
and    %esi,%edx
mov    %dx,(%rax,%rcx,2)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1144 <fl2sh+0x44>
jmp    13a1 <fl2sh+0x2a1>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    139c <fl2sh+0x29c>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0xdc0(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1274 <fl2sh+0x174>
movsd  -0x40(%rbp),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1283 <fl2sh+0x183>
movsd  -0x40(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  0xd78(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12b8 <fl2sh+0x1b8>
movsd  0xd66(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    12e9 <fl2sh+0x1e9>
movsd  0xd50(%rip),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    12e4 <fl2sh+0x1e4>
movsd  0xd3d(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    12e9 <fl2sh+0x1e9>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     132e <fl2sh+0x22e>
movsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
movswl -0x22(%rbp),%esi
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movswl (%rax,%rcx,2),%edx
and    %esi,%edx
mov    %dx,(%rax,%rcx,2)
jmp    1386 <fl2sh+0x286>
movaps 0xceb(%rip),%xmm1
movsd  -0x40(%rbp),%xmm0
pxor   %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
movswl -0x22(%rbp),%esi
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movswl (%rax,%rcx,2),%edx
and    %esi,%edx
mov    %dx,(%rax,%rcx,2)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movswl (%rax,%rcx,2),%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %ax,%dx
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    138b <fl2sh+0x28b>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1222 <fl2sh+0x122>
jmp    13a1 <fl2sh+0x2a1>
mov    -0x30(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<fl2sh>:
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1196 <fl2sh+0x96>
jp     1196 <fl2sh+0x96>
test   %rdi,%rdi
jle    125a <fl2sh+0x15a>
xor    %r8d,%r8d
movss  0xee8(%rip),%xmm0
movss  0xee4(%rip),%xmm1
movsd  0xee8(%rip),%xmm2
movss  0xed8(%rip),%xmm3
movsd  0xee0(%rip),%xmm4
xor    %eax,%eax
jmp    1171 <fl2sh+0x71>
nopl   0x0(%rax)
add    $0x1,%rax
movaps %xmm2,%xmm5
cvttsd2si %xmm5,%r9d
and    %ecx,%r9d
mov    %r9w,(%rdx,%r8,2)
add    $0x1,%r8
cmp    %r8,%rdi
je     1259 <fl2sh+0x159>
movss  (%rsi,%r8,4),%xmm5
mulss  %xmm0,%xmm5
ucomiss %xmm1,%xmm5
ja     1150 <fl2sh+0x50>
ucomiss %xmm5,%xmm3
jbe    1190 <fl2sh+0x90>
add    $0x1,%rax
movaps %xmm4,%xmm5
jmp    1157 <fl2sh+0x57>
xchg   %ax,%ax
cvtss2sd %xmm5,%xmm5
jmp    1157 <fl2sh+0x57>
test   %rdi,%rdi
jle    125a <fl2sh+0x15a>
movapd 0xe59(%rip),%xmm6
xorpd  %xmm0,%xmm6
xor    %r8d,%r8d
movss  0xe5a(%rip),%xmm2
xorps  %xmm3,%xmm3
movsd  0xe5f(%rip),%xmm7
movsd  0xe5f(%rip),%xmm5
movapd 0xe2e(%rip),%xmm9
xorpd  %xmm8,%xmm8
xor    %eax,%eax
jmp    11fe <fl2sh+0xfe>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
movapd %xmm7,%xmm4
cvttsd2si %xmm4,%r9d
and    %ecx,%r9d
mov    %r9w,(%rdx,%r8,2)
add    $0x1,%r8
cmp    %r8,%rdi
je     1259 <fl2sh+0x159>
movss  (%rsi,%r8,4),%xmm4
mulss  %xmm2,%xmm4
ucomiss %xmm3,%xmm4
movapd %xmm6,%xmm1
jb     1215 <fl2sh+0x115>
movapd %xmm0,%xmm1
cvtss2sd %xmm4,%xmm4
addsd  %xmm1,%xmm4
ucomisd %xmm7,%xmm4
ja     11e0 <fl2sh+0xe0>
ucomisd %xmm4,%xmm5
jbe    1240 <fl2sh+0x140>
add    $0x1,%rax
movapd %xmm5,%xmm4
jmp    1247 <fl2sh+0x147>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
ucomisd %xmm8,%xmm4
jae    11e8 <fl2sh+0xe8>
xorpd  %xmm9,%xmm4
cvttsd2si %xmm4,%r9d
and    %ecx,%r9d
neg    %r9d
jmp    11f0 <fl2sh+0xf0>
ret
xor    %eax,%eax
ret

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

LLVM IR (-O2):
define i64 @fl2sh(i64 %0, float* nocapture %1, i16* nocapture writeonly %2, double %3, i32 %4) {
%6 = fcmp oeq double %3, 0.000000e+00
%7 = icmp sgt i64 %0, 0
br i1 %6, label %12, label %8
8: ; preds = %5
br i1 %7, label %9, label %74
9: ; preds = %8
%10 = fneg double %3
%11 = trunc i32 %4 to i16
br label %37
12: ; preds = %5
br i1 %7, label %13, label %74
13: ; preds = %12
%14 = trunc i32 %4 to i16
br label %15
15: ; preds = %13, %29
%16 = phi i64 [ 0, %13 ], [ %35, %29 ]
%17 = phi i64 [ 0, %13 ], [ %30, %29 ]
%18 = getelementptr inbounds float, float* %1, i64 %16
%19 = load float, float* %18
%20 = fmul float %19, 3.276800e+04
%21 = fcmp ogt float %20, 3.276700e+04
br i1 %21, label %22, label %24
22: ; preds = %15
%23 = add nsw i64 %17, 1
br label %29
24: ; preds = %15
%25 = fpext float %20 to double
%26 = fcmp olt float %20, -3.276800e+04
br i1 %26, label %27, label %29
27: ; preds = %24
%28 = add nsw i64 %17, 1
br label %29
29: ; preds = %24, %27, %22
%30 = phi i64 [ %23, %22 ], [ %28, %27 ], [ %17, %24 ]
%31 = phi double [ 3.276700e+04, %22 ], [ -3.276800e+04, %27 ], [ %25, %24 ]
%32 = fptosi double %31 to i16
%33 = getelementptr inbounds i16, i16* %2, i64 %16
%34 = and i16 %32, %14
store i16 %34, i16* %33
%35 = add nuw nsw i64 %16, 1
%36 = icmp eq i64 %35, %0
br i1 %36, label %74, label %15
37: ; preds = %9, %70
%38 = phi i64 [ 0, %9 ], [ %72, %70 ]
%39 = phi i64 [ 0, %9 ], [ %71, %70 ]
%40 = getelementptr inbounds float, float* %1, i64 %38
%41 = load float, float* %40
%42 = fmul float %41, 3.276800e+04
%43 = fpext float %42 to double
%44 = fcmp ult float %42, 0.000000e+00
%45 = select i1 %44, double %10, double %3
%46 = fadd double %45, %43
%47 = fcmp ogt double %46, 3.276700e+04
br i1 %47, label %48, label %50
48: ; preds = %37
%49 = add nsw i64 %39, 1
br label %56
50: ; preds = %37
%51 = fcmp olt double %46, -3.276800e+04
br i1 %51, label %52, label %54
52: ; preds = %50
%53 = add nsw i64 %39, 1
br label %62
54: ; preds = %50
%55 = fcmp ult double %46, 0.000000e+00
br i1 %55, label %62, label %56
56: ; preds = %48, %54
%57 = phi double [ 3.276700e+04, %48 ], [ %46, %54 ]
%58 = phi i64 [ %49, %48 ], [ %39, %54 ]
%59 = fptosi double %57 to i16
%60 = getelementptr inbounds i16, i16* %2, i64 %38
%61 = and i16 %59, %11
store i16 %61, i16* %60
br label %70
62: ; preds = %52, %54
%63 = phi double [ -3.276800e+04, %52 ], [ %46, %54 ]
%64 = phi i64 [ %53, %52 ], [ %39, %54 ]
%65 = fneg double %63
%66 = fptosi double %65 to i16
%67 = getelementptr inbounds i16, i16* %2, i64 %38
%68 = and i16 %66, %11
%69 = sub i16 0, %68
store i16 %69, i16* %67
br label %70
70: ; preds = %56, %62
%71 = phi i64 [ %58, %56 ], [ %64, %62 ]
%72 = add nuw nsw i64 %38, 1
%73 = icmp eq i64 %72, %0
br i1 %73, label %74, label %37
74: ; preds = %70, %29, %8, %12
%75 = phi i64 [ 0, %12 ], [ 0, %8 ], [ %30, %29 ], [ %71, %70 ]
ret i64 %75
}