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

Comparing ray/src/rt/rc3.c (file contents):
Revision 2.3 by greg, Sun Jun 10 05:25:42 2012 UTC vs.
Revision 2.17 by greg, Thu Nov 15 15:26:52 2012 UTC

# Line 6 | Line 6 | static const char RCSid[] = "$Id$";
6   * Controlling process for multiple children
7   */
8  
9 + #include <signal.h>
10   #include "rcontrib.h"
10 #include "platform.h"
11 #include "rtprocess.h"
11   #include "selcall.h"
12  
13 + #define MAXIQ           (int)(PIPE_BUF/(sizeof(FVECT)*2))
14 +
15   /* Modifier contribution queue (results waiting to be output) */
16   typedef struct s_binq {
17 <        int             ndx;            /* index for this entry */
18 <        int             n2add;          /* number left to add */
17 >        RNUMBER         ndx;            /* index for this entry */
18 >        RNUMBER         nadded;         /* accumulated so far */
19          struct s_binq   *next;          /* next in queue */
20          MODCONT         *mca[1];        /* contrib. array (extends struct) */
21   } BINQ;
# Line 22 | Line 23 | typedef struct s_binq {
23   static BINQ     *out_bq = NULL;         /* output bin queue */
24   static BINQ     *free_bq = NULL;        /* free queue entries */
25  
26 < static SUBPROC  kida[MAXPROCESS];       /* child processes */
27 < static FILE     *inq_fp[MAXPROCESS];    /* input streams */
26 > static struct {
27 >        RNUMBER r1;                     /* assigned ray starting index */
28 >        SUBPROC pr;                     /* PID, i/o descriptors */
29 >        FILE    *infp;                  /* file pointer to read from process */
30 >        int     nr;                     /* number of rays to sum (0 if free) */
31 > } kida[MAXPROCESS];             /* our child processes */
32  
33  
34 < /* Get new (empty) bin queue entry */
34 > /* Get new bin queue entry */
35   static BINQ *
36   new_binq()
37   {
38 <        BINQ    *bp = free_bq;
38 >        BINQ    *bp;
39          int     i;
40  
41 <        if (bp != NULL) {               /* something already available? */
41 >        if (free_bq != NULL) {          /* something already available? */
42 >                bp = free_bq;
43                  free_bq = bp->next;
44                  bp->next = NULL;
45 <                bp->n2add = accumulate-1;
45 >                bp->nadded = 0;
46                  return(bp);
47          }
48                                          /* else allocate fresh */
49 <        bp = (BINQ *)malloc(sizeof(BINQ)+(nmods-1)*sizeof(MODCONT *));
49 >        bp = (BINQ *)malloc(sizeof(BINQ) + sizeof(MODCONT *)*(nmods-1));
50          if (bp == NULL)
51                  goto memerr;
52          for (i = nmods; i--; ) {
# Line 53 | Line 59 | new_binq()
59                  /* memset(bp->mca[i]->cbin, 0, sizeof(DCOLOR)*mp->nbins); */
60          }
61          bp->ndx = 0;
62 <        bp->n2add = accumulate-1;
62 >        bp->nadded = 0;
63          bp->next = NULL;
64          return(bp);
65   memerr:
# Line 89 | Line 95 | free_binq(BINQ *bp)
95  
96  
97   /* Add modifier values to accumulation record in queue and clear */
98 < void
98 > static void
99   queue_modifiers()
100   {
101          MODCONT *mpin, *mpout;
# Line 105 | Line 111 | queue_modifiers()
111                          addcolor(mpout->cbin[j], mpin->cbin[j]);
112                  memset(mpin->cbin, 0, sizeof(DCOLOR)*mpin->nbins);
113          }
114 +        out_bq->nadded++;
115   }
116  
117  
118 < /* Sum one modifier record into another (doesn't update n2add) */
118 > /* Sum one modifier record into another (updates nadded) */
119   static void
120   add_modbin(BINQ *dst, BINQ *src)
121   {
# Line 120 | Line 127 | add_modbin(BINQ *dst, BINQ *src)
127                  for (j = mpout->nbins; j--; )
128                          addcolor(mpout->cbin[j], mpin->cbin[j]);
129          }
130 +        dst->nadded += src->nadded;
131   }
132  
133  
# Line 139 | Line 147 | queue_output(BINQ *bp)
147                  }
148                  return;
149          }
150 <        b_last = NULL;                  /* else insert in output queue */
150 >        b_last = NULL;                  /* insert in output queue */
151          for (b_cur = out_bq; b_cur != NULL && b_cur->ndx < bp->ndx;
152                                  b_cur = b_cur->next)
153                  b_last = b_cur;
154 +
155          if (b_last != NULL) {
156                  bp->next = b_cur;
157                  b_last->next = bp;
# Line 154 | Line 163 | queue_output(BINQ *bp)
163                  return;
164          b_cur = out_bq;                 /* else merge accumulation entries */
165          while (b_cur->next != NULL) {
166 <                if (b_cur->n2add <= 0 ||
166 >                if (b_cur->nadded >= accumulate ||
167                                  (b_cur->ndx-1)/accumulate !=
168                                  (b_cur->next->ndx-1)/accumulate) {
169                          b_cur = b_cur->next;
170                          continue;
171                  }
172                  add_modbin(b_cur, b_cur->next);
164                b_cur->n2add--;
173                  b_last = b_cur->next;
174                  b_cur->next = b_last->next;
175                  b_last->next = NULL;
# Line 170 | Line 178 | queue_output(BINQ *bp)
178   }
179  
180  
181 < /* Get current with output queue by producing ready results */
181 > /* Count number of records ready for output */
182   static int
183 < output_catchup()
183 > queue_ready()
184   {
185 +        int     nready = 0;
186 +        BINQ    *bp;
187 +
188 +        for (bp = out_bq; bp != NULL && bp->nadded >= accumulate &&
189 +                                bp->ndx == lastdone+nready*accumulate+1;
190 +                                bp = bp->next)
191 +                ++nready;
192 +
193 +        return(nready);
194 + }
195 +
196 +
197 + /* Catch up with output queue by producing ready results */
198 + static int
199 + output_catchup(int nmax)
200 + {
201          int     nout = 0;
202          BINQ    *bp;
203          int     i;
204 <
205 <        if (accumulate <= 0)            /* just accumulating? */
206 <                return(0);
207 <                                        /* else output ready results */
208 <        while (out_bq != NULL && (out_bq->ndx == lastdone+1) & !out_bq->n2add) {
204 >                                        /* output ready results */
205 >        while (out_bq != NULL && out_bq->nadded >= accumulate
206 >                                && out_bq->ndx == lastdone+1) {
207 >                if ((nmax > 0) & (nout >= nmax))
208 >                        break;
209                  bp = out_bq;                    /* pop off first entry */
210                  out_bq = bp->next;
211                  bp->next = NULL;
# Line 206 | Line 230 | put_zero_record(int ndx)
230          for (i = nmods; i--; )
231                  memset(bp->mca[i]->cbin, 0, sizeof(DCOLOR)*bp->mca[i]->nbins);
232          bp->ndx = ndx;
233 +        bp->nadded = 1;
234          queue_output(bp);
235 <        output_catchup();
235 >        output_catchup(0);
236   }
237  
238  
239 + /* Get results from child process and add to queue */
240 + static void
241 + queue_results(int k)
242 + {
243 +        BINQ    *bq = new_binq();       /* get results holder */
244 +        int     j;
245 +
246 +        bq->ndx = kida[k].r1;
247 +        bq->nadded = kida[k].nr;
248 +                                        /* read from child */
249 +        for (j = 0; j < nmods; j++)
250 +                if (fread(bq->mca[j]->cbin, sizeof(DCOLOR), bq->mca[j]->nbins,
251 +                                        kida[k].infp) != bq->mca[j]->nbins)
252 +                        error(SYSTEM, "read error from render process");
253 +                        
254 +        queue_output(bq);               /* put results in output queue */
255 +        kida[k].nr = 0;                 /* mark child as available */
256 + }
257 +
258 +
259   /* callback to set output spec to NULL (stdout) */
260   static int
261   set_stdout(const LUENT *le, void *p)
# Line 224 | Line 269 | set_stdout(const LUENT *le, void *p)
269   int
270   in_rchild()
271   {
272 < #ifdef _WIN32
273 <        error(WARNING, "multiprocessing unsupported -- running solo");
274 <        nproc = 1;
230 <        return(1);
231 < #else
232 <                                        /* try to fork ourselves */
233 <        while (nchild < nproc) {
234 <                int     p0[2], p1[2];
235 <                int     pid;
236 <                                        /* prepare i/o pipes */
272 >        int     rval;
273 >
274 >        while (nchild < nproc) {        /* fork until target reached */
275                  errno = 0;
276 <                if (pipe(p0) < 0 || pipe(p1) < 0)
277 <                        error(SYSTEM, "pipe() call failed!");
278 <                pid = fork();           /* fork parent process */
279 <                if (pid == 0) {         /* if in child, set up & return true */
280 <                        close(p0[1]); close(p1[0]);
281 <                        dup2(p0[0], 0); close(p0[0]);
282 <                        dup2(p1[1], 1); close(p1[1]);
276 >                rval = open_process(&kida[nchild].pr, NULL);
277 >                if (rval < 0)
278 >                        error(SYSTEM, "open_process() call failed");
279 >                if (rval == 0) {        /* if in child, set up & return true */
280 >                        lu_doall(&modconttab, set_stdout, NULL);
281 >                        lu_done(&ofiletab);
282 >                        while (nchild--) {      /* don't share other pipes */
283 >                                close(kida[nchild].pr.w);
284 >                                fclose(kida[nchild].infp);
285 >                        }
286                          inpfmt = (sizeof(RREAL)==sizeof(double)) ? 'd' : 'f';
287                          outfmt = 'd';
288                          header = 0;
248                        waitflush = xres = 1;
289                          yres = 0;
290                          raysleft = 0;
291 <                        account = accumulate = 1;
292 <                        lu_doall(&modconttab, set_stdout, NULL);
293 <                        nchild = -1;
294 <                        return(1);      /* child return value */
291 >                        if (accumulate == 1) {
292 >                                waitflush = xres = 1;
293 >                                account = accumulate = 1;
294 >                        } else {        /* parent controls accumulation */
295 >                                waitflush = xres = 0;
296 >                                account = accumulate = 0;
297 >                        }
298 >                        return(1);      /* return "true" in child */
299                  }
300 <                if (pid < 0)
301 <                        error(SYSTEM, "fork() call failed!");
302 <                                        /* connect parent's pipes */
303 <                close(p0[0]); close(p1[1]);
304 <                kida[nchild].r = p1[0];
261 <                kida[nchild].w = p0[1];
262 <                kida[nchild].pid = pid;
263 <                kida[nchild].running = -1;
264 <                inq_fp[nchild] = fdopen(p1[0], "rb");
265 <                if (inq_fp[nchild] == NULL)
300 >                if (rval != PIPE_BUF)
301 >                        error(CONSISTENCY, "bad value from open_process()");
302 >                                        /* connect to child's output */
303 >                kida[nchild].infp = fdopen(kida[nchild].pr.r, "rb");
304 >                if (kida[nchild].infp == NULL)
305                          error(SYSTEM, "out of memory in in_rchild()");
306   #ifdef getc_unlocked
307 <                flockfile(inq_fp[nchild]);      /* avoid mutex overhead */
307 >                flockfile(kida[nchild].infp);   /* avoid mutex overhead */
308   #endif
309 <                ++nchild;
309 >                kida[nchild++].nr = 0;  /* mark as available */
310          }
311 <        return(0);                      /* parent return value */
273 < #endif
311 >        return(0);                      /* return "false" in parent */
312   }
313  
314  
315   /* Close child processes */
316   void
317 < end_children()
317 > end_children(int immed)
318   {
319          int     status;
320          
321 <        while (nchild-- > 0) {
322 <                kida[nchild].r = -1;    /* close(-1) error is ignored */
323 <                if ((status = close_process(&kida[nchild])) > 0) {
321 >        while (nchild > 0) {
322 >                nchild--;
323 > #ifdef SIGKILL
324 >                if (immed)              /* error mode -- quick exit */
325 >                        kill(kida[nchild].pr.pid, SIGKILL);
326 > #endif
327 >                if ((status = close_process(&kida[nchild].pr)) > 0 && !immed) {
328                          sprintf(errmsg,
329                                  "rendering process returned bad status (%d)",
330                                          status);
331                          error(WARNING, errmsg);
332                  }
333 <                fclose(inq_fp[nchild]); /* performs actual close() */
333 >                fclose(kida[nchild].infp);
334          }
335   }
336  
337  
338 < /* Wait for the next available child, managing output queue as well */
338 > /* Wait for the next available child, managing output queue simultaneously */
339   static int
340 < next_child_nq(int force_wait)
340 > next_child_nq(int flushing)
341   {
342          static struct timeval   polling;
343 <        struct timeval          *pmode = force_wait | (accumulate <= 0) ?
302 <                                        (struct timeval *)NULL : &polling;
343 >        struct timeval          *pmode;
344          fd_set                  readset, errset;
345 <        int                     i, j, n, nr;
345 >        int                     i, n, nr, nqr;
346  
347 <        if (!force_wait)                /* see if there's one free */
347 >        if (!flushing)                  /* see if there's one free */
348                  for (i = nchild; i--; )
349 <                        if (kida[i].running < 0)
349 >                        if (!kida[i].nr)
350                                  return(i);
351 +
352 +        nqr = queue_ready();            /* choose blocking mode or polling */
353 +        if ((nqr > 0) & !flushing)
354 +                pmode = &polling;
355 +        else
356 +                pmode = NULL;
357 + tryagain:                               /* catch up with output? */
358 +        if (pmode == &polling) {
359 +                if (nqr > nchild)       /* don't get too far behind */
360 +                        nqr -= output_catchup(nqr-nchild);
361 +        } else if (nqr > 0)             /* clear output before blocking */
362 +                nqr -= output_catchup(0);
363                                          /* prepare select() call */
364          FD_ZERO(&readset); FD_ZERO(&errset);
365          n = nr = 0;
366          for (i = nchild; i--; ) {
367 <                if (kida[i].running > 0) {
368 <                        FD_SET(kida[i].r, &readset);
367 >                if (kida[i].nr) {
368 >                        FD_SET(kida[i].pr.r, &readset);
369                          ++nr;
370                  }
371 <                FD_SET(kida[i].r, &errset);
372 <                if (kida[i].r >= n)
373 <                        n = kida[i].r + 1;
371 >                FD_SET(kida[i].pr.r, &errset);
372 >                if (kida[i].pr.r >= n)
373 >                        n = kida[i].pr.r + 1;
374          }
322 tryagain:
323        if (pmode == NULL)              /* catch up in case we block */
324                output_catchup();
375          if (!nr)                        /* nothing to wait for? */
376                  return(-1);
377          if ((nr > 1) | (pmode == &polling)) {
378                  errno = 0;
379 <                i = select(n, &readset, NULL, &errset, pmode);
380 <                if (!i) {
379 >                nr = select(n, &readset, NULL, &errset, pmode);
380 >                if (!nr) {
381                          pmode = NULL;   /* try again, blocking this time */
382                          goto tryagain;
383                  }
384 <                if (i < 0)
384 >                if (nr < 0)
385                          error(SYSTEM, "select() error in next_child_nq()");
386          } else
387                  FD_ZERO(&errset);
388          n = -1;                         /* read results from child(ren) */
389          for (i = nchild; i--; ) {
390 <                BINQ    *bq;
341 <                if (FD_ISSET(kida[i].r, &errset))
390 >                if (FD_ISSET(kida[i].pr.r, &errset))
391                          error(USER, "rendering process died");
392 <                if (!FD_ISSET(kida[i].r, &readset))
393 <                        continue;
345 <                bq = new_binq();        /* get results holder */
346 <                bq->ndx = kida[i].running;
347 <                                        /* read from child */
348 <                for (j = 0; j < nmods; j++) {
349 <                        nr = bq->mca[j]->nbins;
350 <                        if (fread(bq->mca[j]->cbin, sizeof(DCOLOR), nr,
351 <                                                        inq_fp[i]) != nr)
352 <                                error(SYSTEM, "read error from render process");
353 <                }
354 <                queue_output(bq);       /* add results to output queue */
355 <                kida[i].running = -1;   /* mark child as available */
356 <                n = i;
392 >                if (FD_ISSET(kida[i].pr.r, &readset))
393 >                        queue_results(n = i);
394          }
395          return(n);                      /* first available child */
396   }
# Line 363 | Line 400 | tryagain:
400   void
401   parental_loop()
402   {
403 <        static int      ignore_warning_given = 0;
404 <        FVECT           orgdir[2];
405 <        double          d;
406 <        int             i;
403 >        const int       qlimit = (accumulate == 1) ? 1 : MAXIQ-1;
404 >        int             ninq = 0;
405 >        FVECT           orgdir[2*MAXIQ];
406 >        int             i, n;
407                                          /* load rays from stdin & process */
408   #ifdef getc_unlocked
409          flockfile(stdin);               /* avoid lock/unlock overhead */
410   #endif
411 <        while (getvec(orgdir[0]) == 0 && getvec(orgdir[1]) == 0) {
412 <                d = normalize(orgdir[1]);
413 <                                                        /* asking for flush? */
414 <                if ((d == 0.0) & (accumulate != 1)) {
415 <                        if (!ignore_warning_given++)
416 <                                error(WARNING,
417 <                                "dummy ray(s) ignored during accumulation\n");
418 <                        continue;
419 <                }
420 <                if ((d == 0.0) | (lastray+1 < lastray)) {
421 <                        while (next_child_nq(1) >= 0)
422 <                                ;                       /* clear the queue */
423 <                        lastdone = lastray = 0;
424 <                }
425 <                if (d == 0.0) {
426 <                        if ((yres <= 0) | (xres <= 0))
427 <                                waitflush = 1;          /* flush right after */
411 >        while (getvec(orgdir[2*ninq]) == 0 && getvec(orgdir[2*ninq+1]) == 0) {
412 >                if (orgdir[2*ninq+1][0] == 0.0 &&       /* asking for flush? */
413 >                                (orgdir[2*ninq+1][1] == 0.0) &
414 >                                (orgdir[2*ninq+1][2] == 0.0)) {
415 >                        if (ninq) {                     /* preempt our queue */
416 >                                i = next_child_nq(0);
417 >                                n = ninq;
418 >                                memset(orgdir[2*n++], 0, sizeof(FVECT)*2);
419 >                                n *= sizeof(FVECT)*2;
420 >                                if (writebuf(kida[i].pr.w, (char *)orgdir, n) != n)
421 >                                        error(SYSTEM, "pipe write error");
422 >                                kida[i].r1 = lastray+1;
423 >                                lastray += kida[i].nr = ninq;
424 >                                ninq = 0;
425 >                        }
426 >                        if ((yres <= 0) | (xres <= 0) &&
427 >                                        (lastray+1) % accumulate == 0) {
428 >                                while (next_child_nq(1) >= 0)
429 >                                        ;               /* clear the queue */
430 >                                lastdone = lastray %= accumulate;
431 >                                waitflush = 1;          /* flush next */
432 >                        }
433                          put_zero_record(++lastray);
434 <                } else {                                /* else assign ray */
435 <                        i = next_child_nq(0);
436 <                        if (writebuf(kida[i].w, (char *)orgdir,
437 <                                        sizeof(orgdir)) != sizeof(orgdir))
434 >                } else if (++ninq >= qlimit ||
435 >                            lastray/accumulate != (lastray+ninq)/accumulate) {
436 >                        i = next_child_nq(0);           /* manages output */
437 >                        n = ninq;
438 >                        if (accumulate != 1)            /* request flush? */
439 >                                memset(orgdir[2*n++], 0, sizeof(FVECT)*2);
440 >                        n *= sizeof(FVECT)*2;           /* send assignment */
441 >                        if (writebuf(kida[i].pr.w, (char *)orgdir, n) != n)
442                                  error(SYSTEM, "pipe write error");
443 <                        kida[i].running = ++lastray;
443 >                        kida[i].r1 = lastray+1;
444 >                        lastray += kida[i].nr = ninq;   /* mark as busy */
445 >                        if (lastray < lastdone) {       /* RNUMBER wrapped? */
446 >                                while (next_child_nq(1) >= 0)
447 >                                        ;
448 >                                lastray -= ninq;
449 >                                lastdone = lastray %= accumulate;
450 >                        }
451 >                        ninq = 0;
452                  }
453                  if (raysleft && !--raysleft)
454                          break;                          /* preemptive EOI */
455          }
456          while (next_child_nq(1) >= 0)           /* empty results queue */
457                  ;
458 <                                                /* output accumulated record */
459 <        if (accumulate <= 0 || account < accumulate) {
460 <                if (account < accumulate) {
407 <                        error(WARNING, "partial accumulation in final record");
408 <                        accumulate -= account;
409 <                }
410 <                for (i = 0; i < nmods; i++)
411 <                        mod_output(out_bq->mca[i]);
412 <                end_record();
413 <                free_binq(out_bq);
458 >        if (account < accumulate) {
459 >                error(WARNING, "partial accumulation in final record");
460 >                free_binq(out_bq);              /* XXX just ignore it */
461                  out_bq = NULL;
462          }
463 +        free_binq(NULL);                        /* clean up */
464 +        lu_done(&ofiletab);
465          if (raysleft)
466                  error(USER, "unexpected EOF on input");
467 <        free_binq(NULL);                        /* clean up */
467 > }
468 >
469 >
470 > /* Wait for the next available child by monitoring "to" pipes */
471 > static int
472 > next_child_ready()
473 > {
474 >        fd_set                  writeset, errset;
475 >        int                     i, n;
476 >
477 >        for (i = nchild; i--; )         /* see if there's one free first */
478 >                if (!kida[i].nr)
479 >                        return(i);
480 >                                        /* prepare select() call */
481 >        FD_ZERO(&writeset); FD_ZERO(&errset);
482 >        n = 0;
483 >        for (i = nchild; i--; ) {
484 >                FD_SET(kida[i].pr.w, &writeset);
485 >                FD_SET(kida[i].pr.r, &errset);
486 >                if (kida[i].pr.w >= n)
487 >                        n = kida[i].pr.w + 1;
488 >                if (kida[i].pr.r >= n)
489 >                        n = kida[i].pr.r + 1;
490 >        }
491 >        errno = 0;
492 >        n = select(n, NULL, &writeset, &errset, NULL);
493 >        if (n < 0)
494 >                error(SYSTEM, "select() error in next_child_ready()");
495 >        n = -1;                         /* identify waiting child */
496 >        for (i = nchild; i--; ) {
497 >                if (FD_ISSET(kida[i].pr.r, &errset))
498 >                        error(USER, "rendering process died");
499 >                if (FD_ISSET(kida[i].pr.w, &writeset))
500 >                        kida[n = i].nr = 0;
501 >        }
502 >        return(n);                      /* first available child */
503 > }
504 >
505 >
506 > /* Modified parental loop for full accumulation mode (-c 0) */
507 > void
508 > feeder_loop()
509 > {
510 >        static int      ignore_warning_given = 0;
511 >        int             ninq = 0;
512 >        FVECT           orgdir[2*MAXIQ];
513 >        int             i, n;
514 >                                        /* load rays from stdin & process */
515 > #ifdef getc_unlocked
516 >        flockfile(stdin);               /* avoid lock/unlock overhead */
517 > #endif
518 >        while (getvec(orgdir[2*ninq]) == 0 && getvec(orgdir[2*ninq+1]) == 0) {
519 >                if (orgdir[2*ninq+1][0] == 0.0 &&       /* asking for flush? */
520 >                                (orgdir[2*ninq+1][1] == 0.0) &
521 >                                (orgdir[2*ninq+1][2] == 0.0)) {
522 >                        if (!ignore_warning_given++)
523 >                                error(WARNING,
524 >                                "dummy ray(s) ignored during accumulation\n");
525 >                        continue;
526 >                }
527 >                if (++ninq >= MAXIQ) {
528 >                        i = next_child_ready();         /* get eager child */
529 >                        n = sizeof(FVECT)*2 * ninq;     /* give assignment */
530 >                        if (writebuf(kida[i].pr.w, (char *)orgdir, n) != n)
531 >                                error(SYSTEM, "pipe write error");
532 >                        kida[i].r1 = lastray+1;
533 >                        lastray += kida[i].nr = ninq;
534 >                        ninq = 0;
535 >                        if (lastray < lastdone)         /* RNUMBER wrapped? */
536 >                                lastdone = lastray = 0;
537 >                }
538 >                if (raysleft && !--raysleft)
539 >                        break;                          /* preemptive EOI */
540 >        }
541 >        if (ninq) {                             /* polish off input */
542 >                i = next_child_ready();
543 >                n = sizeof(FVECT)*2 * ninq;
544 >                if (writebuf(kida[i].pr.w, (char *)orgdir, n) != n)
545 >                        error(SYSTEM, "pipe write error");
546 >                kida[i].r1 = lastray+1;
547 >                lastray += kida[i].nr = ninq;
548 >                ninq = 0;
549 >        }
550 >        memset(orgdir, 0, sizeof(FVECT)*2);     /* get results */
551 >        for (i = nchild; i--; ) {
552 >                writebuf(kida[i].pr.w, (char *)orgdir, sizeof(FVECT)*2);
553 >                queue_results(i);
554 >        }
555 >        if (recover)                            /* and from before? */
556 >                queue_modifiers();
557 >        end_children(0);                        /* free up file descriptors */
558 >        for (i = 0; i < nmods; i++)
559 >                mod_output(out_bq->mca[i]);     /* output accumulated record */
560 >        end_record();
561 >        free_binq(out_bq);                      /* clean up */
562 >        out_bq = NULL;
563 >        free_binq(NULL);
564          lu_done(&ofiletab);
565 +        if (raysleft)
566 +                error(USER, "unexpected EOF on input");
567   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines