Commit 1eff71bb authored by Tomáš Oberhuber's avatar Tomáš Oberhuber
Browse files

Refactoring operators.

parent 1f41d711

examples/advection/advection.h

+6 −6
Original line number Diff line number Diff line 
#include <TNL/tnlConfig.h>

#include <TNL/Solvers/Solver.h>

#include <TNL/Solvers/BuildConfigTags.h>

#include <TNL/operators/tnlDirichletBoundaryConditions.h>

#include <TNL/operators/tnlNeumannBoundaryConditions.h>

#include <TNL/Operators/DirichletBoundaryConditions.h>

#include <TNL/Operators/NeumannBoundaryConditions.h>

#include <TNL/Functions/Analytic/ConstantFunction.h>

#include "advectionProblem.h"

#include "LaxFridrichs.h"
@@ -86,12 +86,12 @@ class advectionSetter 
             typedef Functions::Analytic::ConstantFunction< Dimensions, Real > ConstantFunction;

             if( boundaryConditionsType == "dirichlet" )

             {

                typedef tnlDirichletBoundaryConditions< MeshType, ConstantFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

                typedef Operators::DirichletBoundaryConditions< MeshType, ConstantFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

                typedef advectionProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

                SolverStarter solverStarter;

                return solverStarter.template run< Problem >( parameters );

             }

             typedef tnlNeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

             typedef Operators::NeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

             typedef advectionProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

             SolverStarter solverStarter;

             return solverStarter.template run< Problem >( parameters );
@@ -99,12 +99,12 @@ class advectionSetter 
          typedef Functions::MeshFunction< MeshType > MeshFunction;

          if( boundaryConditionsType == "dirichlet" )

          {

             typedef tnlDirichletBoundaryConditions< MeshType, MeshFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

             typedef Operators::DirichletBoundaryConditions< MeshType, MeshFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

             typedef advectionProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

             SolverStarter solverStarter;

             return solverStarter.template run< Problem >( parameters );

          }

          typedef tnlNeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

          typedef Operators::NeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

          typedef advectionProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

          SolverStarter solverStarter;

          return solverStarter.template run< Problem >( parameters );

examples/heat-equation/tnl-heat-equation-eoc.h

+6 −6
Original line number Diff line number Diff line
@@ -15,11 +15,11 @@ 
#include <TNL/Solvers/FastBuildConfigTag.h>

#include <TNL/Solvers/BuildConfigTags.h>

#include <TNL/Functions/TestFunction.h>

#include <TNL/operators/diffusion/tnlLinearDiffusion.h>

#include <TNL/operators/diffusion/tnlExactLinearDiffusion.h>

#include <TNL/Operators/diffusion/LinearDiffusion.h>

#include <TNL/Operators/diffusion/ExactLinearDiffusion.h>

#include <TNL/Problems/HeatEquationEocRhs.h>

#include <TNL/Problems/HeatEquationEocProblem.h>

#include <TNL/operators/tnlDirichletBoundaryConditions.h>

#include <TNL/Operators/DirichletBoundaryConditions.h>



using namespace TNL;

