Subversion Repositories gelsvn

Rev

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

Rev 601 Rev 607
Line 160... Line 160...
160 
            Vec3d e(normalize(m.pos(w.vertex()) - vertex));
 160 
            Vec3d e(normalize(m.pos(w.vertex()) - vertex));
161 
            edges.push_back(e);
 161 
            edges.push_back(e);
162 
        }
 162 
        }
163 
        size_t N=edges.size();
 163 
        size_t N=edges.size();
164 
        double angle_sum = 0;
 164 
        double angle_sum = 0;
165 
        for(int 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 
                s_max(-1.0, s_min(1.0, dot(edges[i],edges[(i+1)%N])));
169 
            angle_sum += acos(dot_prod);
 169 
            angle_sum += acos(dot_prod);
170 
        }
 170 
        }