This work was funded by NSF grant #CMS-0403746. We thank ORNL, NERSC and NCSA for providing computer time, and Dave McWilliams of NCSA, in particular, for help in getting the large jobs to work. Computer time at ORNL was provided by CNMS/NTI grant #CNMS2004–028. Computer time at NCSA was provided under proposal #ASC050004. We also thank Xin Yuan at Florida State University for permitting use of his Linux cluster, where our codes were first tested. A.S. gratefully acknowledges his debt to Sri S. S. Baba, but for whose help and inspiration, this work would neither have started nor been accomplished.