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

Comparing ray/src/hd/rhdisp2.c (file contents):
Revision 3.3 by gregl, Thu Nov 20 18:06:35 1997 UTC vs.
Revision 3.33 by greg, Sat Feb 22 02:07:24 2003 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1997 Silicon Graphics, Inc. */
2
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6
4   /*
5   * Holodeck beam tracking for display process
6   */
# Line 12 | Line 9 | static char SCCSid[] = "$SunId$ SGI";
9   #include "rhdisp.h"
10   #include "rhdriver.h"
11  
12 < extern int      *getviewcells();
12 > #ifndef MEYERNG
13 > #define MEYERNG         0.2     /* target mean eye range (rel. to grid) */
14 > #endif
15  
17 typedef struct {
18        int     hd;             /* holodeck section number (-1 if inactive) */
19        int     i[3];           /* voxel index (may be outside section) */
20 } VOXL;                 /* a voxel */
21
22 static VOXL voxel[8] = {        /* current voxel list */
23        {-1}, {-1}, {-1}, {-1},
24        {-1}, {-1}, {-1}, {-1}
25 };
26
16   #define CBEAMBLK        1024    /* cbeam allocation block size */
17  
18 < static struct beamcomp {
30 <        short   wants;          /* flags telling which voxels want us */
31 <        short   hd;             /* holodeck section number */
32 <        int     bi;             /* beam index */
33 < } *cbeam = NULL;        /* current beam list */
34 <
18 > static PACKHEAD *cbeam = NULL;  /* current beam list */
19   static int      ncbeams = 0;    /* number of sorted beams in cbeam */
20   static int      xcbeams = 0;    /* extra (unregistered) beams past ncbeams */
21   static int      maxcbeam = 0;   /* size of cbeam array */
22  
23 < struct cellact {
40 <        short   vi;             /* voxel index */
41 <        short   add;            /* zero means delete */
42 < };              /* action for cell */
23 > VIEWPOINT       cureye;         /* current eye position */
24  
44 struct beamact {
45        struct cellact  ca;     /* cell action */
46        GCOORD  gc;             /* grid coordinate */
47 };              /* beam origin and action */
25  
49
26   int
27   newcbeam()              /* allocate new entry at end of cbeam array */
28   {
# Line 55 | Line 31 | newcbeam()             /* allocate new entry at end of cbeam arra
31          if ((i = ncbeams + xcbeams++) >= maxcbeam) {    /* grow array */
32                  maxcbeam += CBEAMBLK;
33                  if (cbeam == NULL)
34 <                        cbeam = (struct beamcomp *)malloc(
35 <                                        maxcbeam*sizeof(struct beamcomp) );
34 >                        cbeam = (PACKHEAD *)malloc(
35 >                                        maxcbeam*sizeof(PACKHEAD) );
36                  else
37 <                        cbeam = (struct beamcomp *)realloc( (char *)cbeam,
38 <                                        maxcbeam*sizeof(struct beamcomp) );
37 >                        cbeam = (PACKHEAD *)realloc( (char *)cbeam,
38 >                                        maxcbeam*sizeof(PACKHEAD) );
39                  if (cbeam == NULL)
40                          error(SYSTEM, "out of memory in newcbeam");
41          }
# Line 69 | Line 45 | newcbeam()             /* allocate new entry at end of cbeam arra
45  
46   int
47   cbeamcmp(cb1, cb2)      /* compare two cbeam entries for sort: keep orphans */
48 < register struct beamcomp        *cb1, *cb2;
48 > register PACKHEAD       *cb1, *cb2;
49   {
50          register int    c;
51  
# Line 81 | Line 57 | register struct beamcomp       *cb1, *cb2;
57  
58   int
59   cbeamcmp2(cb1, cb2)     /* compare two cbeam entries for sort: no orphans */
60 < register struct beamcomp        *cb1, *cb2;
60 > register PACKHEAD       *cb1, *cb2;
61   {
62          register int    c;
63  
64 <        if (!cb1->wants)                /* put orphans at the end, unsorted */
65 <                return(cb2->wants);
64 >        if (!cb1->nr)                   /* put orphans at the end, unsorted */
65 >                return(cb2->nr);
66 >        if (!cb2->nr)
67 >                return(-1);
68          if ((c = cb1->bi - cb2->bi))    /* sort on beam index first */
69                  return(c);
70          return(cb1->hd - cb2->hd);      /* use hd to resolve matches */
# Line 97 | Line 75 | int
75   findcbeam(hd, bi)       /* find the specified beam in our sorted list */
76   int     hd, bi;
77   {
78 <        struct beamcomp cb;
79 <        register struct beamcomp        *p;
78 >        PACKHEAD        cb;
79 >        register PACKHEAD       *p;
80  
81          if (ncbeams <= 0)
82                  return(-1);
83 <        cb.wants = 0; cb.hd = hd; cb.bi = bi;
84 <        p = (struct beamcomp *)bsearch((char *)&cb, (char *)cbeam, ncbeams,
85 <                        sizeof(struct beamcomp), cbeamcmp);
83 >        cb.hd = hd; cb.bi = bi; cb.nr = cb.nc = 0;
84 >        p = (PACKHEAD *)bsearch((char *)&cb, (char *)cbeam, ncbeams,
85 >                        sizeof(PACKHEAD), cbeamcmp);
86          if (p == NULL)
87                  return(-1);
88          return(p - cbeam);
# Line 130 | Line 108 | int    bi;
108          if (bi < 1 | bi > nbeams(hdlist[hd]))
109                  error(INTERNAL, "illegal beam index in getcbeam");
110          n = newcbeam();         /* allocate and assign */
111 <        cbeam[n].wants = 0; cbeam[n].hd = hd; cbeam[n].bi = bi;
111 >        cbeam[n].hd = hd; cbeam[n].bi = bi; cbeam[n].nr = cbeam[n].nc = 0;
112          return(n);
113   }
114  
# Line 143 | Line 121 | int    adopt;
121          if (!(ncbeams += xcbeams))
122                  return;
123          xcbeams = 0;
124 <        qsort((char *)cbeam, ncbeams, sizeof(struct beamcomp),
124 >        qsort((char *)cbeam, ncbeams, sizeof(PACKHEAD),
125                          adopt ? cbeamcmp : cbeamcmp2);
126          if (adopt)
127                  return;
128          for (i = ncbeams; i--; )        /* identify orphans */
129 <                if (cbeam[i].wants)
129 >                if (cbeam[i].nr)
130                          break;
131          xcbeams = ncbeams - ++i;        /* put orphans after ncbeams */
132          ncbeams = i;
133   }
134  
135  
136 < cbeamadd(bl, n, v, hr, vr)      /* add beams to server list */
137 < register struct beamcomp        *bl;
160 < int     n;
161 < VIEW    *v;
162 < int     hr, vr;
136 > beam_init(fresh)                /* clear beam list for new view(s) */
137 > int     fresh;
138   {
164        register PACKHEAD       *pa;
139          register int    i;
140  
141 <        if (n <= 0)
142 <                return;
143 <        pa = (PACKHEAD *)malloc(n*sizeof(PACKHEAD));
144 <        if (pa == NULL)
145 <                error(SYSTEM, "out of memory in cbeamadd");
146 <        for (i = 0; i < n; i++)
173 <                pa[i].nr = npixels(v, hr, vr,
174 <                                hdlist[pa[i].hd=bl[i].hd],
175 <                                pa[i].bi=bl[i].bi) / 8;
176 <
177 <        serv_request(DR_ADDSET, n*sizeof(PACKHEAD), (char *)pa);
178 <        free((char *)pa);
141 >        if (fresh)                      /* discard old beams? */
142 >                ncbeams = xcbeams = 0;
143 >        else                            /* else clear sample requests */
144 >                for (i = ncbeams+xcbeams; i--; )
145 >                        cbeam[i].nr = 0;
146 >        cureye.rng = 0.;
147   }
148  
149  
150 < cbeamdel(bl, n)         /* delete unwanted beam requests */
151 < register struct beamcomp        *bl;
152 < int     n;
150 > int2 *
151 > beam_view(vn, hr, vr)           /* add beam view (if advisable) */
152 > VIEW    *vn;
153 > int     hr, vr;
154   {
155 <        register PACKHEAD       *pa;
156 <        register int    i;
157 <
158 <        if (n <= 0)
159 <                return;
160 <        pa = (PACKHEAD *)malloc(n*sizeof(PACKHEAD));
161 <        if (pa == NULL)
162 <                error(SYSTEM, "out of memory in cbeamdel");
163 <        for (i = 0; i < n; i++) {
164 <                pa[i].hd = bl[i].hd;
196 <                pa[i].bi = bl[i].bi;
197 <                pa[i].nr = 0;           /* removes any request */
155 >        int2    *slist;
156 >        BEAMLIST        blist;
157 >        double  eravg, d;
158 >        register HOLO   *hp;
159 >        register int    i, n;
160 >                                        /* compute beams for view */
161 >        slist = viewbeams(vn, hr, vr, &blist);
162 >        if (*slist < 0) {
163 >                error(COMMAND, "no sections visible from this view");
164 >                return(NULL);
165          }
166 <        serv_request(DR_DELSET, n*sizeof(PACKHEAD), (char *)pa);
167 <        free((char *)pa);
168 < }
169 <
170 <
171 < int
172 < set_voxels(vl, n)       /* set new voxel array */
206 < VOXL    vl[8];
207 < int     n;
208 < {
209 <        short   wmap[256];
210 <        int     vmap[8];
211 <        int     comn = 0;
212 <        int     no;
213 <        register int    i, j;
214 <                                        /* find common voxels */
215 <        for (j = 0; j < 8 && voxel[j].hd >= 0; j++) {
216 <                vmap[j] = -1;
217 <                for (i = n; i--; )
218 <                        if (!bcmp((char *)(vl+i), (char *)(voxel+j),
219 <                                        sizeof(VOXL))) {
220 <                                vmap[j] = i;
221 <                                comn |= 1<<i;
222 <                                break;
223 <                        }
166 >                                        /* sort current beam list */
167 >        cbeamsort(1);
168 >                                        /* add new beams to list */
169 >        for (i = blist.nb; i--; ) {
170 >                n = getcbeam(blist.bl[i].hd, blist.bl[i].bi);
171 >                if (blist.bl[i].nr > cbeam[n].nr)
172 >                        cbeam[n].nr = blist.bl[i].nr;
173          }
174 <        no = comn ? j : 0;              /* compute flag mapping */
175 <        for (i = 256; i--; ) {
176 <                wmap[i] = 0;
177 <                for (j = no; j--; )
178 <                        if (vmap[j] >= 0 && i & 1<<j)
179 <                                wmap[i] |= 1<<vmap[j];
174 >        free((void *)blist.bl);         /* free list */
175 >        if (MEYERNG <= FTINY)
176 >                return(slist);
177 >                                        /* compute average eye range */
178 >        eravg = 0.;
179 >        for (n = 0; slist[n] >= 0; n++) {
180 >                hp = hdlist[slist[n]];
181 >                eravg +=        MEYERNG/3. / VLEN(hp->wg[0]) +
182 >                                MEYERNG/3. / VLEN(hp->wg[1]) +
183 >                                MEYERNG/3. / VLEN(hp->wg[2]) ;
184          }
185 <                                        /* fix cbeam flags */
186 <        for (i = ncbeams; i--; )
187 <                cbeam[i].wants = wmap[cbeam[i].wants];
188 <                                        /* update our voxel list */
189 <        bcopy((char *)vl, (char *)voxel, n*sizeof(VOXL));
190 <        for (j = n; j < 8; j++)
191 <                voxel[j].hd = -1;
192 <        return(comn);                   /* return bit array of common voxels */
193 < }
241 <
242 <
243 < int
244 < get_voxels(vl, vp)      /* find voxels corresponding to view point */
245 < VOXL    vl[8];
246 < FVECT   vp;
247 < {
248 <        static int      lastn = 0, lastd = -1;
249 <        int     n = 0;
250 <        FVECT   gp;
251 <        double  d;
252 <        int     dist, bestd = 0x7fff;
253 <        VOXL    vox;
254 <        register int    i, j, k;
255 <                                        /* find closest voxels */
256 <        for (i = 0; n < 8 && hdlist[i]; i++) {
257 <                hdgrid(gp, hdlist[i], vp);
258 <                for (j = 0; n < 8 && j < 8; j++) {
259 <                        dist = 0;
260 <                        for (k = 0; k < 3; k++) {
261 <                                d = gp[k] - .5 + (j>>k & 1);
262 <                                if (d < 0.)
263 <                                        dist += -(vox.i[k] = (int)d - 1);
264 <                                else if ((vox.i[k] = d) >= hdlist[i]->grid[k])
265 <                                        dist += vox.i[k] -
266 <                                                        hdlist[i]->grid[k] + 1;
267 <                        }
268 <                        if (dist > bestd)       /* others were closer */
269 <                                continue;
270 <                        if (dist < bestd) {     /* start closer list */
271 <                                n = 0;
272 <                                bestd = dist;
273 <                        }
274 <                        vox.hd = i;             /* add this voxel */
275 <                        copystruct(&vl[n], &vox);
276 <                        n++;
277 <                }
185 >        eravg /= (double)n;
186 >                                        /* add to current eye position */
187 >        if (cureye.rng <= FTINY) {
188 >                VCOPY(cureye.vpt, vn->vp);
189 >                cureye.rng = eravg;
190 >        } else if ((d = sqrt(dist2(vn->vp,cureye.vpt))) + eravg > cureye.rng) {
191 >                for (i = 3; i--; )
192 >                        cureye.vpt[i] = 0.5*(cureye.vpt[i] + vn->vp[i]);
193 >                cureye.rng = 0.5*(cureye.rng + eravg + d);
194          }
195 <                                        /* warn of dangerous moves */
280 <        if (n < lastn && bestd >= lastd)
281 <                error(WARNING, "moving outside holodeck section");
282 <        else if (n > lastn && bestd <= lastd)
283 <                error(WARNING, "moving inside holodeck section");
284 <        lastd = bestd;
285 <        return(lastn = n);
195 >        return(slist);
196   }
197  
198  
199   int
200 < dobeam(gcp, bp)         /* express interest or disintrest in a beam */
201 < GCOORD  *gcp;
292 < register struct beamact *bp;
200 > beam_sync(all)                  /* update beam list on server */
201 > int     all;
202   {
203 <        GCOORD  gc[2];
204 <        int     bi, i;
205 <                                        /* compute beam index */
206 <        copystruct(gc, gcp);
207 <        copystruct(gc+1, &bp->gc);
208 <        if ((bi = hdbindex(hdlist[voxel[bp->ca.vi].hd], gc)) <= 0)
209 <                error(CONSISTENCY, "bad grid coordinate in dobeam");
210 <        if (bp->ca.add) {               /* add it in */
211 <                i = getcbeam(voxel[bp->ca.vi].hd, bi);
303 <                cbeam[i].wants |= 1<<bp->ca.vi; /* say we want it */
304 <                return(i == ncbeams+xcbeams-1); /* return 1 if totally new */
305 <        }
306 <                                        /* else delete it */
307 <        i = findcbeam(voxel[bp->ca.vi].hd, bi);
308 <        if (i >= 0 && cbeam[i].wants & 1<<bp->ca.vi) {
309 <                cbeam[i].wants &= ~(1<<bp->ca.vi);      /* we don't want it */
310 <                return(1);                      /* indicate change */
311 <        }
312 <        return(0);
313 < }
314 <
315 <
316 <
317 < int
318 < docell(gcp, cap)        /* find beams corresponding to cell and voxel */
319 < GCOORD  *gcp;
320 < register struct cellact *cap;
321 < {
322 <        FVECT   org, dir[4];
323 <        FVECT   gp, cv[4], vc;
324 <        struct beamact  bo;
325 <        register int    i;
326 <                                /* compute cell vertices */
327 <        hdcell(cv, hdlist[voxel[cap->vi].hd], gcp);
328 <                                /* compute cell and voxel centers */
329 <        for (i = 0; i < 3; i++) {
330 <                org[i] = 0.5*(cv[0][i] + cv[2][i]);
331 <                gp[i] = voxel[cap->vi].i[i] + 0.5;
332 <        }
333 <        hdworld(vc, hdlist[voxel[cap->vi].hd], gp);
334 <                                /* compute voxel pyramid using vector trick */
335 <        for (i = 0; i < 3; i++) {
336 <                dir[0][i] = vc[i] - cv[0][i];           /* to 3 */
337 <                dir[2][i] = vc[i] - cv[3][i];           /* to 0 */
338 <                if (gcp->w & 1) {       /* watch vertex order! */
339 <                        dir[1][i] = vc[i] - cv[2][i];   /* to 1 */
340 <                        dir[3][i] = vc[i] - cv[1][i];   /* to 2 */
341 <                } else {
342 <                        dir[1][i] = vc[i] - cv[1][i];   /* to 2 */
343 <                        dir[3][i] = vc[i] - cv[2][i];   /* to 1 */
344 <                }
345 <        }
346 <                                /* visit beams for opposite cells */
347 <        copystruct(&bo.ca, cap);
348 <        copystruct(&bo.gc, gcp);
349 <        return(visit_cells(org, dir, hdlist[voxel[cap->vi].hd], dobeam, &bo));
350 < }
351 <
352 <
353 < int
354 < doview(cap, vp)                 /* visit cells for a given view */
355 < struct cellact  *cap;
356 < VIEW    *vp;
357 < {
358 <        FVECT   org, dir[4];
359 <                                        /* compute view pyramid */
360 <        if (vp->type == VT_PAR) goto viewerr;
361 <        if (viewray(org, dir[0], vp, 0., 0.) < -FTINY) goto viewerr;
362 <        if (viewray(org, dir[1], vp, 0., 1.) < -FTINY) goto viewerr;
363 <        if (viewray(org, dir[2], vp, 1., 1.) < -FTINY) goto viewerr;
364 <        if (viewray(org, dir[3], vp, 1., 0.) < -FTINY) goto viewerr;
365 <                                        /* visit cells within pyramid */
366 <        return(visit_cells(org, dir, hdlist[voxel[cap->vi].hd], docell, cap));
367 < viewerr:
368 <        error(INTERNAL, "unusable view in doview");
369 < }
370 <
371 <
372 < mvview(voxi, vold, vnew)        /* move view for a voxel */
373 < int     voxi;
374 < VIEW    *vold, *vnew;
375 < {
376 <        int     netchange = 0;
377 <        int     *ocl, *ncl;
378 <        struct cellact  ca;
379 <        int     ocnt, ncnt;
380 <        int     c;
381 <        register GCOORD *ogcp, *ngcp;
382 <                                /* get old and new cell lists */
383 <        ocl = getviewcells(hdlist[voxel[voxi].hd], vold);
384 <        ncl = getviewcells(hdlist[voxel[voxi].hd], vnew);
385 <        if (ocl != NULL) {
386 <                ocnt = *ocl; ogcp = (GCOORD *)(ocl+1);
203 >                                        /* set new eye position */
204 >        serv_request(DR_VIEWPOINT, sizeof(VIEWPOINT), (char *)&cureye);
205 >                                        /* sort list (put orphans at end) */
206 >        cbeamsort(all < 0);
207 >                                        /* send beam request */
208 >        if (all) {
209 >                if (ncbeams > 0)
210 >                        serv_request(DR_NEWSET,
211 >                                        ncbeams*sizeof(PACKHEAD), cbeam);
212          } else {
213 <                ocnt = 0; ogcp = NULL;
213 >                if (ncbeams+xcbeams > 0)
214 >                        serv_request(DR_ADJSET,
215 >                                (ncbeams+xcbeams)*sizeof(PACKHEAD), cbeam);
216          }
217 <        if (ncl != NULL) {
218 <                ncnt = *ncl; ngcp = (GCOORD *)(ncl+1);
392 <        } else {
393 <                ncnt = 0; ngcp = NULL;
394 <        }
395 <        ca.vi = voxi;           /* add and delete cells */
396 <        while (ocnt > 0 & ncnt > 0)
397 <                if ((c = cellcmp(ogcp, ngcp)) > 0) {
398 <                        ca.add = 1;             /* new cell */
399 <                        netchange += docell(ngcp++, &ca);
400 <                        ncnt--;
401 <                } else if (c < 0) {
402 <                        ca.add = 0;             /* obsolete cell */
403 <                        netchange -= docell(ogcp++, &ca);
404 <                        ocnt--;
405 <                } else {
406 <                        ogcp++; ocnt--;         /* unchanged cell */
407 <                        ngcp++; ncnt--;
408 <                }
409 <                                /* take care of list tails */
410 <        for (ca.add = 1; ncnt > 0; ncnt--)
411 <                netchange += docell(ngcp++, &ca);
412 <        for (ca.add = 0; ocnt > 0; ocnt--)
413 <                netchange -= docell(ogcp++, &ca);
414 <                                /* clean up */
415 <        if (ocl != NULL) free((char *)ocl);
416 <        if (ncl != NULL) free((char *)ncl);
417 <        return(netchange);
217 >        xcbeams = 0;                    /* truncate our list */
218 >        return(ncbeams);
219   }
220  
221  
222 < beam_view(vo, vn)       /* change beam view */
223 < VIEW    *vo, *vn;
222 > gridlines(f)                    /* run through holodeck section grid lines */
223 > int     (*f)();
224   {
225 <        struct cellact  ca;
226 <        VOXL    vlnew[8];
227 <        int     n, comn;
228 <
229 <        if (!vn->type) {                /* clear our beam list */
230 <                set_voxels(vlnew, 0);
231 <                cbeamdel(cbeam, ncbeams);
232 <                ncbeams = 0;
233 <                return;
234 <        }
235 <                                        /* find our new voxels */
236 <        n = get_voxels(vlnew, vn->vp);
237 <                                        /* set the new voxels */
238 <        comn = set_voxels(vlnew, n);
239 <        if (!vo->type)
240 <                comn = 0;
241 <        ca.add = 1;                     /* update our beam list */
242 <        for (ca.vi = n; ca.vi--; )
243 <                if (comn & 1<<ca.vi)    /* change which cells we see */
244 <                        mvview(ca.vi, vo, vn);
245 <                else                    /* else add all new cells */
246 <                        doview(&ca, vn);
247 <                                        /* inform server of new beams */
248 <        cbeamadd(cbeam+ncbeams, xcbeams, vn, odev.hres, odev.vres);
249 <                                        /* sort list to put orphans at end */
250 <        cbeamsort(0);
251 <                                        /* tell server to delete orphans */
252 <        cbeamdel(cbeam+ncbeams, xcbeams);
253 <        xcbeams = 0;                    /* truncate our list */
225 >        register int    hd, w, i;
226 >        int     g0, g1;
227 >        FVECT   wp[2], mov;
228 >        double  d;
229 >                                        /* do each wall on each section */
230 >        for (hd = 0; hdlist[hd] != NULL; hd++)
231 >                for (w = 0; w < 6; w++) {
232 >                        g0 = hdwg0[w];
233 >                        g1 = hdwg1[w];
234 >                        d = 1.0/hdlist[hd]->grid[g0];
235 >                        mov[0] = d * hdlist[hd]->xv[g0][0];
236 >                        mov[1] = d * hdlist[hd]->xv[g0][1];
237 >                        mov[2] = d * hdlist[hd]->xv[g0][2];
238 >                        if (w & 1) {
239 >                                VSUM(wp[0], hdlist[hd]->orig,
240 >                                                hdlist[hd]->xv[w>>1], 1.);
241 >                                VSUM(wp[0], wp[0], mov, 1.);
242 >                        } else
243 >                                VCOPY(wp[0], hdlist[hd]->orig);
244 >                        VSUM(wp[1], wp[0], hdlist[hd]->xv[g1], 1.);
245 >                        for (i = hdlist[hd]->grid[g0]; ; ) {    /* g0 lines */
246 >                                (*f)(wp);
247 >                                if (!--i) break;
248 >                                wp[0][0] += mov[0]; wp[0][1] += mov[1];
249 >                                wp[0][2] += mov[2]; wp[1][0] += mov[0];
250 >                                wp[1][1] += mov[1]; wp[1][2] += mov[2];
251 >                        }
252 >                        d = 1.0/hdlist[hd]->grid[g1];
253 >                        mov[0] = d * hdlist[hd]->xv[g1][0];
254 >                        mov[1] = d * hdlist[hd]->xv[g1][1];
255 >                        mov[2] = d * hdlist[hd]->xv[g1][2];
256 >                        if (w & 1)
257 >                                VSUM(wp[0], hdlist[hd]->orig,
258 >                                                hdlist[hd]->xv[w>>1], 1.);
259 >                        else
260 >                                VSUM(wp[0], hdlist[hd]->orig, mov, 1.);
261 >                        VSUM(wp[1], wp[0], hdlist[hd]->xv[g0], 1.);
262 >                        for (i = hdlist[hd]->grid[g1]; ; ) {    /* g1 lines */
263 >                                (*f)(wp);
264 >                                if (!--i) break;
265 >                                wp[0][0] += mov[0]; wp[0][1] += mov[1];
266 >                                wp[0][2] += mov[2]; wp[1][0] += mov[0];
267 >                                wp[1][1] += mov[1]; wp[1][2] += mov[2];
268 >                        }
269 >                }
270   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines