vectorizer.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_VECTORIZER_H
#define __H2D_VECTORIZER_H

#include "linearizer.h"

namespace Hermes
{
  namespace Hermes2D
  {
    namespace Views
    {
      class HERMES_API Vectorizer : public LinearizerBase
      {
      public:

        Vectorizer();
        ~Vectorizer();

        void process_solution(MeshFunction<double>* xsln, MeshFunction<double>* ysln, int xitem = H2D_FN_VAL_0, int yitem = H2D_FN_VAL_0, double eps = HERMES_EPS_NORMAL);

        int get_num_vertices();
        double4* get_vertices();

        int2* get_dashes();
        int get_num_dashes();

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

        void set_curvature_epsilon(double curvature_epsilon);

        double get_curvature_epsilon();

        void free();

      protected:
        double curvature_epsilon;

        int xitem, component_x, value_type_x;
        int yitem, component_y, value_type_y;

        double4* verts;  
        int2* dashes;

        int dashes_count; 
        int dashes_size; 

        int get_vertex(int p1, int p2, double x, double y, double xvalue, double yvalue);
        int create_vertex(double x, double y, double xvalue, double yvalue);
        void process_dash(int iv1, int iv2);

        int add_vertex();

        void add_dash(int iv1, int iv2);

        void process_triangle(MeshFunction<double>** fns, int iv0, int iv1, int iv2, int level,
          double* xval, double* yval, double* phx, double* phy, int* indices, bool curved);

        void process_quad(MeshFunction<double>** fns, int iv0, int iv1, int iv2, int iv3, int level,
          double* xval, double* yval, double* phx, double* phy, int* indices, bool curved);

        void find_min_max();
      };
    }
  }
}
#endif