float WarpFreq(float fcl, float fcu, float freq, float minFreq, float maxFreq,
               float alpha) {
  if (alpha == 1.0)
    return freq;
  else {
    float scale = 1.0 / alpha;
    float cu = fcu * 2 / (1 + scale);
    float cl = fcl * 2 / (1 + scale);
    float au = (maxFreq - cu * scale) / (maxFreq - cu);
    float al = (cl * scale - minFreq) / (cl - minFreq);
    if (freq > cu)
      return au * (freq - cu) + scale * cu;
    else if (freq < cl)
      return al * (freq - minFreq) + minFreq;
    else
      return scale * freq;
  }
}


x86 assembly (clang, -O0):
<WarpFreq>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
movss  %xmm2,-0x10(%rbp)
movss  %xmm3,-0x14(%rbp)
movss  %xmm4,-0x18(%rbp)
movss  %xmm5,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xecd(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1152 <WarpFreq+0x52>
jp     1152 <WarpFreq+0x52>
movss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    12a6 <WarpFreq+0x1a6>
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe9d(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  0xe94(%rip),%xmm0
mulss  -0xc(%rbp),%xmm0
movss  0xe83(%rip),%xmm1
addss  -0x20(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  0xe71(%rip),%xmm0
mulss  -0x8(%rbp),%xmm0
movss  0xe60(%rip),%xmm1
addss  -0x20(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x18(%rbp),%xmm1
movss  -0x24(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x18(%rbp),%xmm1
subss  -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x28(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x14(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  -0x28(%rbp),%xmm1
subss  -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x10(%rbp),%xmm0
ucomiss -0x24(%rbp),%xmm0
jbe    125e <WarpFreq+0x15e>
movss  -0x2c(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
subss  -0x24(%rbp),%xmm2
movss  -0x20(%rbp),%xmm1
mulss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    12a6 <WarpFreq+0x1a6>
movss  -0x10(%rbp),%xmm1
movss  -0x28(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1297 <WarpFreq+0x197>
movss  -0x30(%rbp),%xmm0
movss  -0x10(%rbp),%xmm2
subss  -0x14(%rbp),%xmm2
movss  -0x14(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    12a6 <WarpFreq+0x1a6>
movss  -0x20(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<WarpFreq>:
ucomiss 0xef9(%rip),%xmm5
jne    110f <WarpFreq+0xf>
jp     110f <WarpFreq+0xf>
movaps %xmm2,%xmm0
ret
movss  0xee9(%rip),%xmm7
movaps %xmm7,%xmm6
divss  %xmm5,%xmm6
addss  %xmm1,%xmm1
addss  %xmm6,%xmm7
divss  %xmm7,%xmm1
ucomiss %xmm1,%xmm2
jbe    1152 <WarpFreq+0x52>
mulss  %xmm1,%xmm6
movaps %xmm4,%xmm0
subss  %xmm6,%xmm0
subss  %xmm1,%xmm4
divss  %xmm4,%xmm0
subss  %xmm1,%xmm2
mulss  %xmm0,%xmm2
addss  %xmm6,%xmm2
movaps %xmm2,%xmm0
ret
addss  %xmm0,%xmm0
divss  %xmm7,%xmm0
ucomiss %xmm2,%xmm0
jbe    117f <WarpFreq+0x7f>
mulss  %xmm0,%xmm6
subss  %xmm3,%xmm6
subss  %xmm3,%xmm0
divss  %xmm0,%xmm6
subss  %xmm3,%xmm2
mulss  %xmm6,%xmm2
addss  %xmm3,%xmm2
movaps %xmm2,%xmm0
ret
mulss  %xmm6,%xmm2
movaps %xmm2,%xmm0
ret

LLVM IR (-O0):
define float @WarpFreq(float %0, float %1, float %2, float %3, float %4, float %5) {
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
store float %0, float* %8
store float %1, float* %9
store float %2, float* %10
store float %3, float* %11
store float %4, float* %12
store float %5, float* %13
%19 = load float, float* %13
%20 = fpext float %19 to double
%21 = fcmp oeq double %20, 1.000000e+00
br i1 %21, label %22, label %24
22: ; preds = %6
%23 = load float, float* %10
store float %23, float* %7
br label %84
24: ; preds = %6
%25 = load float, float* %13
%26 = fpext float %25 to double
%27 = fdiv double 1.000000e+00, %26
%28 = fptrunc double %27 to float
store float %28, float* %14
%29 = load float, float* %9
%30 = fmul float %29, 2.000000e+00
%31 = load float, float* %14
%32 = fadd float 1.000000e+00, %31
%33 = fdiv float %30, %32
store float %33, float* %15
%34 = load float, float* %8
%35 = fmul float %34, 2.000000e+00
%36 = load float, float* %14
%37 = fadd float 1.000000e+00, %36
%38 = fdiv float %35, %37
store float %38, float* %16
%39 = load float, float* %12
%40 = load float, float* %15
%41 = load float, float* %14
%42 = fneg float %40
%43 = call float @llvm.fmuladd.f32(float %42, float %41, float %39)
%44 = load float, float* %12
%45 = load float, float* %15
%46 = fsub float %44, %45
%47 = fdiv float %43, %46
store float %47, float* %17
%48 = load float, float* %16
%49 = load float, float* %14
%50 = load float, float* %11
%51 = fneg float %50
%52 = call float @llvm.fmuladd.f32(float %48, float %49, float %51)
%53 = load float, float* %16
%54 = load float, float* %11
%55 = fsub float %53, %54
%56 = fdiv float %52, %55
store float %56, float* %18
%57 = load float, float* %10
%58 = load float, float* %15
%59 = fcmp ogt float %57, %58
br i1 %59, label %60, label %69
60: ; preds = %24
%61 = load float, float* %17
%62 = load float, float* %10
%63 = load float, float* %15
%64 = fsub float %62, %63
%65 = load float, float* %14
%66 = load float, float* %15
%67 = fmul float %65, %66
%68 = call float @llvm.fmuladd.f32(float %61, float %64, float %67)
store float %68, float* %7
br label %84
69: ; preds = %24
%70 = load float, float* %10
%71 = load float, float* %16
%72 = fcmp olt float %70, %71
br i1 %72, label %73, label %80
73: ; preds = %69
%74 = load float, float* %18
%75 = load float, float* %10
%76 = load float, float* %11
%77 = fsub float %75, %76
%78 = load float, float* %11
%79 = call float @llvm.fmuladd.f32(float %74, float %77, float %78)
store float %79, float* %7
br label %84
80: ; preds = %69
%81 = load float, float* %14
%82 = load float, float* %10
%83 = fmul float %81, %82
store float %83, float* %7
br label %84
84: ; preds = %80, %73, %60, %22
%85 = load float, float* %7
ret float %85
}

LLVM IR (-O2):
define float @WarpFreq(float %0, float %1, float %2, float %3, float %4, float %5) {
%7 = fcmp oeq float %5, 1.000000e+00
br i1 %7, label %35, label %8
8: ; preds = %6
%9 = fdiv float 1.000000e+00, %5
%10 = fmul float %1, 2.000000e+00
%11 = fadd float %9, 1.000000e+00
%12 = fdiv float %10, %11
%13 = fmul float %0, 2.000000e+00
%14 = fdiv float %13, %11
%15 = fneg float %3
%16 = tail call float @llvm.fmuladd.f32(float %14, float %9, float %15)
%17 = fsub float %14, %3
%18 = fdiv float %16, %17
%19 = fcmp olt float %12, %2
br i1 %19, label %20, label %28
20: ; preds = %8
%21 = fneg float %12
%22 = tail call float @llvm.fmuladd.f32(float %21, float %9, float %4)
%23 = fsub float %4, %12
%24 = fdiv float %22, %23
%25 = fsub float %2, %12
%26 = fmul float %9, %12
%27 = tail call float @llvm.fmuladd.f32(float %24, float %25, float %26)
br label %35
28: ; preds = %8
%29 = fcmp ogt float %14, %2
br i1 %29, label %30, label %33
30: ; preds = %28
%31 = fsub float %2, %3
%32 = tail call float @llvm.fmuladd.f32(float %18, float %31, float %3)
br label %35
33: ; preds = %28
%34 = fmul float %9, %2
br label %35
35: ; preds = %20, %30, %33, %6
%36 = phi float [ %2, %6 ], [ %27, %20 ], [ %32, %30 ], [ %34, %33 ]
ret float %36
}