boundarymesh.C

Go to the documentation of this file.

#ifdef HAVE_CONFIG_H
#include <tumble-conf.h>
#endif /* HAVE_CONFIG_H */

#include "boundarymesh.h"

#include "datastore.h"
#include "spline.h"

using namespace std;

BoundaryMesh::BoundaryMesh(PersistantStore& pstore, DataStore *dstore) 
: BoundaryComplex(pstore), beziermesh(NULL), datastore(dstore)
{
}

void BoundaryMesh::set_bezier_mesh(BezierMesh* bez) {
  assert(!beziermesh);
  beziermesh = bez;
}

BoundaryVertex* BoundaryMesh::add_boundary_vertex(const Point2D& p)
{
    ControlPoint cp = datastore->add_cp(p);
    BoundaryVertex *bv = new BoundaryVertex(get_store(), cp);
    add_vertex(bv);
    return bv;
}

BoundaryEdge* BoundaryMesh::add_boundary_edge(BoundaryVertex *v0,BoundaryVertex *v1, const vector<Point2D> &newd, const vector<double> &newk, Movement _fixed, int _color, double _restlength)
{
  BoundaryEdge *be = add_boundary_edge(v0, v1, newd, newk);
  be->set_fixed(_fixed);
  be->set_color(_color);
  be->set_restlength( _restlength );

  return be;
}



BoundaryEdge* BoundaryMesh::add_boundary_edge(BoundaryVertex *v0,
                                              BoundaryVertex *v1,
                                              const vector<Point2D>& newd,
                                              const vector<double>& newk)
{
    /* check for bad values */
    assert(v0); assert(v1);
    assert(newd.size() >= 3);
    assert(newk.size() >= 2);

    /* override endpoints of newd */
    vector<ControlPoint> deboor;
    deboor.push_back(v0->get_control_point());
    for(unsigned i=1; i < newd.size()-1; i++){
        deboor.push_back( datastore->add_cp( newd[ i ] ) );
    }
    deboor.push_back(v1->get_control_point());

    /* Create new edge */
    BoundaryEdge *be = new BoundaryEdge(get_store(), v0, v1, 
                    deboor, newk, datastore);
    add_edge(be);
    return be;
}


BoundaryFace* BoundaryMesh::add_boundary_face(const vector<BoundaryEdge*> &edges,double minangle, int color)
{
  BoundaryFace *bf=new BoundaryFace(get_store(), minangle,color);
  BoundaryEdge **es = new BoundaryEdge*[edges.size()];
  for(unsigned i=0; i < edges.size(); i++) es[i]=edges[i];
  add_face(bf,es,edges.size());
  return bf;
}

void BoundaryMesh::delete_vertex(BoundaryVertex *v)
{
    while(v->num_edges() > 0) {
      delete_edge(v->get_any_edge());
    }
    datastore->rem_cp(v->get_control_point());
    BoundaryComplex::delete_vertex(v);
}


void BoundaryMesh::delete_edge(BoundaryEdge *e)
{
    QBSpline *spline=e->get_spline();
    for(unsigned i = 1; i < spline->b.size() - 1; i++) {
      datastore->rem_cp( spline->b[i] );
    }
    BoundaryComplex::delete_edge( e );
}

void BoundaryMesh::bbox(Point2D &top, Point2D &bot)
{
    Boundary_Edge_Hash_T::iterator i;
    for(i = get_edges_begin(); i != get_edges_end(); ++i){
        double x_max,x_min,y_max,y_min;
        (*i)->get_spline()->bbox(&y_max,&x_min,&y_min,&x_max);
        if(x_max > top.coords[0])  top.coords[0]=x_max;
        if(y_max > top.coords[1])  top.coords[1]=y_max;
        if(x_min < bot.coords[0])  bot.coords[0]=x_min;
        if(y_min < bot.coords[1])  bot.coords[1]=y_min;
    }
    Boundary_Vertex_Hash_T::iterator j;
    for(j = get_vertices_begin(); j!= get_vertices_end(); ++j){
      double x,y;
      x = (*j)->get_control_point()->x(); 
      y = (*j)->get_control_point()->y();
      if(x > top.coords[0])  top.coords[0]=x;
      if(y > top.coords[1])  top.coords[1]=y;
      if(x < bot.coords[0])  bot.coords[0]=x;
      if(y < bot.coords[1])  bot.coords[1]=y;
    }
}

ostream& operator<<(ostream &stream, const BoundaryMesh &bm)
{
    stream<<"BoundaryMesh:"<<endl;
    stream<<"  beziermesh: "<<bm.beziermesh<<endl;
    stream<<"  datastore: "<<bm.datastore<<endl;
    stream<<"  CellComplex: "<<endl;
    stream<<*((BoundaryComplex*)(&bm))<<endl;
    return stream;
}

---