double blend(double a, double b, double k) {
  if (a < b) {
    if (k < a)
      return 1.f;
    if (b < k)
      return 0.f;
    double t = 1.f - (k - a) / (b - a);
    return 3 * t * t - 2 * t * t * t;
  } else {
    if (k > a)
      return 1.f;
    if (b > k)
      return 0.f;
    double t = 1.f - (a - k) / (a - b);
    return 3 * t * t - 2 * t * t * t;
  }
}


x86 assembly (clang, -O0):
<blend>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11e9 <blend+0xe9>
movsd  -0x20(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    114d <blend+0x4d>
movsd  0xecd(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    129e <blend+0x19e>
movsd  -0x18(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    116e <blend+0x6e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    129e <blend+0x19e>
movsd  -0x20(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
divsd  %xmm0,%xmm1
movsd  0xe82(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe69(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  0xe4f(%rip),%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    129e <blend+0x19e>
movsd  -0x20(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    120b <blend+0x10b>
movsd  0xe0f(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    129e <blend+0x19e>
movsd  -0x18(%rbp),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1228 <blend+0x128>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    129e <blend+0x19e>
movsd  -0x10(%rbp),%xmm1
subsd  -0x20(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
divsd  %xmm0,%xmm1
movsd  0xdc8(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xdaf(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  0xd95(%rip),%xmm1
mulsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<blend>:
ucomisd %xmm0,%xmm1
jbe    1130 <blend+0x30>
ucomisd %xmm2,%xmm0
ja     1136 <blend+0x36>
xorpd  %xmm3,%xmm3
ucomisd %xmm1,%xmm2
ja     114a <blend+0x4a>
subsd  %xmm0,%xmm2
subsd  %xmm0,%xmm1
divsd  %xmm1,%xmm2
movsd  0xed6(%rip),%xmm1
subsd  %xmm2,%xmm1
jmp    116b <blend+0x6b>
ucomisd %xmm0,%xmm2
jbe    1140 <blend+0x40>
movsd  0xec2(%rip),%xmm3
jmp    114a <blend+0x4a>
xorpd  %xmm3,%xmm3
ucomisd %xmm2,%xmm1
jbe    114f <blend+0x4f>
movapd %xmm3,%xmm0
ret
movapd %xmm0,%xmm3
subsd  %xmm2,%xmm3
subsd  %xmm1,%xmm0
divsd  %xmm0,%xmm3
movsd  0xe99(%rip),%xmm1
subsd  %xmm3,%xmm1
movsd  0xe95(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0xe91(%rip),%xmm2
mulsd  %xmm1,%xmm2
mulsd  %xmm1,%xmm2
mulsd  %xmm1,%xmm2
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
ret

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

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