void gradientWaveformToShort(const float *gradientsFloat, int points,
                             float maxGradientAmplitudeScanner,
                             short *gradientsShort) {
  float gaussPerCMTo16bit = 32766 / maxGradientAmplitudeScanner;
  int n = 0;
  for (n = 0; n < points - 1; n++)
    gradientsShort[n] = (short)(gradientsFloat[n] * gaussPerCMTo16bit) & 0xfffe;
  gradientsShort[n] = (short)(gradientsFloat[n] * gaussPerCMTo16bit) | 0x0001;
}


x86 assembly (clang, -O0):
<gradientWaveformToShort>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movss  0xee4(%rip),%xmm0
divss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1182 <gradientWaveformToShort+0x82>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x1c(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
cwtl
and    $0xfffe,%eax
mov    %ax,%dx
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1134 <gradientWaveformToShort+0x34>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x1c(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
cwtl
or     $0x1,%eax
mov    %ax,%dx
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
pop    %rbp
ret

x86 assembly (clang, -O2):
<gradientWaveformToShort>:
movss  0xef8(%rip),%xmm1
divss  %xmm0,%xmm1
cmp    $0x2,%esi
jl     1123 <gradientWaveformToShort+0x23>
add    $0xffffffff,%esi
mov    %esi,%r9d
cmp    $0x8,%esi
jae    112b <gradientWaveformToShort+0x2b>
xor    %ecx,%ecx
jmp    1220 <gradientWaveformToShort+0x120>
xor    %r9d,%r9d
jmp    123f <gradientWaveformToShort+0x13f>
mov    %r9d,%ecx
and    $0xfffffff8,%ecx
movaps %xmm1,%xmm0
shufps $0x0,%xmm1,%xmm0
lea    -0x8(%rcx),%rax
mov    %rax,%r8
shr    $0x3,%r8
add    $0x1,%r8
test   %rax,%rax
je     1252 <gradientWaveformToShort+0x152>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
xor    %esi,%esi
movdqa 0xeaf(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups (%rdi,%rsi,4),%xmm3
movups 0x10(%rdi,%rsi,4),%xmm4
mulps  %xmm0,%xmm3
mulps  %xmm0,%xmm4
cvttps2dq %xmm3,%xmm3
packssdw %xmm3,%xmm3
cvttps2dq %xmm4,%xmm4
packssdw %xmm4,%xmm4
pand   %xmm2,%xmm3
pand   %xmm2,%xmm4
punpcklqdq %xmm4,%xmm3
movdqu %xmm3,(%rdx,%rsi,2)
movups 0x20(%rdi,%rsi,4),%xmm3
movups 0x30(%rdi,%rsi,4),%xmm4
mulps  %xmm0,%xmm3
mulps  %xmm0,%xmm4
cvttps2dq %xmm3,%xmm3
packssdw %xmm3,%xmm3
cvttps2dq %xmm4,%xmm4
packssdw %xmm4,%xmm4
pand   %xmm2,%xmm3
pand   %xmm2,%xmm4
punpcklqdq %xmm4,%xmm3
movdqu %xmm3,0x10(%rdx,%rsi,2)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rax
jne    1170 <gradientWaveformToShort+0x70>
test   $0x1,%r8b
je     121a <gradientWaveformToShort+0x11a>
movups (%rdi,%rsi,4),%xmm2
movups 0x10(%rdi,%rsi,4),%xmm3
mulps  %xmm0,%xmm2
mulps  %xmm0,%xmm3
cvttps2dq %xmm2,%xmm0
packssdw %xmm0,%xmm0
cvttps2dq %xmm3,%xmm2
packssdw %xmm2,%xmm2
movdqa 0xe07(%rip),%xmm3
pand   %xmm3,%xmm0
pand   %xmm3,%xmm2
punpcklqdq %xmm2,%xmm0
movdqu %xmm0,(%rdx,%rsi,2)
cmp    %r9,%rcx
je     123f <gradientWaveformToShort+0x13f>
nop
movss  (%rdi,%rcx,4),%xmm0
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
and    $0xfffe,%eax
mov    %ax,(%rdx,%rcx,2)
add    $0x1,%rcx
cmp    %rcx,%r9
jne    1220 <gradientWaveformToShort+0x120>
mulss  (%rdi,%r9,4),%xmm1
cvttss2si %xmm1,%eax
or     $0x1,%eax
mov    %ax,(%rdx,%r9,2)
ret
xor    %esi,%esi
test   $0x1,%r8b
jne    11e2 <gradientWaveformToShort+0xe2>
jmp    121a <gradientWaveformToShort+0x11a>

LLVM IR (-O0):
define void @gradientWaveformToShort(float* %0, i32 %1, float %2, i16* %3) {
%5 = alloca float*
%6 = alloca i32
%7 = alloca float
%8 = alloca i16*
%9 = alloca float
%10 = alloca i32
store float* %0, float** %5
store i32 %1, i32* %6
store float %2, float* %7
store i16* %3, i16** %8
%11 = load float, float* %7
%12 = fdiv float 3.276600e+04, %11
store float %12, float* %9
store i32 0, i32* %10
store i32 0, i32* %10
br label %13
13: ; preds = %34, %4
%14 = load i32, i32* %10
%15 = load i32, i32* %6
%16 = sub nsw i32 %15, 1
%17 = icmp slt i32 %14, %16
br i1 %17, label %18, label %37
18: ; preds = %13
%19 = load float*, float** %5
%20 = load i32, i32* %10
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds float, float* %19, i64 %21
%23 = load float, float* %22
%24 = load float, float* %9
%25 = fmul float %23, %24
%26 = fptosi float %25 to i16
%27 = sext i16 %26 to i32
%28 = and i32 %27, 65534
%29 = trunc i32 %28 to i16
%30 = load i16*, i16** %8
%31 = load i32, i32* %10
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds i16, i16* %30, i64 %32
store i16 %29, i16* %33
br label %34
34: ; preds = %18
%35 = load i32, i32* %10
%36 = add nsw i32 %35, 1
store i32 %36, i32* %10
br label %13
37: ; preds = %13
%38 = load float*, float** %5
%39 = load i32, i32* %10
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
%42 = load float, float* %41
%43 = load float, float* %9
%44 = fmul float %42, %43
%45 = fptosi float %44 to i16
%46 = sext i16 %45 to i32
%47 = or i32 %46, 1
%48 = trunc i32 %47 to i16
%49 = load i16*, i16** %8
%50 = load i32, i32* %10
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i16, i16* %49, i64 %51
store i16 %48, i16* %52
ret void
}

LLVM IR (-O2):
define void @gradientWaveformToShort(float* nocapture %0, i32 %1, float %2, i16* nocapture writeonly %3) {
%5 = fdiv float 3.276600e+04, %2
%6 = add i32 %1, -1
%7 = icmp sgt i32 %1, 1
br i1 %7, label %8, label %99
8: ; preds = %4
%9 = zext i32 %6 to i64
%10 = icmp ult i32 %6, 8
br i1 %10, label %85, label %11
11: ; preds = %8
%12 = and i64 %9, 4294967288
%13 = insertelement <4 x float> poison, float %5, i64 0
%14 = shufflevector <4 x float> %13, <4 x float> poison, <4 x i32> zeroinitializer
%15 = insertelement <4 x float> poison, float %5, i64 0
%16 = shufflevector <4 x float> %15, <4 x float> poison, <4 x i32> zeroinitializer
%17 = add nsw i64 %12, -8
%18 = lshr exact i64 %17, 3
%19 = add nuw nsw i64 %18, 1
%20 = and i64 %19, 1
%21 = icmp eq i64 %17, 0
br i1 %21, label %63, label %22
22: ; preds = %11
%23 = and i64 %19, 4611686018427387902
br label %24
24: ; preds = %24, %22
%25 = phi i64 [ 0, %22 ], [ %60, %24 ]
%26 = phi i64 [ 0, %22 ], [ %61, %24 ]
%27 = getelementptr inbounds float, float* %0, i64 %25
%28 = bitcast float* %27 to <4 x float>*
%29 = load <4 x float>, <4 x float>* %28
%30 = getelementptr inbounds float, float* %27, i64 4
%31 = bitcast float* %30 to <4 x float>*
%32 = load <4 x float>, <4 x float>* %31
%33 = fmul <4 x float> %14, %29
%34 = fmul <4 x float> %16, %32
%35 = fptosi <4 x float> %33 to <4 x i16>
%36 = fptosi <4 x float> %34 to <4 x i16>
%37 = and <4 x i16> %35, <i16 -2, i16 -2, i16 -2, i16 -2>
%38 = and <4 x i16> %36, <i16 -2, i16 -2, i16 -2, i16 -2>
%39 = getelementptr inbounds i16, i16* %3, i64 %25
%40 = bitcast i16* %39 to <4 x i16>*
store <4 x i16> %37, <4 x i16>* %40
%41 = getelementptr inbounds i16, i16* %39, i64 4
%42 = bitcast i16* %41 to <4 x i16>*
store <4 x i16> %38, <4 x i16>* %42
%43 = or i64 %25, 8
%44 = getelementptr inbounds float, float* %0, i64 %43
%45 = bitcast float* %44 to <4 x float>*
%46 = load <4 x float>, <4 x float>* %45
%47 = getelementptr inbounds float, float* %44, i64 4
%48 = bitcast float* %47 to <4 x float>*
%49 = load <4 x float>, <4 x float>* %48
%50 = fmul <4 x float> %14, %46
%51 = fmul <4 x float> %16, %49
%52 = fptosi <4 x float> %50 to <4 x i16>
%53 = fptosi <4 x float> %51 to <4 x i16>
%54 = and <4 x i16> %52, <i16 -2, i16 -2, i16 -2, i16 -2>
%55 = and <4 x i16> %53, <i16 -2, i16 -2, i16 -2, i16 -2>
%56 = getelementptr inbounds i16, i16* %3, i64 %43
%57 = bitcast i16* %56 to <4 x i16>*
store <4 x i16> %54, <4 x i16>* %57
%58 = getelementptr inbounds i16, i16* %56, i64 4
%59 = bitcast i16* %58 to <4 x i16>*
store <4 x i16> %55, <4 x i16>* %59
%60 = add nuw i64 %25, 16
%61 = add i64 %26, 2
%62 = icmp eq i64 %61, %23
br i1 %62, label %63, label %24
63: ; preds = %24, %11
%64 = phi i64 [ 0, %11 ], [ %60, %24 ]
%65 = icmp eq i64 %20, 0
br i1 %65, label %83, label %66
66: ; preds = %63
%67 = getelementptr inbounds float, float* %0, i64 %64
%68 = bitcast float* %67 to <4 x float>*
%69 = load <4 x float>, <4 x float>* %68
%70 = getelementptr inbounds float, float* %67, i64 4
%71 = bitcast float* %70 to <4 x float>*
%72 = load <4 x float>, <4 x float>* %71
%73 = fmul <4 x float> %14, %69
%74 = fmul <4 x float> %16, %72
%75 = fptosi <4 x float> %73 to <4 x i16>
%76 = fptosi <4 x float> %74 to <4 x i16>
%77 = and <4 x i16> %75, <i16 -2, i16 -2, i16 -2, i16 -2>
%78 = and <4 x i16> %76, <i16 -2, i16 -2, i16 -2, i16 -2>
%79 = getelementptr inbounds i16, i16* %3, i64 %64
%80 = bitcast i16* %79 to <4 x i16>*
store <4 x i16> %77, <4 x i16>* %80
%81 = getelementptr inbounds i16, i16* %79, i64 4
%82 = bitcast i16* %81 to <4 x i16>*
store <4 x i16> %78, <4 x i16>* %82
br label %83
83: ; preds = %63, %66
%84 = icmp eq i64 %12, %9
br i1 %84, label %97, label %85
85: ; preds = %8, %83
%86 = phi i64 [ 0, %8 ], [ %12, %83 ]
br label %87
87: ; preds = %85, %87
%88 = phi i64 [ %95, %87 ], [ %86, %85 ]
%89 = getelementptr inbounds float, float* %0, i64 %88
%90 = load float, float* %89
%91 = fmul float %5, %90
%92 = fptosi float %91 to i16
%93 = and i16 %92, -2
%94 = getelementptr inbounds i16, i16* %3, i64 %88
store i16 %93, i16* %94
%95 = add nuw nsw i64 %88, 1
%96 = icmp eq i64 %95, %9
br i1 %96, label %97, label %87
97: ; preds = %87, %83
%98 = zext i32 %6 to i64
br label %99
99: ; preds = %4, %97
%100 = phi i64 [ %98, %97 ], [ 0, %4 ]
%101 = getelementptr inbounds float, float* %0, i64 %100
%102 = load float, float* %101
%103 = fmul float %5, %102
%104 = fptosi float %103 to i16
%105 = or i16 %104, 1
%106 = getelementptr inbounds i16, i16* %3, i64 %100
store i16 %105, i16* %106
ret void
}