  A fast, parallel algorithm for distant-dependent calculation and simulation
of crystal properties is presented along with speedup results and methods of
application. An illustrative example is used to compute the Lennard-Jones
lattice constants up to 32 significant figures for $4\leq p\leq30$ in the
simple cubic, face-centered cubic, body-centered cubic, hexagonal-close-pack,
and diamond lattices. In most cases, the known precision of these constants is
more than doubled, and in some cases, corrected from previously published
figures. The tools and strategies to make this computation possible are
detailed along with application to other potentials, including those that model
defects.
