Subversion Repositories gelsvn

Rev

Rev 526 | Rev 543 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 526 Rev 529
Line 264... Line 264...
264 
            curvature_tensors[*v] = curvature_tensor_from_edges(m, *v);
 264 
            curvature_tensors[*v] = curvature_tensor_from_edges(m, *v);
265 
    }
 265 
    }
266 
 
 266 
 
267 
    void smooth_curvature_tensors(const Manifold& m, VertexAttributeVector<Mat3x3d>& curvature_tensors)
 267 
    void smooth_curvature_tensors(const Manifold& m, VertexAttributeVector<Mat3x3d>& curvature_tensors)
268 
    {
 268 
    {
269 
//        assert(curvature_tensors.size() == m.no_vertices(false));
269 
        assert(curvature_tensors.size() == m.total_vertices());
270 
        VertexAttributeVector<Mat3x3d> tmp_curvature_tensors(m.total_vertices());
 270 
        VertexAttributeVector<Mat3x3d> tmp_curvature_tensors(m.total_vertices());
271 
        double tmp_area;
 271 
        double tmp_area;
272 
 
 272 
 
273 
        for(VertexIDIterator v = m.vertices_begin(); v != m.vertices_end(); ++v){
 273 
        for(VertexIDIterator v = m.vertices_begin(); v != m.vertices_end(); ++v){
274 
            if(boundary(m, *v))
 274 
            if(boundary(m, *v))
Line 343... Line 343...
343 
                                const VertexAttributeVector<Mat3x3d>& curvature_tensors,
 343 
                                const VertexAttributeVector<Mat3x3d>& curvature_tensors,
344 
                                VertexAttributeVector<Vec3d>& min_curv_direction,
 344 
                                VertexAttributeVector<Vec3d>& min_curv_direction,
345 
                                VertexAttributeVector<Vec3d>& max_curv_direction,
 345 
                                VertexAttributeVector<Vec3d>& max_curv_direction,
346 
                                VertexAttributeVector<double>& curvature)
 346 
                                VertexAttributeVector<double>& curvature)
347 
    {
 347 
    {
348 
//        assert(curvature_tensors.size() == m.no_vertices(false));
348 
        assert(curvature_tensors.size() == m.total_vertices());
349 
 
 349 
 
350 
        double max_val = -1e30;
 350 
        double max_val = -1e30;
351 
        for(VertexIDIterator v = m.vertices_begin(); v != m.vertices_end(); ++v){
 351 
        for(VertexIDIterator v = m.vertices_begin(); v != m.vertices_end(); ++v){
352 
            Mat3x3d C,Q,L;
 352 
            Mat3x3d C,Q,L;
353 
            C = curvature_tensors[*v];
 353 
            C = curvature_tensors[*v];