WARNING: this code has not yet been cleaned up it. It repeats the computational experiment from Section 4.1 of the paper on a slightly smaller scale.