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

Comparing ray/src/cv/mgf2meta.c (file contents):
Revision 2.3 by greg, Tue Apr 11 13:32:47 1995 UTC vs.
Revision 2.13 by greg, Thu Apr 21 03:11:55 2022 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1995 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Convert MGF (Materials and Geometry Format) to Metafile 2-d graphics
6   */
7  
8   #include <stdio.h>
9 + #include <stdlib.h>
10 + #include <string.h>
11   #include <math.h>
13 – #include "mgflib/parser.h"
12  
13 < #define MX(v)   (int)(((1<<14)-1)*(v)[(proj_axis+1)%3])
14 < #define MY(v)   (int)(((1<<14)-1)*(v)[(proj_axis+2)%3])
13 > #include "meta.h"
14 > #include "random.h"
15 > #include "mgf_parser.h"
16 > #include "plocate.h" /* XXX shouldn't this rather be in rtmath.h? */
17  
18 < int     r_face();
18 > #define MSIZE   ((1<<14)-1)
19 > #define MX(v)   (int)(MSIZE*(v)[(proj_axis+1)%3])
20 > #define MY(v)   (int)(MSIZE*(v)[(proj_axis+2)%3])
21 >
22   int     proj_axis;
23   double  limit[3][2];
24   int     layer;
25 + long    rthresh = 1;
26  
27   extern int      mg_nqcdivs;
28  
29 + static int r_face(int ac, char **av);
30 + static void newlayer(void);
31 + static int doline(int v1x, int v1y, int v2x, int v2y);
32  
33 < main(argc, argv)                /* convert files to stdout */
34 < int     argc;
35 < char    *argv[];
33 >
34 > int
35 > main(           /* convert files to stdout */
36 >        int     argc,
37 >        char    *argv[]
38 > )
39   {
40          int     i;
41                                  /* initialize dispatch table */
# Line 36 | Line 46 | char   *argv[];
46          mg_nqcdivs = 3;         /* reduce object subdivision */
47          mg_init();              /* initialize the parser */
48                                          /* get arguments */
49 +        if (argc > 9 && !strcmp(argv[1], "-t")) {
50 +                rthresh = atof(argv[2])*MSIZE + 0.5;
51 +                rthresh *= rthresh;
52 +                argv += 2;
53 +                argc -= 2;
54 +        }
55          if (argc < 8 || (proj_axis = argv[1][0]-'x') < 0 || proj_axis > 2)
56                  goto userr;
57          limit[0][0] = atof(argv[2]); limit[0][1] = atof(argv[3]);
# Line 49 | Line 65 | char   *argv[];
65                  for (i = 8; i < argc; i++) {
66                          if (mg_load(argv[i]) != MG_OK)
67                                  exit(1);
68 <                        if (++layer >= 16) {
53 <                                mendpage();
54 <                                layer = 0;
55 <                        }
68 >                        newlayer();
69                  }
70          mendpage();                     /* print page */
71          mdone();                        /* close output */
72          exit(0);
73   userr:
74 <        fprintf(stderr, "Usage: %s {x|y|z} xmin xmax ymin ymax zmin zmax [file.mgf] ..\n",
75 <                        argv[0]);
74 >        fputs("Usage: mgf2meta [-t thresh] {x|y|z} xmin xmax ymin ymax zmin zmax [file.mgf] ..\n",
75 >                        stderr);
76          exit(1);
77   }
78  
79  
80   int
81 < r_face(ac, av)                  /* convert a face */
82 < int     ac;
83 < char    **av;
81 > r_face(                 /* convert a face */
82 >        int     ac,
83 >        char    **av
84 > )
85   {
86          static FVECT    bbmin = {0,0,0}, bbmax = {1,1,1};
87          register int    i, j;
88          register C_VERTEX       *cv;
89          FVECT   v1, v2, vo;
76 –        int     newline = 1;
90  
91          if (ac < 4)
92                  return(MG_EARGC);
# Line 91 | Line 104 | char   **av;
104                          v2[j] = (v2[j] - limit[j][0])/(limit[j][1]-limit[j][0]);
105                  VCOPY(v1, vo);
106                  VCOPY(vo, v2);
107 <                if (clip(v1, v2, bbmin, bbmax)) {
108 <                        mline(MX(v1), MY(v1), layer/4, 0, layer%4);
96 <                        mdraw(MX(v2), MY(v2));
97 <                }
107 >                if (clip(v1, v2, bbmin, bbmax))
108 >                        doline(MX(v1), MY(v1), MX(v2), MY(v2));
109          }
110          return(MG_OK);
111 + }
112 +
113 +
114 + #define  HTBLSIZ        16381           /* prime hash table size */
115 +
116 + short   hshtab[HTBLSIZ][4];             /* done line segments */
117 +
118 + #define  hash(mx1,my1,mx2,my2)  ((long)(mx1)<<15 ^ (long)(my1)<<10 ^ \
119 +                                        (long)(mx2)<<5 ^ (long)(my2))
120 +
121 +
122 + void
123 + newlayer(void)                          /* start a new layer */
124 + {
125 +        (void)memset((char *)hshtab, '\0', sizeof(hshtab));
126 +        if (++layer >= 16) {
127 +                mendpage();
128 +                layer = 0;
129 +        }
130 + }
131 +
132 +
133 + int
134 + doline(         /* draw line conditionally */
135 +        int     v1x,
136 +        int v1y,
137 +        int v2x,
138 +        int v2y
139 + )
140 + {
141 +        register int    h;
142 +
143 +        if (v1x > v2x || (v1x == v2x && v1y > v2y)) {   /* sort endpoints */
144 +                h=v1x; v1x=v2x; v2x=h;
145 +                h=v1y; v1y=v2y; v2y=h;
146 +        }
147 +        h = hash(v1x, v1y, v2x, v2y) % HTBLSIZ;
148 +        if (hshtab[h][0] == v1x && hshtab[h][1] == v1y &&
149 +                        hshtab[h][2] == v2x && hshtab[h][3] == v2y)
150 +                return(0);
151 +        hshtab[h][0] = v1x; hshtab[h][1] = v1y;
152 +        hshtab[h][2] = v2x; hshtab[h][3] = v2y;
153 +        if ((long)(v2x-v1x)*(v2x-v1x) + (long)(v2y-v1y)*(v2y-v1y)
154 +                        <= irandom(rthresh))
155 +                return(0);
156 +        mline(v1x, v1y, layer/4, 0, layer%4);
157 +        mdraw(v2x, v2y);
158 +        return(1);
159   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines