include/shape/intersector_parametric.hpp

Go to the documentation of this file.

#ifndef SHAPE_INTERSECTOR_PARAMETRIC_HPP
#define SHAPE_INTERSECTOR_PARAMETRIC_HPP

# include <shape/ssi/ssi_dsearch.hpp>
# include <shape/ssi/ssiqtsl.hpp>

#define TMPL template<class C, class V>

namespace mmx { namespace shape {

template<class C,class V= default_env>
struct intersector_parametric {
    typedef typename use<ssi_def,C,V>::ParametricSurface Input;
    typedef typename use<ssi_def,C,V>::Mesh Output;
    typedef typename use<ssi_def,C,V>::Mesh::Point Point;

    intersector_parametric(): u_spl(200), v_spl(200) {
        input1 = NULL;
        input2 = NULL;
        output = NULL;
    }

    intersector_parametric(Input *s): u_spl(200), v_spl(200), input1(s), input2(s) {
        output = NULL;
    }

    intersector_parametric(Input *s1, Input* s2): u_spl(200), v_spl(200), input1(s1), input2(s2) {
        output = NULL;
    }
    void set_sample(int n) { u_spl=n; v_spl=n;}
    void set_sample(int n, int m) { u_spl=n; v_spl=m;}

    void set_input (Input* s) { input1 = s; input2 = s; }

    void set_input (Input* s1, Input* s2) { input1 = s1; input2 = s2; }
    void set_input1(Input* s) { input1 = s; }
    void set_input2(Input* s) { input2 = s; }

    Output* get_output() {return output;}

    void run ();

private:
    void run_inter (void);
    void run_self  (void);
    int u_spl, v_spl;
    Input  * input1, * input2;
    Output * output;
};

TMPL void intersector_parametric<C,V>::run(void) {
    if( input1 == input2 )
        this->run_self();
    else
        this->run_inter();
}
TMPL void intersector_parametric<C,V>::run_self(void) {
    mmx::ssi::dsearch<C,V> r(input1,u_spl,v_spl);

    output = new Output;

    C x,y,z;
    for ( unsigned i = 0; i < r.nbcurves; i++ )
        for ( unsigned k = 0; k < r.sizes[i]; k++ )
        {
            //std::cout<< r.lefts[i][k][0] <<" "<< r.lefts[i][k][1]<<std::endl;
            input1->eval(x,y,z, r.lefts[i][k][0],r.lefts[i][k][1]);
            output->push_back_vertex(new Point(x,y,z));
            Point* p= new Point;input1->eval(p->x(),p->y(),p->z(), r.lefts[i][k][0],r.lefts[i][k][1]);
            //std::cout<<"Eval "<<p->x()<<" "<<p->y()<< " "<<p->z()<<std::endl;
            //output->push_back_vertex(p);

        };

    int np=0;
    for ( unsigned i = 0; i < r.nbcurves; i++ )
        for (unsigned k = 0; k < r.sizes[i]; k+=2 )
        {
            output->push_back_edge(np,np+1);
            np+=2;
        };

}

TMPL void intersector_parametric<C,V>::run_inter(void) {
    std::cerr<<"Intersect:"<< u_spl <<std::endl;

    ssiqtsl<C,V> r(input1, input2, u_spl);

    output = new Output;
    int nbv =r.spcs.size();

    for ( int i = 0; i < nbv; i ++) {
        Point* p = new Point(r.spcs[i][0], r.spcs[i][1], r.spcs[i][2]);
        output->push_back_vertex(p);
    }

    for ( int i = 0; i < nbv; i+=2) {
        output->push_back_edge(i, i+1);
    }

}

//=================================================================
}
              }

#undef TMPL

#endif // SHAPE_INTERSECTOR_PARAMETRIC_HPP

---