home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 October / usenetsourcesnewsgroupsinfomagicoctober1994disk2.iso / misc / volume23 / rayshade / patch01 next >
Text File  |  1991-10-19  |  44KB  |  1,467 lines
  1. Newsgroups: comp.sources.misc
  2. From: Craig Kolb <rayshade-request@cs.princeton.edu>
  3. Subject:  v23i075:  rayshade - A raytracing package for UNIX, Patch01
  4. Message-ID: <1991Oct19.024232.13698@sparky.imd.sterling.com>
  5. X-Md4-Signature: 6348dbe45469af21d17e7a255c013923
  6. Date: Sat, 19 Oct 1991 02:42:32 GMT
  7. Approved: kent@sparky.imd.sterling.com
  8.  
  9. Submitted-by: Craig Kolb <rayshade-request@cs.princeton.edu>
  10. Posting-number: Volume 23, Issue 75
  11. Archive-name: rayshade/patch01
  12. Environment: UNIX
  13. Patch-To: rayshade: Volume 21, Issue 3-22
  14.  
  15. System: rayshade version 4.0
  16. Patch #: 1
  17. Priority: HIGH
  18.  
  19. Description:
  20.  
  21.     Added "transform" support.
  22.  
  23.     New "window" and "crop" support.  Window uses pixel coordinates
  24.     to define a subwindow, crop uses normalized coordinates and
  25.     crops the subwindow.
  26.  
  27.     Rayview supports spotlights and transparency (given appropriate
  28.     hardware.)
  29.  
  30.     Faster bbox transform, a la Jim Arvo.
  31.  
  32.     Fixed view initialization weirdness that lead to bugs
  33.     when rendering animated stereo pairs.
  34.  
  35.     Don't declare free or exit void if stdlib.h is used.
  36.  
  37.     Fixed problems with coincident bbox faces.
  38.  
  39.     Removed a number of floating point comparisons
  40.  
  41.     Corrected help message.
  42.  
  43.     Potential roundoff problem in triangle dPdU code.
  44.  
  45.     CoordSysTransform did not detect up-Z == -1.
  46.  
  47.     Defined MAKE macro in Makefiles for braindamaged machines.
  48.     Note that "make -n" on such machines will almost certainly do
  49.     the wrong thing.
  50.  
  51.     Fixed various typos.
  52.  
  53.  
  54. Fix:    From rn, say "| patch -p -N -d DIR", where DIR is your rayshade source
  55.     directory.  Outside of rn, say "cd DIR; patch -p -N <thisarticle".
  56.     If you don't have the patch program, apply the following by hand,
  57.     or get patch (version 2.0, latest patchlevel).
  58.  
  59.     After patching:
  60.  
  61.         *** DO NOTHING--INSTALL PATCH #2 FIRST ***
  62.  
  63.     If patch indicates that patchlevel is the wrong version, you may need
  64.     to apply one or more previous patches, or the patch may already
  65.     have been applied.  See the patchlevel.h file to find out what has or
  66.     has not been applied.  In any event, don't continue with the patch.
  67.  
  68.     If you are missing previous patches (hah!), they can be obtained
  69.     through anonymous ftp from weedeater.math.yale.edu.  Questions
  70.     and problems should be addressed to:
  71.  
  72.         rayshade-request@cs.princeton.edu
  73.  
  74.  
  75. Index: patchlevel.h
  76. Prereq: 0
  77. 1c1
  78. < #define PATCHLEVEL 0
  79. ---
  80. > #define PATCHLEVEL 1
  81.  
  82. Index: rayview/glmethods.c
  83. Prereq: 4.0
  84. *** rayview/glmethods.c.old    1991/07/17 17:38:43
  85. --- rayview/glmethods.c    1991/10/10 22:45:27
  86. ***************
  87. *** 14,22 ****
  88.    * There is no warranty or other guarantee of fitness of this software
  89.    * for any purpose.  It is provided solely "as is".
  90.    *
  91. !  * $Id: glmethods.c,v 4.0 91/07/17 17:38:43 kolb Exp Locker: kolb $
  92.    *
  93.    * $Log:    glmethods.c,v $
  94.    * Revision 4.0  91/07/17  17:38:43  kolb
  95.    * Initial version
  96.    * 
  97. --- 14,28 ----
  98.    * There is no warranty or other guarantee of fitness of this software
  99.    * for any purpose.  It is provided solely "as is".
  100.    *
  101. !  * $Id: glmethods.c,v 4.0.1.2 91/10/10 22:45:27 cek Exp $
  102.    *
  103.    * $Log:    glmethods.c,v $
  104. +  * Revision 4.0.1.2  91/10/10  22:45:27  cek
  105. +  * patch1: Added spotlight support.
  106. +  * 
  107. +  * Revision 4.0.1.1  1991/09/29  15:28:36  cek
  108. +  * patch1: Added support for transparency.
  109. +  *
  110.    * Revision 4.0  91/07/17  17:38:43  kolb
  111.    * Initial version
  112.    * 
  113. ***************
  114. *** 117,123 ****
  115.       PlaneMethodRegister(GLPlaneDraw);
  116.       PolygonMethodRegister(GLPolygonDraw);
  117.       SphereMethodRegister(GLSphereDraw);
  118. -     /* Have to write sphere routine for RS6000... */
  119.       /*TorusMethodRegister(GLTorusDraw);*/
  120.       TriangleMethodRegister(GLTriangleDraw);
  121.   
  122. --- 123,128 ----
  123. ***************
  124. *** 133,138 ****
  125. --- 138,152 ----
  126.       short val;
  127.       float tmp;
  128.   
  129. +     /*
  130. +      * We're only sampling the scene once, so we need
  131. +      * not do lots of work to determine exactly what
  132. +      * animated transformations are doing...
  133. +      */
  134. +     Options.samples = 1;
  135. +     SamplingSetOptions(Options.samples, Options.gaussian,
  136. +                Options.filterwidth);
  137. + 
  138.       DrawInit();
  139.       qdevice(ESCKEY);
  140.       qdevice(SPACEKEY);
  141. ***************
  142. *** 197,203 ****
  143. --- 211,219 ----
  144.           doublebuffer();
  145.       }
  146.   
  147. + 
  148.       gconfig();
  149. +     blendfunction(BF_SA, BF_MSA);
  150.       /*
  151.        * Initialize viewing matrix.
  152.        */
  153. ***************
  154. *** 329,339 ****
  155.                    DIFFUSE, 0, 0, 0,
  156.                    SPECULAR, 0, 0, 0,
  157.                    SHININESS, 0,
  158.                    LMNULL};
  159.   static float    *amb = &surfprops[1],
  160.           *diff = &surfprops[5],
  161.           *spec = &surfprops[9],
  162. !         *shine = &surfprops[13];
  163.   
  164.   static void
  165.   GLPushSurface(surf)
  166. --- 345,357 ----
  167.                    DIFFUSE, 0, 0, 0,
  168.                    SPECULAR, 0, 0, 0,
  169.                    SHININESS, 0,
  170. +                  ALPHA, 1,
  171.                    LMNULL};
  172.   static float    *amb = &surfprops[1],
  173.           *diff = &surfprops[5],
  174.           *spec = &surfprops[9],
  175. !         *shine = &surfprops[13],
  176. !         *alpha = &surfprops[15];
  177.   
  178.   static void
  179.   GLPushSurface(surf)
  180. ***************
  181. *** 357,363 ****
  182.           spec[0] = surf->spec.r; spec[1] = surf->spec.g;
  183.               spec[2] = surf->spec.b;
  184.           shine[0] = surf->srexp;
  185. !         lmdef(DEFMATERIAL, cursurf, 15, surfprops);
  186.           lastsurf = surf;
  187.       }
  188.       lmbind(MATERIAL, cursurf);
  189. --- 375,383 ----
  190.           spec[0] = surf->spec.r; spec[1] = surf->spec.g;
  191.               spec[2] = surf->spec.b;
  192.           shine[0] = surf->srexp;
  193. !         *alpha = 1. - surf->transp;
  194. !         lmdef(DEFMATERIAL, cursurf, sizeof(surfprops)/sizeof(float),
  195. !             surfprops);
  196.           lastsurf = surf;
  197.       }
  198.       lmbind(MATERIAL, cursurf);
  199. ***************
  200. *** 630,639 ****
  201.   
  202.   float lightprops[] =    {POSITION, 0., 0., 0., 0.,
  203.                LCOLOR, 0, 0, 0,
  204.                LMNULL};
  205.   
  206.   float     *lpos = &lightprops[1],
  207. !     *lcolor = &lightprops[6];
  208.   
  209.   float lmodel[] =    {AMBIENT, 1., 1., 1.,
  210.                ATTENUATION, 1., 0.,
  211. --- 650,664 ----
  212.   
  213.   float lightprops[] =    {POSITION, 0., 0., 0., 0.,
  214.                LCOLOR, 0, 0, 0,
  215. +              SPOTDIRECTION, 0., 0., 0.,
  216. +              SPOTLIGHT, 0., 180.,
  217.                LMNULL};
  218.   
  219.   float     *lpos = &lightprops[1],
  220. !     *lcolor = &lightprops[6],
  221. !     *spotdir = &lightprops[10],
  222. !     *spotexp = &lightprops[14],
  223. !     *spotspread = &lightprops[15];
  224.   
  225.   float lmodel[] =    {AMBIENT, 1., 1., 1.,
  226.                ATTENUATION, 1., 0.,
  227. ***************
  228. *** 676,686 ****
  229.       lmodel[2] = Options.ambient.g;
  230.       lmodel[3] = Options.ambient.b;
  231.   
  232. ! #ifdef sgi
  233. !     lmdef(DEFLMODEL, 1, 14, lmodel);
  234. ! #else
  235. !     lmdef(DEFLMODEL, 1, 10, lmodel);
  236. ! #endif
  237.       lmbind(LMODEL, 1);
  238.   }
  239.   
  240. --- 701,707 ----
  241.       lmodel[2] = Options.ambient.g;
  242.       lmodel[3] = Options.ambient.b;
  243.   
  244. !     lmdef(DEFLMODEL, 1, sizeof(lmodel) / sizeof(float), lmodel);
  245.       lmbind(LMODEL, 1);
  246.   }
  247.   
  248. ***************
  249. *** 704,710 ****
  250.       lpos[1] = ext->pos.y;
  251.       lpos[2] = ext->pos.z;
  252.       lpos[3] = 1.;
  253. !     lmdef(DEFLIGHT, curlight++, 10, lightprops);
  254.   }
  255.   
  256.   
  257. --- 725,732 ----
  258.       lpos[1] = ext->pos.y;
  259.       lpos[2] = ext->pos.z;
  260.       lpos[3] = 1.;
  261. !     lmdef(DEFLIGHT, curlight++, sizeof(lightprops)/sizeof(float),
  262. !         lightprops);
  263.   }
  264.   
  265.   
  266. ***************
  267. *** 716,722 ****
  268.       lpos[1] = inf->dir.y;
  269.       lpos[2] = inf->dir.z;
  270.       lpos[3] = 0.;
  271. !     lmdef(DEFLIGHT, curlight++, 10, lightprops);
  272.   }
  273.   
  274.   static void
  275. --- 738,745 ----
  276.       lpos[1] = inf->dir.y;
  277.       lpos[2] = inf->dir.z;
  278.       lpos[3] = 0.;
  279. !     lmdef(DEFLIGHT, curlight++, sizeof(lightprops)/sizeof(float),
  280. !         lightprops);
  281.   }
  282.   
  283.   static void
  284. ***************
  285. *** 727,733 ****
  286.       lpos[1] = pt->pos.y;
  287.       lpos[2] = pt->pos.z;
  288.       lpos[3] = 1.;
  289. !     lmdef(DEFLIGHT, curlight++, 10, lightprops);
  290.   }
  291.   
  292.   static void
  293. --- 750,757 ----
  294.       lpos[1] = pt->pos.y;
  295.       lpos[2] = pt->pos.z;
  296.       lpos[3] = 1.;
  297. !     lmdef(DEFLIGHT, curlight++, sizeof(lightprops) / sizeof(float),
  298. !             lightprops);
  299.   }
  300.   
  301.   static void
  302. ***************
  303. *** 734,739 ****
  304. --- 758,776 ----
  305.   GLSpotLight(spot)
  306.   Spotlight *spot;
  307.   {
  308. +     lpos[0] = spot->pos.x;
  309. +     lpos[1] = spot->pos.y;
  310. +     lpos[2] = spot->pos.z;
  311. +     lpos[3] = 1.;
  312. +     spotdir[0] = spot->pos.x - spot->dir.x;
  313. +     spotdir[1] = spot->pos.y - spot->dir.y;
  314. +     spotdir[2] = spot->pos.z - spot->dir.z;
  315. +     *spotexp = spot->coef;
  316. +     *spotspread = spot->falloff;
  317. +     lmdef(DEFLIGHT, curlight++, sizeof(lightprops) / sizeof(float),
  318. +             lightprops);
  319. +     *spotspread = 180.;
  320. +     *spotexp = 1.;
  321.   }
  322.   
  323.   static float boxfaces[6][4][3] = {
  324.  
  325. Index: etc/rsconvert/yacc.y
  326. Prereq: 4.0
  327. *** etc/rsconvert/yacc.y.old    1991/07/17 17:11:26
  328. --- etc/rsconvert/yacc.y    1991/09/28 20:44:09
  329. ***************
  330. *** 11,17 ****
  331.   /*                                        */
  332.   /* There is no warranty or other guarantee of fitness of this software       */
  333.   /* for any purpose.  It is provided solely "as is".               */
  334. ! /* $Id: yacc.y,v 4.0 91/07/17 17:11:26 kolb Exp Locker: kolb $ */
  335.   %{
  336.   #include <stdio.h>
  337.   #include "libcommon/common.h"
  338. --- 11,17 ----
  339.   /*                                        */
  340.   /* There is no warranty or other guarantee of fitness of this software       */
  341.   /* for any purpose.  It is provided solely "as is".               */
  342. ! /* $Id: yacc.y,v 4.0.1.1 91/09/28 20:44:09 cek Exp Locker: cek $ */
  343.   %{
  344.   #include <stdio.h>
  345.   #include "libcommon/common.h"
  346. ***************
  347. *** 399,404 ****
  348. --- 399,405 ----
  349.               if ($10) {
  350.                   WriteString("\ttranslu");
  351.                   WriteFloat($10);
  352. +                 WriteString("1 1 1");
  353.                   WriteFloat($11);
  354.               }
  355.           }
  356.  
  357. Index: libshade/yacc.y
  358. Prereq: 4.0
  359. *** libshade/yacc.y.old    1991/07/17 17:39:39
  360. --- libshade/yacc.y    1991/10/08 20:40:35
  361. ***************
  362. *** 11,17 ****
  363.   /*                                        */
  364.   /* There is no warranty or other guarantee of fitness of this software       */
  365.   /* for any purpose.  It is provided solely "as is".               */
  366. ! /* $Id: yacc.y,v 4.0 91/07/17 17:39:39 kolb Exp Locker: kolb $ */
  367.   %{
  368.   #include "rayshade.h"
  369.   
  370. --- 11,17 ----
  371.   /*                                        */
  372.   /* There is no warranty or other guarantee of fitness of this software       */
  373.   /* for any purpose.  It is provided solely "as is".               */
  374. ! /* $Id: yacc.y,v 4.0.1.2 91/10/08 20:40:35 cek Exp Locker: cek $ */
  375.   %{
  376.   #include "rayshade.h"
  377.   
  378. ***************
  379. *** 72,77 ****
  380. --- 72,78 ----
  381.   #include "libcommon/rotate.h"
  382.   #include "libcommon/scale.h"
  383.   #include "libcommon/translate.h"
  384. + #include "libcommon/xform.h"
  385.   
  386.   Geom *NewAggregate();
  387.   char yyfilename[BUFSIZ];            /* Input filename */
  388. ***************
  389. *** 126,132 ****
  390.   %token tTEXTURE tCHECKER tWOOD tCONTRAST tCUTOFF tCLOUD
  391.   %token tAMBIENT tDIFFUSE tREFLECT tTRANSP tSPECULAR tSPECPOW
  392.   %token tINDEX tATMOSPHERE tNOSHADOW tAREA tTRANSLU tTORUS
  393. ! %token tEYESEP tSHADOWTRANSP tREPORT tVERBOSE tQUIET tWINDOW tSTRIPE
  394.   %token tMAP tUV tSPHERICAL tCYLINDRICAL tPLANAR
  395.   %token tIMAGE tSMOOTH tCOMPONENT tTEXTSURF tRANGE tTILE tSTARTTIME tFRAMELENGTH
  396.   %token tNAME tFILTER tGAUSS tBODY tSAMPLE tEXTINCT tWINDY tMOUNT
  397. --- 127,133 ----
  398.   %token tTEXTURE tCHECKER tWOOD tCONTRAST tCUTOFF tCLOUD
  399.   %token tAMBIENT tDIFFUSE tREFLECT tTRANSP tSPECULAR tSPECPOW
  400.   %token tINDEX tATMOSPHERE tNOSHADOW tAREA tTRANSLU tTORUS
  401. ! %token tEYESEP tSHADOWTRANSP tREPORT tVERBOSE tQUIET tWINDOW tCROP tSTRIPE
  402.   %token tMAP tUV tSPHERICAL tCYLINDRICAL tPLANAR
  403.   %token tIMAGE tSMOOTH tCOMPONENT tTEXTSURF tRANGE tTILE tSTARTTIME tFRAMELENGTH
  404.   %token tNAME tFILTER tGAUSS tBODY tSAMPLE tEXTINCT tWINDY tMOUNT
  405. ***************
  406. *** 166,171 ****
  407. --- 167,173 ----
  408.           | Fov
  409.           | Screen
  410.           | Window
  411. +         | Crop
  412.           | Report
  413.           | Aperture
  414.           | Focaldist
  415. ***************
  416. *** 489,494 ****
  417. --- 491,533 ----
  418.               TransRotateSetTheta($$, $5);
  419.               if (!$$->animated)
  420.                   TransPropagate($$);
  421. +         }
  422. +         | tTRANSFORM    AnimExpr AnimExpr AnimExpr
  423. +                 AnimExpr AnimExpr AnimExpr
  424. +                 AnimExpr AnimExpr AnimExpr
  425. +         {
  426. +             $$ = TransXformCreate();
  427. +             TransXformSetX0($$, $2);
  428. +             TransXformSetY0($$, $3);
  429. +             TransXformSetZ0($$, $4);
  430. +             TransXformSetX1($$, $5);
  431. +             TransXformSetY1($$, $6);
  432. +             TransXformSetZ1($$, $7);
  433. +             TransXformSetX2($$, $8);
  434. +             TransXformSetY2($$, $9);
  435. +             TransXformSetZ2($$, $10);
  436. +             if (!$$->animated)
  437. +                 TransPropagate($$);
  438. +         }
  439. +         | tTRANSFORM    AnimExpr AnimExpr AnimExpr
  440. +                 AnimExpr AnimExpr AnimExpr
  441. +                 AnimExpr AnimExpr AnimExpr
  442. +                 AnimExpr AnimExpr AnimExpr
  443. +         {
  444. +             TransXformSetX0($$, $2);
  445. +             TransXformSetY0($$, $3);
  446. +             TransXformSetZ0($$, $4);
  447. +             TransXformSetX1($$, $5);
  448. +             TransXformSetY1($$, $6);
  449. +             TransXformSetZ1($$, $7);
  450. +             TransXformSetX2($$, $8);
  451. +             TransXformSetY2($$, $9);
  452. +             TransXformSetZ2($$, $10);
  453. +             TransXformSetXt($$, $11);
  454. +             TransXformSetYt($$, $12);
  455. +             TransXformSetZt($$, $13);
  456. +             if (!$$->animated)
  457. +                 TransPropagate($$);
  458.           };
  459.   Eyep        : tEYEP Vector Transforms
  460.           {
  461. ***************
  462. *** 602,616 ****
  463.               }
  464.           }
  465.           ;
  466. ! Window        : tWINDOW Expr Expr Expr Expr
  467.           {
  468.               if (!Options.window_set) {
  469. !                 Options.xmin = $2;
  470. !                 Options.xmax = $3;
  471. !                 Options.ymin = $4;
  472. !                 Options.ymax = $5;
  473.               }
  474.           }
  475.   Report        : tREPORT Verbose Quiet IExpr Filename
  476.           {
  477.               if (!Options.freq_set)
  478. --- 641,666 ----
  479.               }
  480.           }
  481.           ;
  482. ! Window        : tWINDOW IExpr IExpr IExpr IExpr
  483.           {
  484.               if (!Options.window_set) {
  485. !                 Options.window[LOW][X] = $2;
  486. !                 Options.window[HIGH][X] = $3;
  487. !                 Options.window[LOW][Y] = $4;
  488. !                 Options.window[HIGH][Y] = $5;
  489.               }
  490.           }
  491. +         ;
  492. + Crop        : tCROP Expr Expr Expr Expr
  493. +         {
  494. +             if (!Options.crop_set) {
  495. +                 Options.crop[LOW][X] = $2;
  496. +                 Options.crop[HIGH][X] = $3;
  497. +                 Options.crop[LOW][Y] = $4;
  498. +                 Options.crop[HIGH][Y] = $5;
  499. +             }
  500. +         }
  501. +         ;
  502.   Report        : tREPORT Verbose Quiet IExpr Filename
  503.           {
  504.               if (!Options.freq_set)
  505. ***************
  506. *** 667,673 ****
  507.           ;
  508.   Shadowtransp    : tSHADOWTRANSP
  509.           {
  510. !             Options.shadowtransp = TRUE;
  511.           }
  512.           ;
  513.   Light        : LightType
  514. --- 717,723 ----
  515.           ;
  516.   Shadowtransp    : tSHADOWTRANSP
  517.           {
  518. !             Options.shadowtransp = !Options.shadowtransp;
  519.           }
  520.           ;
  521.   Light        : LightType
  522.  
  523. Index: libray/libobj/grid.c
  524. Prereq: 4.0
  525. *** libray/libobj/grid.c.old    1991/07/17 14:38:02
  526. --- libray/libobj/grid.c    1991/10/04 15:55:37
  527. ***************
  528. *** 13,21 ****
  529.    * There is no warranty or other guarantee of fitness of this software
  530.    * for any purpose.  It is provided solely "as is".
  531.    *
  532. !  * $Id: grid.c,v 4.0 91/07/17 14:38:02 kolb Exp Locker: kolb $
  533.    *
  534.    * $Log:    grid.c,v $
  535.    * Revision 4.0  91/07/17  14:38:02  kolb
  536.    * Initial version.
  537.    * 
  538. --- 13,24 ----
  539.    * There is no warranty or other guarantee of fitness of this software
  540.    * for any purpose.  It is provided solely "as is".
  541.    *
  542. !  * $Id: grid.c,v 4.0.1.1 91/10/04 15:55:37 cek Exp Locker: cek $
  543.    *
  544.    * $Log:    grid.c,v $
  545. +  * Revision 4.0.1.1  91/10/04  15:55:37  cek
  546. +  * patch1: Removed straight floating point comparisons.
  547. +  * 
  548.    * Revision 4.0  91/07/17  14:38:02  kolb
  549.    * Initial version.
  550.    * 
  551. ***************
  552. *** 120,132 ****
  553.       x = x2voxel(grid, curpos.x);
  554.       if (x == grid->xsize)
  555.           x--;
  556. !     if (ray->dir.x < 0.) {
  557.           tMaxX = offset + (voxel2x(grid, x) - curpos.x) / ray->dir.x;
  558.           tDeltaX = grid->voxsize[X] / - ray->dir.x;
  559.           stepX = outX = -1;
  560.           raybounds[LOW][X] = &np.x;
  561.           raybounds[HIGH][X] = &curpos.x;
  562. !     } else if (ray->dir.x > 0.) {
  563.           tMaxX = offset + (voxel2x(grid, x+1) - curpos.x) / ray->dir.x;
  564.           tDeltaX = grid->voxsize[X] / ray->dir.x;
  565.           stepX = 1;
  566. --- 123,140 ----
  567.       x = x2voxel(grid, curpos.x);
  568.       if (x == grid->xsize)
  569.           x--;
  570. !     if (fabs(ray->dir.x) < EPSILON) {
  571. !         tMaxX = FAR_AWAY;
  572. !         raybounds[LOW][X] = &curpos.x;
  573. !         raybounds[HIGH][X] = &np.x;
  574. !         tDeltaX = 0.;
  575. !     } else if (ray->dir.x < 0.) {
  576.           tMaxX = offset + (voxel2x(grid, x) - curpos.x) / ray->dir.x;
  577.           tDeltaX = grid->voxsize[X] / - ray->dir.x;
  578.           stepX = outX = -1;
  579.           raybounds[LOW][X] = &np.x;
  580.           raybounds[HIGH][X] = &curpos.x;
  581. !     } else {
  582.           tMaxX = offset + (voxel2x(grid, x+1) - curpos.x) / ray->dir.x;
  583.           tDeltaX = grid->voxsize[X] / ray->dir.x;
  584.           stepX = 1;
  585. ***************
  586. *** 133,155 ****
  587.           outX = grid->xsize;
  588.           raybounds[LOW][X] = &curpos.x;
  589.           raybounds[HIGH][X] = &np.x;
  590. -     } else {
  591. -         tMaxX = FAR_AWAY;
  592. -         raybounds[LOW][X] = &curpos.x;
  593. -         raybounds[HIGH][X] = &np.x;
  594. -         tDeltaX = 0.;
  595.       }
  596.   
  597.       y = y2voxel(grid, curpos.y);
  598.       if (y == grid->ysize)
  599.           y--;
  600. !     if (ray->dir.y < 0.) {
  601.           tMaxY = offset + (voxel2y(grid, y) - curpos.y) / ray->dir.y;
  602.           tDeltaY = grid->voxsize[Y] / - ray->dir.y;
  603.           stepY = outY = -1;
  604.           raybounds[LOW][Y] = &np.y;
  605.           raybounds[HIGH][Y] = &curpos.y;
  606. !     } else if (ray->dir.y > 0.) {
  607.           tMaxY = offset + (voxel2y(grid, y+1) - curpos.y) / ray->dir.y;
  608.           tDeltaY = grid->voxsize[Y] / ray->dir.y;
  609.           stepY = 1;
  610. --- 141,164 ----
  611.           outX = grid->xsize;
  612.           raybounds[LOW][X] = &curpos.x;
  613.           raybounds[HIGH][X] = &np.x;
  614.       }
  615.   
  616.       y = y2voxel(grid, curpos.y);
  617.       if (y == grid->ysize)
  618.           y--;
  619. ! 
  620. !     if (fabs(ray->dir.y) < EPSILON) {
  621. !         tMaxY = FAR_AWAY;
  622. !         raybounds[LOW][Y] = &curpos.y;
  623. !         raybounds[HIGH][Y] = &np.y;
  624. !         tDeltaY = 0.;
  625. !     } else if (ray->dir.y < 0.) {
  626.           tMaxY = offset + (voxel2y(grid, y) - curpos.y) / ray->dir.y;
  627.           tDeltaY = grid->voxsize[Y] / - ray->dir.y;
  628.           stepY = outY = -1;
  629.           raybounds[LOW][Y] = &np.y;
  630.           raybounds[HIGH][Y] = &curpos.y;
  631. !     } else {
  632.           tMaxY = offset + (voxel2y(grid, y+1) - curpos.y) / ray->dir.y;
  633.           tDeltaY = grid->voxsize[Y] / ray->dir.y;
  634.           stepY = 1;
  635. ***************
  636. *** 156,178 ****
  637.           outY = grid->ysize;
  638.           raybounds[LOW][Y] = &curpos.y;
  639.           raybounds[HIGH][Y] = &np.y;
  640. -     } else {
  641. -         tMaxY = FAR_AWAY;
  642. -         raybounds[LOW][Y] = &curpos.y;
  643. -         raybounds[HIGH][Y] = &np.y;
  644. -         tDeltaY = 0.;
  645.       }
  646.   
  647.       z = z2voxel(grid, curpos.z);
  648.       if (z == grid->zsize)
  649.           z--;
  650. !     if (ray->dir.z < 0.) {
  651.           tMaxZ = offset + (voxel2z(grid, z) - curpos.z) / ray->dir.z;
  652.           tDeltaZ = grid->voxsize[Z] / - ray->dir.z;
  653.           stepZ = outZ = -1;
  654.           raybounds[LOW][Z] = &np.z;
  655.           raybounds[HIGH][Z] = &curpos.z;
  656. !     } else if (ray->dir.z > 0.) {
  657.           tMaxZ = offset + (voxel2z(grid, z+1) - curpos.z) / ray->dir.z;
  658.           tDeltaZ = grid->voxsize[Z] / ray->dir.z;
  659.           stepZ = 1;
  660. --- 165,187 ----
  661.           outY = grid->ysize;
  662.           raybounds[LOW][Y] = &curpos.y;
  663.           raybounds[HIGH][Y] = &np.y;
  664.       }
  665.   
  666.       z = z2voxel(grid, curpos.z);
  667.       if (z == grid->zsize)
  668.           z--;
  669. !     if (fabs(ray->dir.z) < EPSILON) {
  670. !         tMaxZ = FAR_AWAY;
  671. !         raybounds[LOW][Z] = &curpos.z;
  672. !         raybounds[HIGH][Z] = &np.z;
  673. !         tDeltaZ = 0.;
  674. !     } else if (ray->dir.z < 0.) {
  675.           tMaxZ = offset + (voxel2z(grid, z) - curpos.z) / ray->dir.z;
  676.           tDeltaZ = grid->voxsize[Z] / - ray->dir.z;
  677.           stepZ = outZ = -1;
  678.           raybounds[LOW][Z] = &np.z;
  679.           raybounds[HIGH][Z] = &curpos.z;
  680. !     } else {
  681.           tMaxZ = offset + (voxel2z(grid, z+1) - curpos.z) / ray->dir.z;
  682.           tDeltaZ = grid->voxsize[Z] / ray->dir.z;
  683.           stepZ = 1;
  684. ***************
  685. *** 179,189 ****
  686.           outZ = grid->zsize;
  687.           raybounds[LOW][Z] = &curpos.z;
  688.           raybounds[HIGH][Z] = &np.z;
  689. -     } else {
  690. -         tMaxZ = FAR_AWAY;
  691. -         raybounds[LOW][Z] = &curpos.z;
  692. -         raybounds[HIGH][Z] = &np.z;
  693. -         tDeltaZ = 0.;
  694.       }
  695.   
  696.       VecScale(tDeltaX, ray->dir, &pDeltaX);
  697. --- 188,193 ----
  698.  
  699. Index: libshade/viewing.c
  700. Prereq: 4.0
  701. *** libshade/viewing.c.old    1991/07/17 14:48:18
  702. --- libshade/viewing.c    1991/09/29 15:53:09
  703. ***************
  704. *** 13,21 ****
  705.    * There is no warranty or other guarantee of fitness of this software
  706.    * for any purpose.  It is provided solely "as is".
  707.    *
  708. !  * $Id: viewing.c,v 4.0 91/07/17 14:48:18 kolb Exp Locker: kolb $
  709.    *
  710.    * $Log:    viewing.c,v $
  711.    * Revision 4.0  91/07/17  14:48:18  kolb
  712.    * Initial version.
  713.    * 
  714. --- 13,24 ----
  715.    * There is no warranty or other guarantee of fitness of this software
  716.    * for any purpose.  It is provided solely "as is".
  717.    *
  718. !  * $Id: viewing.c,v 4.0.1.1 91/09/29 15:53:09 cek Exp Locker: cek $
  719.    *
  720.    * $Log:    viewing.c,v $
  721. +  * Revision 4.0.1.1  91/09/29  15:53:09  cek
  722. +  * patch1: Added support for window and crop commands.
  723. +  * 
  724.    * Revision 4.0  91/07/17  14:48:18  kolb
  725.    * Initial version.
  726.    * 
  727. ***************
  728. *** 71,76 ****
  729. --- 74,85 ----
  730.       VecScale(magnitude, Screen.scrni, &Screen.scrnx);
  731.       magnitude = 2.*Camera.lookdist * tan(deg2rad(0.5*Camera.vfov)) /
  732.                   Screen.yres;
  733. + #ifndef URT
  734. +     /*
  735. +      * If using "generic" file format, render top-to-bottom (yick).
  736. +      */
  737. +     magnitude *= -1;
  738. + #endif
  739.       VecScale(magnitude, Screen.scrnj, &Screen.scrny);
  740.   
  741.       Screen.firstray.x -= 0.5*Screen.yres*Screen.scrny.x +
  742. ***************
  743. *** 119,124 ****
  744. --- 128,138 ----
  745.   void
  746.   ViewingSetup()
  747.   {
  748. + #define SWAP(a,b)    (tmp = (a), (a) = (b), (b) = tmp)
  749. + 
  750. +     Float tmp;
  751. +     int xwidth, ywidth;
  752. + 
  753.       if (Options.stereo && Options.eyesep == UNSET)
  754.           RLerror(RL_PANIC,
  755.               "No eye separation specified for stereo rendering.\n");
  756. ***************
  757. *** 137,147 ****
  758.       if (Screen.yres == UNSET)
  759.           Screen.yres = YRESOLUTION;
  760.   
  761.   #ifdef URT
  762.       /*
  763. !      * If using the URT, we can use the RLE file header to
  764. !      * determine window size.  Screen size (Screen.xres, Screen.yres)
  765. !      * is determined as usual (from command line or input file).
  766.        */
  767.       if (Options.appending) {
  768.           /*
  769. --- 151,205 ----
  770.       if (Screen.yres == UNSET)
  771.           Screen.yres = YRESOLUTION;
  772.   
  773. +     /*
  774. +      * The window to be rendered is defined by applying
  775. +      * the crop window to the sub window.  The subwindow
  776. +      * is defined using pixel numbers, and must be within
  777. +      * [0, xres -1],[0, yres -1].  The default is the entire
  778. +      * screen.  The crop window is specified using normalized
  779. +      * coordinates.
  780. +      */
  781. + 
  782. +     if (!Options.window_set) {
  783. +         /* If no window set, set equal to entire screen. */    
  784. +         Options.window[LOW][X] = Options.window[LOW][Y] = 0;
  785. +         Options.window[HIGH][X] = Screen.xres -1;
  786. +         Options.window[HIGH][Y] = Screen.yres -1;
  787. +     }
  788. + 
  789. +     /* Truncate crop window to legal limits. */
  790. +     if (Options.crop[LOW][X] > Options.crop[HIGH][X])
  791. +         SWAP(Options.crop[LOW][X], Options.crop[HIGH][X]);
  792. +     if (Options.crop[LOW][Y] > Options.crop[HIGH][Y])
  793. +         SWAP(Options.crop[LOW][Y], Options.crop[HIGH][Y]);
  794. +     if (Options.crop[LOW][X] < 0.) Options.crop[LOW][X] = 0.;
  795. +     if (Options.crop[LOW][Y] < 0.) Options.crop[LOW][Y] = 0.;
  796. +     if (Options.crop[HIGH][X] > 1.) Options.crop[HIGH][X] = 1.;
  797. +     if (Options.crop[HIGH][Y] > 1.) Options.crop[HIGH][Y] = 1.;
  798. + 
  799. +     xwidth = Options.window[HIGH][X] - Options.window[LOW][X];
  800. +     ywidth = Options.window[HIGH][Y] - Options.window[LOW][Y];
  801. + 
  802. +     /* Compute x and y extents of window to be renered. */
  803. +     Screen.minx = (int)(Options.window[LOW][X] +
  804. +             Options.crop[LOW][X] * xwidth);
  805. +     Screen.maxx = (int)(Options.window[LOW][X] +
  806. +             Options.crop[HIGH][X] * xwidth);
  807. +     Screen.miny = (int)(Options.window[LOW][Y] +
  808. +             Options.crop[LOW][Y] * ywidth);
  809. +     Screen.maxy = (int)(Options.window[LOW][Y] +
  810. +             Options.crop[HIGH][Y] * ywidth);
  811. + 
  812.   #ifdef URT
  813.       /*
  814. !      * If using the URT, we should use the RLE file header to
  815. !      * determine cropped window size.  Screen size
  816. !      * (Screen.xres, Screen.yres) is determined from command
  817. !      * line or input file, as usual.
  818. !      *
  819. !      * If the cropped window computed in PictureSetWindow()
  820. !      * is not equal the cropped window computed above,
  821. !      * a warning message is issued.
  822.        */
  823.       if (Options.appending) {
  824.           /*
  825. ***************
  826. *** 148,165 ****
  827.            * Read image header to determine window size.
  828.            */
  829.           PictureSetWindow();
  830. -     } else {
  831. - #endif
  832. -     if (Options.xmin < 0.) Options.xmin = 0.;
  833. -     if (Options.ymin < 0.) Options.ymin = 0.;
  834. -     if (Options.xmax > 1.) Options.xmax = 1.;
  835. -     if (Options.xmax > 1.) Options.xmax = 1.;
  836. - 
  837. -     Screen.minx = (int)(Options.xmin * Screen.xres);
  838. -     Screen.maxx = (int)(Options.xmax * Screen.xres -1);
  839. -     Screen.miny = (int)(Options.ymin * Screen.yres);
  840. -     Screen.maxy = (int)(Options.ymax * Screen.yres -1);
  841. - #ifdef URT
  842.       }
  843.   #endif
  844.   
  845. --- 206,211 ----
  846.  
  847. Index: README
  848. *** README.old    1991/07/17 14:27:21
  849. --- README    1991/10/10 18:16:26
  850. ***************
  851. *** 35,41 ****
  852.   conversion utility, rsconvert.  See Doc/Changes for a list of some of
  853.   the more important changes.
  854.   
  855. ! Rayshade v4.0 has been tested on several different UNIX-based computers,
  856.   including: SGI 4D, IBM RS6000, Sun Sparcstation 1, Sun 3, DECstation,
  857.   Apollo DN10000, NeXT.  If your machine has a C compiler, enough memory
  858.   (at least 2Mb), and runs something resembling UNIX, rayshade should be fairly
  859. --- 35,41 ----
  860.   conversion utility, rsconvert.  See Doc/Changes for a list of some of
  861.   the more important changes.
  862.   
  863. ! Rayshade v4.0 Beta has been tested on several different UNIX-based computers,
  864.   including: SGI 4D, IBM RS6000, Sun Sparcstation 1, Sun 3, DECstation,
  865.   Apollo DN10000, NeXT.  If your machine has a C compiler, enough memory
  866.   (at least 2Mb), and runs something resembling UNIX, rayshade should be fairly
  867. ***************
  868. *** 90,96 ****
  869.   Standard Procedural Databases through awk using the etc/nff2shade.awk
  870.   script, and then feed the results to rayshade.  Note that for several
  871.   of the databases, you will need to tweak the data a bit to achieve
  872. ! optimal runtimes.  See spd2shade.awk for details.
  873.   
  874.   Rayshade was not designed to be, nor is it, the be-all and end-all of
  875.   ray tracers.  It was meant to be reasonably fast, portable, and easy to
  876. --- 90,96 ----
  877.   Standard Procedural Databases through awk using the etc/nff2shade.awk
  878.   script, and then feed the results to rayshade.  Note that for several
  879.   of the databases, you will need to tweak the data a bit to achieve
  880. ! optimal runtimes.  See nff2shade.awk for details.
  881.   
  882.   Rayshade was not designed to be, nor is it, the be-all and end-all of
  883.   ray tracers.  It was meant to be reasonably fast, portable, and easy to
  884. ***************
  885. *** 136,146 ****
  886.   for the 'net tracer' conversations and for the inspiration to do something
  887.   to clean up rayshade.
  888.   
  889. ! Thanks to Marc Andreessen, Ray Bellis, Dominique Boisvert, William Bouma,
  890. ! Allen Braunsdorf, Jeff Butterworth, Nick Carriero, Nancy Everson, Tom Friedel,
  891. ! Robert Funchess, David Gelernter, Mike Gigante, Ed Herderick, John Knuston,
  892. ! Raphael Manfredi, Lee Moore, Dietmar Saupe, Brian Wyvill,
  893. ! and everybody else for their bugfixes, suggestions, input files,
  894.   encouragement, and other feedback.
  895.   
  896.   Thanks also to Mark Podlipec for providing the blob and torus objects.
  897. --- 136,146 ----
  898.   for the 'net tracer' conversations and for the inspiration to do something
  899.   to clean up rayshade.
  900.   
  901. ! Thanks to Marc Andreessen, Tony Apodaca, Ray Bellis, Dominique Boisvert,
  902. ! William Bouma, Allen Braunsdorf, Jeff Butterworth, Nick Carriero,
  903. ! Nancy Everson, Tom Friedel, Robert Funchess, David Gelernter, Mike Gigante,
  904. ! Ed Herderick, John Knuston, Raphael Manfredi, Lee Moore, Dietmar Saupe,
  905. ! Brian Wyvill, and everybody else for their bugfixes, suggestions, input files,
  906.   encouragement, and other feedback.
  907.   
  908.   Thanks also to Mark Podlipec for providing the blob and torus objects.
  909. ***************
  910. *** 156,179 ****
  911.   particularly braindamaged section of code, please let us know.  We can be
  912.   most easily reached by email.
  913.   
  914. ! Once you start using rayshade, send a message to
  915. ! rayshade@weedeater.math.yale.edu and you'll be added to the rayshade
  916. ! mailing list.
  917.   
  918.   Craig Kolb
  919.   Rod Bogart
  920. ! 17 July 1991
  921.   rayshade@weedeater.math.yale.edu
  922.   
  923.   Craig Kolb 
  924. - Department of Mathematics
  925. - Yale University
  926. - 10 Hillhouse Ave.
  927. - Box 2155 Yale Station
  928. - New Haven, CT  06520-2155
  929. - (Until September, 1991)
  930. - 
  931.   Computer Science Department
  932.   Princeton University
  933.   Princeton, NJ 08544
  934. - (September, 1991 onwards)
  935. --- 156,172 ----
  936.   particularly braindamaged section of code, please let us know.  We can be
  937.   most easily reached by email.
  938.   
  939. ! Once you start using rayshade, and you'd like to be kept abreast of
  940. ! changes and extensions, send a message to rayshade-request@cs.princeton.edu,
  941. ! and we'll add you to the rayshade mailing list.
  942.   
  943.   Craig Kolb
  944.   Rod Bogart
  945. ! 9 October 1991
  946.   rayshade@weedeater.math.yale.edu
  947.   
  948.   Craig Kolb 
  949.   Computer Science Department
  950.   Princeton University
  951. + 35 Olden St.
  952.   Princeton, NJ 08544
  953.  
  954. Index: libshade/options.c
  955. Prereq: 4.0
  956. *** libshade/options.c.old    1991/07/17 14:46:47
  957. --- libshade/options.c    1991/10/05 18:23:15
  958. ***************
  959. *** 13,21 ****
  960.    * There is no warranty or other guarantee of fitness of this software
  961.    * for any purpose.  It is provided solely "as is".
  962.    *
  963. !  * $Id: options.c,v 4.0 91/07/17 14:46:47 kolb Exp Locker: kolb $
  964.    *
  965.    * $Log:    options.c,v $
  966.    * Revision 4.0  91/07/17  14:46:47  kolb
  967.    * Initial version.
  968.    * 
  969. --- 13,25 ----
  970.    * There is no warranty or other guarantee of fitness of this software
  971.    * for any purpose.  It is provided solely "as is".
  972.    *
  973. !  * $Id: options.c,v 4.0.1.1 91/10/05 18:23:15 cek Exp Locker: cek $
  974.    *
  975.    * $Log:    options.c,v $
  976. +  * Revision 4.0.1.1  91/10/05  18:23:15  cek
  977. +  * patch1: Changes for new window/crop options.
  978. +  * patch1: Corrected help message.
  979. +  * 
  980.    * Revision 4.0  91/07/17  14:46:47  kolb
  981.    * Initial version.
  982.    * 
  983. ***************
  984. *** 210,239 ****
  985.                   OpenStatsFile();
  986.                   argv++; argc--;
  987.                   break;
  988. - #ifdef URT
  989.               case 'W':
  990. !                 Options.xmin = atof(argv[1]);
  991. !                 Options.xmax = atof(argv[2]);
  992. !                 Options.ymin = atof(argv[3]);
  993. !                 Options.ymax = atof(argv[4]);
  994.                   Options.window_set = TRUE;
  995.                   argv += 4; argc -= 4;
  996.                   break;
  997. - #endif
  998. - #ifdef LINDA
  999.               case 'X':
  1000. !                 Options.workers = atoi(argv[1]);
  1001. !                 if (Options.workers<0 || Options.workers>17) {
  1002. !                     RLerror(RL_PANIC, "%d workers?\n",
  1003. !                             Options.workers);
  1004. !                 }
  1005. !                 argv++; argc--;
  1006. !                 break;
  1007. !             case 'w':
  1008. !                 Options.verbose_worker =
  1009. !                     !Options.verbose_worker;
  1010.                   break;
  1011. - #endif
  1012.               default:
  1013.                   RLerror(RL_PANIC,"Bad argument: %s\n",argv[0]);
  1014.           }
  1015. --- 214,235 ----
  1016.                   OpenStatsFile();
  1017.                   argv++; argc--;
  1018.                   break;
  1019.               case 'W':
  1020. !                 Options.window[LOW][X] = atoi(argv[1]);
  1021. !                 Options.window[HIGH][X] = atoi(argv[3]);
  1022. !                 Options.window[LOW][Y] = atoi(argv[3]);
  1023. !                 Options.window[HIGH][Y] = atoi(argv[4]);
  1024.                   Options.window_set = TRUE;
  1025.                   argv += 4; argc -= 4;
  1026.                   break;
  1027.               case 'X':
  1028. !                 Options.crop[LOW][X] = atof(argv[1]);
  1029. !                 Options.crop[HIGH][X] = atof(argv[2]);
  1030. !                 Options.crop[LOW][Y] = atof(argv[3]);
  1031. !                 Options.crop[HIGH][Y] = atof(argv[4]);
  1032. !                 Options.crop_set = TRUE;
  1033. !                 argv += 4; argc -= 4;
  1034.                   break;
  1035.               default:
  1036.                   RLerror(RL_PANIC,"Bad argument: %s\n",argv[0]);
  1037.           }
  1038. ***************
  1039. *** 304,312 ****
  1040. --- 300,314 ----
  1041.       fprintf(stderr,"\t-e \t\t(Write exponential RLE file.)\n");
  1042.   #endif
  1043.       fprintf(stderr,"\t-F freq\t\t(Set frequency of status report.)\n");
  1044. +     fprintf(stderr,"\t-f \t\t(Flip all triangle normals.)\n");
  1045.       fprintf(stderr,"\t-G gamma\t(Use given gamma correction exponent.)\n");
  1046. +     fprintf(stderr,"\t-g \t\t(Use Gaussian pixel filter.)\n");
  1047.       fprintf(stderr,"\t-h \t\t(Print this message.)\n");
  1048. +     fprintf(stderr,"\t-j \t\t(Toggle jittered sampling.)\n");
  1049.       fprintf(stderr,"\t-l \t\t(Render image for left eye view.)\n");
  1050. + #ifdef URT
  1051. +     fprintf(stderr,"\t-m \t\t(Output sample map in alpha channel.)\n");
  1052. + #endif
  1053.       fprintf(stderr,"\t-N number\t(Render given number of frames.)\n");
  1054.       fprintf(stderr,"\t-n \t\t(Do not render shadows.)\n");
  1055.       fprintf(stderr,"\t-O outfile \t(Set output file name.)\n");
  1056. ***************
  1057. *** 321,331 ****
  1058.       fprintf(stderr,"\t-T r g b\t(Set contrast threshold (0. - 1.).)\n");
  1059.       fprintf(stderr,"\t-V filename \t(Write verbose output to filename.)\n");
  1060.       fprintf(stderr,"\t-v \t\t(Verbose output.)\n");
  1061. ! #ifdef URT
  1062. !     fprintf(stderr,"\t-W x y x y \t(Render subwindow.)\n");
  1063. ! #endif
  1064. ! #ifdef LINDA
  1065. !     fprintf(stderr,"\t-X workers\t(Number of workers.)\n");
  1066. !     fprintf(stderr,"\t-w \t\t(Verbose worker output.)\n");
  1067. ! #endif
  1068.   }
  1069. --- 323,328 ----
  1070.       fprintf(stderr,"\t-T r g b\t(Set contrast threshold (0. - 1.).)\n");
  1071.       fprintf(stderr,"\t-V filename \t(Write verbose output to filename.)\n");
  1072.       fprintf(stderr,"\t-v \t\t(Verbose output.)\n");
  1073. !     fprintf(stderr,"\t-W x x y y \t(Render subwindow.)\n");
  1074. !     fprintf(stderr,"\t-X l r b t \t(Crop window.)\n");
  1075.   }
  1076.  
  1077. Index: libray/libobj/bounds.c
  1078. Prereq: 4.0
  1079. *** libray/libobj/bounds.c.old    1991/07/17 14:36:21
  1080. --- libray/libobj/bounds.c    1991/10/05 18:17:29
  1081. ***************
  1082. *** 13,21 ****
  1083.    * There is no warranty or other guarantee of fitness of this software
  1084.    * for any purpose.  It is provided solely "as is".
  1085.    *
  1086. !  * $Id: bounds.c,v 4.0 91/07/17 14:36:21 kolb Exp Locker: kolb $
  1087.    *
  1088.    * $Log:    bounds.c,v $
  1089.    * Revision 4.0  91/07/17  14:36:21  kolb
  1090.    * Initial version.
  1091.    * 
  1092. --- 13,27 ----
  1093.    * There is no warranty or other guarantee of fitness of this software
  1094.    * for any purpose.  It is provided solely "as is".
  1095.    *
  1096. !  * $Id: bounds.c,v 4.0.1.2 91/10/05 18:17:29 cek Exp Locker: cek $
  1097.    *
  1098.    * $Log:    bounds.c,v $
  1099. +  * Revision 4.0.1.2  91/10/05  18:17:29  cek
  1100. +  * patch1: Faster box transform, a la Jim Arvo.
  1101. +  * 
  1102. +  * Revision 4.0.1.1  1991/09/29  15:42:05  cek
  1103. +  * patch1: Justified floating point compare...
  1104. +  *
  1105.    * Revision 4.0  91/07/17  14:36:21  kolb
  1106.    * Initial version.
  1107.    * 
  1108. ***************
  1109. *** 134,139 ****
  1110. --- 140,149 ----
  1111.       } else if (pos < bounds[LOW][Z] || pos > bounds[HIGH][Z])
  1112.           return FALSE;
  1113.   
  1114. +     /*
  1115. +      * If tmin == mindist, then there was no "near"
  1116. +      * intersection farther than EPSILON away.
  1117. +      */
  1118.       if (tmin == mindist) {
  1119.           if (tmax < *maxdist) {
  1120.               *maxdist = tmax;
  1121. ***************
  1122. *** 157,165 ****
  1123.   RSMatrix *trans;
  1124.   Float objbounds[2][3];
  1125.   {
  1126. !     Vector tmp;
  1127. !     Float bounds[2][3];
  1128. !     int x, y, z;
  1129.   
  1130.       /*
  1131.        * Can't (and shouldn't) do anything with unbounded objects.
  1132. --- 167,174 ----
  1133.   RSMatrix *trans;
  1134.   Float objbounds[2][3];
  1135.   {
  1136. !     Float bounds[2][3], a, b;
  1137. !     int i, j;
  1138.   
  1139.       /*
  1140.        * Can't (and shouldn't) do anything with unbounded objects.
  1141. ***************
  1142. *** 166,200 ****
  1143.        */
  1144.       if (objbounds[LOW][X] > objbounds[HIGH][X])
  1145.           return;
  1146. -         
  1147. -     BoundsInit(bounds);
  1148.   
  1149. !     /*
  1150. !      * Find bounding box of transformed corners of bounding box.
  1151. !      */
  1152. !     for (x = 0 ; x < 2; x++) {
  1153. !         for (y = 0; y < 2; y++) {
  1154. !             for (z = 0; z < 2; z++) {
  1155. !                 tmp.x = objbounds[x][X];
  1156. !                 tmp.y = objbounds[y][Y];
  1157. !                 tmp.z = objbounds[z][Z];
  1158. !                 PointTransform(&tmp, trans);
  1159. !                 if (tmp.x < bounds[LOW][X])
  1160. !                     bounds[LOW][X] = tmp.x;
  1161. !                 if (tmp.x > bounds[HIGH][X])
  1162. !                     bounds[HIGH][X] = tmp.x;
  1163. !                 if (tmp.y < bounds[LOW][Y])
  1164. !                     bounds[LOW][Y] = tmp.y;
  1165. !                 if (tmp.y > bounds[HIGH][Y])
  1166. !                     bounds[HIGH][Y] = tmp.y;
  1167. !                 if (tmp.z < bounds[LOW][Z])
  1168. !                     bounds[LOW][Z] = tmp.z;
  1169. !                 if (tmp.z > bounds[HIGH][Z])
  1170. !                     bounds[HIGH][Z] = tmp.z;
  1171.               }
  1172.           }
  1173.       }
  1174. - 
  1175.       BoundsCopy(bounds, objbounds);
  1176.   }
  1177.   
  1178. --- 175,198 ----
  1179.        */
  1180.       if (objbounds[LOW][X] > objbounds[HIGH][X])
  1181.           return;
  1182.   
  1183. !     bounds[LOW][X] = bounds[HIGH][X] = trans->translate.x;
  1184. !     bounds[LOW][Y] = bounds[HIGH][Y] = trans->translate.y;
  1185. !     bounds[LOW][Z] = bounds[HIGH][Z] = trans->translate.z;
  1186. ! 
  1187. !     for (i = 0; i < 3; i++) {
  1188. !         for (j = 0; j < 3; j++) {
  1189. !             a = trans->matrix[j][i] * objbounds[LOW][j];
  1190. !             b = trans->matrix[j][i] * objbounds[HIGH][j];
  1191. !             if (a < b) {
  1192. !                 bounds[LOW][i] += a;
  1193. !                 bounds[HIGH][i] += b;
  1194. !             } else {
  1195. !                 bounds[LOW][i] += b;
  1196. !                 bounds[HIGH][i] += a;
  1197.               }
  1198.           }
  1199.       }
  1200.       BoundsCopy(bounds, objbounds);
  1201.   }
  1202.   
  1203.  
  1204. Index: Doc/quickref.txt
  1205. *** Doc/quickref.txt.old    1991/07/17 14:25:15
  1206. --- Doc/quickref.txt    1991/10/05 18:00:35
  1207. ***************
  1208. *** 40,46 ****
  1209.   -r             Right eye view         -S samples     Use Samples^2 samples
  1210.   -s             Toggle shadow caching  -T r g b       Contrast threshold
  1211.   -V filename    Verbose file output    -v             Verbose output
  1212. ! -W lx ly hx hy Render subwindow
  1213.   -------------------------------------------------------------------------------
  1214.   
  1215.   File: /* Input file consists of...*/
  1216. --- 40,46 ----
  1217.   -r             Right eye view         -S samples     Use Samples^2 samples
  1218.   -s             Toggle shadow caching  -T r g b       Contrast threshold
  1219.   -V filename    Verbose file output    -v             Verbose output
  1220. ! -W lx hx ly hy Render subwindow       -X l r b t     Crop window
  1221.   -------------------------------------------------------------------------------
  1222.   
  1223.   File: /* Input file consists of...*/
  1224. ***************
  1225. *** 70,76 ****
  1226.           shutter Speed           /* Shutter speed (0 --> no blur) */
  1227.           framelength Length      /* Length of a singelf frame (1) */
  1228.           screen Xsize Ysize      /* Screen size */
  1229. !         window Xmin Ymin Xmax Ymax /* Window (0 0 xsize-1 ysize-1) */
  1230.           eyesep Separation       /* eye separation (0) */
  1231.   
  1232.   SurfDef: /* Give a name to a set of surface attributes. */
  1233. --- 70,77 ----
  1234.           shutter Speed           /* Shutter speed (0 --> no blur) */
  1235.           framelength Length      /* Length of a singelf frame (1) */
  1236.           screen Xsize Ysize      /* Screen size */
  1237. !         window Xmin Xmax Ymin Ymax /* Window (0 xsize-1 0 ysize-1) */
  1238. !         crop left right bot top /* Crop window (0 1 0 1) */
  1239.           eyesep Separation       /* eye separation (0) */
  1240.   
  1241.   SurfDef: /* Give a name to a set of surface attributes. */
  1242. ***************
  1243. *** 199,205 ****
  1244.   Mapping:
  1245.           map uv
  1246.           map cylindrical [Xorigin Yorigin Zorigin Xup Yup Zup Xu Yu Zu] 
  1247. !         map linear      [Xorigin Yorigin Zorigin Xv  Yv  Zv  Xu Yu Zu]
  1248.           map spherical   [Xorigin Yorigin Zorigin Xup Yup Zup Xu Yu Zu]
  1249.   
  1250.   Light:
  1251. --- 200,206 ----
  1252.   Mapping:
  1253.           map uv
  1254.           map cylindrical [Xorigin Yorigin Zorigin Xup Yup Zup Xu Yu Zu] 
  1255. !         map planar      [Xorigin Yorigin Zorigin Xv  Yv  Zv  Xu Yu Zu]
  1256.           map spherical   [Xorigin Yorigin Zorigin Xup Yup Zup Xu Yu Zu]
  1257.   
  1258.   Light:
  1259.  
  1260. Index: libshade/symtab.c
  1261. Prereq: 4.0
  1262. *** libshade/symtab.c.old    1991/07/17 14:48:02
  1263. --- libshade/symtab.c    1991/10/05 18:24:29
  1264. ***************
  1265. *** 13,21 ****
  1266.    * There is no warranty or other guarantee of fitness of this software
  1267.    * for any purpose.  It is provided solely "as is".
  1268.    *
  1269. !  * $Id: symtab.c,v 4.0 91/07/17 14:48:02 kolb Exp Locker: kolb $
  1270.    *
  1271.    * $Log:    symtab.c,v $
  1272.    * Revision 4.0  91/07/17  14:48:02  kolb
  1273.    * Initial version.
  1274.    * 
  1275. --- 13,24 ----
  1276.    * There is no warranty or other guarantee of fitness of this software
  1277.    * for any purpose.  It is provided solely "as is".
  1278.    *
  1279. !  * $Id: symtab.c,v 4.0.1.1 91/10/05 18:24:29 cek Exp Locker: cek $
  1280.    *
  1281.    * $Log:    symtab.c,v $
  1282. +  * Revision 4.0.1.1  91/10/05  18:24:29  cek
  1283. +  * patch1: Added casts for builtin functions.
  1284. +  * 
  1285.    * Revision 4.0  91/07/17  14:48:02  kolb
  1286.    * Initial version.
  1287.    * 
  1288. ***************
  1289. *** 29,42 ****
  1290.       {"pi", 3.141592, NULL, FLOAT_EXPR, FALSE, 0},
  1291.       {"dtor", 0.017453, NULL, FLOAT_EXPR, FALSE, 0},
  1292.       {"rtod", 57.29578, NULL, FLOAT_EXPR, FALSE, 0},
  1293. !     {"cos", 0.0, cos, BUILTIN_EXPR, FALSE, 1},
  1294. !     {"sin", 0.0, sin, BUILTIN_EXPR, FALSE, 1},
  1295. !     {"tan", 0.0, tan, BUILTIN_EXPR, FALSE, 1},
  1296. !     {"sqrt", 0.0, sqrt, BUILTIN_EXPR, FALSE, 1},
  1297. !     {"acos", 0.0, acos, BUILTIN_EXPR, FALSE, 1},
  1298. !     {"asin", 0.0, asin, BUILTIN_EXPR, FALSE, 1},
  1299. !     {"atan", 0.0, atan, BUILTIN_EXPR, FALSE, 1},
  1300. !     {"hypot", 0.0, hypot, BUILTIN_EXPR, FALSE, 2},
  1301.       {"time", 0.0, NULL, FLOAT_EXPR, TRUE, 0},
  1302.       {"frame", 0.0, NULL, FLOAT_EXPR, TRUE, 0},
  1303.       {"linear", 0.0, LinearTime, BUILTIN_EXPR, TRUE, 4},
  1304. --- 32,45 ----
  1305.       {"pi", 3.141592, NULL, FLOAT_EXPR, FALSE, 0},
  1306.       {"dtor", 0.017453, NULL, FLOAT_EXPR, FALSE, 0},
  1307.       {"rtod", 57.29578, NULL, FLOAT_EXPR, FALSE, 0},
  1308. !     {"cos", 0.0, (Float (*)())cos, BUILTIN_EXPR, FALSE, 1},
  1309. !     {"sin", 0.0, (Float (*)())sin, BUILTIN_EXPR, FALSE, 1},
  1310. !     {"tan", 0.0, (Float (*)())tan, BUILTIN_EXPR, FALSE, 1},
  1311. !     {"sqrt", 0.0, (Float (*)())sqrt, BUILTIN_EXPR, FALSE, 1},
  1312. !     {"acos", 0.0, (Float (*)())acos, BUILTIN_EXPR, FALSE, 1},
  1313. !     {"asin", 0.0, (Float (*)())asin, BUILTIN_EXPR, FALSE, 1},
  1314. !     {"atan", 0.0, (Float (*)())atan, BUILTIN_EXPR, FALSE, 1},
  1315. !     {"hypot", 0.0, (Float (*)())hypot, BUILTIN_EXPR, FALSE, 2},
  1316.       {"time", 0.0, NULL, FLOAT_EXPR, TRUE, 0},
  1317.       {"frame", 0.0, NULL, FLOAT_EXPR, TRUE, 0},
  1318.       {"linear", 0.0, LinearTime, BUILTIN_EXPR, TRUE, 4},
  1319.  
  1320. Index: Doc/Guide/running.tex
  1321. *** Doc/Guide/running.tex.old    1991/07/17 15:46:02
  1322. --- Doc/Guide/running.tex    1991/10/05 18:02:29
  1323. ***************
  1324. *** 153,171 ****
  1325.   the windows into a final image.
  1326.   
  1327.   \begin{defkey}{window}{{\em minx maxx miny maxy}}
  1328. !     Render the image in the given screen subwindow,
  1329. !     specified in normalized units.
  1330.   \end{defkey}
  1331.   The window must be properly
  1332.   contained within the screen, i.e., {\em minx} and {\em miny} must
  1333.   be greater than or equal to zero, while {\em maxx} and
  1334. ! {\em maxy} must be less than or equal to one.
  1335. ! The {\em mtv} image file format does not support windows.  The
  1336. ! Utah Raster tool
  1337.   {\em rlecomp\/} is useful for reconstructing the full image from
  1338.   sub-images.
  1339.   By default, the window
  1340. ! is equal to ($0., 1., 0, 1$), or the entire screen.
  1341.   
  1342.   {\em Gamma correction} may also be applied to the three output color
  1343.   channels.  See Appendix A for more details.
  1344. --- 153,184 ----
  1345.   the windows into a final image.
  1346.   
  1347.   \begin{defkey}{window}{{\em minx maxx miny maxy}}
  1348. !     Render the image in the given window.
  1349.   \end{defkey}
  1350.   The window must be properly
  1351.   contained within the screen, i.e., {\em minx} and {\em miny} must
  1352.   be greater than or equal to zero, while {\em maxx} and
  1353. ! {\em maxy} must be less than {\em xsize} and {\em ysize}, respectively.
  1354. ! The Utah Raster tool
  1355.   {\em rlecomp\/} is useful for reconstructing the full image from
  1356.   sub-images.
  1357.   By default, the window
  1358. ! is equivalent to the entire screen.
  1359. ! 
  1360. ! It is also convenient to be able to render a small portion
  1361. ! of the window by specifying a subregion using normalized coordinates.
  1362. ! 
  1363. ! \begin{defkey}{crop}{{\em left right bottom top}}
  1364. !     Crop the rendering window.
  1365. ! \end{defkey}
  1366. ! The rendering window is cropped by rendering the screen
  1367. ! area that falls within $minx + left(maxx - minx)$ and
  1368. ! $minx + right(maxx - minx)$ in the $X$ direction, and similarly
  1369. ! for the $Y$ direction.
  1370. ! {\em Left} and {\em bottom} must be greater than or equal to zero.
  1371. ! {\em Right} and {\em top} must be less than or equal to one.
  1372. ! If {\em left} is greater than {\em right}, the two values are
  1373. ! swapped, and similarly for {\em bottom} and {\em top}.
  1374.   
  1375.   {\em Gamma correction} may also be applied to the three output color
  1376.   channels.  See Appendix A for more details.
  1377.  
  1378. Index: libray/libcommon/rotate.c
  1379. Prereq: 4.0
  1380. *** libray/libcommon/rotate.c.old    1991/07/17 14:31:18
  1381. --- libray/libcommon/rotate.c    1991/09/29 15:35:38
  1382. ***************
  1383. *** 13,21 ****
  1384.    * There is no warranty or other guarantee of fitness of this software
  1385.    * for any purpose.  It is provided solely "as is".
  1386.    *
  1387. !  * $Id: rotate.c,v 4.0 91/07/17 14:31:18 kolb Exp Locker: kolb $
  1388.    *
  1389.    * $Log:    rotate.c,v $
  1390.    * Revision 4.0  91/07/17  14:31:18  kolb
  1391.    * Initial version.
  1392.    * 
  1393. --- 13,24 ----
  1394.    * There is no warranty or other guarantee of fitness of this software
  1395.    * for any purpose.  It is provided solely "as is".
  1396.    *
  1397. !  * $Id: rotate.c,v 4.0.1.1 91/09/29 15:35:38 cek Exp Locker: cek $
  1398.    *
  1399.    * $Log:    rotate.c,v $
  1400. +  * Revision 4.0.1.1  91/09/29  15:35:38  cek
  1401. +  * patch1: Added comments.
  1402. +  * 
  1403.    * Revision 4.0  91/07/17  14:31:18  kolb
  1404.    * Initial version.
  1405.    * 
  1406. ***************
  1407. *** 26,31 ****
  1408. --- 29,37 ----
  1409.   TransMethods *iRotateMethods;
  1410.   void RotationMatrix();
  1411.   
  1412. + /*
  1413. +  * Create and return reference to Rotate structure.
  1414. +  */
  1415.   Rotate *
  1416.   RotateCreate()
  1417.   {
  1418. ***************
  1419. *** 37,42 ****
  1420. --- 43,52 ----
  1421.       return res;
  1422.   }
  1423.   
  1424. + /*
  1425. +  * Return a pointer to collection of methods for the
  1426. +  * Rotate transformation.
  1427. +  */
  1428.   TransMethods *
  1429.   RotateMethods()
  1430.   {
  1431. ***************
  1432. *** 48,53 ****
  1433. --- 58,68 ----
  1434.       return iRotateMethods;    
  1435.   }
  1436.   
  1437. + /*
  1438. +  * Given a Rotate structure and forward and inverse transformations,
  1439. +  * propagate the information in the Rotate structure to the
  1440. +  * transformations.
  1441. +  */
  1442.   void
  1443.   RotatePropagate(rotate, trans, itrans)
  1444.   Rotate *rotate;
  1445. ***************
  1446. *** 62,67 ****
  1447. --- 77,86 ----
  1448.       MatrixInvert(trans, itrans);
  1449.   }
  1450.   
  1451. + /*
  1452. +  * Initialize a rotation matrix given an axis of rotation and an
  1453. +  * angle.  Right-handed rotation is applied.
  1454. +  */
  1455.   void
  1456.   RotationMatrix(x, y, z, theta, trans)
  1457.   Float x, y, z, theta;
  1458.  
  1459. *** End of Patch 1 ***
  1460.  
  1461. exit 0 # Just in case...
  1462. -- 
  1463. Kent Landfield                   INTERNET: kent@sparky.IMD.Sterling.COM
  1464. Sterling Software, IMD           UUCP:     uunet!sparky!kent
  1465. Phone:    (402) 291-8300         FAX:      (402) 291-4362
  1466. Please send comp.sources.misc-related mail to kent@uunet.uu.net.
  1467.