using namespace TNL::Problems;
@@ -58,12 +58,12 @@ class heatEquationSetter 
   static bool run( const Config::ParameterContainer& parameters )

   {

      enum { Dimensions = MeshType::meshDimensions };

      typedef tnlLinearDiffusion< MeshType, Real, Index > ApproximateOperator;

      typedef tnlExactLinearDiffusion< Dimensions > ExactOperator;

      typedef Operators::LinearDiffusion< MeshType, Real, Index > ApproximateOperator;

      typedef Operators::ExactLinearDiffusion< Dimensions > ExactOperator;

      typedef Functions::TestFunction< MeshType::meshDimensions, Real, Device > TestFunction;

      typedef HeatEquationEocRhs< ExactOperator, TestFunction > RightHandSide;

      typedef Vectors::StaticVector < MeshType::meshDimensions, Real > Vertex;

      typedef tnlDirichletBoundaryConditions< MeshType, TestFunction, Dimensions, Real, Index > BoundaryConditions;

      typedef Operators::DirichletBoundaryConditions< MeshType, TestFunction, Dimensions, Real, Index > BoundaryConditions;

      typedef HeatEquationEocProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Solver;

      SolverStarter solverStarter;

      return solverStarter.template run< Solver >( parameters );

examples/heat-equation/tnl-heat-equation.h

+10 −10
Original line number Diff line number Diff line
@@ -14,9 +14,9 @@ 
#include <TNL/Solvers/Solver.h>

#include <TNL/Solvers/FastBuildConfigTag.h>

#include <TNL/Solvers/BuildConfigTags.h>

#include <TNL/operators/diffusion/tnlLinearDiffusion.h>

#include <TNL/operators/tnlDirichletBoundaryConditions.h>

#include <TNL/operators/tnlNeumannBoundaryConditions.h>

#include <TNL/Operators/diffusion/LinearDiffusion.h>

#include <TNL/Operators/DirichletBoundaryConditions.h>

#include <TNL/Operators/NeumannBoundaryConditions.h>

#include <TNL/Functions/Analytic/ConstantFunction.h>

#include <TNL/Functions/MeshFunction.h>

#include <TNL/Problems/HeatEquationProblem.h>
@@ -41,8 +41,8 @@ class heatEquationConfig 


         typedef tnlGrid< 1, double, Devices::Host, int > Mesh;

         typedef Functions::MeshFunction< Mesh > MeshFunction;

         tnlDirichletBoundaryConditions< Mesh, MeshFunction >::configSetup( config );

         tnlDirichletBoundaryConditions< Mesh, Functions::Analytic::ConstantFunction< 1 > >::configSetup( config );

         Operators::DirichletBoundaryConditions< Mesh, MeshFunction >::configSetup( config );

         Operators::DirichletBoundaryConditions< Mesh, Functions::Analytic::ConstantFunction< 1 > >::configSetup( config );

         config.addEntry< String >( "boundary-conditions-file", "File with the values of the boundary conditions.", "boundary.tnl" );

         config.addEntry< double >( "boundary-conditions-constant", "This sets a value in case of the constant boundary conditions." );

         config.addEntry< double >( "right-hand-side-constant", "This sets a constant value for the right-hand side.", 0.0 );
@@ -69,7 +69,7 @@ class heatEquationSetter 
   static bool run( const Config::ParameterContainer& parameters )

   {

      enum { Dimensions = MeshType::meshDimensions };

      typedef tnlLinearDiffusion< MeshType, Real, Index > ApproximateOperator;

      typedef Operators::LinearDiffusion< MeshType, Real, Index > ApproximateOperator;

      typedef Functions::Analytic::ConstantFunction< Dimensions, Real > RightHandSide;

      typedef Vectors::StaticVector < MeshType::meshDimensions, Real > Vertex;
@@ -79,12 +79,12 @@ class heatEquationSetter 
         typedef Functions::Analytic::ConstantFunction< Dimensions, Real > ConstantFunction;

         if( boundaryConditionsType == "dirichlet" )

         {

            typedef tnlDirichletBoundaryConditions< MeshType, ConstantFunction > BoundaryConditions;

            typedef Operators::DirichletBoundaryConditions< MeshType, ConstantFunction > BoundaryConditions;

            typedef HeatEquationProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

            SolverStarter solverStarter;

            return solverStarter.template run< Problem >( parameters );

         }

         typedef tnlNeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

         typedef Operators::NeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

         typedef HeatEquationProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

         SolverStarter solverStarter;

         return solverStarter.template run< Problem >( parameters );
@@ -92,12 +92,12 @@ class heatEquationSetter 
      typedef Functions::MeshFunction< MeshType > MeshFunction;

      if( boundaryConditionsType == "dirichlet" )

      {

         typedef tnlDirichletBoundaryConditions< MeshType, MeshFunction > BoundaryConditions;

         typedef Operators::DirichletBoundaryConditions< MeshType, MeshFunction > BoundaryConditions;

         typedef HeatEquationProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

         SolverStarter solverStarter;

         return solverStarter.template run< Problem >( parameters );

      }

      typedef tnlNeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

      typedef Operators::NeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

      typedef HeatEquationProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

      SolverStarter solverStarter;

      return solverStarter.template run< Problem >( parameters );

examples/inviscid-flow/1d/euler.h

+6 −6
Original line number Diff line number Diff line 
#include <TNL/tnlConfig.h>

#include <TNL/Solvers/Solver.h>

#include <TNL/Solvers/BuildConfigTags.h>

#include <TNL/operators/tnlDirichletBoundaryConditions.h>

#include <TNL/operators/tnlNeumannBoundaryConditions.h>

#include <TNL/Operators/DirichletBoundaryConditions.h>

#include <TNL/Operators/NeumannBoundaryConditions.h>

#include <TNL/Functions/Analytic/ConstantFunction.h>

#include "eulerProblem.h"

#include "LaxFridrichs.h"
@@ -82,12 +82,12 @@ class eulerSetter 
             typedef Functions::Analytic::ConstantFunction< Dimensions, Real > ConstantFunction;

             if( boundaryConditionsType == "dirichlet" )

             {

                typedef tnlDirichletBoundaryConditions< MeshType, ConstantFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

                typedef Operators::DirichletBoundaryConditions< MeshType, ConstantFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

                typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

                SolverStarter solverStarter;

                return solverStarter.template run< Problem >( parameters );

             }

             typedef tnlNeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

             typedef Operators::NeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

             typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

             SolverStarter solverStarter;

             return solverStarter.template run< Problem >( parameters );
@@ -95,12 +95,12 @@ class eulerSetter 
          typedef Functions::MeshFunction< MeshType > MeshFunction;

          if( boundaryConditionsType == "dirichlet" )

          {

             typedef tnlDirichletBoundaryConditions< MeshType, MeshFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

             typedef Operators::DirichletBoundaryConditions< MeshType, MeshFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

             typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

             SolverStarter solverStarter;

             return solverStarter.template run< Problem >( parameters );

          }

          typedef tnlNeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

          typedef Operators::NeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

          typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

          SolverStarter solverStarter;

          return solverStarter.template run< Problem >( parameters );

examples/inviscid-flow/2d/euler.h

+6 −6
Original line number Diff line number Diff line 
#include <TNL/tnlConfig.h>

#include <TNL/Solvers/Solver.h>

#include <TNL/Solvers/BuildConfigTags.h>

#include <TNL/operators/tnlDirichletBoundaryConditions.h>

#include <TNL/operators/tnlNeumannBoundaryConditions.h>

#include <TNL/Operators/DirichletBoundaryConditions.h>

#include <TNL/Operators/NeumannBoundaryConditions.h>

#include <TNL/Functions/Analytic/ConstantFunction.h>

#include "eulerProblem.h"

#include "LaxFridrichs.h"
@@ -83,12 +83,12 @@ class eulerSetter 
             typedef Functions::Analytic::ConstantFunction< Dimensions, Real > ConstantFunction;

             if( boundaryConditionsType == "dirichlet" )

             {

                typedef tnlDirichletBoundaryConditions< MeshType, ConstantFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

                typedef Operators::DirichletBoundaryConditions< MeshType, ConstantFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

                typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

                SolverStarter solverStarter;

                return solverStarter.template run< Problem >( parameters );

             }

             typedef tnlNeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

             typedef Operators::NeumannBoundaryConditions< MeshType, ConstantFunction, Real, Index > BoundaryConditions;

             typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

             SolverStarter solverStarter;

             return solverStarter.template run< Problem >( parameters );
@@ -96,12 +96,12 @@ class eulerSetter 
          typedef Functions::MeshFunction< MeshType > MeshFunction;

          if( boundaryConditionsType == "dirichlet" )

          {

             typedef tnlDirichletBoundaryConditions< MeshType, MeshFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

             typedef Operators::DirichletBoundaryConditions< MeshType, MeshFunction, MeshType::getMeshDimensions(), Real, Index > BoundaryConditions;

             typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

             SolverStarter solverStarter;

             return solverStarter.template run< Problem >( parameters );

          }

          typedef tnlNeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

          typedef Operators::NeumannBoundaryConditions< MeshType, MeshFunction, Real, Index > BoundaryConditions;

          typedef eulerProblem< MeshType, BoundaryConditions, RightHandSide, ApproximateOperator > Problem;

          SolverStarter solverStarter;

          return solverStarter.template run< Problem >( parameters );