#ifndef SCALAR_NODAL_SCHEME_HPP
#define SCALAR_NODAL_SCHEME_HPP

#include <algebra/CRSMatrix.hpp>
#include <algebra/CRSMatrixDescriptor.hpp>
#include <algebra/LeastSquareSolver.hpp>
#include <algebra/LinearSolver.hpp>
#include <algebra/TinyVector.hpp>
#include <algebra/Vector.hpp>
#include <language/utils/InterpolateItemValue.hpp>
#include <mesh/Connectivity.hpp>
#include <mesh/DualConnectivityManager.hpp>
#include <mesh/DualMeshManager.hpp>
#include <mesh/Mesh.hpp>
#include <mesh/MeshData.hpp>
#include <mesh/MeshDataManager.hpp>
#include <mesh/MeshFaceBoundary.hpp>
#include <mesh/MeshNodeBoundary.hpp>
#include <mesh/PrimalToDiamondDualConnectivityDataMapper.hpp>
#include <mesh/SubItemValuePerItem.hpp>
#include <scheme/DirichletBoundaryConditionDescriptor.hpp>
#include <scheme/FourierBoundaryConditionDescriptor.hpp>
#include <scheme/IDiscreteFunction.hpp>
#include <scheme/NeumannBoundaryConditionDescriptor.hpp>
#include <scheme/SymmetryBoundaryConditionDescriptor.hpp>

class ScalarNodalSchemeHandler
{
 private:
  class IScalarNodalScheme;

  template <size_t Dimension>
  class ScalarNodalScheme;

  template <size_t Dimension>
  class InterpolationWeightsManager;

 public:
  std::unique_ptr<IScalarNodalScheme> m_scheme;

  std::shared_ptr<const IDiscreteFunction> solution() const;

  ScalarNodalSchemeHandler(const std::shared_ptr<const IDiscreteFunction>& cell_k,
                           const std::shared_ptr<const IDiscreteFunction>& f,
                           const std::vector<std::shared_ptr<const IBoundaryConditionDescriptor>>& bc_descriptor_list,
                           const std::shared_ptr<const IDiscreteFunction>& P,
                           const double& dt);

  ~ScalarNodalSchemeHandler();
};

#endif   // SCALAR_NODAL_SCHEME_HPP