4fdd4c2adc
git-svn-id: http://www.openmesh.org/svnrepo/OpenMesh/trunk@362 fdac6126-5c0c-442c-9429-916003d36597
137 lines
5.5 KiB
C++
137 lines
5.5 KiB
C++
/*===========================================================================*\
|
|
* *
|
|
* OpenMesh *
|
|
* Copyright (C) 2001-2011 by Computer Graphics Group, RWTH Aachen *
|
|
* www.openmesh.org *
|
|
* *
|
|
*---------------------------------------------------------------------------*
|
|
* This file is part of OpenMesh. *
|
|
* *
|
|
* OpenMesh is free software: you can redistribute it and/or modify *
|
|
* it under the terms of the GNU Lesser General Public License as *
|
|
* published by the Free Software Foundation, either version 3 of *
|
|
* the License, or (at your option) any later version with the *
|
|
* following exceptions: *
|
|
* *
|
|
* If other files instantiate templates or use macros *
|
|
* or inline functions from this file, or you compile this file and *
|
|
* link it with other files to produce an executable, this file does *
|
|
* not by itself cause the resulting executable to be covered by the *
|
|
* GNU Lesser General Public License. This exception does not however *
|
|
* invalidate any other reasons why the executable file might be *
|
|
* covered by the GNU Lesser General Public License. *
|
|
* *
|
|
* OpenMesh 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 Lesser General Public License for more details. *
|
|
* *
|
|
* You should have received a copy of the GNU LesserGeneral Public *
|
|
* License along with OpenMesh. If not, *
|
|
* see <http://www.gnu.org/licenses/>. *
|
|
* *
|
|
\*===========================================================================*/
|
|
|
|
/*===========================================================================*\
|
|
* *
|
|
* $Revision$ *
|
|
* $Date$ *
|
|
* *
|
|
\*===========================================================================*/
|
|
|
|
/** \file ModQuadricT.cc
|
|
Bodies of template member function.
|
|
*/
|
|
|
|
//=============================================================================
|
|
//
|
|
// CLASS ModQuadric - IMPLEMENTATION
|
|
//
|
|
//=============================================================================
|
|
|
|
#define OPENMESH_DECIMATER_MODQUADRIC_CC
|
|
|
|
//== INCLUDES =================================================================
|
|
|
|
#include <OpenMesh/Tools/Decimater/ModQuadricT.hh>
|
|
|
|
|
|
//== NAMESPACE ===============================================================
|
|
|
|
namespace OpenMesh { // BEGIN_NS_OPENMESH
|
|
namespace Decimater { // BEGIN_NS_DECIMATER
|
|
|
|
|
|
//== IMPLEMENTATION ==========================================================
|
|
|
|
|
|
template<class DecimaterType>
|
|
void
|
|
ModQuadricT<DecimaterType>::
|
|
initialize()
|
|
{
|
|
using Geometry::Quadricd;
|
|
// alloc quadrics
|
|
if (!quadrics_.is_valid())
|
|
Base::mesh().add_property( quadrics_ );
|
|
|
|
// clear quadrics
|
|
typename Mesh::VertexIter v_it = Base::mesh().vertices_begin(),
|
|
v_end = Base::mesh().vertices_end();
|
|
|
|
for (; v_it != v_end; ++v_it)
|
|
Base::mesh().property(quadrics_, v_it).clear();
|
|
|
|
// calc (normal weighted) quadric
|
|
typename Mesh::FaceIter f_it = Base::mesh().faces_begin(),
|
|
f_end = Base::mesh().faces_end();
|
|
|
|
typename Mesh::FaceVertexIter fv_it;
|
|
typename Mesh::VertexHandle vh0, vh1, vh2;
|
|
typedef Vec3d Vec3;
|
|
double a,b,c,d, area;
|
|
|
|
for (; f_it != f_end; ++f_it)
|
|
{
|
|
fv_it = Base::mesh().fv_iter(f_it.handle());
|
|
vh0 = fv_it.handle(); ++fv_it;
|
|
vh1 = fv_it.handle(); ++fv_it;
|
|
vh2 = fv_it.handle();
|
|
|
|
Vec3 v0, v1, v2;
|
|
{
|
|
using namespace OpenMesh;
|
|
|
|
v0 = vector_cast<Vec3>(Base::mesh().point(vh0));
|
|
v1 = vector_cast<Vec3>(Base::mesh().point(vh1));
|
|
v2 = vector_cast<Vec3>(Base::mesh().point(vh2));
|
|
}
|
|
|
|
Vec3 n = (v1-v0) % (v2-v0);
|
|
area = n.norm();
|
|
if (area > FLT_MIN)
|
|
{
|
|
n /= area;
|
|
area *= 0.5;
|
|
}
|
|
|
|
a = n[0];
|
|
b = n[1];
|
|
c = n[2];
|
|
d = -(vector_cast<Vec3>(Base::mesh().point(vh0))|n);
|
|
|
|
Quadricd q(a, b, c, d);
|
|
q *= area;
|
|
|
|
Base::mesh().property(quadrics_, vh0) += q;
|
|
Base::mesh().property(quadrics_, vh1) += q;
|
|
Base::mesh().property(quadrics_, vh2) += q;
|
|
}
|
|
}
|
|
|
|
|
|
//=============================================================================
|
|
} // END_NS_DECIMATER
|
|
} // END_NS_OPENMESH
|
|
//=============================================================================
|