I am trying to impose homogeneous Neumann boundary condition for discontinuous Galerkin method. But the int1d(…) command for the Neumann boundary condition is not working for discontinuous Galerkin case. It’s not showing any error message but also not imposing this boundary condition in the output solution. When I use this command for continuous Galerkin method, it works well.
So, how can I impose Neumann boundary condition for solving a problem using discontinuous Galerkin method ?