Export method deltaT

dune/fem-dg/solver/dg.hh

@@ -134,6 +134,8 @@ namespace Fem
       tp_.provideTimeStepEstimate( dt );
     }
 
+    double deltaT() const { return tp_.deltaT(); }
+
   protected:
     const DiscreteFunctionSpaceType&      space_;
 

pydemo/euler/shock_tube.py

@@ -68,6 +68,7 @@ def useODESolver():
     while t < 0.4:
         u_h_n.assign(u_h)
         operator(u_h_n, u_h)
+        dt = operator.deltaT()
         t += dt
         if t > saveTime:
             count += 1

python/dune/femdg/_operators.py

@@ -125,4 +125,7 @@ def createFemDGSolver(Model, space):
     # add method activated to inspect limited cells.
     setTimeStepSize = Method('setTimeStepSize', '&'+typeName+'::setTimeStepSize')
 
-    return load(includes, typeName, constructor, setTimeStepSize).Operator( space, advModel, diffModel )
+    # add method to obtain time step size
+    deltaT = Method('deltaT', '&'+typeName+'::deltaT')
+
+    return load(includes, typeName, constructor, setTimeStepSize, deltaT ).Operator( space, advModel, diffModel )