Debugging

Set S=[32; 32; 32] and R=diag([8 8 8]) in ``setup.m'' and run ``setup.m''. You should find a basis set compression ratio of about 18.3, a little larger than the theoretical value. (We still have relatively few points so that the number of points actually falling within a sphere isn't exactly the volume of the sphere.)