refinementcontainer.h

Go to the documentation of this file.

#ifndef _REFINEMENTCONTAINER_H_
#define _REFINEMENTCONTAINER_H_

#include <deal.II/lac/vector.h>
#include "dopetypes.h"
#include <limits>

namespace DOpE
{
  class RefinementContainer
  {
    public:
      RefinementContainer(DOpEtypes::RefinementType ref_type =
          DOpEtypes::RefinementType::global);
      virtual
      ~RefinementContainer()
      {
      }

      virtual const dealii::Vector<float>&
      GetLocalErrorIndicators() const;
      virtual double
      GetTopFraction() const;
      virtual double
      GetBottomFraction() const;
      virtual unsigned int
      GetMaxNElements() const;
      virtual double
      GetConvergenceOrder() const;

      DOpEtypes::RefinementType
      GetRefType() const;

      bool
      UsesCoarsening() const;
    protected:
      bool _coarsening;
    private:
      const dealii::Vector<float> _dummy;
      const DOpEtypes::RefinementType _ref_type;

  };

  /***************************************************************/

  class LocalRefinement : public RefinementContainer
  {
    public:
      virtual
      ~LocalRefinement()
      {
      }

      virtual const dealii::Vector<float>&
      GetLocalErrorIndicators() const;

    protected:
      LocalRefinement(const dealii::Vector<float>&,
          DOpEtypes::RefinementType ref_type);
    private:
      LocalRefinement();
      const dealii::Vector<float>& _indicators;
  };

  /***************************************************************/

  class RefineFixedFraction : public LocalRefinement
  {
    public:
      RefineFixedFraction(const dealii::Vector<float>& indicators,
          double top_fraction = 0.1, double bottom_fraction = 0.0,
          const unsigned int max_n_elements =
              std::numeric_limits<unsigned int>::max());

      virtual
      ~ RefineFixedFraction()
      {
      }

      virtual double
      GetTopFraction() const;
      virtual double
      GetBottomFraction() const;
    private:
      const double _top_fraction, _bottom_fraction;
      const unsigned int _max_n_elements;
  };

  /***************************************************************/

  class RefineFixedNumber : public LocalRefinement
  {
    public:
      RefineFixedNumber(const dealii::Vector<float>& indicators,
          double top_fraction = 0.1, double bottom_fraction = 0.0,
          const unsigned int max_n_elements =
              std::numeric_limits<unsigned int>::max());

      virtual
      ~ RefineFixedNumber()
      {
      }

      virtual double
      GetTopFraction() const;
      virtual double
      GetBottomFraction() const;

    private:
      const double _top_fraction, _bottom_fraction;
      const unsigned int _max_n_elements;
  };

  /***************************************************************/

  class RefineOptimized : public LocalRefinement
  {
    public:
      RefineOptimized(const dealii::Vector<float>& indicators,
          double convergence_order = 2.);

      virtual
      ~ RefineOptimized()
      {
      }

      virtual double
      GetConvergenceOrder() const;

    private:
      const double _convergence_order;
  };

}

#endif /* _RefinementContainer_H_ */