#ifndef HeatEquationBenchmarkRHS_H_
#define HeatEquationBenchmarkRHS_H_
#include <TNL/Functions/Domain.h>
template< typename Mesh, typename Real >class HeatEquationBenchmarkRhs
  : public Functions::Domain< Mesh::getMeshDimension(), Functions::MeshDomain > 
 {
   public:
      using MeshType = Mesh;
      using RealType = Real;

      bool setup( const Config::ParameterContainer& parameters,
                  const String& prefix = "" )
      {
         return true;
      }

      template< typename MeshEntity >
      __cuda_callable__
      Real operator()( const MeshEntity& entity,
                       const Real& time = 0.0 ) const
      {
         //typedef typename MeshEntity::MeshType::PointType PointType;
         //PointType v = entity.getCenter();
         return 0.0;
      }
};

#endif /* HeatEquationBenchmarkRHS_H_ */
