  Within the next decade the Laser Interferometer Space Antenna (LISA) is due
to be launched, providing the opportunity to extract physics from stellar
objects and systems, such as \textit{Extreme Mass Ratio Inspirals}, (EMRIs)
otherwise undetectable to ground based interferometers and Pulsar Timing Arrays
(PTA). Unlike previous sources detected by the currently available
observational methods, these sources can \textit{only} be simulated using an
accurate computation of the gravitational self-force. Whereas the field has
seen outstanding progress in the frequency domain, metric reconstruction and
self-force calculations are still an open challenge in the time domain. Such
computations would not only further corroborate frequency domain calculations
and models, but also allow for full self-consistent evolution of the orbit
under the effect of the self-force. Given we have \textit{a priori} information
about the local structure of the discontinuity at the particle, we will show
how to construct discontinuous spatial and temporal discretisations by
operating on discontinuous Lagrange and Hermite interpolation formulae and
hence recover higher order accuracy. In this work we demonstrate how this
technique in conjunction with well-suited gauge choice (hyperboloidal slicing)
and numerical (discontinuous collocation with time symmetric) methods can
provide a relatively simple method of lines numerical algorithm to the problem.
This is the first of a series of papers studying the behaviour of a
point-particle prescribing circular geodesic motion in Schwarzschild in the
\textit{time domain}. In this work we describe the numerical machinery
necessary for these computations and show not only our work is capable of
highly accurate flux radiation measurements but it also shows suitability for
evaluation of the necessary field and it's derivatives at the particle limit.
