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.
Commit : 401a43b0dfc43b29bebfdfa96e20a189733d8316
Issue
The volume integral is carried out using the
volume_integral_fluxfunction inwave_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.