Also pass jacobian values to boundaryFlux.

dune/fem-dg/operator/dg/discretemodelcommon.hh

@@ -286,7 +286,7 @@ namespace Dune {
         else
         {
           // returns advection wave speed
-          return model_.boundaryFlux( left, left.values()[uVar], gLeft );
+          return model_.boundaryFlux( left, left.values()[uVar], left.jacobians()[uVar], gLeft );
         }
       }
       else