Commit 2ff4b2bf authored by Tomáš Oberhuber's avatar Tomáš Oberhuber
Browse files

Fixing the Lax-Fridrichs scheme.

parent 89b43c94
Loading
Loading
Loading
Loading
+16 −3
Original line number Diff line number Diff line
@@ -33,6 +33,19 @@ tnlLaxFridrichs< tnlGrid< 2, Real, Device, Index, GridGeometry >,
{
}

template< typename Real,
          typename Device,
          typename Index,
          typename PressureGradient,
          template< int, typename, typename, typename > class GridGeometry >
tnlString tnlLaxFridrichs< tnlGrid< 2, Real, Device, Index, GridGeometry >,
                           PressureGradient > :: getTypeStatic()
{
   return tnlString( "tnlLaxFridrichs< " ) +
          tnlGrid< 2, Real, Device, Index, GridGeometry > :: getTypeStatic() + ", " +
          PressureGradient :: getTypeStatic() + " >";
}

template< typename Real,
          typename Device,
          typename Index,
+4 −2
Original line number Diff line number Diff line
@@ -23,8 +23,8 @@
#include <mesh/tnlIdenticalGridGeometry.h>
#include <schemes/gradient/tnlCentralFDMGradient.h>

template< typename MeshType,
          typename PressureGradient = tnlCentralFDMGradient >
template< typename Mesh,
          typename PressureGradient = tnlCentralFDMGradient< Mesh > >
class tnlLaxFridrichs
{
};
@@ -47,6 +47,8 @@ class tnlLaxFridrichs< tnlGrid< 2, Real, Device, Index, GridGeometry >, Pressure

   tnlLaxFridrichs();

   static tnlString getTypeStatic();

   void getExplicitRhs( const IndexType centralVolume,
                        RealType& rho_t,
                        RealType& rho_u1_t,