/Users/mourrain/Devel/mmx/shape/include/shape/semialgebraic_domain2d_axel.hpp

Go to the documentation of this file.

#pragma once


# include  <shape/graph.hpp>
# include  <shape/viewer_axel.hpp>
# include  <shape/cell2d.hpp>

# define TMPL template<class T>
# define AXEL viewer<axel,K>
# define BoundingBox (*T)::BoundingBox

//#define CHECK if ( ! v[i]->is_intersected() || v[i]->is_active()  ) continue
//#define CHECK if ( ! v[i]->is_intersected() ) continue
//#define CHECK if ( ! v[i]->is_active()  ) continue
#define CHECK if (false) {}

namespace mmx {
  
  namespace shape {
    
    // Function to be passed to the DFS procedure
    template<class T>
    void write_ok(gNode<T>* v)
    {
      std::cout<<"Ok"<<std::endl;
    }
  }//namespace shape
  
  
  namespace shape {
    
    template<class K, class T> AXEL&
    print_boxes(AXEL& out, Graph<T> g) {
      
      
      //Seq<T> vertices;
      double xmin,xmax,ymin,ymax;
      Seq<T> v;
      
      g.dfs(v);
        //g.dfs(& shape::write_ok);
        //g.vertex_list(v);
      
      unsigned n=0;
      
      
      for (unsigned i=0; i<v.size(); i++)
      {
        CHECK;
        n+=4;
      }
      
      out<<" <curve type=\"mesh\" color=\"255 255 0\" >\n<vect>\nVECT\n";
      out<<n   <<" "
         <<2*n <<" "
         <<n   <<"\n";
      
      for(unsigned i=0; i<n; i++) out<<"2 ";//
      out<<"\n";
      
      for(unsigned i=0; i<n; i++) out<<"1 ";//
      out<<"\n";
      
                        //print edges
      for (unsigned i=0; i<v.size(); i++)
      {
        CHECK;
        
        xmin= v[i]->boundingBox().xmin();
        xmax= v[i]->boundingBox().xmax();
        ymin= v[i]->boundingBox().ymin();
        ymax= v[i]->boundingBox().ymax();
        
        //std::cout <<xmin  <<" "<<xmax  <<" "<<ymin  <<" "<<ymax <<"\n";
        
        out <<xmin  <<" "<<ymin <<" 0 "
            <<xmax  <<" "<<ymin <<" 0 "         
            <<"\n";
        
          out <<xmax  <<" "<<ymin  <<" 0 "
              <<xmax  <<" "<<ymax  <<" 0 "              
              <<"\n";
          
          out <<xmax  <<" "<<ymax  <<" 0 "
              <<xmin  <<" "<<ymax  <<" 0 "              
              <<"\n";
          
          out <<xmin  <<" "<<ymax  <<" 0 "
              <<xmin  <<" "<<ymin  <<" 0 "              
              <<"\n";
      }
      
      
      for(unsigned i=0; i<n; i++) 
        out<< "0.314 0.979 1 1\n";
      
      out<<" </vect>\n </curve>\n";
      
      return out;
    }
    
  }
}

# undef TMPL
# undef CHECK
# undef BoundingBox
# undef AXEL