  An algorithm for Monte Carlo simulations is proposed in which the parameter
controlling the strength of the transition becomes a dynamical variable and in
which efficient transitions are achieved by cluster steps. It allows to avoid
the strongly suppressed tunneling between the phases by travelling easily via
the second order region. Numerical results for the Potts model are presented
which demonstrate the advantages of the method.
