mesh.hpp

Go to the documentation of this file.

// Copyright (c) 2024 Andrea Enrico Franzoni (andreaenrico.franzoni@gmail.com)
//
// This file is part of PPCKO
//
// Permission is hereby granted, free of charge, to any person obtaining a copy
// of PPCKO and associated documentation files (the PPCKO software), to deal
// PPCKO without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of PPCKO, and to permit persons to whom PPCKO is
// furnished to do so, subject to the following conditions:
//
// PPCKO IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
// OUT OF OR IN CONNECTION WITH PPCKO OR THE USE OR OTHER DEALINGS IN
// PPCKO.
 
#ifndef _HH_MESH_HH
#define _HH_MESH_HH
#include "domain.hpp"
#include "meshGenerators.hpp"
#include <functional>
#include <vector>

 
namespace Geometry
{
class Mesh1D
{
public:
  Mesh1D() = default;

  Mesh1D(Domain1D const &d, unsigned int const &n);

  Mesh1D(Geometry::OneDMeshGenerator const &gf)
    : myDomain{gf.getDomain()}, myNodes{gf()} {};


  void reset(OneDMeshGenerator const &mg);

  unsigned int
  numNodes() const
  {
    return myNodes.size();
  }

  double
  operator[](int i) const
  {
    return myNodes[i];
  }

  std::vector<double>
  nodes() const
    {
      return myNodes;
    }

  std::vector<double>::iterator
  begin()
  {
    return myNodes.begin();
  }
  std::vector<double>::const_iterator
  cbegin() const
  {
    return myNodes.cbegin();
  }
  std::vector<double>::iterator
  end()
  {
    return myNodes.end();
  }
  std::vector<double>::const_iterator
  cend() const
  {
    return myNodes.cend();
  }

  Domain1D
  domain() const
  {
    return myDomain;
  }

  double hmin() const;
  double hmax() const;
 
private:
  Domain1D            myDomain;
  std::vector<double> myNodes;
};
 
} // namespace Geometry
#endif