A quick test with C++'s <random> library agrees. For n = 100,000,000 points, the rejection method beat out the improved spherical by 7.8%. Using the normal distribution was even slower.
XYZ: 16.418057 seconds
Gauss: 23.415560 seconds
Spherical: 17.705553 seconds