slab_tria.hh

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 - 2023 by the ideal.II authors
//
// This file is part of the ideal.II library.
//
// The ideal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 3.0 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of ideal.II.
//
// ---------------------------------------------------------------------
 
#ifndef INCLUDE_IDEAL_II_GRID_SLAB_TRIA_HH_
#define INCLUDE_IDEAL_II_GRID_SLAB_TRIA_HH_
 
#include <deal.II/grid/tria.h>
 
#include <list>
#include <memory>
 
namespace idealii::slab
{
  template <int dim>
  class Triangulation
  {
  public:
    Triangulation(std::shared_ptr<dealii::Triangulation<dim>> space_tria,
                  double                                      startpoint,
                  double                                      endpoint);
 
    Triangulation(std::shared_ptr<dealii::Triangulation<dim>> space_tria,
                  std::vector<double>                         step_sizes,
                  double                                      startpoint,
                  double                                      endpoint);
 
    Triangulation(const Triangulation &other);
    std::shared_ptr<dealii::Triangulation<dim>>
    spatial();
 
    std::shared_ptr<dealii::Triangulation<1>>
    temporal();
 
    double
    startpoint();
 
    double
    endpoint();
 
    void
    update_temporal_triangulation(std::vector<double> step_sizes,
                                  double              startpoint,
                                  double              endpoint);
 
  private:
    std::shared_ptr<dealii::Triangulation<dim>> _spatial_tria;
    std::shared_ptr<dealii::Triangulation<1>>   _temporal_tria;
    double                                      _startpoint;
    double                                      _endpoint;
  };
 
  template <int dim>
  using TriaIterator = typename std::list<Triangulation<dim>>::iterator;
} // namespace idealii::slab
 
#endif /* INCLUDE_IDEAL_II_GRID_SLAB_TRIA_HH_ */