  We adapt the VMR (valleys-to-mountains reflections) algorithm, originally
devised by us for simulations of SOS models, to the BCSOS model. It is the
first time that a cluster algorithm is used for a model with constraints. The
performance of this new algorithm is studied in detail in both phases of the
model, including a finite size scaling analysis of the autocorrelations.
