We have not considered that approach, rather (as most people do) the U parameters are viewed as fitting parameters. If you are looking at semiconductors, I would furthermore not recommend the Hubbard method, but rather look at DFT+1/2 or maybe meta-GGA (TB09). And in fact, these days even HSE06 is computationally reasonable, at least with LCAO basis sets.