orderizer.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_ORDERIZER_H
#define __H2D_ORDERIZER_H

#include "../space/space.h"
#include "global.h"
#ifndef NOGLUT
#include <pthread.h>
#endif
#include "../quadrature/quad_all.h"
#include "../mesh/traverse.h"

namespace Hermes
{
  namespace Hermes2D
  {
    namespace Views
    {
      class HERMES_API Orderizer
      {
      public:

        Orderizer();
        ~Orderizer();

        template<typename Scalar>
        void save_orders_vtk(SpaceSharedPtr<Scalar> space, const char* file_name);

        template<typename Scalar>
        void save_markers_vtk(SpaceSharedPtr<Scalar> space, const char* file_name);

        template<typename Scalar>
        void save_mesh_vtk(SpaceSharedPtr<Scalar> space, const char* file_name);

        void calc_vertices_aabb(double* min_x, double* max_x,
          double* min_y, double* max_y) const;

        template<typename Scalar>
        void process_space(SpaceSharedPtr<Scalar> space, bool show_edge_orders = false);

        int get_labels(int*& lvert, char**& ltext, double2*& lbox) const;

        int get_num_vertices();
        double3* get_vertices();

        void lock_data() const;
        void unlock_data() const;

        int3* get_triangles();
        int* get_triangle_markers();
        int get_num_triangles();

        int2* get_edges();
        int* get_edge_markers();
        int get_num_edges();

        void free();
      protected:
#ifndef NOGLUT
        mutable pthread_mutex_t data_mutex;
#endif
        int2* edges;
        int* edge_markers;
        void add_edge(int iv1, int iv2, int marker);

        void reallocate(MeshSharedPtr mesh);

        char  buffer[1000];
        char* labels[11][11];

        double3* verts;
        int  label_size, label_count;
        int* lvert;
        char** ltext;
        double2* lbox;

        int3* tris;
        int* tri_markers;

        int vertex_count, triangle_count, edges_count;
        int vertex_size, triangle_size, edges_size;

        void add_triangle(int iv0, int iv1, int iv2, int marker);

        static void calc_aabb(double* x, double* y, int stride, int num, double* min_x, double* max_x, double* min_y, double* max_y);

        int add_vertex();

        void make_vert(int & index, double x, double y, double val);
      };
    }
  }
}
#endif