Enable Limtier again.

dune/fem-dg/models/modelwrapper.hh

@@ -374,10 +374,6 @@ namespace Fem
                      RangeType& jump) const
     {
       jump = AdditionalType :: jump( it, x, uLeft, uRight );
-      // take pressure as shock detection values
-      //const RangeFieldType pl = pressure( uLeft );
-      //const RangeFieldType pr = pressure( uRight );
-      //jump  = (pl-pr)/(0.5*(pl+pr));
     }
 
     // calculate jump between left and right value

dune/fem-dg/solver/dg.hh

@@ -30,7 +30,7 @@ namespace Fem
              class Additional,
              Solver::Enum solverId             = Solver::Enum::fem,
              Formulation::Enum formId          = Formulation::Enum::primal,
-             AdvectionLimiter::Enum limiterId  = AdvectionLimiter::Enum::unlimited,
+             AdvectionLimiter::Enum limiterId  = AdvectionLimiter::Enum::limited,
              AdvectionFlux::Enum advFluxId     = AdvectionFlux::Enum::llf,
              DiffusionFlux::Enum diffFluxId    = DiffusionFlux::Enum::primal
              >

pydemo/euler/shock_tube.py

@@ -11,7 +11,7 @@ from euler import Model,Sod
 
 parameter.append({"fem.verboserank": -1})
 
-grid = structuredGrid([-1, 0], [1, 0.1], [200, 10])
+grid = structuredGrid([-1, 0], [1, 0.1], [20, 5])
 dimR      = 4
 t = 0
 dt = 1e-5
@@ -58,7 +58,7 @@ def useGalerkinOp():
 def useODESolver():
     global count, t, dt, saveTime
     spaceName = "dgonb"
-    polOrder = 0
+    polOrder = 2
     space = create.space(spaceName, grid, order=polOrder, dimrange=dimR)
     u_h   = space.interpolate(Sod(), name='u_h')
     u_h_n = u_h.copy(name="previous")
@@ -71,6 +71,7 @@ def useODESolver():
         u_h_n.assign(u_h)
         operator(u_h_n, u_h)
         dt = operator.deltaT()
+        print('dt = ', dt )
         t += dt
         if t > saveTime:
             count += 1