bug fix

dune/fem-dg/models/modelwrapper.hh

@@ -417,7 +417,7 @@ namespace Fem
     }
 
   protected:
-    ProblemType impl_;
+    const ProblemType &impl_;
     ModifiedRangeType modified_;
   };
 

pydemo/euler/shock_tube.py

@@ -4,19 +4,15 @@ from dune.fem import parameter
 import dune.create as create
 from dune.models.elliptic.formfiles import loadModels
 from llf import NumFlux
-from ufl import *
 from dune.femdg import createFemDGSolver
+from ufl import *
 
-exec(open("euler.ufl").read())
+from euler import Model,Sod
 
 parameter.append({"fem.verboserank": -1})
 
 grid = structuredGrid([-1, 0], [1, 0.1], [200, 10])
 dimR      = 4
-UL = as_vector( Model.toCons([1,0,0,1]) )
-UR = as_vector( Model.toCons([0.125,0,0,0.1]) )
-x = SpatialCoordinate(space.cell())
-initial = conditional(x[0]<0,UL,UR)
 t = 0
 saveStep = 0.01
 saveTime = saveStep
@@ -32,7 +28,7 @@ def useGalerkinOp():
 
     n = FacetNormal(space.cell())
 
-    u_h   = space.interpolate(initial, name='u_h')
+    u_h   = space.interpolate(Sod(), name='u_h')
     u_h_n = u_h.copy(name="previous")
 
     fullModel = inner( Model.F_c(u), grad(v) ) * dx -\
@@ -61,9 +57,9 @@ def useODESolver():
     spaceName = "dgonb"
     polOrder = 0
     space = create.space(spaceName, grid, order=polOrder, dimrange=dimR)
-    u_h   = space.interpolate(initial, name='u_h')
+    u_h   = space.interpolate(Sod(), name='u_h')
     u_h_n = u_h.copy(name="previous")
-    operator = createFemDGSolver( Compressible2DEuler, space )
+    operator = createFemDGSolver( Model, space )
 
     start = time.time()
     grid.writeVTK('sod', pointdata=[u_h], number=count)

python/dune/femdg/_operators.py

@@ -114,7 +114,7 @@ def createFemDGSolver(Model, space):
                                'const '+advModelType + ' &advectionModel',
                                'const '+diffModelType + ' &diffusionModel'
                               ],
-                              ['return new ' + typeName + '(space);'],
+                              ['return new ' + typeName + '(space, advectionModel, diffusionModel);'],
                               ['"space"_a',
                                '"advectionModel"_a',
                                '"diffusionModel"_a',