Volume integral not computed for large number of elements

[Balavarun5]

Commit : 401a43b

Issue

The volume integral is carried out using the volume_integral_flux function in wave_equation.py.
When the integration is carried out using either the Gauss-Lobatto or the Gaussian quadrature, the code after the else condition is executed.

This gives errors while dealing with a large number of elements.