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);

      virtual
      ~ RefineFixedFraction()
      {
      }

      virtual double
      GetTopFraction() const;
      virtual double
      GetBottomFraction() const;
    private:
      const double top_fraction_, bottom_fraction_;
  };

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

  class RefineFixedNumber : public LocalRefinement
  {
    public:
      RefineFixedNumber(const dealii::Vector<float>& indicators,
                        double top_fraction = 0.1, double bottom_fraction = 0.0);
      virtual
      ~ RefineFixedNumber()
      {
      }

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

    private:
      const double top_fraction_, bottom_fraction_;
  };

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

  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_ */