Subversion Repositories gelsvn

Rev

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

Rev 601 Rev 607
Line 97... Line 97...
97
        }
97
        }
98
 
98
 
99
        Vec3d norm(0);
99
        Vec3d norm(0);
100
        N = one_ring.size();
100
        N = one_ring.size();
101
        Vec3d p0 = m.pos(v);
101
        Vec3d p0 = m.pos(v);
102
        for(int i = 0; i < N; ++i){
102
        for(size_t i = 0; i < N; ++i){
103
            Vec3d p1 = one_ring[i];
103
            Vec3d p1 = one_ring[i];
104
            Vec3d p2 = one_ring[(i+1) % N];
104
            Vec3d p2 = one_ring[(i+1) % N];
105
            Vec3d e0 = normalize(p1 - p0);
105
            Vec3d e0 = normalize(p1 - p0);
106
            Vec3d e1 = normalize(p2 - p0);
106
            Vec3d e1 = normalize(p2 - p0);
107
            norm += cross(e0, e1) * acos(dot(e0, e1));
107
            norm += cross(e0, e1) * acos(dot(e0, e1));