mesh.cpp

/*:
**     mesh.cpp - mesh ops
**
**   Copyright (C) 2004 Luiz Velho.
*/

#include "a48.h"
#include <map>

using namespace A48;

Mesh::Mesh(Surface* s) : surf_(s)
{
  int np, nf, *fcs;
  surf()->base_mesh(&np, &fcs, &nf);
  set_base_mesh(np, fcs, nf);
  if (!is_triquad())
    make_triquad();
  new_front();
}

Mesh::~Mesh()
{
  for (VertexIter p = vc_.begin(); p != vc_.end(); p++) 
    delete *p;
  vc_.clear();
  for (EdgeIter e = ec_.begin(); e != ec_.end(); e++)
    delete *e;
  ec_.clear();
  for (FaceIter f = fc_.begin(); f != fc_.end(); f++)
    delete *f;
  fc_.clear();
}

Face* Mesh::add_face(Hedge *e0, Hedge *e1, Hedge *e2)
{
  Face *f = new Face(e0, e1, e2);
  add_face(f);
  return f;
}

Hedge* Mesh::add_edge(Vertex *v0, Vertex *v1)
{
  Edge *e = new Edge(v0, v1);
  add_edge(e);
  return e->hedge(0);
}

Vertex* Mesh::add_vertex(void)
{
  Vertex *v = surf()->new_vertex();
  add_vertex(v);
  return v;
}

void Mesh::link_mesh()
{
  for (FaceIter f = fc_.begin(); f != fc_.end(); f++)
    (*f)->link_star_verts();
  for (EdgeIter e = ec_.begin(); e != ec_.end(); e++) {
    Hedge* h = (*e)->hedge(0);
    if (h->is_bdry())
      h->dst()->set_star(h);
  }
}

Hedge* Mesh::get_hedge(int i0, int i1, Vertex* verts[], HedgeMap *hedges)
{
  bool mate = false;
  if (i0 > i1) {
    std::swap<int>(i0, i1);
    mate = true;
  }
  Hedge* he; HedgeMap::iterator ei  = hedges->find(Ipair(i0,i1));
  if (ei == hedges->end())
    (*hedges)[Ipair(i0,i1)] = he = add_edge(verts[i0], verts[i1]);
  else
    he = (*ei).second;
  return (mate)? he->mate() : he ;
}

void Mesh::put_face(int i0, int i1, int i2, Vertex* verts[], HedgeMap* hedges)
{
  Hedge* e0 = get_hedge(i1, i2, verts, hedges);
  Hedge* e1 = get_hedge(i2, i0, verts, hedges);
  Hedge* e2 = get_hedge(i0, i1, verts, hedges);
  add_face(e0, e1, e2);
}

void Mesh::set_base_mesh(int np, int* tris, int nf)
{
  Vertex** verts = new Vertex*[np];
  for (int i=0; i < np; i++) {
    verts[i] = surf()->new_vertex();
    surf()->sample(i, verts[i]);
    add_vertex(verts[i]);
  }
  HedgeMap hedges;
  for (int i=0; i < nf; i++, tris += 3) 
    put_face(tris[2], tris[0],tris[1], verts, &hedges);
  link_mesh();
  delete(verts);
}

---