compute_rdf extension

Present version has an upper bound on the number of atom pairs; need an extension.