home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / unix / volume21 / rayshade / patch5 < prev    next >
Encoding:
Text File  |  1990-04-11  |  49.3 KB  |  1,965 lines
  1. System: rayshade version 3.0
  2. Patch #: 5
  3. Priority: HIGH
  4. Subject: Patch #5
  5.  
  6. Description:
  7.  
  8.     Rayshade assumed that a freed pointer could be safely
  9.     dereferenced up until the next call to malloc().  It
  10.     is now assumed that free() is destructive.
  11.  
  12.     Pointers were being freed twice in several places.
  13.  
  14.     The test for edge-on intersection with planar primitives
  15.     was not particularly robust.  In particular, near edge-on
  16.     intersections with textured planes caused TLOSS errors under
  17.     certain circumstances.
  18.  
  19.     HUGE was too huge for some machines.  FAR_AWAY is used instead.
  20.  
  21.     The matrices to transform cylinders and cones to their canonical
  22.     forms were applied as if they were user-specified.  This caused
  23.     cylinders and cones to be always textured in "canonical"
  24.     space.  The matrices are now stored and applied locally to the
  25.     cylinder/cone routines.
  26.  
  27.     Simplified cone and cylinder intersection routines by making
  28.     canonical primitives truly canonical.
  29.  
  30.     It was possible to divide by zero when rendering cones.
  31.  
  32.     The validity of colormap entires was not checked.
  33.  
  34.     times() is (void) on some SYSV systems and (long) on others.
  35.     Removed checking of its return value altogether.
  36.  
  37.     Various pieces of lint removed throughout the code.
  38.  
  39.     Replaced outdated code in primitive tutorial.
  40.  
  41.  
  42. Fix:    From rn, say "| patch -p -N -d DIR", where DIR is your top-level
  43.     rayshade source directory.  Outside of rn, say
  44.     "cd DIR; patch -p -N <thisarticle".
  45.     If you don't have the patch program, apply the following by hand,
  46.     or get patch (version 2.0, latest patchlevel).
  47.  
  48.     After patching:
  49.         make depend
  50.         make
  51.         make install
  52.  
  53.     If patch indicates that patchlevel is the wrong version, you may need
  54.     to apply one or more previous patches, or the patch may already
  55.     have been applied.  See the patchlevel.h file to find out what has or
  56.     has not been applied.  In any event, don't continue with the patch.
  57.  
  58.     If you are missing previous patches they can be obtained via anonymous
  59.     ftp from weedeater.math.yale.edu (130.132.23.17).
  60.  
  61.     Craig Kolb
  62.     kolb@yale.edu
  63.  
  64.  
  65. Index: src/patchlevel.h
  66. Prereq: 4
  67. 1c1
  68. < #define PATCHLEVEL 4
  69. ---
  70. > #define PATCHLEVEL 5
  71.  
  72. Index: Changes
  73. *** Changes.old    Tue Apr 10 16:50:08 1990
  74. --- Changes    Tue Apr 10 16:50:09 1990
  75. ***************
  76. *** 1,3 ****
  77. --- 1,15 ----
  78. + Patch #5:
  79. + ---------
  80. + Changes to reflect that free() may be destructive.
  81. + Fixed several other free()-related bugs.
  82. + Cylinders and cones were always textured in canonical cylinder/cone space.
  83. + Fixed divide by zero bug in cylinder code.
  84. + Optimizations in cylinder and cone code.
  85. + HUGE in bounds.c was too huge for some machines.  Replaced by FAR_AWAY.
  86. + Improved handling of near edge-on intersection with triangles and planes.
  87. + Added check for invalid colormap entries.
  88. + Lint removal.
  89. + 
  90.   Patch #4:
  91.   ---------
  92.   Added gloss texture.
  93.  
  94. Index: src/bounds.c
  95. *** src/bounds.c.old    Tue Apr 10 16:47:54 1990
  96. --- src/bounds.c    Tue Apr 10 16:47:54 1990
  97. ***************
  98. *** 18,26 ****
  99.    * name of the person performing the modification, the date of modification,
  100.    * and the reason for such modification.
  101.    *
  102. !  * $Id: bounds.c,v 3.0 89/10/27 02:05:46 craig Exp $
  103.    *
  104.    * $Log:    bounds.c,v $
  105.    * Revision 3.0  89/10/27  02:05:46  craig
  106.    * Baseline for first official release.
  107.    * 
  108. --- 18,29 ----
  109.    * name of the person performing the modification, the date of modification,
  110.    * and the reason for such modification.
  111.    *
  112. !  * $Id: bounds.c,v 3.0.1.1 90/04/04 19:03:47 craig Exp $
  113.    *
  114.    * $Log:    bounds.c,v $
  115. +  * Revision 3.0.1.1  90/04/04  19:03:47  craig
  116. +  * patch5: Substituded FAR_AWAY for HUGE in bounds initialization routines.
  117. +  * 
  118.    * Revision 3.0  89/10/27  02:05:46  craig
  119.    * Baseline for first official release.
  120.    * 
  121. ***************
  122. *** 31,40 ****
  123.   #include "typedefs.h"
  124.   #include "funcdefs.h"
  125.   
  126. - #ifndef HUGE
  127. - #define HUGE    1.701411e38
  128. - #endif
  129. - 
  130.   /*
  131.    * Ray-bounding box intersection test.
  132.    */
  133. --- 34,39 ----
  134. ***************
  135. *** 188,195 ****
  136.   init_bounds(bounds)
  137.   double bounds[2][3];
  138.   {
  139. !     bounds[LOW][X] = bounds[LOW][Y] = bounds[LOW][Z] = HUGE;
  140. !     bounds[HIGH][X] = bounds[HIGH][Y] = bounds[HIGH][Z] = -HUGE;
  141.   }
  142.   
  143.   /*
  144. --- 187,194 ----
  145.   init_bounds(bounds)
  146.   double bounds[2][3];
  147.   {
  148. !     bounds[LOW][X] = bounds[LOW][Y] = bounds[LOW][Z] = FAR_AWAY;
  149. !     bounds[HIGH][X] = bounds[HIGH][Y] = bounds[HIGH][Z] = -FAR_AWAY;
  150.   }
  151.   
  152.   /*
  153.  
  154. Index: src/cone.c
  155. *** src/cone.c.old    Tue Apr 10 16:47:59 1990
  156. --- src/cone.c    Tue Apr 10 16:48:00 1990
  157. ***************
  158. *** 18,26 ****
  159.    * name of the person performing the modification, the date of modification,
  160.    * and the reason for such modification.
  161.    *
  162. !  * $Id: cone.c,v 3.0.1.3 90/02/12 13:20:58 craig Exp $
  163.    *
  164.    * $Log:    cone.c,v $
  165.    * Revision 3.0.1.3  90/02/12  13:20:58  craig
  166.    * patch4: Changes to avoid rotation about null axis.
  167.    * 
  168. --- 18,33 ----
  169.    * name of the person performing the modification, the date of modification,
  170.    * and the reason for such modification.
  171.    *
  172. !  * $Id: cone.c,v 3.0.1.5 90/04/09 14:30:08 craig Exp $
  173.    *
  174.    * $Log:    cone.c,v $
  175. +  * Revision 3.0.1.5  90/04/09  14:30:08  craig
  176. +  * patch5: Transformation information now stored locally.
  177. +  * patch5: Canonical cone now truly canonical.
  178. +  * 
  179. +  * Revision 3.0.1.4  90/04/04  14:51:25  craig
  180. +  * patch5: Fixed divide by zero problem in intcone().
  181. +  * 
  182.    * Revision 3.0.1.3  90/02/12  13:20:58  craig
  183.    * patch4: Changes to avoid rotation about null axis.
  184.    * 
  185. ***************
  186. *** 41,56 ****
  187.   #include "constants.h"
  188.   
  189.   Object *
  190. ! makcone(surf, cent, ax, br, ar, trans)
  191.   char *surf;
  192.   Vector *cent, *ax;
  193.   double br, ar;
  194. - TransInfo *trans;
  195.   {
  196.       Cone *cone;
  197.       Primitive *prim;
  198.       Object *newobj;
  199. !     double len, dtmp;
  200.       Vector axis, base, tmp;
  201.   
  202.       prim = mallocprim();
  203. --- 48,62 ----
  204.   #include "constants.h"
  205.   
  206.   Object *
  207. ! makcone(surf, cent, ax, br, ar)
  208.   char *surf;
  209.   Vector *cent, *ax;
  210.   double br, ar;
  211.   {
  212.       Cone *cone;
  213.       Primitive *prim;
  214.       Object *newobj;
  215. !     double tantheta, lprime, tlen, len, dtmp;
  216.       Vector axis, base, tmp;
  217.   
  218.       prim = mallocprim();
  219. ***************
  220. *** 61,78 ****
  221.       prim->objpnt.p_cone = cone;
  222.   
  223.       /*
  224. !      * Cones are defined by a basepoint, an apex point, and
  225. !      * base and apex radii.  The cone is stored as
  226. !      * the origin of the cone, the change in radius per
  227. !      * unit distance from the origin of the cone, the maximum z-value
  228. !      * of the cone, and "start_pos",
  229. !      * the distance along the axis from the cone origin where
  230. !      * the first endcap appears (where the passed "basepoint"
  231. !      * appears).
  232.        *
  233. !      * The intcone() routine intersects a ray with a cone aligned
  234. !      * along the Z axis.  Thus, we must define a transformation
  235. !      * matrix which will transform an axis-aligned cone to the desired.
  236.        */
  237.   
  238.       /*
  239. --- 67,78 ----
  240.       prim->objpnt.p_cone = cone;
  241.   
  242.       /*
  243. !      * To render a cone, we transform the desired cone into
  244. !      * a canonical, Z-axis aligned, unit length, unit radius
  245. !      * at the apex cone.
  246.        *
  247. !      * Here, we construct the transformation matrix to map
  248. !      * from cone<-->world space.
  249.        */
  250.   
  251.       /*
  252. ***************
  253. *** 93,99 ****
  254.            * If the base and apex radii are equal, then we
  255.            * can treat the cone as a cylinder.
  256.            */
  257. !         return makcyl(surf, cent, ax, br, trans);
  258.       }
  259.       /*
  260.        * Find the axis and axis length.
  261. --- 93,99 ----
  262.            * If the base and apex radii are equal, then we
  263.            * can treat the cone as a cylinder.
  264.            */
  265. !         return makcyl(surf, cent, ax, br);
  266.       }
  267.       /*
  268.        * Find the axis and axis length.
  269. ***************
  270. *** 107,135 ****
  271.           free((char *)newobj);
  272.           return (Object *)0;
  273.       }
  274. -     cone->apex_rad = ar;
  275.       /*
  276.        * "tantheta" is the change in radius per unit length along
  277.        * the cone axis.
  278.        */
  279. !     cone->tantheta = (ar - br) / len;
  280.       /*
  281. !      * Start pos defines the point along the axis where the first
  282.        * endcap should be placed.
  283.        */
  284. !     cone->start_pos = br / cone->tantheta;
  285.       /*
  286.        * Find the true base (origin) of the cone.
  287.        */
  288. !     scalar_prod(-cone->start_pos, axis, &base);
  289.       vecadd(base, *cent, &base);
  290.       /*
  291. !      * The apex endcap is placed cone->len units from the cone
  292. !      * origin.
  293.        */
  294. !     cone->end_pos = cone->start_pos + len;
  295.       /*
  296. !      * Calculate rotation matrix to map from world space to cone space.
  297.        */
  298.       if (abs(axis.z) == 1.) {
  299.           tmp.x = 1;
  300. --- 107,138 ----
  301.           free((char *)newobj);
  302.           return (Object *)0;
  303.       }
  304.       /*
  305.        * "tantheta" is the change in radius per unit length along
  306.        * the cone axis.
  307.        */
  308. !     tantheta = (ar - br) / len;
  309.       /*
  310. !      * lprime defines the distance along the axis where the first
  311.        * endcap should be placed.
  312.        */
  313. !     lprime = br / tantheta;
  314.       /*
  315.        * Find the true base (origin) of the cone.
  316.        */
  317. !     scalar_prod(-lprime, axis, &base);
  318.       vecadd(base, *cent, &base);
  319.       /*
  320. !      * tlen is the total length of the cone.
  321. !      */
  322. !     tlen = lprime + len;
  323. !     /*
  324. !      * start_dist is the distance from the origin of the canonical
  325. !      * cone at which the first endcap appears.
  326.        */
  327. !     cone->start_dist = lprime / tlen;
  328.       /*
  329. !      * Calculate transformation to map from cone space to world space.
  330.        */
  331.       if (abs(axis.z) == 1.) {
  332.           tmp.x = 1;
  333. ***************
  334. *** 139,147 ****
  335.           tmp.y = -axis.x;
  336.           tmp.z = 0.;
  337.       }
  338. !     RS_rotate(trans, &tmp, acos(axis.z));
  339. !     RS_translate(trans, &base);
  340. !     cone->tantheta *= cone->tantheta;
  341.   
  342.       return newobj;
  343.   }
  344. --- 142,152 ----
  345.           tmp.y = -axis.x;
  346.           tmp.z = 0.;
  347.       }
  348. !     init_trans(&cone->trans.obj2world);
  349. !     RS_scale(&cone->trans.obj2world, ar, ar, tlen);
  350. !     RS_rotate(&cone->trans.obj2world, &tmp, acos(axis.z));
  351. !     RS_translate(&cone->trans.obj2world, &base);
  352. !     invert_trans(&cone->trans.world2obj, &cone->trans.obj2world);
  353.   
  354.       return newobj;
  355.   }
  356. ***************
  357. *** 156,162 ****
  358.   Primitive       *obj;
  359.   {
  360.       double t1, t2, a, b, c, disc, zpos, et1, et2;
  361. !     double x, y;
  362.       extern unsigned long primtests[];
  363.       Cone *cone;
  364.   
  365. --- 161,170 ----
  366.   Primitive       *obj;
  367.   {
  368.       double t1, t2, a, b, c, disc, zpos, et1, et2;
  369. !     double x, y, distfact;
  370. !     Ray newray;
  371. !     Vector nray, npos;
  372. !     extern double TransformRay();
  373.       extern unsigned long primtests[];
  374.       Cone *cone;
  375.   
  376. ***************
  377. *** 164,183 ****
  378.       cone = obj->objpnt.p_cone;
  379.   
  380.       /*
  381. !      * Recall that 'tantheta' is really tantheta^2.
  382.        */
  383. !     a = ray->x * ray->x + ray->y * ray->y - ray->z*ray->z*cone->tantheta;
  384. !     b = ray->x * pos->x + ray->y * pos->y - cone->tantheta*ray->z*pos->z;
  385. !     c = pos->x*pos->x + pos->y*pos->y - cone->tantheta*pos->z*pos->z;
  386.   
  387.       if (equal(a, 0.)) {
  388.           /*
  389.            * Only one intersection point...
  390.            */
  391.           t1 = -c / b;
  392. !         zpos = pos->z + t1 * ray->z;
  393. !         if (t1 < EPSILON || zpos < cone->start_pos ||
  394. !             zpos > cone->end_pos)
  395.               t1 = FAR_AWAY;
  396.           t2 = FAR_AWAY;
  397.       } else {
  398. --- 172,196 ----
  399.       cone = obj->objpnt.p_cone;
  400.   
  401.       /*
  402. !      * Transform ray from world to cone space.
  403.        */
  404. !     newray.pos = *pos;
  405. !     newray.dir = *ray;
  406. !     distfact = TransformRay(&newray, &cone->trans.world2obj);
  407. !     nray = newray.dir;
  408. !     npos = newray.pos;
  409.   
  410. +     a = nray.x * nray.x + nray.y * nray.y - nray.z*nray.z;
  411. +     b = nray.x * npos.x + nray.y * npos.y - nray.z*npos.z;
  412. +     c = npos.x*npos.x + npos.y*npos.y - npos.z*npos.z;
  413. + 
  414.       if (equal(a, 0.)) {
  415.           /*
  416.            * Only one intersection point...
  417.            */
  418.           t1 = -c / b;
  419. !         zpos = npos.z + t1 * nray.z;
  420. !         if (t1 < EPSILON || zpos < cone->start_dist || zpos > 1.)
  421.               t1 = FAR_AWAY;
  422.           t2 = FAR_AWAY;
  423.       } else {
  424. ***************
  425. *** 190,220 ****
  426.           /*
  427.            * Clip intersection points.
  428.            */
  429. !         zpos = pos->z + t1 * ray->z;
  430. !         if (t1 < EPSILON || zpos < cone->start_pos ||
  431. !             zpos > cone->end_pos)
  432.               t1 = FAR_AWAY;
  433. !         zpos = pos->z + t2 * ray->z;
  434. !         if (t2 < EPSILON || zpos < cone->start_pos ||
  435. !             zpos > cone->end_pos)
  436.               t2 = FAR_AWAY;
  437.       }
  438.       /*
  439.        * Find t for both endcaps.
  440.        */
  441. !     et1 = (cone->start_pos - pos->z) / ray->z;
  442. !     x = pos->x + et1 * ray->x;
  443. !     y = pos->y + et1 * ray->y;
  444. !     if (x*x + y*y > cone->start_pos*cone->start_pos*cone->tantheta)
  445.           et1 = FAR_AWAY;
  446. !     et2 = (cone->end_pos - pos->z) / ray->z;
  447. !     x = pos->x + et2 * ray->x;
  448. !     y = pos->y + et2 * ray->y;
  449. !     if (x*x + y*y > cone->end_pos*cone->end_pos*cone->tantheta)
  450.           et2 = FAR_AWAY;
  451.   
  452.       t1 = min(t1, min(t2, min(et1, et2)));
  453. !     return (t1 == FAR_AWAY ? 0. : t1);
  454.   }
  455.   
  456.   /*
  457. --- 203,237 ----
  458.           /*
  459.            * Clip intersection points.
  460.            */
  461. !         zpos = npos.z + t1 * nray.z;
  462. !         if (t1 < EPSILON || zpos < cone->start_dist || zpos > 1.)
  463.               t1 = FAR_AWAY;
  464. !         zpos = npos.z + t2 * nray.z;
  465. !         if (t2 < EPSILON || zpos < cone->start_dist || zpos > 1.)
  466.               t2 = FAR_AWAY;
  467.       }
  468. + 
  469. +     if (equal(nray.z, 0.)) {
  470. +         t1 = min(t1, t2);
  471. +         return (t1 == FAR_AWAY ? 0. : t1 / distfact);
  472. +     }
  473. + 
  474.       /*
  475.        * Find t for both endcaps.
  476.        */
  477. !     et1 = (cone->start_dist - npos.z) / nray.z;
  478. !     x = npos.x + et1 * nray.x;
  479. !     y = npos.y + et1 * nray.y;
  480. !     if (x*x + y*y > cone->start_dist*cone->start_dist)
  481.           et1 = FAR_AWAY;
  482. !     et2 = (1. - npos.z) / nray.z;
  483. !     x = npos.x + et2 * nray.x;
  484. !     y = npos.y + et2 * nray.y;
  485. !     if (x*x + y*y > 1.)
  486.           et2 = FAR_AWAY;
  487.   
  488.       t1 = min(t1, min(t2, min(et1, et2)));
  489. !     return (t1 == FAR_AWAY ? 0. : t1 / distfact);
  490.   }
  491.   
  492.   /*
  493. ***************
  494. *** 225,237 ****
  495.   Primitive *obj;
  496.   {
  497.       Cone *cone;
  498.   
  499.       cone = obj->objpnt.p_cone;
  500.   
  501. !     if (equal(pos->z, cone->start_pos)) {
  502.           nrm->x = nrm->y = 0.;
  503.           nrm->z = -1.;
  504. !     } else if (equal(pos->z, cone->end_pos)) {
  505.           nrm->x = nrm->y = 0.;
  506.           nrm->z = 1.;
  507.       } else {
  508. --- 242,261 ----
  509.   Primitive *obj;
  510.   {
  511.       Cone *cone;
  512. +     Vector npos;
  513.   
  514.       cone = obj->objpnt.p_cone;
  515.   
  516. !     /*
  517. !      * Transform intersection point to cone space.
  518. !      */
  519. !     npos = *pos;
  520. !     transform_point(&npos, &cone->trans.world2obj);
  521. !     
  522. !     if (equal(npos.z, cone->start_dist)) {
  523.           nrm->x = nrm->y = 0.;
  524.           nrm->z = -1.;
  525. !     } else if (equal(npos.z, 1)) {
  526.           nrm->x = nrm->y = 0.;
  527.           nrm->z = 1.;
  528.       } else {
  529. ***************
  530. *** 239,248 ****
  531.            * The following is equal to
  532.            * (pos X (0, 0, 1)) X pos
  533.            */
  534. !         nrm->x = pos->x * pos->z;
  535. !         nrm->y = pos->y * pos->z;
  536. !         nrm->z = -pos->x * pos->x - pos->y * pos->y;
  537.       }
  538.   }
  539.   
  540.   /*
  541. --- 263,276 ----
  542.            * The following is equal to
  543.            * (pos X (0, 0, 1)) X pos
  544.            */
  545. !         nrm->x = npos.x * npos.z;
  546. !         nrm->y = npos.y * npos.z;
  547. !         nrm->z = -npos.x * npos.x - npos.y * npos.y;
  548.       }
  549. +     /*
  550. +      * Transform normal back to world space.
  551. +      */
  552. +     TransformNormal(nrm, &cone->trans.world2obj);
  553.   }
  554.   
  555.   /*
  556. ***************
  557. *** 256,263 ****
  558.   
  559.       cone = o->objpnt.p_cone;
  560.   
  561. !     bounds[LOW][X] = bounds[LOW][Y] = -cone->apex_rad;
  562. !     bounds[HIGH][X] = bounds[HIGH][Y] = cone->apex_rad;
  563. !     bounds[LOW][Z] = cone->start_pos;
  564. !     bounds[HIGH][Z] = cone->end_pos;
  565.   }
  566. --- 284,295 ----
  567.   
  568.       cone = o->objpnt.p_cone;
  569.   
  570. !     bounds[LOW][X] = bounds[LOW][Y] = -1;
  571. !     bounds[HIGH][X] = bounds[HIGH][Y] = 1;
  572. !     bounds[LOW][Z] = cone->start_dist;
  573. !     bounds[HIGH][Z] = 1;
  574. !     /*
  575. !      * Transform bounding box to world space.
  576. !      */
  577. !     transform_bounds(&cone->trans.obj2world, bounds);
  578.   }
  579.  
  580. Index: src/cylinder.c
  581. *** src/cylinder.c.old    Tue Apr 10 16:48:04 1990
  582. --- src/cylinder.c    Tue Apr 10 16:48:05 1990
  583. ***************
  584. *** 18,26 ****
  585.    * name of the person performing the modification, the date of modification,
  586.    * and the reason for such modification.
  587.    *
  588. !  * $Id: cylinder.c,v 3.0.1.4 90/03/07 22:43:04 craig Exp $
  589.    *
  590.    * $Log:    cylinder.c,v $
  591.    * Revision 3.0.1.4  90/03/07  22:43:04  craig
  592.    * patch4: Fixed typo in previous fix.
  593.    * 
  594. --- 18,30 ----
  595.    * name of the person performing the modification, the date of modification,
  596.    * and the reason for such modification.
  597.    *
  598. !  * $Id: cylinder.c,v 3.0.1.5 90/04/09 14:29:43 craig Exp $
  599.    *
  600.    * $Log:    cylinder.c,v $
  601. +  * Revision 3.0.1.5  90/04/09  14:29:43  craig
  602. +  * patch5: Transformation information now stored locally.
  603. +  * patch5: Canonical cylinder now truly canonical.
  604. +  * 
  605.    * Revision 3.0.1.4  90/03/07  22:43:04  craig
  606.    * patch4: Fixed typo in previous fix.
  607.    * 
  608. ***************
  609. *** 44,54 ****
  610.   #include "constants.h"
  611.   
  612.   Object *
  613. ! makcyl(surf, cent, ax, r, trans)
  614.   char *surf;
  615.   Vector *cent, *ax;
  616.   double r;
  617. - TransInfo *trans;
  618.   {
  619.       Cylinder *cyl;
  620.       Primitive *prim;
  621. --- 48,57 ----
  622.   #include "constants.h"
  623.   
  624.   Object *
  625. ! makcyl(surf, cent, ax, r)
  626.   char *surf;
  627.   Vector *cent, *ax;
  628.   double r;
  629.   {
  630.       Cylinder *cyl;
  631.       Primitive *prim;
  632. ***************
  633. *** 80,89 ****
  634.           return (Object *)0;
  635.       }
  636.   
  637. -     cyl->rad = r*r;
  638. -     cyl->len = len;
  639.       /*
  640. !      * Define matrix to transform from axis-aligned to desired cylinder.
  641.        */
  642.       if (abs(axis.z) == 1.) {
  643.           dir.x = 1.;
  644. --- 83,91 ----
  645.           return (Object *)0;
  646.       }
  647.   
  648.       /*
  649. !      * Define matrix to transform from Z axis aligned unit cylinder
  650. !      * to desired cylinder.
  651.        */
  652.       if (abs(axis.z) == 1.) {
  653.           dir.x = 1.;
  654. ***************
  655. *** 93,101 ****
  656.           dir.y = -axis.x;
  657.           dir.z = 0.;
  658.       }
  659. !     RS_rotate(trans, &dir, acos(axis.z));
  660. !     RS_translate(trans, cent);
  661. ! 
  662.       return newobj;
  663.   }
  664.   
  665. --- 95,105 ----
  666.           dir.y = -axis.x;
  667.           dir.z = 0.;
  668.       }
  669. !     init_trans(&cyl->trans.obj2world);
  670. !     RS_scale(&cyl->trans.obj2world, r, r, len);
  671. !     RS_rotate(&cyl->trans.obj2world, &dir, acos(axis.z));
  672. !     RS_translate(&cyl->trans.obj2world, cent);
  673. !     invert_trans(&cyl->trans.world2obj, &cyl->trans.obj2world);
  674.       return newobj;
  675.   }
  676.   
  677. ***************
  678. *** 108,131 ****
  679.   Primitive       *obj;
  680.   {
  681.       double t1, t2, a, b, c, zpos1, zpos2, et1, et2, x, y, disc;
  682.       extern unsigned long primtests[];
  683.       Cylinder *cyl;
  684.   
  685.       primtests[CYL]++;
  686.       cyl = obj->objpnt.p_cylinder;
  687.   
  688. !     a = ray->x * ray->x + ray->y * ray->y;
  689. !     c = pos->x*pos->x + pos->y*pos->y - cyl->rad;
  690.   
  691. !     if (a < EPSILON*EPSILON) {    /* |ray->z| == 1. */
  692.           if(c < EPSILON*EPSILON)    /* Within endcap */
  693.               /* Wrong if origin is inside cylinder. */
  694. !             return min(-pos->z / ray->z,
  695. !                   (cyl->len - pos->z) / ray->z);
  696.           return 0.;
  697.       }
  698.   
  699. !     b = ray->x * pos->x + ray->y * pos->y;
  700.       disc = b*b - a*c;
  701.       if(disc < 0.)
  702.           return 0.;
  703. --- 112,148 ----
  704.   Primitive       *obj;
  705.   {
  706.       double t1, t2, a, b, c, zpos1, zpos2, et1, et2, x, y, disc;
  707. +     double distfact;
  708.       extern unsigned long primtests[];
  709. +     extern double TransformRay();
  710. +     Ray newray;
  711. +     Vector nray, npos;
  712.       Cylinder *cyl;
  713.   
  714.       primtests[CYL]++;
  715.       cyl = obj->objpnt.p_cylinder;
  716.   
  717. !     /*
  718. !      * Transform ray into canonical cylinder space.
  719. !      */
  720. !     newray.dir = *ray;
  721. !     newray.pos = *pos;
  722. !     distfact = TransformRay(&newray, &cyl->trans.world2obj);
  723. !     nray = newray.dir;
  724. !     npos = newray.pos;
  725. ! 
  726. !     a = nray.x * nray.x + nray.y * nray.y;
  727. !     c = npos.x*npos.x + npos.y*npos.y - 1;
  728.   
  729. !     if (a < EPSILON*EPSILON) {    /* |nray.z| == 1. */
  730.           if(c < EPSILON*EPSILON)    /* Within endcap */
  731.               /* Wrong if origin is inside cylinder. */
  732. !             return min(-npos.z / nray.z,
  733. !                   (1. - npos.z) / nray.z) / distfact;
  734.           return 0.;
  735.       }
  736.   
  737. !     b = nray.x * npos.x + nray.y * npos.y;
  738.       disc = b*b - a*c;
  739.       if(disc < 0.)
  740.           return 0.;
  741. ***************
  742. *** 134,142 ****
  743.       t2 = (-b - disc) / a;
  744.       if(t1 < EPSILON && t2 < EPSILON)
  745.           return 0.;
  746. !     zpos1 = pos->z + t1 * ray->z;
  747. !     zpos2 = pos->z + t2 * ray->z;
  748. !     if ((zpos1 > cyl->len && zpos2 > cyl->len) ||
  749.           (zpos1 < 0. && zpos2 < 0.))
  750.           return 0.;
  751.       if (t1 < EPSILON)
  752. --- 151,159 ----
  753.       t2 = (-b - disc) / a;
  754.       if(t1 < EPSILON && t2 < EPSILON)
  755.           return 0.;
  756. !     zpos1 = npos.z + t1 * nray.z;
  757. !     zpos2 = npos.z + t2 * nray.z;
  758. !     if ((zpos1 > 1 && zpos2 > 1) ||
  759.           (zpos1 < 0. && zpos2 < 0.))
  760.           return 0.;
  761.       if (t1 < EPSILON)
  762. ***************
  763. *** 143,176 ****
  764.           t1 = FAR_AWAY;
  765.       if (t2 < EPSILON)
  766.           t2 = FAR_AWAY;
  767. -     if (t1 == FAR_AWAY && t2 == FAR_AWAY)
  768. -         return 0.;
  769.       /*
  770.        * Don't bother checking endcaps if both intersection points
  771.        * are on the cylinder.
  772.        */
  773. !     if ((zpos1 > 0. && zpos1 < cyl->len && zpos2 > 0. && zpos2 < cyl->len))
  774. !         return min(t1, t2);
  775.       /*
  776.         * It's possible to get rid of the ray-disc intersection tests
  777.        * (by looking at t1, t2 and zpos1, zpos), but the code gets messy.
  778.        */
  779. !     if (zpos1 < 0. || zpos1 > cyl->len)
  780.           t1 = FAR_AWAY;
  781. !     if (zpos2 < 0. || zpos2 > cyl->len)
  782.           t2 = FAR_AWAY;
  783. !     et1 = -pos->z / ray->z;
  784. !     x = pos->x + et1 * ray->x;
  785. !     y = pos->y + et1 * ray->y;
  786. !     if (x*x + y*y > cyl->rad)
  787.           et1 = FAR_AWAY;
  788. !     et2 = (cyl->len - pos->z) / ray->z;
  789. !     x = pos->x + et2 * ray->x;
  790. !     y = pos->y + et2 * ray->y;
  791. !     if (x*x + y*y > cyl->rad)
  792.           et2 = FAR_AWAY;
  793.       t1 = min(t1, min(t2, min(et1, et2)));
  794. !     return (t1 == FAR_AWAY ? 0. : t1);
  795.   }
  796.   
  797.   nrmcyl(pos, obj, nrm)
  798. --- 160,191 ----
  799.           t1 = FAR_AWAY;
  800.       if (t2 < EPSILON)
  801.           t2 = FAR_AWAY;
  802.       /*
  803.        * Don't bother checking endcaps if both intersection points
  804.        * are on the cylinder.
  805.        */
  806. !     if ((zpos1 > 0. && zpos1 < 1. && zpos2 > 0. && zpos2 < 1.))
  807. !         return min(t1, t2) / distfact;
  808.       /*
  809.         * It's possible to get rid of the ray-disc intersection tests
  810.        * (by looking at t1, t2 and zpos1, zpos), but the code gets messy.
  811.        */
  812. !     if (zpos1 < 0. || zpos1 > 1.)
  813.           t1 = FAR_AWAY;
  814. !     if (zpos2 < 0. || zpos2 > 1.)
  815.           t2 = FAR_AWAY;
  816. !     et1 = -npos.z / nray.z;
  817. !     x = npos.x + et1 * nray.x;
  818. !     y = npos.y + et1 * nray.y;
  819. !     if (x*x + y*y > 1.)
  820.           et1 = FAR_AWAY;
  821. !     et2 = (1. - npos.z) / nray.z;
  822. !     x = npos.x + et2 * nray.x;
  823. !     y = npos.y + et2 * nray.y;
  824. !     if (x*x + y*y > 1.)
  825.           et2 = FAR_AWAY;
  826.       t1 = min(t1, min(t2, min(et1, et2)));
  827. !     return (t1 == FAR_AWAY ? 0. : t1 / distfact);
  828.   }
  829.   
  830.   nrmcyl(pos, obj, nrm)
  831. ***************
  832. *** 178,190 ****
  833.   Primitive *obj;
  834.   {
  835.       Cylinder *cyl;
  836.       double dist;
  837.   
  838.       cyl = obj->objpnt.p_cylinder;
  839.   
  840. !     dist = pos->x*pos->x + pos->y*pos->y;
  841. !     if (dist+EPSILON < cyl->rad) {
  842. !         if (equal(pos->z,0.)) {
  843.               /*
  844.                * Hit on lower endcap.
  845.                */
  846. --- 193,212 ----
  847.   Primitive *obj;
  848.   {
  849.       Cylinder *cyl;
  850. +     Vector npos;
  851.       double dist;
  852.   
  853.       cyl = obj->objpnt.p_cylinder;
  854.   
  855. !     /*
  856. !      * Transform position into cylinder space.
  857. !      */
  858. !     npos = *pos;
  859. !     transform_point(&npos, &cyl->trans.world2obj);
  860. ! 
  861. !     dist = npos.x*npos.x + npos.y*npos.y;
  862. !     if (dist+EPSILON < 1) {
  863. !         if (equal(npos.z,0.)) {
  864.               /*
  865.                * Hit on lower endcap.
  866.                */
  867. ***************
  868. *** 198,208 ****
  869.               nrm->z = 1.;
  870.           }
  871.       } else {    /* Hit along cylinder. */
  872. !         nrm->x = pos->x;
  873. !         nrm->y = pos->y;
  874.           nrm->z = 0.;
  875.           /* Will be normalized by ShadeRay(). */
  876.       }
  877.   }
  878.   
  879.   cylextent(o, bounds)
  880. --- 220,235 ----
  881.               nrm->z = 1.;
  882.           }
  883.       } else {    /* Hit along cylinder. */
  884. !         nrm->x = npos.x;
  885. !         nrm->y = npos.y;
  886.           nrm->z = 0.;
  887.           /* Will be normalized by ShadeRay(). */
  888.       }
  889. + 
  890. +     /*
  891. +      * Tranform normal back to world space.
  892. +      */
  893. +     TransformNormal(nrm, &cyl->trans.world2obj);
  894.   }
  895.   
  896.   cylextent(o, bounds)
  897. ***************
  898. *** 210,221 ****
  899.   double bounds[2][3];
  900.   {
  901.       Cylinder *cyl;
  902. !     double r;
  903.       cyl = o->objpnt.p_cylinder;
  904.   
  905. !     r = sqrt(cyl->rad);
  906. !     bounds[LOW][X] = bounds[LOW][Y] = -r;
  907. !     bounds[HIGH][X] = bounds[HIGH][Y] = r;
  908.       bounds[LOW][Z] = 0.;
  909. !     bounds[HIGH][Z] = cyl->len;
  910.   }
  911. --- 237,251 ----
  912.   double bounds[2][3];
  913.   {
  914.       Cylinder *cyl;
  915. ! 
  916.       cyl = o->objpnt.p_cylinder;
  917.   
  918. !     bounds[LOW][X] = bounds[LOW][Y] = -1;
  919. !     bounds[HIGH][X] = bounds[HIGH][Y] = 1;
  920.       bounds[LOW][Z] = 0.;
  921. !     bounds[HIGH][Z] = 1;
  922. !     /*
  923. !      * Transform bounding box to world space.
  924. !      */
  925. !     transform_bounds(&cyl->trans.obj2world, bounds);
  926.   }
  927.  
  928. Index: src/hf.c
  929. *** src/hf.c.old    Tue Apr 10 16:48:12 1990
  930. --- src/hf.c    Tue Apr 10 16:48:13 1990
  931. ***************
  932. *** 18,26 ****
  933.    * name of the person performing the modification, the date of modification,
  934.    * and the reason for such modification.
  935.    *
  936. !  * $Id: hf.c,v 3.0.1.3 90/02/12 13:30:14 craig Exp $
  937.    *
  938.    * $Log:    hf.c,v $
  939.    * Revision 3.0.1.3  90/02/12  13:30:14  craig
  940.    * patch4: Fixed bug involving computation of hf->lsize[i]
  941.    * 
  942. --- 18,29 ----
  943.    * name of the person performing the modification, the date of modification,
  944.    * and the reason for such modification.
  945.    *
  946. !  * $Id: hf.c,v 3.0.1.4 90/04/04 14:52:30 craig Exp $
  947.    *
  948.    * $Log:    hf.c,v $
  949. +  * Revision 3.0.1.4  90/04/04  14:52:30  craig
  950. +  * patch5: Lint removal.
  951. +  * 
  952.    * Revision 3.0.1.3  90/02/12  13:30:14  craig
  953.    * patch4: Fixed bug involving computation of hf->lsize[i]
  954.    * 
  955. ***************
  956. *** 118,124 ****
  957.       /*
  958.        * Get HF size.
  959.        */
  960. !     if (fread(&hf->size, sizeof(int), 1, fp) == 0)
  961.           yyerror("Cannot read height field size.");
  962.   
  963.       hf->data = (float **)share_malloc(hf->size * sizeof(float *));
  964. --- 121,127 ----
  965.       /*
  966.        * Get HF size.
  967.        */
  968. !     if (fread((char *)&hf->size, sizeof(int), 1, fp) == 0)
  969.           yyerror("Cannot read height field size.");
  970.   
  971.       hf->data = (float **)share_malloc(hf->size * sizeof(float *));
  972. ***************
  973. *** 127,133 ****
  974.           /*
  975.            * Read in row of HF data.
  976.            */
  977. !         if (fread(hf->data[i],sizeof(float),hf->size,fp) != hf->size)
  978.               yyerror("Not enough heightfield data.");
  979.           for (j = 0; j < hf->size; j++) {
  980.               val = hf->data[i][j];
  981. --- 130,136 ----
  982.           /*
  983.            * Read in row of HF data.
  984.            */
  985. !         if (fread((char *)hf->data[i],sizeof(float),hf->size,fp) != hf->size)
  986.               yyerror("Not enough heightfield data.");
  987.           for (j = 0; j < hf->size; j++) {
  988.               val = hf->data[i][j];
  989. ***************
  990. *** 145,151 ****
  991.                   hf->minz = val;
  992.           }
  993.       }
  994. !     fclose(fp);
  995.       /*
  996.        * Allocate levels of grid.  hf->levels = log base BestSize of hf->size
  997.        */
  998. --- 148,154 ----
  999.                   hf->minz = val;
  1000.           }
  1001.       }
  1002. !     (void)fclose(fp);
  1003.       /*
  1004.        * Allocate levels of grid.  hf->levels = log base BestSize of hf->size
  1005.        */
  1006. ***************
  1007. *** 154,160 ****
  1008.               ;
  1009.       hf->levels++;
  1010.       hf->qsize = CACHESIZE;
  1011. !     hf->q = (hfTri **)Calloc(hf->qsize, sizeof(hfTri *));
  1012.       hf->qtail = 0;
  1013.   
  1014.       hf->lsize = (int *)share_malloc(hf->levels * sizeof(int));
  1015. --- 157,163 ----
  1016.               ;
  1017.       hf->levels++;
  1018.       hf->qsize = CACHESIZE;
  1019. !     hf->q = (hfTri **)Calloc((unsigned)hf->qsize, sizeof(hfTri *));
  1020.       hf->qtail = 0;
  1021.   
  1022.       hf->lsize = (int *)share_malloc(hf->levels * sizeof(int));
  1023. ***************
  1024. *** 563,568 ****
  1025. --- 566,572 ----
  1026.   /*
  1027.    * Compute normal to height field.
  1028.    */
  1029. + /*ARGSUSED*/
  1030.   nrmhf(pos, prim, nrm)
  1031.   Vector *pos, *nrm;
  1032.   Primitive *prim;
  1033.  
  1034. Index: src/input_yacc.y
  1035. *** src/input_yacc.y.old    Tue Apr 10 16:48:20 1990
  1036. --- src/input_yacc.y    Tue Apr 10 16:48:23 1990
  1037. ***************
  1038. *** 17,23 ****
  1039.   /* name of the person performing the modification, the date of modification,*/
  1040.   /* and the reason for such modification.                    */
  1041.   /*                                        */
  1042. ! /* $Id: input_yacc.y,v 3.0.1.8 90/03/07 21:26:42 craig Exp $ */
  1043.   %{
  1044.   #include <stdio.h>
  1045.   #include "constants.h"
  1046. --- 17,23 ----
  1047.   /* name of the person performing the modification, the date of modification,*/
  1048.   /* and the reason for such modification.                    */
  1049.   /*                                        */
  1050. ! /* $Id: input_yacc.y,v 3.0.1.9 90/04/04 18:58:50 craig Exp $ */
  1051.   %{
  1052.   #include <stdio.h>
  1053.   #include "constants.h"
  1054. ***************
  1055. *** 196,203 ****
  1056.           {
  1057.               CurObj->data->texture = CurText;
  1058.               CurText = (Texture *)0;
  1059.               free((char *)CurObj);
  1060. !             CurObj = CurObj->next;
  1061.           }
  1062.           ;
  1063.   Objectdef    : Startdef Objdefs tENDDEF
  1064. --- 196,204 ----
  1065.           {
  1066.               CurObj->data->texture = CurText;
  1067.               CurText = (Texture *)0;
  1068. +             ListTmp = CurObj->next;
  1069.               free((char *)CurObj);
  1070. !             CurObj = ListTmp;
  1071.           }
  1072.           ;
  1073.   Objectdef    : Startdef Objdefs tENDDEF
  1074.  
  1075. Index: src/main.c
  1076. *** src/main.c.old    Tue Apr 10 16:48:29 1990
  1077. --- src/main.c    Tue Apr 10 16:48:30 1990
  1078. ***************
  1079. *** 1,5 ****
  1080.   char rcsid[] =
  1081. !     "$Id: main.c,v 3.0.1.3 90/03/07 21:33:09 craig Exp $";
  1082.   /*
  1083.    * main.c
  1084.    *
  1085. --- 1,5 ----
  1086.   char rcsid[] =
  1087. !     "$Id: main.c,v 3.0.1.4 90/04/04 14:50:44 craig Exp $";
  1088.   /*
  1089.    * main.c
  1090.    *
  1091. ***************
  1092. *** 21,26 ****
  1093. --- 21,29 ----
  1094.    * and the reason for such modification.
  1095.    *
  1096.    * $Log:    main.c,v $
  1097. +  * Revision 3.0.1.4  90/04/04  14:50:44  craig
  1098. +  * patch5: Removed check of return value of times().
  1099. +  * 
  1100.    * Revision 3.0.1.3  90/03/07  21:33:09  craig
  1101.    * patch4: Return value of times() is now checked.
  1102.    * patch4: Workers is externed iff LINDA is defined.
  1103. ***************
  1104. *** 179,188 ****
  1105.   double *usertime, *systime;
  1106.   {
  1107.       struct tms time;
  1108. -     long times();
  1109.   
  1110. !     if (times(&time) < 0)
  1111. !         RSwarning("Cannot get cpu time?!\n");
  1112.       *usertime = (double)time.tms_utime / (double)HZ;
  1113.       *systime = (double)time.tms_stime / (double)HZ;
  1114.   }
  1115. --- 182,189 ----
  1116.   double *usertime, *systime;
  1117.   {
  1118.       struct tms time;
  1119.   
  1120. !     (void)times(&time);
  1121.       *usertime = (double)time.tms_utime / (double)HZ;
  1122.       *systime = (double)time.tms_stime / (double)HZ;
  1123.   }
  1124.  
  1125. Index: src/memory.c
  1126. *** src/memory.c.old    Tue Apr 10 16:48:33 1990
  1127. --- src/memory.c    Tue Apr 10 16:48:33 1990
  1128. ***************
  1129. *** 18,26 ****
  1130.    * name of the person performing the modification, the date of modification,
  1131.    * and the reason for such modification.
  1132.    *
  1133. !  * $Id: memory.c,v 3.0.1.1 89/12/06 16:33:59 craig Exp $
  1134.    *
  1135.    * $Log:    memory.c,v $
  1136.    * Revision 3.0.1.1  89/12/06  16:33:59  craig
  1137.    * patch2: Added calls to new error/warning routines.
  1138.    * 
  1139. --- 18,29 ----
  1140.    * name of the person performing the modification, the date of modification,
  1141.    * and the reason for such modification.
  1142.    *
  1143. !  * $Id: memory.c,v 3.0.1.2 90/04/04 18:56:08 craig Exp $
  1144.    *
  1145.    * $Log:    memory.c,v $
  1146. +  * Revision 3.0.1.2  90/04/04  18:56:08  craig
  1147. +  * patch5: Lint removal.
  1148. +  * 
  1149.    * Revision 3.0.1.1  89/12/06  16:33:59  craig
  1150.    * patch2: Added calls to new error/warning routines.
  1151.    * 
  1152. ***************
  1153. *** 29,34 ****
  1154. --- 32,40 ----
  1155.    * 
  1156.    */
  1157.   #include <stdio.h>
  1158. + #ifdef SYSV
  1159. + #include <memory.h>
  1160. + #endif
  1161.   #include "typedefs.h"
  1162.   #include "funcdefs.h"
  1163.   
  1164. ***************
  1165. *** 56,64 ****
  1166.   
  1167.       res = Malloc(nelem*elen);
  1168.   #ifdef SYSV
  1169. !     memset(res, (char)0, nelem*elen);
  1170.   #else
  1171. !     bzero(res, nelem*elen);
  1172.   #endif
  1173.       return res;
  1174.   }
  1175. --- 62,70 ----
  1176.   
  1177.       res = Malloc(nelem*elen);
  1178.   #ifdef SYSV
  1179. !     (void)memset(res, (char)0, (int)(nelem*elen));
  1180.   #else
  1181. !     bzero(res, (int)nelem*elen);
  1182.   #endif
  1183.       return res;
  1184.   }
  1185.  
  1186. Index: src/noise.c
  1187. *** src/noise.c.old    Tue Apr 10 16:48:36 1990
  1188. --- src/noise.c    Tue Apr 10 16:48:36 1990
  1189. ***************
  1190. *** 18,26 ****
  1191.    * name of the person performing the modification, the date of modification,
  1192.    * and the reason for such modification.
  1193.    *
  1194. !  * $Id: noise.c,v 3.0 89/10/27 02:05:57 craig Exp $
  1195.    *
  1196.    * $Log:    noise.c,v $
  1197.    * Revision 3.0  89/10/27  02:05:57  craig
  1198.    * Baseline for first official release.
  1199.    * 
  1200. --- 18,29 ----
  1201.    * name of the person performing the modification, the date of modification,
  1202.    * and the reason for such modification.
  1203.    *
  1204. !  * $Id: noise.c,v 3.0.1.1 90/04/04 18:59:22 craig Exp $
  1205.    *
  1206.    * $Log:    noise.c,v $
  1207. +  * Revision 3.0.1.1  90/04/04  18:59:22  craig
  1208. +  * patch5: Lint removal.
  1209. +  * 
  1210.    * Revision 3.0  89/10/27  02:05:57  craig
  1211.    * Baseline for first official release.
  1212.    * 
  1213. ***************
  1214. *** 100,106 ****
  1215.       (void)srandom(0);
  1216.   #endif
  1217.   
  1218. !     hashTable = (short int *) malloc(4096*sizeof(short int));
  1219.       for (i = 0; i < 4096; i++)
  1220.           hashTable[i] = i;
  1221.       for (i = 4095; i > 0; i--) {
  1222. --- 103,109 ----
  1223.       (void)srandom(0);
  1224.   #endif
  1225.   
  1226. !     hashTable = (short *) Malloc(4096*sizeof(short int));
  1227.       for (i = 0; i < 4096; i++)
  1228.           hashTable[i] = i;
  1229.       for (i = 4095; i > 0; i--) {
  1230.  
  1231. Index: src/outputp.c
  1232. *** src/outputp.c.old    Tue Apr 10 16:48:40 1990
  1233. --- src/outputp.c    Tue Apr 10 16:48:40 1990
  1234. ***************
  1235. *** 18,26 ****
  1236.    * name of the person performing the modification, the date of modification,
  1237.    * and the reason for such modification.
  1238.    *
  1239. !  * $Id: outputp.c,v 3.0.1.2 90/02/12 16:09:05 craig Exp $
  1240.    *
  1241.    * $Log:    outputp.c,v $
  1242.    * Revision 3.0.1.2  90/02/12  16:09:05  craig
  1243.    * patch4: Fixed bug with duplicate headers in mtv files when Appending.
  1244.    * 
  1245. --- 18,29 ----
  1246.    * name of the person performing the modification, the date of modification,
  1247.    * and the reason for such modification.
  1248.    *
  1249. !  * $Id: outputp.c,v 3.0.1.3 90/04/04 14:53:04 craig Exp $
  1250.    *
  1251.    * $Log:    outputp.c,v $
  1252. +  * Revision 3.0.1.3  90/04/04  14:53:04  craig
  1253. +  * patch5: Lint removal.
  1254. +  * 
  1255.    * Revision 3.0.1.2  90/02/12  16:09:05  craig
  1256.    * patch4: Fixed bug with duplicate headers in mtv files when Appending.
  1257.    * 
  1258. ***************
  1259. *** 101,107 ****
  1260.       buffer[1]=(unsigned char *)Malloc(sizeof(unsigned char)*(unsigned)Xres);
  1261.       buffer[2]=(unsigned char *)Malloc(sizeof(unsigned char)*(unsigned)Xres);
  1262.       if (Appending)
  1263. !         fseek(imgfile, 0L, 2);        /* Go to end of file */
  1264.   }
  1265.   
  1266.   
  1267. --- 104,110 ----
  1268.       buffer[1]=(unsigned char *)Malloc(sizeof(unsigned char)*(unsigned)Xres);
  1269.       buffer[2]=(unsigned char *)Malloc(sizeof(unsigned char)*(unsigned)Xres);
  1270.       if (Appending)
  1271. !         (void)fseek(imgfile, 0L, 2);    /* Go to end of file */
  1272.   }
  1273.   
  1274.   
  1275. ***************
  1276. *** 134,140 ****
  1277.   endpic()
  1278.   {
  1279.       sv_puteof(&sv_globals);
  1280. !     fclose(imgfile);
  1281.   }
  1282.   
  1283.   #else /* NORLE */
  1284. --- 137,143 ----
  1285.   endpic()
  1286.   {
  1287.       sv_puteof(&sv_globals);
  1288. !     (void)fclose(imgfile);
  1289.   }
  1290.   
  1291.   #else /* NORLE */
  1292. ***************
  1293. *** 166,174 ****
  1294.       extern int Xres;
  1295.   
  1296.       for (i = 0; i < Xres; i++) {
  1297. !         fputc(correct(255.*buf[i].r), imgfile);
  1298. !         fputc(correct(255.*buf[i].g), imgfile);
  1299. !         fputc(correct(255.*buf[i].b), imgfile);
  1300.       }
  1301.       fflush(imgfile);
  1302.   }
  1303. --- 169,177 ----
  1304.       extern int Xres;
  1305.   
  1306.       for (i = 0; i < Xres; i++) {
  1307. !         (void)fputc(correct(255.*buf[i].r), imgfile);
  1308. !         (void)fputc(correct(255.*buf[i].g), imgfile);
  1309. !         (void)fputc(correct(255.*buf[i].b), imgfile);
  1310.       }
  1311.       fflush(imgfile);
  1312.   }
  1313. ***************
  1314. *** 175,181 ****
  1315.   
  1316.   endpic()
  1317.   {
  1318. !     fclose(imgfile);
  1319.   }
  1320.   #endif /* NORLE */
  1321.   
  1322. --- 178,184 ----
  1323.   
  1324.   endpic()
  1325.   {
  1326. !     (void)fclose(imgfile);
  1327.   }
  1328.   #endif /* NORLE */
  1329.   
  1330.  
  1331. Index: src/plane.c
  1332. *** src/plane.c.old    Tue Apr 10 16:48:43 1990
  1333. --- src/plane.c    Tue Apr 10 16:48:44 1990
  1334. ***************
  1335. *** 18,26 ****
  1336.    * name of the person performing the modification, the date of modification,
  1337.    * and the reason for such modification.
  1338.    *
  1339. !  * $Id: plane.c,v 3.0.1.1 89/12/06 16:34:03 craig Exp $
  1340.    *
  1341.    * $Log:    plane.c,v $
  1342.    * Revision 3.0.1.1  89/12/06  16:34:03  craig
  1343.    * patch2: Added calls to new error/warning routines.
  1344.    * 
  1345. --- 18,29 ----
  1346.    * name of the person performing the modification, the date of modification,
  1347.    * and the reason for such modification.
  1348.    *
  1349. !  * $Id: plane.c,v 3.0.1.2 90/04/04 19:02:42 craig Exp $
  1350.    *
  1351.    * $Log:    plane.c,v $
  1352. +  * Revision 3.0.1.2  90/04/04  19:02:42  craig
  1353. +  * patch5: Test for edge-on intersection more robust.
  1354. +  * 
  1355.    * Revision 3.0.1.1  89/12/06  16:34:03  craig
  1356.    * patch2: Added calls to new error/warning routines.
  1357.    * 
  1358. ***************
  1359. *** 77,83 ****
  1360.       plane = obj->objpnt.p_plane;
  1361.   
  1362.       denom = dotp(&plane->norm, ray);
  1363. !     if (denom == 0.)
  1364.           return 0.;
  1365.       dist = (plane->d - dotp(&plane->norm, pos)) / denom;
  1366.       return (dist > FAR_AWAY ? 0. : dist);
  1367. --- 80,86 ----
  1368.       plane = obj->objpnt.p_plane;
  1369.   
  1370.       denom = dotp(&plane->norm, ray);
  1371. !     if (equal(denom, 0.))
  1372.           return 0.;
  1373.       dist = (plane->d - dotp(&plane->norm, pos)) / denom;
  1374.       return (dist > FAR_AWAY ? 0. : dist);
  1375.  
  1376. Index: src/poly.c
  1377. *** src/poly.c.old    Tue Apr 10 16:48:47 1990
  1378. --- src/poly.c    Tue Apr 10 16:48:47 1990
  1379. ***************
  1380. *** 18,26 ****
  1381.    * name of the person performing the modification, the date of modification,
  1382.    * and the reason for such modification.
  1383.    *
  1384. !  * $Id: poly.c,v 3.0.1.1 89/12/06 16:33:38 craig Exp $
  1385.    *
  1386.    * $Log:    poly.c,v $
  1387.    * Revision 3.0.1.1  89/12/06  16:33:38  craig
  1388.    * patch2: Added calls to new error/warning routines.
  1389.    * 
  1390. --- 18,29 ----
  1391.    * name of the person performing the modification, the date of modification,
  1392.    * and the reason for such modification.
  1393.    *
  1394. !  * $Id: poly.c,v 3.0.1.2 90/04/04 19:01:19 craig Exp $
  1395.    *
  1396.    * $Log:    poly.c,v $
  1397. +  * Revision 3.0.1.2  90/04/04  19:01:19  craig
  1398. +  * patch5: Corrected free()-related problems.
  1399. +  * 
  1400.    * Revision 3.0.1.1  89/12/06  16:33:38  craig
  1401.    * patch2: Added calls to new error/warning routines.
  1402.    * 
  1403. ***************
  1404. *** 49,55 ****
  1405.       Object *newobj;
  1406.       double indexval;
  1407.       Vector edge1, edge2, anorm;
  1408. !     PointList *cur;
  1409.       int i;
  1410.       extern int TrashBadPoly;
  1411.   
  1412. --- 52,58 ----
  1413.       Object *newobj;
  1414.       double indexval;
  1415.       Vector edge1, edge2, anorm;
  1416. !     PointList *cur, *pltmp;
  1417.       int i;
  1418.       extern int TrashBadPoly;
  1419.   
  1420. ***************
  1421. *** 71,81 ****
  1422.        * to worry about doing so.
  1423.        */
  1424.       i = npoints -1;
  1425. !     for(cur = plist;cur;cur = cur->next) {
  1426.           poly->points[i--] = cur->vec;
  1427.           free((char *)cur);
  1428.       }
  1429. -     free((char *)plist);
  1430.   
  1431.       /*
  1432.        * Find normal to polygon.  Check all edges before giving
  1433. --- 74,84 ----
  1434.        * to worry about doing so.
  1435.        */
  1436.       i = npoints -1;
  1437. !     for(cur = plist; cur != (PointList *)0; cur = pltmp) {
  1438.           poly->points[i--] = cur->vec;
  1439. +         pltmp = cur->next;
  1440.           free((char *)cur);
  1441.       }
  1442.   
  1443.       /*
  1444.        * Find normal to polygon.  Check all edges before giving
  1445.  
  1446. Index: src/primobj.h
  1447. *** src/primobj.h.old    Tue Apr 10 16:48:50 1990
  1448. --- src/primobj.h    Tue Apr 10 16:48:50 1990
  1449. ***************
  1450. *** 18,26 ****
  1451.    * name of the person performing the modification, the date of modification,
  1452.    * and the reason for such modification.
  1453.    *
  1454. !  * $Id: primobj.h,v 3.0 89/10/27 02:06:00 craig Exp $
  1455.    *
  1456.    * $Log:    primobj.h,v $
  1457.    * Revision 3.0  89/10/27  02:06:00  craig
  1458.    * Baseline for first official release.
  1459.    * 
  1460. --- 18,29 ----
  1461.    * name of the person performing the modification, the date of modification,
  1462.    * and the reason for such modification.
  1463.    *
  1464. !  * $Id: primobj.h,v 3.0.1.1 90/04/09 13:39:59 craig Exp $
  1465.    *
  1466.    * $Log:    primobj.h,v $
  1467. +  * Revision 3.0.1.1  90/04/09  13:39:59  craig
  1468. +  * patch5: New Cylinder and Cone types.
  1469. +  * 
  1470.    * Revision 3.0  89/10/27  02:06:00  craig
  1471.    * Baseline for first official release.
  1472.    * 
  1473. ***************
  1474. *** 30,41 ****
  1475.    * Cylinder
  1476.    */
  1477.   typedef struct cylinder    {
  1478. !     double rad, len;        /* Radius and length. */
  1479.   } Cylinder;
  1480.   
  1481.   typedef struct cone {
  1482. !     double start_pos, end_pos, apex_rad;    /* End cap Z positions */
  1483. !     double tantheta;            /* "slope" */
  1484.   } Cone;
  1485.   
  1486.   /*
  1487. --- 33,44 ----
  1488.    * Cylinder
  1489.    */
  1490.   typedef struct cylinder    {
  1491. !     Trans trans;
  1492.   } Cylinder;
  1493.   
  1494.   typedef struct cone {
  1495. !     double start_dist;
  1496. !     Trans trans;
  1497.   } Cone;
  1498.   
  1499.   /*
  1500.  
  1501. Index: src/setup.c
  1502. *** src/setup.c.old    Tue Apr 10 16:48:55 1990
  1503. --- src/setup.c    Tue Apr 10 16:48:56 1990
  1504. ***************
  1505. *** 18,26 ****
  1506.    * name of the person performing the modification, the date of modification,
  1507.    * and the reason for such modification.
  1508.    *
  1509. !  * $Id: setup.c,v 3.0.1.2 90/02/12 13:26:49 craig Exp $
  1510.    *
  1511.    * $Log:    setup.c,v $
  1512.    * Revision 3.0.1.2  90/02/12  13:26:49  craig
  1513.    * patch4: Added calls to new error/warning routines.
  1514.    * patch4: Attempting to render empty scene now results in an error message.
  1515. --- 18,29 ----
  1516.    * name of the person performing the modification, the date of modification,
  1517.    * and the reason for such modification.
  1518.    *
  1519. !  * $Id: setup.c,v 3.0.1.3 90/04/04 19:00:19 craig Exp $
  1520.    *
  1521.    * $Log:    setup.c,v $
  1522. +  * Revision 3.0.1.3  90/04/04  19:00:19  craig
  1523. +  * patch5: Lint removal.
  1524. +  * 
  1525.    * Revision 3.0.1.2  90/02/12  13:26:49  craig
  1526.    * patch4: Added calls to new error/warning routines.
  1527.    * patch4: Attempting to render empty scene now results in an error message.
  1528. ***************
  1529. *** 178,184 ****
  1530.        * transformation.
  1531.        */
  1532.       for (i = 0; i < nlight; i++) {
  1533. !         light[i].cache = (Primitive **)Calloc(maxlevel+1,
  1534.                       sizeof(Primitive *));
  1535.           light[i].trans = (TransInfo *)Malloc((maxlevel+1)*
  1536.                       sizeof(TransInfo));
  1537. --- 181,187 ----
  1538.        * transformation.
  1539.        */
  1540.       for (i = 0; i < nlight; i++) {
  1541. !         light[i].cache = (Primitive **)Calloc((unsigned)maxlevel+1,
  1542.                       sizeof(Primitive *));
  1543.           light[i].trans = (TransInfo *)Malloc((maxlevel+1)*
  1544.                       sizeof(TransInfo));
  1545.  
  1546. Index: src/shade.c
  1547. *** src/shade.c.old    Tue Apr 10 16:49:00 1990
  1548. --- src/shade.c    Tue Apr 10 16:49:02 1990
  1549. ***************
  1550. *** 18,26 ****
  1551.    * name of the person performing the modification, the date of modification,
  1552.    * and the reason for such modification.
  1553.    *
  1554. !  * $Id: shade.c,v 3.0.1.2 90/03/07 21:31:21 craig Exp $
  1555.    *
  1556.    * $Log:    shade.c,v $
  1557.    * Revision 3.0.1.2  90/03/07  21:31:21  craig
  1558.    * patch4: Replaced call to LightCoordSys() with call to VectCoordSys().
  1559.    * 
  1560. --- 18,29 ----
  1561.    * name of the person performing the modification, the date of modification,
  1562.    * and the reason for such modification.
  1563.    *
  1564. !  * $Id: shade.c,v 3.0.1.3 90/04/04 19:04:19 craig Exp $
  1565.    *
  1566.    * $Log:    shade.c,v $
  1567. +  * Revision 3.0.1.3  90/04/04  19:04:19  craig
  1568. +  * patch5: Corrected free()-related problems.
  1569. +  * 
  1570.    * Revision 3.0.1.2  90/03/07  21:31:21  craig
  1571.    * patch4: Replaced call to LightCoordSys() with call to VectCoordSys().
  1572.    * 
  1573. ***************
  1574. *** 180,186 ****
  1575.                    * of two refracting surfaces. Skip it.
  1576.                    */
  1577.                   return;
  1578. -             free((char *)NewRay.media);
  1579.               NewRay.media = NewRay.media->next;
  1580.               if (refract(&NewRay.dir, surf->kref,
  1581.                   NewRay.media ? NewRay.media->surf->kref :
  1582. --- 183,188 ----
  1583.  
  1584. Index: src/texture.c
  1585. *** src/texture.c.old    Tue Apr 10 16:49:06 1990
  1586. --- src/texture.c    Tue Apr 10 16:49:08 1990
  1587. ***************
  1588. *** 18,26 ****
  1589.    * name of the person performing the modification, the date of modification,
  1590.    * and the reason for such modification.
  1591.    *
  1592. !  * $Id: texture.c,v 3.0.1.3 90/03/07 21:30:22 craig Exp $
  1593.    *
  1594.    * $Log:    texture.c,v $
  1595.    * Revision 3.0.1.3  90/03/07  21:30:22  craig
  1596.    * patch4: Fixed thresholding bug in fBmText().
  1597.    * patch4: Colormap is now set to NULL in new_texture().
  1598. --- 18,29 ----
  1599.    * name of the person performing the modification, the date of modification,
  1600.    * and the reason for such modification.
  1601.    *
  1602. !  * $Id: texture.c,v 3.0.1.4 90/04/04 17:22:39 craig Exp $
  1603.    *
  1604.    * $Log:    texture.c,v $
  1605. +  * Revision 3.0.1.4  90/04/04  17:22:39  craig
  1606. +  * patch5: Added checking for invalid colormap entries.
  1607. +  * 
  1608.    * Revision 3.0.1.3  90/03/07  21:30:22  craig
  1609.    * patch4: Fixed thresholding bug in fBmText().
  1610.    * patch4: Colormap is now set to NULL in new_texture().
  1611. ***************
  1612. *** 114,119 ****
  1613. --- 117,123 ----
  1614.       return text;
  1615.   }
  1616.   
  1617. + /*ARGSUSED*/
  1618.   WoodText(text, pos, norm, surf)
  1619.   Texture *text;
  1620.   Vector *pos, *norm;
  1621. ***************
  1622. *** 120,126 ****
  1623.   Surface *surf;
  1624.   {
  1625.       double red, grn, blu;
  1626. !     double chaos, midBrown, brownLayer, greenLayer;
  1627.       double perturb, brownPerturb, greenPerturb, grnPerturb;
  1628.       double t;
  1629.   
  1630. --- 124,130 ----
  1631.   Surface *surf;
  1632.   {
  1633.       double red, grn, blu;
  1634. !     double chaos, brownLayer, greenLayer;
  1635.       double perturb, brownPerturb, greenPerturb, grnPerturb;
  1636.       double t;
  1637.   
  1638. ***************
  1639. *** 265,270 ****
  1640. --- 269,275 ----
  1641.   /*
  1642.    * Apply "blotch" texture.
  1643.    */
  1644. + /*ARGSUSED*/
  1645.   BlotchText(text, pos, norm, surf)
  1646.   Texture *text;
  1647.   Vector *pos, *norm;
  1648. ***************
  1649. *** 282,287 ****
  1650. --- 287,293 ----
  1651.       }
  1652.   }
  1653.   
  1654. + /*ARGSUSED*/
  1655.   fBmText(text, pos, norm, surf)
  1656.   Texture *text;
  1657.   Vector *pos, *norm;
  1658. ***************
  1659. *** 314,319 ****
  1660. --- 320,326 ----
  1661.   /*
  1662.    * Apply a "checker" texture.
  1663.    */
  1664. + /*ARGSUSED*/
  1665.   CheckerText(text, pos, norm, surf)
  1666.   Texture *text;
  1667.   Vector *pos, *norm;
  1668. ***************
  1669. *** 330,335 ****
  1670. --- 337,343 ----
  1671.       /* else surface stays the same. */
  1672.   }
  1673.   
  1674. + /*ARGSUSED*/
  1675.   fBmBumpText(text, pos, norm, surf)
  1676.   Texture *text;
  1677.   Vector *pos, *norm;
  1678. ***************
  1679. *** 348,353 ****
  1680. --- 356,362 ----
  1681.   /*
  1682.    * Apply a "bump" texture.
  1683.    */
  1684. + /*ARGSUSED*/
  1685.   BumpText(text, pos, norm, surf)
  1686.   Texture *text;
  1687.   Vector *pos, *norm;
  1688. ***************
  1689. *** 362,367 ****
  1690. --- 371,377 ----
  1691.       (void)normalize(norm);
  1692.   }
  1693.   
  1694. + /*ARGSUSED*/
  1695.   MarbleText(text, pos, norm, surf)
  1696.   Texture *text;
  1697.   Vector *pos, *norm;
  1698. ***************
  1699. *** 385,390 ****
  1700. --- 395,401 ----
  1701.       }
  1702.   }
  1703.   
  1704. + /*ARGSUSED*/
  1705.   GlossText(text, pos, norm, surf)
  1706.   Texture *text;
  1707.   Vector *pos, *norm;
  1708. ***************
  1709. *** 420,426 ****
  1710.       FILE *fp;
  1711.       Color *map;
  1712.       char buf[BUFSIZ];
  1713. !     int i;
  1714.   
  1715.       fp = fopen(filename, "r");
  1716.       if (fp == (FILE *)NULL)
  1717. --- 431,437 ----
  1718.       FILE *fp;
  1719.       Color *map;
  1720.       char buf[BUFSIZ];
  1721. !     int i, num;
  1722.   
  1723.       fp = fopen(filename, "r");
  1724.       if (fp == (FILE *)NULL)
  1725. ***************
  1726. *** 429,436 ****
  1727.       map = (Color *)Calloc(256, sizeof(Color));
  1728.   
  1729.       for (i = 0; fgets(buf,BUFSIZ,fp) != NULL && i < 256; i++) {
  1730. !         sscanf(buf,"%lf %lf %lf",&map[i].r, &map[i].g, &map[i].b);
  1731.           ScaleColor(1. / 255., map[i], &map[i]);
  1732.       }
  1733.       return map;
  1734.   }
  1735. --- 440,452 ----
  1736.       map = (Color *)Calloc(256, sizeof(Color));
  1737.   
  1738.       for (i = 0; fgets(buf,BUFSIZ,fp) != NULL && i < 256; i++) {
  1739. !         num = sscanf(buf,"%lf %lf %lf",&map[i].r, &map[i].g, &map[i].b);
  1740. !         if (num != 3) {
  1741. !             RSerror("%s, line %d:  Bad color specification.\n",
  1742. !                 filename, i+1);
  1743. !         }
  1744.           ScaleColor(1. / 255., map[i], &map[i]);
  1745.       }
  1746. +     (void)fclose(fp);
  1747.       return map;
  1748.   }
  1749.  
  1750. Index: src/triangle.c
  1751. *** src/triangle.c.old    Tue Apr 10 16:49:11 1990
  1752. --- src/triangle.c    Tue Apr 10 16:49:13 1990
  1753. ***************
  1754. *** 18,26 ****
  1755.    * name of the person performing the modification, the date of modification,
  1756.    * and the reason for such modification.
  1757.    *
  1758. !  * $Id: triangle.c,v 3.0.1.1 89/12/06 16:33:32 craig Exp $
  1759.    *
  1760.    * $Log:    triangle.c,v $
  1761.    * Revision 3.0.1.1  89/12/06  16:33:32  craig
  1762.    * patch2: Added calls to new error/warning routines.
  1763.    * 
  1764. --- 18,30 ----
  1765.    * name of the person performing the modification, the date of modification,
  1766.    * and the reason for such modification.
  1767.    *
  1768. !  * $Id: triangle.c,v 3.0.1.2 90/04/04 18:55:45 craig Exp $
  1769.    *
  1770.    * $Log:    triangle.c,v $
  1771. +  * Revision 3.0.1.2  90/04/04  18:55:45  craig
  1772. +  * patch5: Test for edge-on intersection is more robust.
  1773. +  * patch5: Lint removal.
  1774. +  * 
  1775.    * Revision 3.0.1.1  89/12/06  16:33:32  craig
  1776.    * patch2: Added calls to new error/warning routines.
  1777.    * 
  1778. ***************
  1779. *** 146,152 ****
  1780.        * Plane intersection.
  1781.        */
  1782.       k = dotp(&triangle->nrm, ray);
  1783. !     if (k == 0.)
  1784.           return 0.;
  1785.       s = (triangle->d - dotp(&triangle->nrm, pos)) / k;
  1786.       if (s <= 0.)
  1787. --- 150,156 ----
  1788.        * Plane intersection.
  1789.        */
  1790.       k = dotp(&triangle->nrm, ray);
  1791. !     if (equal(k, 0.))
  1792.           return 0.;
  1793.       s = (triangle->d - dotp(&triangle->nrm, pos)) / k;
  1794.       if (s <= 0.)
  1795. ***************
  1796. *** 211,216 ****
  1797. --- 215,221 ----
  1798.       return s;
  1799.   }
  1800.   
  1801. + /*ARGSUSED*/
  1802.   nrmtri(pos, obj, nrm)
  1803.   Vector           *pos, *nrm;
  1804.   Primitive       *obj;
  1805.  
  1806. Index: src/viewing.c
  1807. *** src/viewing.c.old    Tue Apr 10 16:49:17 1990
  1808. --- src/viewing.c    Tue Apr 10 16:49:17 1990
  1809. ***************
  1810. *** 18,26 ****
  1811.    * name of the person performing the modification, the date of modification,
  1812.    * and the reason for such modification.
  1813.    *
  1814. !  * $Id: viewing.c,v 3.0.1.3 89/12/02 14:38:37 craig Exp $
  1815.    *
  1816.    * $Log:    viewing.c,v $
  1817.    * Revision 3.0.1.3  89/12/02  14:38:37  craig
  1818.    * patch2: Added depth of field code, courtesy of Rodney G. Bogart.
  1819.    * 
  1820. --- 18,29 ----
  1821.    * name of the person performing the modification, the date of modification,
  1822.    * and the reason for such modification.
  1823.    *
  1824. !  * $Id: viewing.c,v 3.0.1.4 90/04/04 18:59:54 craig Exp $
  1825.    *
  1826.    * $Log:    viewing.c,v $
  1827. +  * Revision 3.0.1.4  90/04/04  18:59:54  craig
  1828. +  * patch5: Lint removal.
  1829. +  * 
  1830.    * Revision 3.0.1.3  89/12/02  14:38:37  craig
  1831.    * patch2: Added depth of field code, courtesy of Rodney G. Bogart.
  1832.    * 
  1833. ***************
  1834. *** 126,132 ****
  1835.       vecadd(aperture_inc, eyep, &(ray->pos));
  1836.       scalar_prod(focaldist, ray->dir, &(ray->dir));
  1837.       vecsub(ray->dir, aperture_inc, &(ray->dir));
  1838. !     normalize(&ray->dir);
  1839.   }
  1840.   
  1841.   /*
  1842. --- 129,135 ----
  1843.       vecadd(aperture_inc, eyep, &(ray->pos));
  1844.       scalar_prod(focaldist, ray->dir, &(ray->dir));
  1845.       vecsub(ray->dir, aperture_inc, &(ray->dir));
  1846. !     (void)normalize(&ray->dir);
  1847.   }
  1848.   
  1849.   /*
  1850.  
  1851. Index: src/voxels.c
  1852. *** src/voxels.c.old    Tue Apr 10 16:49:21 1990
  1853. --- src/voxels.c    Tue Apr 10 16:49:21 1990
  1854. ***************
  1855. *** 18,26 ****
  1856.    * name of the person performing the modification, the date of modification,
  1857.    * and the reason for such modification.
  1858.    *
  1859. !  * $Id: voxels.c,v 3.0.1.2 90/02/12 13:19:50 craig Exp $
  1860.    *
  1861.    * $Log:    voxels.c,v $
  1862.    * Revision 3.0.1.2  90/02/12  13:19:50  craig
  1863.    * patch4: Added reporting of total number of primitives.
  1864.    * 
  1865. --- 18,29 ----
  1866.    * name of the person performing the modification, the date of modification,
  1867.    * and the reason for such modification.
  1868.    *
  1869. !  * $Id: voxels.c,v 3.0.1.3 90/04/04 19:05:23 craig Exp $
  1870.    *
  1871.    * $Log:    voxels.c,v $
  1872. +  * Revision 3.0.1.3  90/04/04  19:05:23  craig
  1873. +  * patch5: Changes to reflect that free() may be destructive.
  1874. +  * 
  1875.    * Revision 3.0.1.2  90/02/12  13:19:50  craig
  1876.    * patch4: Added reporting of total number of primitives.
  1877.    * 
  1878. ***************
  1879. *** 152,158 ****
  1880.   {
  1881.       Grid *grid;
  1882.       Object *otmp;
  1883. !     ObjList *ltmp;
  1884.       int x, y, i;
  1885.       extern ObjList *find_bounds();
  1886.   
  1887. --- 155,161 ----
  1888.   {
  1889.       Grid *grid;
  1890.       Object *otmp;
  1891. !     ObjList *ltmp, *nltmp;
  1892.       int x, y, i;
  1893.       extern ObjList *find_bounds();
  1894.   
  1895. ***************
  1896. *** 190,199 ****
  1897.       /*
  1898.        * obj->data now holds a linked list of bounded objects.
  1899.        */
  1900. !     for(ltmp = (ObjList *)obj->data; ltmp; ltmp = ltmp->next) {
  1901.           otmp = ltmp->data;
  1902.           engrid(otmp, grid);
  1903. !         free(ltmp);
  1904.       }
  1905.       obj->type = GRID;
  1906.       obj->data = (char *)grid;
  1907. --- 193,203 ----
  1908.       /*
  1909.        * obj->data now holds a linked list of bounded objects.
  1910.        */
  1911. !     for(ltmp = (ObjList *)obj->data; ltmp != (ObjList *)0; ltmp = nltmp) {
  1912.           otmp = ltmp->data;
  1913.           engrid(otmp, grid);
  1914. !         nltmp = ltmp->next;
  1915. !         free((char *)ltmp);
  1916.       }
  1917.       obj->type = GRID;
  1918.       obj->data = (char *)grid;
  1919.  
  1920. Index: doc/primitive.ms
  1921. *** doc/primitive.ms.old    Tue Apr 10 17:22:56 1990
  1922. --- doc/primitive.ms    Tue Apr 10 17:22:57 1990
  1923. ***************
  1924. *** 3,11 ****
  1925.   .\" Brief tutorial on adding new primitives to rayshade.
  1926.   .\" Craig Kolb 10/89
  1927.   .\"
  1928. ! .\" $Id: primitive.ms,v 3.0 89/10/24 18:07:30 craig Exp $
  1929.   .\"
  1930.   .\" $Log:    primitive.ms,v $
  1931.   .\" Revision 3.0  89/10/24  18:07:30  craig
  1932.   .\" Baseline for first official release.
  1933.   .\" 
  1934. --- 3,14 ----
  1935.   .\" Brief tutorial on adding new primitives to rayshade.
  1936.   .\" Craig Kolb 10/89
  1937.   .\"
  1938. ! .\" $Id: primitive.ms,v 3.0.1.1 90/04/10 17:22:36 craig Exp $
  1939.   .\"
  1940.   .\" $Log:    primitive.ms,v $
  1941. + .\" Revision 3.0.1.1  90/04/10  17:22:36  craig
  1942. + .\" patch5: addscaledvec() no longer returns a structure...
  1943. + .\" 
  1944.   .\" Revision 3.0  89/10/24  18:07:30  craig
  1945.   .\" Baseline for first official release.
  1946.   .\" 
  1947. ***************
  1948. *** 505,511 ****
  1949.           /*
  1950.            *  Find difference between point of intersection and center of disc.
  1951.            */
  1952. !         vecsub(addscaledvec(*pos, dist, *ray), disc->center, &hit);
  1953.           /*
  1954.            * If hit point is <= disc->radius from center, we've hit the disc.
  1955.            */
  1956. --- 508,515 ----
  1957.           /*
  1958.            *  Find difference between point of intersection and center of disc.
  1959.            */
  1960. !         addscaledvec(*pos, dist, *ray, &hit);
  1961. !         vecsub(hit, disc->center, &hit);
  1962.           /*
  1963.            * If hit point is <= disc->radius from center, we've hit the disc.
  1964.            */
  1965.