QuazarTech / DG_Maxwell

Solution of Maxwell equations using Discontinuous Galerkin method.
GNU General Public License v3.0
9 stars 5 forks source link

An error floor for 1D advection #54

Open Balavarun5 opened 7 years ago

Balavarun5 commented 7 years ago

Commit : 5a0b7437f411cce7eb85156f44c13a788c2a54e1

Issue

The L1 norm of an advected 1D wave has an error floor of 1e-8 as shown in the image below. advection_n_quad_lgl_plus_1

After a certain N_LGL, the error doesn't decay as N_LGL ** (- N_LGL)

This error could be a result of the Issue discussed in #27 where the volume integral term has a non trivial error. The plot of L1 norm of volume integral term vs N_LGL at the time of raising this issue is shown below vol_int_quad_lgl_plus_1