vertex.cpp

/***************************************************************************
 *   Copyright (C) 2007 by Pablo Diaz-Gutierrez   *
 *   pablo@ics.uci.edu   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU Library General Public License as       *
 *   published by the Free Software Foundation; either version 2 of the    *
 *   License, or (at your option) any later version.                       *
 *                                                                         *
 *   This program 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 Library General Public     *
 *   License along with this program; if not, write to the                 *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

#ifdef _WIN32
#pragma warning(disable:4786)
#endif

#include <cassert>
#include <iostream>
#include "vertex.h"
#include "halfedge.h"
#include "face.h"

using namespace std;
using namespace Geometry;
using namespace HE;


Vector3Df Vertex::normal() const
{
        const_edge_circulator e = begin();
        const_edge_circulator sentinel = e;

        Vector3Df N(0,0,0);
        do
        {
                assert(this != (*e)->dst());
                Vector3Df v((*e)->tangent());
                ++e;
                
                // Do not use the normal of faces which are actually holes
                if ((*e)->face()->hole())
                        continue;

                Vector3Df v2((*e)->tangent());
                float a = v.angle(v2);
                Vector3Df n = v2 ^ v;
                N += a*n;
        } while (e != sentinel) ;

        N.normalize();
        return N;
}


ostream& HE::operator<< (ostream& out, const Vertex& v)
{
        //return out << '[' << v.index() << ": " << v.position() << ", e" << v.edge() << ']';
        return out << '[' << v.index() << ": " << v.position() << ']';
}

int Vertex::valence() const
{
        const_edge_circulator e = begin();
        const_edge_circulator sentinel = e;

        int n(0);
        do
        {
                n++;
                ++e;
        } while (e != sentinel) ;

        return n;
}


HalfEdge* Vertex::edgeTo(int v)
{
  edge_circulator e = begin();
  edge_circulator sentinel = e;

  do
  {
    if (v == (*e)->dst()->index())
      return *e;
    ++e;
  } while (e != sentinel) ;

  return 0;
}

const HalfEdge* Vertex::edgeTo(int v) const
{
  const_edge_circulator e = begin();
  const_edge_circulator sentinel = e;

  do
  {
    if (v == (*e)->dst()->index())
      return *e;
    ++e;
  } while (e != sentinel) ;

  return 0;
}


void Vertex::edge_circulator::operator++()
{
        _here = _here->twin()->next();
}

void Vertex::const_edge_circulator::operator++()
{
        _here = _here->twin()->next();
}


bool Vertex::isBoundary() const
{
  const_edge_circulator e = begin();
  const_edge_circulator sentinel = e;

  do
  {
    if ((*e)->isBoundary())
      return true;
    ++e;
  } while (e != sentinel) ;

  return false;  
}

---