For a second-order phase transition the critical energy range of interest is larger than the energy range covered by a canonical Monte Carlo simulation at the critical temperature. Such an extended energy range can be covered by performing a Wang-Landau recursion for the spectral density followed by a multicanonical simulation with fixed weights. But in the conventional approach one loses the advantage due to cluster algorithms. A cluster version of the Wang-Landau recursion together with a subsequent multibondic simulation improves for 2D and 3D Ising models the efficiency of the conventional Wang-Landau or multicanonical approach by power laws in the lattice size. In our simulations real gains in CPU time reach 2 orders of magnitude.