/******************************************************************************* * CGoGN: Combinatorial and Geometric modeling with Generic N-dimensional Maps * * version 0.1 * * Copyright (C) 2009-2012, IGG Team, LSIIT, University of Strasbourg * * * * This library 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 2.1 of the License, or (at your * * option) any later version. * * * * This library 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 Lesser General Public License * * along with this library; if not, write to the Free Software Foundation, * * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. * * * * Web site: http://cgogn.unistra.fr/ * * Contact information: cgogn@unistra.fr * * * *******************************************************************************/ #include "Topology/generic/traversor/traversorCell.h" #include "Topology/generic/traversor/traversor2.h" #include "Algo/Geometry/basic.h" namespace CGoGN { namespace Algo { namespace Surface { namespace Geometry { template ATTR_TYPE computeLaplacianTopoVertex( typename PFP::MAP& map, Dart d, const VertexAttribute& attr) { ATTR_TYPE l(0) ; ATTR_TYPE value = attr[d] ; unsigned int wSum = 0 ; Traversor2VE t(map, d) ; for(Dart it = t.begin(); it != t.end(); it = t.next()) { l += attr[map.phi1(it)] - value ; ++wSum ; } l /= wSum ; return l ; } template ATTR_TYPE computeLaplacianCotanVertex( typename PFP::MAP& map, Dart d, const EdgeAttribute& edgeWeight, const VertexAttribute& vertexArea, const VertexAttribute& attr) { typedef typename PFP::REAL REAL; ATTR_TYPE l(0) ; REAL vArea = vertexArea[d] ; ATTR_TYPE value = attr[d] ; REAL wSum = 0 ; Traversor2VE t(map, d) ; for(Dart it = t.begin(); it != t.end(); it = t.next()) { REAL w = edgeWeight[it] / vArea ; l += (attr[map.phi1(it)] - value) * w ; wSum += w ; } l /= wSum ; return l ; } template void computeLaplacianTopoVertices( typename PFP::MAP& map, const VertexAttribute& attr, VertexAttribute& laplacian) { TraversorV t(map) ; for(Dart d = t.begin(); d != t.end(); d = t.next()) laplacian[d] = computeLaplacianTopoVertex(map, d, attr) ; } template void computeLaplacianCotanVertices( typename PFP::MAP& map, const EdgeAttribute& edgeWeight, const VertexAttribute& vertexArea, const VertexAttribute& attr, VertexAttribute& laplacian) { TraversorV t(map) ; for(Dart d = t.begin(); d != t.end(); d = t.next()) laplacian[d] = computeLaplacianCotanVertex(map, d, edgeWeight, vertexArea, attr) ; } template typename PFP::REAL computeCotanWeightEdge( typename PFP::MAP& map, Dart d, const VertexAttribute& position) { typedef typename PFP::VEC3 VEC3; typedef typename PFP::REAL REAL; if(map.isBoundaryEdge(d)) { const VEC3& p1 = position[d] ; const VEC3& p2 = position[map.phi1(d)] ; const VEC3& p3 = position[map.phi_1(d)] ; REAL cot_alpha = 1 / tan(Geom::angle(p1 - p3, p2 - p3)) ; return 0.5 * cot_alpha ; } else { const VEC3& p1 = position[d] ; const VEC3& p2 = position[map.phi1(d)] ; const VEC3& p3 = position[map.phi_1(d)] ; const VEC3& p4 = position[map.phi_1(map.phi2(d))] ; REAL cot_alpha = 1 / tan(Geom::angle(p1 - p3, p2 - p3)) ; REAL cot_beta = 1 / tan(Geom::angle(p2 - p4, p1 - p4)) ; return 0.5 * ( cot_alpha + cot_beta ) ; } } template void computeCotanWeightEdges( typename PFP::MAP& map, const VertexAttribute& position, EdgeAttribute& edgeWeight) { TraversorE t(map) ; for(Dart d = t.begin(); d != t.end(); d = t.next()) edgeWeight[d] = computeCotanWeightEdge(map, d, position) ; } } // namespace Geometry } // namespace Surface namespace Volume { namespace Geometry { template ATTR_TYPE computeLaplacianTopoVertex( typename PFP::MAP& map, Dart d, const VertexAttribute& attr) { ATTR_TYPE l(0) ; ATTR_TYPE value = attr[d] ; unsigned int wSum = 0 ; Traversor3VE t(map, d) ; for(Dart it = t.begin(); it != t.end(); it = t.next()) { l += attr[map.phi1(it)] - value ; ++wSum ; } l /= wSum ; return l ; } template void computeLaplacianTopoVertices( typename PFP::MAP& map, const VertexAttribute& attr, VertexAttribute& laplacian) { TraversorV t(map) ; for(Dart d = t.begin(); d != t.end(); d = t.next()) laplacian[d] = computeLaplacianTopoVertex(map, d, attr) ; } } // namespace Geometry } // namespace Volume } // namespace Algo } // namespace CGoGN