The authors consider the automatic generation of parallel templates as basis for the implementation of n/sup 2/-problem solvers on parallel systems. The automatic template generator relies on the embedding of a virtual systolic loop on the given machine's network and optimizes hyper-systolic data communication by means of simulated annealing techniques. They have generated a parallel template to efficiently program an n-body code on the parallel SIMD system Quadrics.