precalc.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_PRECALC_H
#define __H2D_PRECALC_H

#include "../function/function.h"
#include "../shapeset/shapeset.h"

namespace Hermes
{
  namespace Hermes2D
  {
    enum SpaceType;
    class HERMES_API PrecalcShapeset : public Function<double>
    {
    public:
      SpaceType get_space_type() const;

      PrecalcShapeset(Shapeset* shapeset);

      PrecalcShapeset(PrecalcShapeset* master_pss);

      virtual ~PrecalcShapeset();

      virtual void set_active_element(Element* e);

      void set_active_shape(int index);

    private:
      virtual void set_quad_2d(Quad2D* quad_2d);

      virtual void free();

      virtual void handle_overflow_idx();


      int get_active_shape() const;

      Shapeset* get_shapeset() const;

      void set_master_transform();

      virtual int get_edge_fn_order(int edge);

      virtual void push_transform(int son);

      virtual void pop_transform();

      Shapeset* shapeset;

      LightArray<std::map<uint64_t, LightArray<Node*>*>*> tables;

      int index;

      int max_index[H2D_NUM_MODES];

      PrecalcShapeset* master_pss;

      bool is_slave() const;

      virtual void precalculate(int order, int mask);

      void update_max_index();

      void force_transform(uint64_t sub_idx, Trf* ctm);

      friend class RefMap;
      template<typename T> friend class KellyTypeAdapt;
      template<typename T> friend class Adapt;
      template<typename T> friend class Func;
      template<typename T> friend class Solution;
      template<typename T> friend class DiscontinuousFunc;
      template<typename T> friend class DiscreteProblem;
      template<typename T> friend class NeighborSearch;
      friend class CurvMap;
    };
  }
}
#endif