home *** CD-ROM | disk | FTP | other *** search
/ Amiga Tools 3 / Amiga Tools 3.iso / grafik / raytracing / rayshade-4.0.6.3 / fixes / fix024 / libshade / shade.c.diff < prev    next >
Encoding:
Text File  |  1994-08-09  |  5.8 KB  |  188 lines
  1. *** shade.c    Mon Oct 11 14:59:56 1993
  2. --- shade.c.frac    Mon Oct 11 14:46:03 1993
  3. ***************
  4. *** 27,32 ****
  5. --- 27,33 ----
  6.   #include "libsurf/atmosphere.h"
  7.   #include "options.h"
  8.   #include "stats.h"
  9. + #include "viewing.h"
  10.   
  11.   Medium    TopMedium;
  12.   Atmosphere *AtmosEffects;
  13. ***************
  14. *** 119,125 ****
  15.        */
  16.       for (lp = Lights; lp; lp = lp->next)
  17.           LightRay(lp, pos, nrm, gnrm, smooth, &refl, surf,
  18. !                 ray->depth, ray->sample, ray->time, color);
  19.   
  20.       if (ray->depth >= Options.maxdepth)
  21.           /*
  22. --- 120,126 ----
  23.        */
  24.       for (lp = Lights; lp; lp = lp->next)
  25.           LightRay(lp, pos, nrm, gnrm, smooth, &refl, surf,
  26. !                 ray, color);
  27.   
  28.       if (ray->depth >= Options.maxdepth)
  29.           /*
  30. ***************
  31. *** 162,185 ****
  32.    * Lighting calculations
  33.    */
  34.   static void
  35. ! LightRay(lp, pos, norm, gnorm, smooth, reflect, surf, depth, samp, time, color)
  36.   Light *lp;            /* Light source */
  37.   Vector *pos, *norm, *gnorm;    /* hit pos, shade norm, geo norm */
  38.   int smooth;            /* true if shade and geo norm differ */
  39.   Vector *reflect;        /* reflection direction */
  40.   Surface *surf;            /* surface characteristics */
  41. ! int depth, samp;        /* ray depth, sample # */
  42. ! Float time;
  43.   Color *color;            /* resulting color */
  44.   {
  45.       Color lcolor;
  46.       Ray newray;
  47. !     Float costheta, cosalpha, dist;
  48.   
  49.       newray.pos = *pos;
  50. !     newray.depth = depth;
  51. !     newray.sample = samp;
  52. !     newray.time = time; 
  53.       newray.media = (Medium *)NULL;    
  54.   
  55.       LightDirection(lp, pos, &newray.dir, &dist);
  56. --- 163,186 ----
  57.    * Lighting calculations
  58.    */
  59.   static void
  60. ! LightRay(lp, pos, norm, gnorm, smooth, reflect, surf, ray, color)
  61.   Light *lp;            /* Light source */
  62.   Vector *pos, *norm, *gnorm;    /* hit pos, shade norm, geo norm */
  63.   int smooth;            /* true if shade and geo norm differ */
  64.   Vector *reflect;        /* reflection direction */
  65.   Surface *surf;            /* surface characteristics */
  66. ! Ray *ray;
  67.   Color *color;            /* resulting color */
  68.   {
  69.       Color lcolor;
  70.       Ray newray;
  71. !     Float costheta, cosalpha, cosphi, dist;
  72.   
  73. +         newray.type = SHADOW_RAY;
  74.       newray.pos = *pos;
  75. !     newray.depth = ray->depth;
  76. !     newray.sample = ray->sample;
  77. !     newray.time = ray->time; 
  78.       newray.media = (Medium *)NULL;    
  79.   
  80.       LightDirection(lp, pos, &newray.dir, &dist);
  81. ***************
  82. *** 186,191 ****
  83. --- 187,202 ----
  84.   
  85.       costheta = dotp(&newray.dir, norm);
  86.   
  87. + /* ray viewpoint is the position at the light source */
  88. +         VecAddScaled(newray.pos, dist, newray.dir, &newray.viewpoint);
  89. + /* shadow ray "sees" same area at intersection position as parent ray */
  90. +         cosphi = -dotp(&ray->dir, norm);
  91. +         if (cosphi < EPSILON) cosphi = EPSILON;
  92. +         newray.width = PixelSize(ray, VecDist(&newray.pos, &ray->pos)) / cosphi * costheta;
  93. + /* shadow rays become smaller as they approach the light source */
  94. +         newray.viewdist = -dist;   
  95. +         newray.stretch = 1.;
  96. + 
  97.       if (smooth) {
  98.           cosalpha = dotp(&newray.dir, gnorm); 
  99.           /*
  100. ***************
  101. *** 278,287 ****
  102.   {
  103.       int total_int_refl = FALSE;
  104.       Ray NewRay;
  105. !     Float dist;
  106.       Color newcol;
  107.       HitList hittmp;        /* Geom intersection record */
  108.   
  109.       NewRay.pos = *pos;        /* Origin == hit point */
  110.       NewRay.media = ray->media;    /* Media == old media */
  111.       NewRay.sample = ray->sample;
  112. --- 289,299 ----
  113.   {
  114.       int total_int_refl = FALSE;
  115.       Ray NewRay;
  116. !     Float dist, cosi, cosr;
  117.       Color newcol;
  118.       HitList hittmp;        /* Geom intersection record */
  119.   
  120. +         NewRay.type = REFRACT_RAY;
  121.       NewRay.pos = *pos;        /* Origin == hit point */
  122.       NewRay.media = ray->media;    /* Media == old media */
  123.       NewRay.sample = ray->sample;
  124. ***************
  125. *** 322,327 ****
  126. --- 334,360 ----
  127.        */
  128.   
  129.       if (!total_int_refl) {
  130. + /* good for planar refraction only */
  131. + /* cosi = cos angle(incident ray, surface normal)
  132. +  * cosr = cos angle(refracted ray, surface normal) */
  133. +                 cosi = -dotp(&ray->dir, norm);                
  134. +         if (cosi < EPSILON) cosi = EPSILON;
  135. +                 cosr = -dotp(&NewRay.dir, norm);                
  136. + /*
  137. +  *                       nr   cosr  2
  138. +  * NewRay.viewdist = y = -- ( ---- )  * x 
  139. +  *                       ni   cosi
  140. +  * where x = distance newray origin to oldray viewpoint.
  141. +  */
  142. +                 NewRay.viewdist = (NewRay.media ? NewRay.media->index : TopMedium.index) / (ray->media ? ray->media->index : TopMedium.index) * 
  143. +                                 (cosr * cosr) / (cosi * cosi) * 
  144. +                                 (ray->viewdist < 0. ? -1. : 1.) * VecDist(&ray->viewpoint, &NewRay.pos);
  145. + /* NewRay.viewpoint = NewRay.pos - y * NewRay.dir */
  146. +                 VecAddScaled(NewRay.pos, -NewRay.viewdist, NewRay.dir, &NewRay.viewpoint);
  147. + /* both rays "see" same area at intersection point */
  148. +                 NewRay.width = PixelSize(ray, VecDist(&ray->pos, &NewRay.pos)) * cosr / cosi;
  149. +                 NewRay.stretch = 1.;
  150. + 
  151.           Stats.RefractRays++;
  152.           hittmp.nodes = 0;
  153.           dist = FAR_AWAY;
  154. ***************
  155. *** 360,365 ****
  156. --- 393,399 ----
  157.       Color newcol;
  158.       Float dist;
  159.   
  160. +         NewRay.type = REFLECT_RAY;
  161.       NewRay.pos = *pos;        /* Origin == hit point */
  162.       NewRay.dir = *dir;        /* Direction == reflection */
  163.       NewRay.media = ray->media;    /* Medium == old medium */
  164. ***************
  165. *** 366,371 ****
  166. --- 400,413 ----
  167.       NewRay.sample = ray->sample;
  168.       NewRay.time = ray->time;
  169.       NewRay.depth = ray->depth + 1;
  170. + 
  171. + /* good for planar reflections only */
  172. + /* new viewpoint is mirrored old one. width at intersection point is the same */
  173. +         NewRay.viewdist = (ray->viewdist < 0. ? -1. : 1.) * VecDist(&ray->viewpoint, &NewRay.pos);
  174. +         NewRay.width = PixelSize(ray, VecDist(&ray->pos, &NewRay.pos));
  175. +         VecAddScaled(NewRay.pos, -NewRay.viewdist, NewRay.dir, &NewRay.viewpoint);
  176. +         NewRay.stretch = 1.;
  177. + 
  178.       Stats.ReflectRays++;
  179.       hittmp.nodes = 0;
  180.       dist = FAR_AWAY;
  181. ***************
  182. *** 374,376 ****
  183. --- 416,419 ----
  184.       ColorMultiply(newcol, *intens, &newcol);
  185.       ColorAdd(*color, newcol, color);
  186.   }
  187. + 
  188.