Dominic Kempf · Dominic Kempf
--- a/tutorial07/src/numericalflux.hh 0 → 100644
+++ b/tutorial07/src/numericalflux.hh 0 → 100644
+  using DF = RF;
+
+  VariableFluxVectorSplitting (const MODEL& model_)
+    : model(model_)
+  {
+  }
+
+  template<typename E, typename X>
+  void numericalFlux(const E& inside, const X& x_inside,
+                     const E& outside, const X& x_outside,
+                     const Dune::FieldVector<DF,dim> n_F,
+                     const Dune::FieldVector<RF,m>& u_s,
+                     const Dune::FieldVector<RF,m>& u_n,Dune::FieldVector<RF,m>& f) const
+  {
+    // check for discontinuity
+    if (model.problem.c(inside,x_inside)==model.problem.c(outside,x_outside))
--- a/tutorial07/src/shallowwaterproblem.hh 0 → 100644
+++ b/tutorial07/src/shallowwaterproblem.hh 0 → 100644
+#ifndef ACOUSTICS_RIEMANNPROBLEM
--- a/tutorial07/src/driver.hh 0 → 100644
+++ b/tutorial07/src/driver.hh 0 → 100644
+  using GFS = Dune::PDELab::PowerGridFunctionSpace<GFSDG,m,VBE,OrderingTag>;
+  GFS gfs(gfsdg);
+
+  typedef typename GFS::template ConstraintsContainer<RF>::Type C;
+  C cg;
+  gfs.update(); // initializing the gfs
+  std::cout << "degrees of freedom: " << gfs.globalSize() << std::endl;
+
+  // Make instationary grid operator
+  using LOP = Dune::PDELab::DGLinearHyperbolicSpatialOperator<NUMFLUX,FEMDG>;
+  LOP lop(numflux);
+  using TLOP = Dune::PDELab::DGLinearHyperbolicTemporalOperator<NUMFLUX,FEMDG>;
+  TLOP tlop(numflux);
+
+  using MBE = Dune::PDELab::ISTL::BCRSMatrixBackend<>;
+  MBE mbe(5); // Maximal number of nonzeroes per row
--- a/tutorial07/src/linearacoustics.hh 0 → 100644
+++ b/tutorial07/src/linearacoustics.hh 0 → 100644
+    \param c speed of sound
+    \param n unit outer normal vector
+    \param RT matrix to be filled
+*/
+
+template<int dim, typename PROBLEM>
+class Model ;
+
+template<typename PROBLEM>
+class Model<2,PROBLEM>
+{
+public:
+  static constexpr int dim = 2;   // space dimension
+  static constexpr int m = dim+1; // system size
+  static constexpr int mplus = 1; // number of positive eigenvectors
+  static constexpr int mminus = 1; // number of negative eigenvectors