Commit 9907cc47 authored by Tomáš Oberhuber's avatar Tomáš Oberhuber
Browse files

Adding vector field to tnl-view. 

Adding Gnuplot writer for vector field.
parent 3d213169

examples/inviscid-flow/RiemannProblemInitialCondition.h

+4 −3
Original line number Diff line number Diff line
@@ -35,9 +35,9 @@ class RiemannProblemInitialCondition 
      

      RiemannProblemInitialCondition()

         : discontinuityPlacement( 0.5 ),

           leftDensity( 1.0 ), rightDensity( 0.0 ),

           leftVelocity( 1.0 ), rightVelocity( 0.0 ),

           leftPressure( 1.0e5 ), rightPressure( 0.0 ),

           leftDensity( 1.0 ), rightDensity( 1.0 ),

           leftVelocity( -2.0 ), rightVelocity( 2.0 ),

           leftPressure( 0.4 ), rightPressure( 0.4 ),

           gamma( 1.67 ){}



      static void configSetup( Config::ConfigDescription& config,
@@ -169,6 +169,7 @@ class RiemannProblemInitialCondition 
         initialCondition->getOperator().setPositiveValue( leftDensity );

         initialCondition->getOperator().setNegativeValue( rightDensity );

         evaluator.evaluate( conservativeVariables.getDensity(), initialCondition );

         conservativeVariables.getDensity()->write( "density.gplt", "gnuplot" );

         

         /****

          * Momentum

examples/inviscid-flow/eulerProblem_impl.h

+3 −1
Original line number Diff line number Diff line
@@ -135,9 +135,11 @@ setInitialCondition( const Config::ParameterContainer& parameters, 
      if( ! initialCondition.setup( parameters ) )

         return false;

      initialCondition.setInitialCondition( conservativeVariables );

   }

      return true;

   }

   std::cerr << "Uknown initial condition " << initialConditionType << std::endl;

   return false;

}



template< typename Mesh,

          typename BoundaryCondition,

examples/inviscid-flow/run-euler

+18 −5
Original line number Diff line number Diff line
@@ -10,7 +10,20 @@ tnl-grid-setup --dimensions 2 \ 


tnl-init --test-function sin-wave \

         --output-file init.tnl

./euler --time-discretisation explicit \

tnl-euler-2d --discontinuity-placement-0 0.3 \

             --discontinuity-placement-1 0.3 \

             --discontinuity-placement-2 0.3 \

             --left-density 1.0 \

             --right-density 1.0 \

             --left-velocity-0 -2.0 \

             --left-velocity-1 -2.0 \

             --left-velocity-2 -2.0 \

             --right-velocity-0 2.0 \

             --right-velocity-1 2.0 \

             --right-velocity-2 2.0 \

             --left-pressure 0.4 \

             --right-pressure 0.4 \

             --time-discretisation explicit \

             --boundary-conditions-constant 0 \

             --discrete-solver merson \

             --snapshot-period 0.01 \

src/TNL/Functions/CMakeLists.txt

+5 −1
Original line number Diff line number Diff line
@@ -14,7 +14,11 @@ SET( headers Domain.h 
             MeshFunctionVTKWriter_impl.h             

             OperatorFunction.h

             TestFunction.h             

             TestFunction_impl.h )

             TestFunction_impl.h

             VectorField.h

             VectorFieldGnuplotWriter.h

             VectorFieldGnuplotWriter_impl.h

 )



SET( CURRENT_DIR ${CMAKE_SOURCE_DIR}/src/TNL/Functions )

set( common_SOURCES

src/TNL/Functions/MeshFunction.h

+1 −1
Original line number Diff line number Diff line
@@ -92,7 +92,7 @@ class MeshFunction : 
      

      const MeshPointer& getMeshPointer() const;

      

      __cuda_callable__ IndexType getDofs( const MeshPointer& meshPointer ) const;

      __cuda_callable__ static IndexType getDofs( const MeshPointer& meshPointer );

      

      __cuda_callable__ const VectorType& getData() const;