The primary unknowns are the wave_scalar s and its first time
derivative wave_fscalar
(
as TOCHNOG only solves first order in time equations, the first time
derivative of s also becomes primary unknown in order to turn this second
order in time equation into a set of first order in time equations).
Further notation: x space coordinate, t time and c speed of sound.