[bugfix][ShallowWater] disable user defined flux, because it segfaults.

dune/fem-dg/operator/fluxes/advection/python.hh

@@ -83,9 +83,10 @@ namespace Fem
      */
     DGAdvectionFluxPythonUserDefine (const ModelImp& modelImp,
                                      const ParameterType& parameter )
-      : BaseType( *(new ModelType(modelImp)), parameter )
+      //: BaseType( *(new ModelType(modelImp)), parameter )
+      : BaseType( modelImp, parameter )
     {
-      modelPtr_.reset( &this->model_ );
+      //modelPtr_.reset( &this->model_ );
     }
 
     /**
@@ -111,7 +112,7 @@ namespace Fem
     }
 
   protected:
-    std::unique_ptr< const ModelType > modelPtr_;
+    //std::unique_ptr< const ModelType > modelPtr_;
   };
 
 } // end namespace Fem

pydemo/shallowWater/main.py

+from dune.fem import threading
+
+threading.use = 8
+
 from dune.fem import parameter
 from dune.femdg.testing import run
 

pydemo/shallowWater/shallowwater.py

@@ -40,11 +40,11 @@ def ShallowWater(topo,g):
         # boundary = {range(1,5): lambda t,x,u,n: Model.F_c(t,x,u)*n}
         boundary = {range(1,5): lambda t,x,u: u}
         topography   = topo
-        def NumericalF_c(model,clsName,includes):
-            print("NumericalFlux_c")
-            from dune.generator.importclass import load
-            from dune.generator import path
-            return load(clsName,[path(__file__)+"llf.hh"]+includes,model)
+        #def NumericalF_c(model,clsName,includes):
+        #    print("NumericalFlux_c")
+        #    from dune.generator.importclass import load
+        #    from dune.generator import path
+        #    return load(clsName,[path(__file__)+"llfwb.hh"]+includes,model)
 
     return Model