Subversion Repositories gelsvn

Rev

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

Rev 572 Rev 585
Line 158... Line 158...
158 
        vector<Vec3d> edges;
 158 
        vector<Vec3d> edges;
159 
        for(HalfEdgeWalker w = m.halfedgewalker(v); !w.full_circle(); w = w.circulate_vertex_cw()){
 159 
        for(HalfEdgeWalker w = m.halfedgewalker(v); !w.full_circle(); w = w.circulate_vertex_cw()){
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 
        int 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(int 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])));
Line 225... Line 225...
225 
 
 225 
 
226 
        vector<Vec3d> points;
 226 
        vector<Vec3d> points;
227 
        for(HalfEdgeWalker w = m.halfedgewalker(v); !w.full_circle(); w = w.circulate_vertex_cw())
 227 
        for(HalfEdgeWalker w = m.halfedgewalker(v); !w.full_circle(); w = w.circulate_vertex_cw())
228 
            points.push_back(Vec3d(m.pos(w.vertex())));
 228 
            points.push_back(Vec3d(m.pos(w.vertex())));
229 
 
 229 
 
230 
        int N = points.size();
 230 
        size_t N = points.size();
231 
 
 231 
 
232 
        CVector b(N);
 232 
        CVector b(N);
233 
        // Compute the matrix of parameter values
233 
        // Compute the matrix of parameter values
234 
        CMatrix PMat(N, 3);
 234 
        CMatrix PMat(N, 3);
235 
        for(int i = 0; i < N; ++i){
 235 
        for(int i = 0; i < N; ++i){