Cluster mean field plus density matrix renormalization theory for the Bose Hubbard models. (1st July 2022)