tfb_discrete_cosine_transform function

ComputesY = g(X) = DCT(X), where DCT type is indicated by the type arg