run_stage function

Run a stage of the PMwG sampler