[bugfix][ModelCaller] Need at least size 1 for jacobian vectors.

dune/fem-dg/pass/modelcaller.hh

@@ -84,8 +84,10 @@ namespace Dune
         time_( 0. ),
         discreteFunctions_( FilterType::apply( argument ) ),
         localFunctionsInside_( Dune::DereferenceTuple< DiscreteFunctionPointerTupleType >::apply( discreteFunctions_ ) ),
-        localFunctionsOutside_( Dune::DereferenceTuple< DiscreteFunctionPointerTupleType >::apply( discreteFunctions_ ) )
-      {}
+        localFunctionsOutside_( Dune::DereferenceTuple< DiscreteFunctionPointerTupleType >::apply( discreteFunctions_ ) ),
+        jacobians_( 1 ) // need at least size 1 for temporary variables
+      {
+      }
 
       // return true, if discrete model has flux
       bool hasFlux () const { return discreteModel().hasFlux(); }
@@ -309,12 +311,15 @@ namespace Dune
 
       CDGDiscreteModelCaller ( ArgumentType &argument, DiscreteModelType &discreteModel )
       : BaseType( argument, discreteModel )
+        , jacobiansInside_( 1 )
+        , jacobiansOutside_( 1 )
 #ifndef NDEBUG
         , quadInnerId_( 0 )
         , quadOuterId_( 0 )
         , quadId_( 0 )
 #endif
-      {}
+      {
+      }
 
       void setEntity ( const EntityType &entity, const VolumeQuadratureType &quadrature )
       {

dune/fem-dg/pass/threadpass.hh

@@ -679,6 +679,7 @@ namespace Fem
         // create NB checker, noSkip of intersections
         NBChecker< skipNone > nbChecker( iterators_, thread );
 
+        assert( arg_ );
         const Iterator endit = iterators_.end();
         for (Iterator it = iterators_.begin(); it != endit; ++it)
         {
@@ -687,7 +688,6 @@ namespace Fem
           myPass.applyLocal( entity, nbChecker );
         }
 
-        assert( arg_ );
         // dest can also be null pointer
         // when the operator is evaluated only
         // for evaluation of the estimators