julia> x = rand(1000);
julia> y1 = @btime FastTransforms.cheb2ultra($x, 2);
  24.224 ms (2 allocations: 7.98 KiB)
julia> y2 = @btime ApproxFunBase.mul_coefficients(Conversion(Chebyshev(), Ultraspherical(2)), $x);
  45.568 μs (13 allocations: 63.27 KiB)
julia> y1 ≈ y2
true