quad.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_QUAD_H
#define __H2D_QUAD_H

#include "../global.h"
namespace Hermes
{
  namespace Hermes2D
  {
    enum GIP2DIndices {
      H2D_GIP2D_X = 0, 
      H2D_GIP2D_Y = 1, 
      H2D_GIP2D_W = 2  
    };

    static int H2D_GIP1D_X = 0;
    static int H2D_GIP1D_W = 1;

    const int g_max_quad = 24;
    const int g_max_tri = 20;

    class HERMES_API Quad1D
    {
    public:

      inline double2* get_points(int order) const { return tables[order]; }
      inline int get_num_points(int order) const { return np[order]; };

      inline int get_max_order() const { return max_order; }
      inline double get_ref_vertex(int n) const { return ref_vert[n]; }

    protected:

      double2** tables;
      int* np;

      double ref_vert[H2D_NUM_MODES];
      int max_order;

      virtual void dummy_fn() = 0; // to prevent this class from being instantiated
    };

    class HERMES_API Quad2D
    {
    public:
      inline int get_num_points(int order, ElementMode2D mode)  const { assert(order < num_tables[mode]); return np[mode][order]; };
      inline double3* get_points(int order, ElementMode2D mode) const { assert(order < num_tables[mode]); return tables[mode][order]; }
      inline int get_edge_points(int edge, int order, ElementMode2D mode) {assert(order < num_tables[mode]);  return  max_order[mode]+1 + (3*(1-mode) + 4*mode)*order + edge;}

      inline int get_max_order(ElementMode2D mode) const { return max_order[mode]; }
      inline int get_safe_max_order(ElementMode2D mode) const { return safe_max_order[mode]; }
      inline int get_num_tables(ElementMode2D mode) const { return num_tables[mode]; }

      inline double2* get_ref_vertex(int n, ElementMode2D mode) { return &ref_vert[mode][n]; }

    protected:
      double3*** tables;
      int** np;

      int num_tables[2];
      int max_order[2], safe_max_order[2];
      int max_edge_order;

      double2 ref_vert[2][H2D_MAX_NUMBER_VERTICES];
    };
  }
}
#endif