Subversion Repositories gelsvn

Rev

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

Rev 631 Rev 632
Line 163... Line 163...
163 
        size_t N=edges.size();
 163 
        size_t N=edges.size();
164 
        double angle_sum = 0;
 164 
        double angle_sum = 0;
165 
        for(size_t i = 0; i < N; ++i)
 165 
        for(size_t i = 0; i < N; ++i)
166 
        {
 166 
        {
167 
            double dot_prod =
167 
            double dot_prod =
168 
                s_max(-1.0, s_min(1.0, dot(edges[i],edges[(i+1)%N])));
 168 
                std::max(-1.0, std::min(1.0, dot(edges[i],edges[(i+1)%N])));
169 
            angle_sum += acos(dot_prod);
 169 
            angle_sum += acos(dot_prod);
170 
        }
 170 
        }
171 
        return (2*M_PI - angle_sum)/voronoi_area(m, v);
 171 
        return (2*M_PI - angle_sum)/voronoi_area(m, v);
172 
 
 172 
 
173 
    }
 173 
    }