Thanks for your reply. I just took a rough look at this paper. It seems the algorithm is not suitable for my case. My design is a high speed modem(above 10M symbol rate). And the system frequency is only around 50M.
So my modem can't use so many multiplication. As the generated new frequency may alias back.