element_to_refine.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_ELEMENT_TO_REFINE_H
#define __H2D_ELEMENT_TO_REFINE_H

#include "global.h"

namespace Hermes
{
  namespace Hermes2D
  {
    enum RefinementType {
      H2D_REFINEMENT_P = 0,
      H2D_REFINEMENT_H = 1,
      H2D_REFINEMENT_H_ANISO_H = 2,
      H2D_REFINEMENT_H_ANISO_V = 3
    };


    class HERMES_API ElementToRefine {
    public:
      ElementToRefine();


      ElementToRefine(int id, unsigned short comp, bool space_only = false);

      ElementToRefine& operator=(const ElementToRefine& orig);

      bool valid;
      bool space_only;

      int id;
      unsigned short comp;
      RefinementType split;
      unsigned short refinement_polynomial_order[H2D_MAX_ELEMENT_SONS];
      unsigned short best_refinement_polynomial_order_type[4][H2D_MAX_ELEMENT_SONS];
      double errors[H2D_MAX_ELEMENT_SONS];


      unsigned short get_num_sons() const;


      static void copy_orders(unsigned short* dest, const unsigned short* src);
      static void copy_errors(double* dest, const double* src);
    private:
      bool refinement_polynomial_order_changed[H2D_MAX_ELEMENT_SONS];
      template<typename T> friend class Adapt;
      template<typename T, typename S> friend class AdaptSolver;
    };
  }
}
#endif