home *** CD-ROM | disk | FTP | other *** search
/ C/C++ Users Group Library 1996 July / C-C++ Users Group Library July 1996.iso / listings / v_13_08 / dwyer / trig.c < prev   
Encoding:
C/C++ Source or Header  |  1995-07-10  |  3.1 KB  |  111 lines
  1. /* Listing 2        Fixed Point Trig Functions */
  2.  
  3.  
  4. #include <stdio.h>
  5. #include <stdlib.h>
  6.  
  7.  
  8.  
  9. #define DEGREES_90            0x0400
  10. #define DEGREES_180           0x0800
  11.  
  12.  
  13.  
  14. typedef short Trig;     /*  describes a 2.14 fixed point number */
  15. typedef short Angle;    /*  describes a 12.4 fixed point number */
  16.  
  17.  
  18.  
  19. static Trig cosineTable[ DEGREES_90 + 1 ] =
  20. /*
  21.     This example uses a circle divided uniformly into 4096 "degrees".  The
  22.     cosine values for the first quadrant are stored here.  Only the first
  23.     quadrant is necessary because the lookup routine is able to take advantage
  24.     of cosine's symmetry.
  25.     
  26.     This table was generated by another small program using floating point.
  27.     All it needed to do was compute the cosine every (1/4096)(2*pi) radians
  28.     and convert it into a 2.14 fixed point number.  The values stored here
  29.     are truncated, but some applications may benefit from rounding, instead.
  30.     
  31.     Note that this table contains one more entry than there are divisions in
  32.     a quadrant.  This is to ensure that all angles in the range [0,90] are
  33.     accounted for.
  34. */
  35. {
  36.    0x4000, 0x3fff, 0x3fff, 0x3fff, 0x3fff, 0x3fff, 0x3fff, 0x3fff,
  37.    0x3ffe, 0x3ffe, 0x3ffe, 0x3ffd, 0x3ffd, 0x3ffc, 0x3ffc, 0x3ffb,
  38.  
  39.    ... [full table available on code disk -mb]
  40.  
  41.    0x0192, 0x0178, 0x015f, 0x0146, 0x012d, 0x0114, 0x00fb, 0x00e2,
  42.    0x00c9, 0x00af, 0x0096, 0x007d, 0x0064, 0x004b, 0x0032, 0x0019,
  43.    0x0000
  44. };
  45.  
  46.  
  47.  
  48. Trig cosine( Angle angle )
  49. /*
  50.     This routine returns the cosine of angle, returning a fixed point value in
  51.     the range [-1, 1]. Since this is performed with a simple table lookup, the
  52.     Trig type returned may have any format desired.  This example treats it as
  53.     a 2.14 fixed point number.
  54.     
  55.     The Angle type describes a 12.4 fixed point number.  This divides the unit
  56.     circle into 4096 integral "degrees".  Since every angle may be represented
  57.     in 16 bits, overflow is impossible: the implicit modulus of limited-
  58.     precision integer math guarantees this.
  59. */
  60. {
  61.    Trig value;
  62.  
  63.    /*  Convert angle from 12.4 fixed point to integer. */
  64.    angle >>= 4;
  65.  
  66.    /*  Since cos(x) == cos(-x), combine both cases. */
  67.    if( 0 > angle )
  68.       angle = -angle;
  69.  
  70.    /*  
  71.           Check which quadrant contains angle.  This isn't
  72.           strictly necessary, but by doing so this routine
  73.           is able to exploit the symmetry of the cosine
  74.           function, reducing the size of the lookup table
  75.           by three quarters.
  76.    */
  77.    if( DEGREES_90 < angle )
  78.    {
  79.       /*
  80.             If not in first quadrant, adjust lookup index
  81.           by noting that cos(x) = -cos(180 - x).
  82.       */
  83.       angle = DEGREES_180 - angle;
  84.       value = -cosineTable[ angle ];
  85.    }
  86.    else
  87.       value = cosineTable[ angle ];
  88.  
  89.    return( value );
  90. }
  91.  
  92.  
  93.  
  94. Trig sine( Angle angle )
  95. /*
  96.     This routine returns the sine of angle.  Like the routine above, it
  97.     returns a Trig value, which represents a 2.14 fixed point number.
  98. */
  99. {
  100.    /*
  101.           Make use of the fact that sin(x) = cos(x - 90).
  102.           Note that the constant must be converted to 12.4
  103.        fixed point before it may be subtracted.
  104.    */
  105.    angle -= angle - (DEGREES_90 << 4);
  106.  
  107.    return( cosine( angle ) );
  108. }
  109.  
  110.  
  111.