Numerical Implementation

The only involved part of this thorn arise in working out where the sources are located (if at all) on each local grid for a multiprocessor run. The source terms are not numerically evolved, but are calculated exactly, based on the physical time and their orbital velocity.

A routine is scheduled to run after the homogeneous equation for the scalar field has been evolved, and simply updates the value of the scalar field by adding on the source contribution.