added case for DG displacemnet plot

src/elasticity/assemble/IElasticity.cpp

@@ -170,6 +170,9 @@ void IElasticity::PrintPressureIteration(const Vector &u) const {
 }
 
 void IElasticity::PlotIteration(const Vectors &mechValues, int step) const {
+  if (mechValues[0].GetDisc().DiscName().find("DG") != std::string::npos) {
+    THROW("TODO: Override in DGElasticity as in PlotDisplacement!")
+  }
   step = step < 0 ? Step() : step;
 
   if (step == 0) PlotBoundary(mechValues[0], mechValues[0]);

src/elasticity/assemble/IElasticity.hpp

@@ -139,7 +139,7 @@ public:
 
   virtual void GetInvariant(const Vector &u, Vector &iota4, Vector &iota) const { return GetInvariant(u,iota4); }
 
-  void PlotDisplacement(const Vector &u, int step, const string &varname) const;
+  virtual void PlotDisplacement(const Vector &u, int step, const string &varname) const;
 
   /**
    *   Creates a new vector with the discretization of target (coarse vector) and projects fine