In this work, we investigate numerically a one-dimensional wave equation in generalized form. The system considers the presence of constant damping and functional anomalous diffusion of the Riesz type. Reaction terms are also considered, in such way that the mathematical model can be presented in variational form when damping is not present. As opposed to previous efforts available in the literature, the reaction terms are not only functions of the solution. Instead, we consider the presence of smooth functions that depend on fractional derivatives of the solution function. Using a finite-difference approach, we propose a numerical scheme to approximate the solutions of the fractional wave equation. Along with this integrator, we propose discrete forms of the local and the total energy operators. In a first stage, we show rigorously that the energy properties of the continuous system are mimicked by our discrete methodology. In particular, we prove that the discrete system is dissipative (respectively, conservative) when damping is present (respectively, absent), in agreement with the continuous model. The theoretical numerical analysis of this system is more complicated in light of the presence of the functional form of the anomalous diffusion. To solve this problem, some novel technical lemmas are proved and used to establish the stability and the quadratic convergence of the scheme. Finally, we provide some computer simulations to show the capability of the scheme to conserve/dissipate the energy. Various fractional problems with functional forms of the anomalous diffusion of the solution are considered to that effect.