Remark 2 In problem (REF ) and accordingly (REF ), for simplicity, we only consider constant \(\omega \)  and assume that the solution \(u(\cdot \,,t)\)  satisfies homogeneous boundary condition. One requires slight modifications of the discrete bilinear and linear functions for the cases of heterogeneous propagation speed and non-homogeneous boundary conditions [1]}, [2]}.

