selector.h

// This file is part of Hermes2D.
//
// Hermes2D is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 2 of the License, or
// (at your option) any later version.
//
// Hermes2D is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with Hermes2D.  If not, see <http://www.gnu.org/licenses/>.

#ifndef __H2D_REFINEMENT_SELECTOR_H
#define __H2D_REFINEMENT_SELECTOR_H

#ifndef _MSC_VER
#include "../mesh/refinement_type.h"

namespace Hermes
{
  namespace Hermes2D
  {
    class ElementToRefine;
    struct Element;
    template<typename Scalar> class Solution;
  }
}
#else
#include "../mesh/element_to_refine.h"
#endif
#include "../mesh/mesh.h"

namespace Hermes
{
  namespace Hermes2D
  {
#define H2DRS_DEFAULT_ORDER -1 
#define H2DRS_MAX_ORDER 9 


    namespace RefinementSelectors {

      template<typename Scalar>
      class HERMES_API Selector : public Hermes::Mixins::Loggable, public Hermes::Mixins::TimeMeasurable
      {
      public:
        virtual ~Selector() {};
      protected:
        const int max_order; 


        Selector(int max_order = H2DRS_DEFAULT_ORDER) : max_order(max_order), isAClone(false) {};

        virtual Selector<Scalar>* clone() = 0;


        virtual bool select_refinement(Element* element, int quad_order, Solution<Scalar>* rsln, ElementToRefine& refinement) = 0;


        virtual void generate_shared_mesh_orders(const Element* element, const int orig_quad_order, const int refinement, int tgt_quad_orders[H2D_MAX_ELEMENT_SONS], const int* suggested_quad_orders) = 0;

        template<typename T> friend class Adapt;
        template<typename T> friend class KellyTypeAdapt;
        
      protected:
        bool isAClone;
      };

      template<typename Scalar>
      class HERMES_API HOnlySelector : public Selector<Scalar> {
      public:
        HOnlySelector() : Selector<Scalar>() {};

        virtual Selector<Scalar>* clone();

      protected:

        virtual bool select_refinement(Element* element, int quad_order, Solution<Scalar>* rsln, ElementToRefine& refinement);


        virtual void generate_shared_mesh_orders(const Element* element, const int orig_quad_order, const int refinement, int tgt_quad_orders[H2D_MAX_ELEMENT_SONS], const int* suggested_quad_orders);
        template<typename T> friend class Adapt;
        template<typename T> friend class KellyTypeAdapt;
      };

      template<typename Scalar>
      class HERMES_API POnlySelector : public Selector<Scalar> {
        const int order_h_inc; 
        const int order_v_inc; 
      public:

        POnlySelector(int max_order, int order_h_inc, int order_v_inc);

        virtual Selector<Scalar>* clone();

      protected:

        virtual bool select_refinement(Element* element, int quad_order, Solution<Scalar>* rsln, ElementToRefine& refinement);


        virtual void generate_shared_mesh_orders(const Element* element, const int orig_quad_order, const int refinement, int tgt_quad_orders[H2D_MAX_ELEMENT_SONS], const int* suggested_quad_orders);
        template<typename T> friend class Adapt;
        template<typename T> friend class KellyTypeAdapt;
      };
    }
  }
}
#endif