[bugfix][Models] Use RangeFieldType for pressure and temperature.

dune/fem-dg/models/defaultmodel.hh

@@ -156,6 +156,7 @@ namespace Fem
     typedef FunctionSpaceType                            DFunctionSpaceType;
     typedef FunctionSpaceType                            RFunctionSpaceType;
     typedef typename Traits::DomainType                  DomainType;
+    typedef typename FunctionSpaceType::RangeFieldType   RangeFieldType;
     typedef typename Traits::RangeType                   RangeType;
     typedef typename Traits::GradientType                GradientType;
     typedef typename Traits::FluxRangeType               FluxRangeType;
@@ -222,7 +223,7 @@ namespace Fem
      *  \brief compute pressure and temperature for compressible flow
      */
     template <class State>
-    void pressureTemperature( const State& u, double& p, double& T ) const
+    void pressureTemperature( const State& u, RangeFieldType& p, RangeFieldType& T ) const
     {
       p = 0.0;
       T = 0.0;

dune/fem-dg/models/modelwrapper.hh

@@ -250,7 +250,7 @@ namespace Fem
     }
 
     template <class State> // double[dim+2]
-    void pressureTemperature( const State& u, double& p, double& T ) const
+    void pressureTemperature( const State& u, RangeFieldType& p, RangeFieldType& T ) const
     {
       if constexpr ( hasAdvection )
       {