home *** CD-ROM | disk | FTP | other *** search
/ Otherware / Otherware_1_SB_Development.iso / mac / developm / source / povsrc.sit / SOURCE / VECTOR.H < prev    next >
Encoding:
C/C++ Source or Header  |  1992-07-03  |  2.7 KB  |  69 lines
  1. /****************************************************************************
  2. *                vector.h
  3. *
  4. *  This module contains macros to perform operations on vectors.
  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.  
  24. /* Misc. Vector Math Macro Definitions */
  25.  
  26. extern DBL VTemp;
  27.  
  28. /* Vector Add */
  29. #define VAdd(a, b, c) {(a).x=(b).x+(c).x;(a).y=(b).y+(c).y;(a).z=(b).z+(c).z;}
  30.  
  31. /* Vector Subtract */
  32. #define VSub(a, b, c) {(a).x=(b).x-(c).x;(a).y=(b).y-(c).y;(a).z=(b).z-(c).z;}
  33.  
  34. /* Scale - Multiply Vector by a Scalar */
  35. #define VScale(a, b, k) {(a).x=(b).x*(k);(a).y=(b).y*(k);(a).z=(b).z*(k);}
  36.  
  37. /* Inverse Scale - Divide Vector by a Scalar */
  38. #define VInverseScale(a, b, k) {(a).x=(b).x/(k);(a).y=(b).y/(k);(a).z=(b).z/(k);}
  39.  
  40. /* Dot Product - Gives Scalar angle (a) between two vectors (b) and (c) */
  41. #define VDot(a, b, c) {a=(b).x*(c).x+(b).y*(c).y+(b).z*(c).z;}
  42.  
  43. /* Cross Product - returns Vector (a) = (b) x (c) 
  44.    WARNING:  a must be different from b and c.*/
  45. #define VCross(a,b,c) {(a).x=(b).y*(c).z-(b).z*(c).y; \
  46. (a).y=(b).z*(c).x-(b).x*(c).z; \
  47.                        (a).z=(b).x*(c).y-(b).y*(c).x;}
  48.  
  49. /* Evaluate - returns Vector (a) = Multiply Vector (b) by Vector (c) */
  50. #define VEvaluate(a, b, c) {(a).x=(b).x*(c).x;(a).y=(b).y*(c).y;(a).z=(b).z*(c).z;}
  51.  
  52. /* Square a Vector */
  53. #define    VSqr(a) {(a).x*(a).x;(a).y*(a).y;(a).z*(a).z;}
  54.  
  55. /* Simple Scalar Square Macro */
  56. #define    Sqr(a)    ((a)*(a))
  57.  
  58. /* Square a Vector (b) and Assign to another Vector (a) */
  59. #define VSquareTerms(a, b) {(a).x=(b).x*(b).x;(a).y=(b).y*(b).y;(a).z=(b).z*(b).z;}
  60.  
  61. /* Vector Length - returs Scalar Euclidean Length (a) of Vector (b) */
  62. #define VLength(a, b) {a=sqrt((b).x*(b).x+(b).y*(b).y+(b).z*(b).z);}
  63.  
  64. /* Normalize a Vector - returns a vector (length of 1) that points at (b) */
  65. #define VNormalize(a,b) {VTemp=sqrt((b).x*(b).x+(b).y*(b).y+(b).z*(b).z);(a).x=(b).x/VTemp;(a).y=(b).y/VTemp;(a).z=(b).z/VTemp;}
  66.  
  67. /* Compute a Vector (a) Halfway Between Two Given Vectors (b) and (c) */
  68. #define VHalf(a, b, c) {(a).x=0.5*((b).x+(c).x);(a).y=0.5*((b).y+(c).y);(a).z=0.5*((b).z+(c).z);}
  69.