ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/ot/wfconv.c
(Generate patch)

Comparing ray/src/ot/wfconv.c (file contents):
Revision 2.10 by greg, Fri Nov 8 16:49:04 2013 UTC vs.
Revision 2.11 by greg, Fri Jan 24 01:26:44 2014 UTC

# Line 9 | Line 9 | static const char RCSid[] = "$Id$";
9   #include "copyright.h"
10   #include "standard.h"
11   #include "cvmesh.h"
12 + #include "triangulate.h"
13   #include <ctype.h>
14  
15   typedef int     VNDX[3];        /* vertex index (point,map,normal) */
# Line 244 | Line 245 | cvtndx(                                /* convert vertex string to index */
245          return(1);
246   }
247  
248 + /* determine dominant axis for triangle */
249 + static int
250 + dominant_axis(char *v1, char *v2, char *v3)
251 + {
252 +        int     v1i = atoi(v1), v2i = atoi(v2), v3i = atoi(v3);
253 +        FVECT   e1, e2, vn;
254 +        int     i, imax;
255  
256 +        VSUB(e1, vlist[v2i], vlist[v1i]);
257 +        VSUB(e2, vlist[v3i], vlist[v2i]);
258 +        VCROSS(vn, e1, e2);
259 +        for (i = imax = 2; i--; )
260 +                if (vn[i]*vn[i] > vn[imax]*vn[imax])
261 +                        imax = i;
262 +        return(vn[imax]*vn[imax] > FTINY ? imax : -1);
263 + }
264 +
265 + /* callback for triangle output from polygon */
266   static int
267 + tri_out(const Vert2_list *tp, int a, int b, int c)
268 + {
269 +        return( puttri( ((char **)tp->p)[a],
270 +                        ((char **)tp->p)[b],
271 +                        ((char **)tp->p)[c] ) );
272 + }
273 +
274 + static int
275   putface(                                /* put out an N-sided polygon */
276          int     ac,
277          char    **av
278   )
279   {
280 <        char            *cp;
281 <        int     i;
280 >        Vert2_list      *poly = polyAlloc(ac);
281 >        int             i, ax, ay;
282  
283 <        while (ac > 3) {                /* break into triangles */
284 <                if (!puttri(av[0], av[1], av[2]))
283 >        if (poly == NULL)
284 >                return(0);
285 >        poly->p = (void *)av;
286 >        for (i = ac-3; i >= 0; i--)     /* identify dominant axis */
287 >                if ((ax = dominant_axis(av[i], av[i+1], av[i+2])) >= 0)
288 >                        break;
289 >        if (ax < 0)
290 >                return(0);
291 >        if (++ax >= 3) ax = 0;
292 >        ay = ax;
293 >        if (++ay >= 3) ay = 0;
294 >        for (i = 0; i < ac; i++) {      /* convert to 2-D polygon */
295 >                VNDX    vi;
296 >                if (!cvtndx(vi, av[i])) {
297 >                        error(WARNING, "bad vertex reference");
298 >                        polyFree(poly);
299                          return(0);
300 <                ac--;                   /* remove vertex & rotate */
301 <                cp = av[0];
302 <                for (i = 0; i < ac-1; i++)
263 <                        av[i] = av[i+2];
264 <                av[i] = cp;
300 >                }
301 >                poly->v[i].mX = vlist[vi[0]][ax];
302 >                poly->v[i].mY = vlist[vi[0]][ay];
303          }
304 <        return(puttri(av[0], av[1], av[2]));
304 >                                        /* break into triangles & output */
305 >        if (!polyTriangulate(poly, &tri_out))
306 >                return(0);
307 >        polyFree(poly);
308 >        return(1);
309   }
310  
311  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines