WIP: Fix/ldg implementation
LDG was not properly implemented, there was only a single numerical flux instead of computing left and right fluxes. Also, the boundary condition is now implemented as u_ext = - u_int + 2*u_Dirichelet instead of u_ext = u_Dirichlet.
-
1D test Neumann -
LDG convergence study -
update LDGNS implementation
Edited by Giacomo Castiglioni