home *** CD-ROM | disk | FTP | other *** search
/ Otherware / Otherware_1_SB_Development.iso / mac / developm / source / povsrc.sit / SOURCE / POINT.C < prev    next >
Encoding:
C/C++ Source or Header  |  1992-07-03  |  4.8 KB  |  166 lines
  1. /****************************************************************************
  2. *                point.c
  3. *
  4. *  This module implements the point & spot light source primitive.
  5. *
  6. *  from Persistence of Vision Raytracer 
  7. *  Copyright 1992 Persistence of Vision Team
  8. *---------------------------------------------------------------------------
  9. *  Copying, distribution and legal info is in the file povlegal.doc which
  10. *  should be distributed with this file. If povlegal.doc is not available
  11. *  or for more info please contact:
  12. *
  13. *       Drew Wells [POV-Team Leader] 
  14. *       CIS: 73767,1244  Internet: 73767.1244@compuserve.com
  15. *       Phone: (213) 254-4041
  16. * 
  17. * This program is based on the popular DKB raytracer version 2.12.
  18. * DKBTrace was originally written by David K. Buck.
  19. * DKBTrace Ver 2.0-2.12 were written by David K. Buck & Aaron A. Collins.
  20. *
  21. *****************************************************************************/
  22.  
  23. #include "frame.h"
  24. #include "vector.h"
  25. #include "povproto.h"
  26.  
  27. METHODS Point_Methods =
  28. { Object_Intersect, All_Point_Intersections,
  29.    Inside_Point, NULL,
  30.    Copy_Point,
  31.    Translate_Point, Rotate_Point,
  32.    Scale_Point, Invert_Point};
  33.  
  34. static DBL cubic_spline PARAMS(( DBL low,DBL high,DBL pos));
  35. extern LIGHT_SHAPE *Get_Light_Source_Shape();
  36.  
  37. int All_Point_Intersections (Object, Ray, Depth_Queue)
  38. OBJECT *Object;
  39. RAY *Ray;
  40. PRIOQ *Depth_Queue;
  41. {
  42.  
  43.    return(FALSE);
  44. }
  45.  
  46.  
  47. int Inside_Point (Test_Point, Object)
  48. VECTOR *Test_Point;
  49. OBJECT *Object;
  50. {
  51.    return(FALSE);
  52. }
  53.  
  54.  
  55. void *Copy_Point (Object)
  56. OBJECT *Object;
  57. {
  58.    LIGHT_SHAPE *New_Shape;
  59.  
  60.    New_Shape = Get_Light_Source_Shape ();
  61.    *New_Shape = *((LIGHT_SHAPE *) Object);
  62.    New_Shape -> Next_Object = NULL;
  63.  
  64.    if (New_Shape->Shape_Texture != NULL)
  65.       New_Shape->Shape_Texture = Copy_Texture (New_Shape->Shape_Texture);
  66.  
  67.    return (New_Shape);
  68. }
  69.  
  70. void Translate_Point (Object, Vector)
  71. OBJECT *Object;
  72. VECTOR *Vector;
  73. {
  74.    VAdd (((LIGHT_SHAPE *) Object)->Center, ((LIGHT_SHAPE *) Object)->Center, *Vector);
  75.    VAdd (((LIGHT_SHAPE *) Object)->Points_At,
  76.       ((LIGHT_SHAPE *) Object)->Points_At, *Vector);
  77. }
  78.  
  79. void Rotate_Point (Object, Vector)
  80. OBJECT *Object;
  81. VECTOR *Vector;
  82. {
  83.    TRANSFORMATION Transformation;
  84.    Get_Rotation_Transformation (&Transformation, Vector);
  85.    MTransformVector (&((LIGHT_SHAPE *) Object)->Center,
  86.       &((LIGHT_SHAPE *) Object)->Center, &Transformation);
  87.    MTransformVector (&((LIGHT_SHAPE *) Object)->Points_At,
  88.       &((LIGHT_SHAPE *) Object)->Points_At, &Transformation);
  89. }
  90.  
  91. void Scale_Point (Object, Vector)
  92. OBJECT *Object;
  93. VECTOR *Vector;
  94. {
  95.    TRANSFORMATION Transformation;
  96.    Get_Scaling_Transformation (&Transformation, Vector);
  97.    MTransformVector (&((LIGHT_SHAPE *) Object)->Center,
  98.       &((LIGHT_SHAPE *) Object)->Center, &Transformation);
  99.    MTransformVector (&((LIGHT_SHAPE *) Object)->Points_At,
  100.       &((LIGHT_SHAPE *) Object)->Points_At, &Transformation);
  101.    Scale_Texture (&((LIGHT_SHAPE *) Object)->Shape_Texture, Vector);
  102. }
  103.  
  104. void Invert_Point (Object)
  105. OBJECT *Object;
  106. {
  107.    ((LIGHT_SHAPE *) Object)->Inverted ^= TRUE;
  108. }
  109.  
  110.  
  111. /* Cubic spline that has tangents of slope 0 at x == low and at x == high.
  112.    For a given value "pos" between low and high the spline value is returned */
  113. static DBL cubic_spline(low, high, pos)
  114. DBL low, high, pos;
  115. {
  116.    /* Check to see if the position is within the proper boundaries */
  117.    if (pos < low)
  118.       return 0.0;
  119.    else if (pos > high)
  120.       return 1.0;
  121.    if (high == low)
  122.       return 0.0;
  123.  
  124.    /* Normalize to the interval 0->1 */
  125.    pos = (pos - low) / (high - low);
  126.  
  127.    /* See where it is on the cubic curve */
  128.    return (3 - 2 * pos) * pos * pos;
  129. }
  130.  
  131. DBL Attenuate_Light (Light_Source, Light_Source_Ray)
  132. LIGHT_SHAPE *Light_Source;
  133. RAY *Light_Source_Ray;
  134. {
  135.    DBL Len,costheta;
  136.    DBL Attenuation = 1.0;
  137.    VECTOR Spot_Direction;
  138.  
  139.    /* If this is a spotlight then attenuate based on the incidence angle */
  140.    if (Light_Source->Type == SPOT_LIGHT_TYPE) {
  141.       VSub(Spot_Direction, Light_Source->Points_At, Light_Source->Center);
  142.       VLength(Len, Spot_Direction);
  143.       if (Len > 0.0) {
  144.          VInverseScale(Spot_Direction, Spot_Direction, Len);
  145.          VDot(costheta, Light_Source_Ray->Direction, Spot_Direction);
  146.          costheta *= -1.0;
  147.          if (costheta > 0.0) {
  148.             Attenuation = pow(costheta, Light_Source->Coeff);
  149.             /* If there is a soft falloff region associated with the light then
  150.                do an interpolation of values between the hot center and the
  151.                direction at which light falls to nothing. */
  152.             if (Light_Source->Radius > 0.0)
  153.                Attenuation *= cubic_spline(Light_Source->Falloff,
  154.                   Light_Source->Radius,
  155.                   costheta);
  156.             /* printf("Atten: %lg\n", Attenuation); */
  157.          }
  158.          else
  159.             Attenuation = 0.0;
  160.       }
  161.       else
  162.          Attenuation = 0.0;
  163.    }    
  164.    return(Attenuation);
  165. }    
  166.