float ct_poly_blep(float t, const float dt) {
  if (t < dt) {
    t /= dt;
    return t + t - t * t - 1.f;
  } else if (t > 1.f - dt) {
    t = (t - 1.0f) / dt;
    return t * t + t + t + 1.;
  } else {
    return 0.0f;
  }
}


x86 assembly (clang, -O0):
<ct_poly_blep>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x8(%rbp)
movss  %xmm1,-0xc(%rbp)
movss  -0x8(%rbp),%xmm1
movss  -0xc(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1173 <ct_poly_blep+0x73>
movss  -0xc(%rbp),%xmm1
movss  -0x8(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x8(%rbp)
movss  -0x8(%rbp),%xmm1
addss  -0x8(%rbp),%xmm1
movss  -0x8(%rbp),%xmm0
movss  -0x8(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  0xe9b(%rip),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11eb <ct_poly_blep+0xeb>
movss  -0x8(%rbp),%xmm0
movss  0xe80(%rip),%xmm1
subss  -0xc(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    11e3 <ct_poly_blep+0xe3>
movss  -0x8(%rbp),%xmm0
movss  0xe65(%rip),%xmm1
subss  %xmm1,%xmm0
divss  -0xc(%rbp),%xmm0
movss  %xmm0,-0x8(%rbp)
movss  -0x8(%rbp),%xmm0
movss  -0x8(%rbp),%xmm2
movss  -0x8(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
addss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe37(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11eb <ct_poly_blep+0xeb>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ct_poly_blep>:
ucomiss %xmm0,%xmm1
jbe    1124 <ct_poly_blep+0x24>
divss  %xmm1,%xmm0
movaps %xmm0,%xmm2
addss  %xmm0,%xmm2
mulss  %xmm0,%xmm0
subss  %xmm0,%xmm2
addss  0xee4(%rip),%xmm2
movaps %xmm2,%xmm0
ret
movss  0xed4(%rip),%xmm3
subss  %xmm1,%xmm3
xorps  %xmm2,%xmm2
ucomiss %xmm3,%xmm0
jbe    115b <ct_poly_blep+0x5b>
addss  0xec4(%rip),%xmm0
divss  %xmm1,%xmm0
movaps %xmm0,%xmm2
mulss  %xmm0,%xmm2
addss  %xmm0,%xmm2
addss  %xmm0,%xmm2
addss  0xea5(%rip),%xmm2
movaps %xmm2,%xmm0
ret

LLVM IR (-O0):
define float @ct_poly_blep(float %0, float %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca float
store float %0, float* %4
store float %1, float* %5
%6 = load float, float* %4
%7 = load float, float* %5
%8 = fcmp olt float %6, %7
br i1 %8, label %9, label %21
9: ; preds = %2
%10 = load float, float* %5
%11 = load float, float* %4
%12 = fdiv float %11, %10
store float %12, float* %4
%13 = load float, float* %4
%14 = load float, float* %4
%15 = fadd float %13, %14
%16 = load float, float* %4
%17 = load float, float* %4
%18 = fneg float %16
%19 = call float @llvm.fmuladd.f32(float %18, float %17, float %15)
%20 = fsub float %19, 1.000000e+00
store float %20, float* %3
br label %41
21: ; preds = %2
%22 = load float, float* %4
%23 = load float, float* %5
%24 = fsub float 1.000000e+00, %23
%25 = fcmp ogt float %22, %24
br i1 %25, label %26, label %40
26: ; preds = %21
%27 = load float, float* %4
%28 = fsub float %27, 1.000000e+00
%29 = load float, float* %5
%30 = fdiv float %28, %29
store float %30, float* %4
%31 = load float, float* %4
%32 = load float, float* %4
%33 = load float, float* %4
%34 = call float @llvm.fmuladd.f32(float %31, float %32, float %33)
%35 = load float, float* %4
%36 = fadd float %34, %35
%37 = fpext float %36 to double
%38 = fadd double %37, 1.000000e+00
%39 = fptrunc double %38 to float
store float %39, float* %3
br label %41
40: ; preds = %21
store float 0.000000e+00, float* %3
br label %41
41: ; preds = %40, %26, %9
%42 = load float, float* %3
ret float %42
}

LLVM IR (-O2):
define float @ct_poly_blep(float %0, float %1) {
%3 = fcmp olt float %0, %1
br i1 %3, label %4, label %10
4: ; preds = %2
%5 = fdiv float %0, %1
%6 = fadd float %5, %5
%7 = fneg float %5
%8 = tail call float @llvm.fmuladd.f32(float %7, float %5, float %6)
%9 = fadd float %8, -1.000000e+00
br label %19
10: ; preds = %2
%11 = fsub float 1.000000e+00, %1
%12 = fcmp olt float %11, %0
br i1 %12, label %13, label %19
13: ; preds = %10
%14 = fadd float %0, -1.000000e+00
%15 = fdiv float %14, %1
%16 = tail call float @llvm.fmuladd.f32(float %15, float %15, float %15)
%17 = fadd float %15, %16
%18 = fadd float %17, 1.000000e+00
br label %19
19: ; preds = %10, %13, %4
%20 = phi float [ %9, %4 ], [ %18, %13 ], [ 0.000000e+00, %10 ]
ret float %20
}