sys/src/libdraw/icossin.c

   1 #include        <u.h>
   2 #include        <libc.h>
   3 #include        <draw.h>
   4
   5 /*
   6  * Integer sine and cosine for integral degree argument.
   7  * Tables computed by (sin,cos)(PI*d/180).
   8  */
   9 static short sinus[91] = {
  10         0,      /* 0 */
  11         18,     /* 1 */
  12         36,     /* 2 */
  13         54,     /* 3 */
  14         71,     /* 4 */
  15         89,     /* 5 */
  16         107,    /* 6 */
  17         125,    /* 7 */
  18         143,    /* 8 */
  19         160,    /* 9 */
  20         178,    /* 10 */
  21         195,    /* 11 */
  22         213,    /* 12 */
  23         230,    /* 13 */
  24         248,    /* 14 */
  25         265,    /* 15 */
  26         282,    /* 16 */
  27         299,    /* 17 */
  28         316,    /* 18 */
  29         333,    /* 19 */
  30         350,    /* 20 */
  31         367,    /* 21 */
  32         384,    /* 22 */
  33         400,    /* 23 */
  34         416,    /* 24 */
  35         433,    /* 25 */
  36         449,    /* 26 */
  37         465,    /* 27 */
  38         481,    /* 28 */
  39         496,    /* 29 */
  40         512,    /* 30 */
  41         527,    /* 31 */
  42         543,    /* 32 */
  43         558,    /* 33 */
  44         573,    /* 34 */
  45         587,    /* 35 */
  46         602,    /* 36 */
  47         616,    /* 37 */
  48         630,    /* 38 */
  49         644,    /* 39 */
  50         658,    /* 40 */
  51         672,    /* 41 */
  52         685,    /* 42 */
  53         698,    /* 43 */
  54         711,    /* 44 */
  55         724,    /* 45 */
  56         737,    /* 46 */
  57         749,    /* 47 */
  58         761,    /* 48 */
  59         773,    /* 49 */
  60         784,    /* 50 */
  61         796,    /* 51 */
  62         807,    /* 52 */
  63         818,    /* 53 */
  64         828,    /* 54 */
  65         839,    /* 55 */
  66         849,    /* 56 */
  67         859,    /* 57 */
  68         868,    /* 58 */
  69         878,    /* 59 */
  70         887,    /* 60 */
  71         896,    /* 61 */
  72         904,    /* 62 */
  73         912,    /* 63 */
  74         920,    /* 64 */
  75         928,    /* 65 */
  76         935,    /* 66 */
  77         943,    /* 67 */
  78         949,    /* 68 */
  79         956,    /* 69 */
  80         962,    /* 70 */
  81         968,    /* 71 */
  82         974,    /* 72 */
  83         979,    /* 73 */
  84         984,    /* 74 */
  85         989,    /* 75 */
  86         994,    /* 76 */
  87         998,    /* 77 */
  88         1002,   /* 78 */
  89         1005,   /* 79 */
  90         1008,   /* 80 */
  91         1011,   /* 81 */
  92         1014,   /* 82 */
  93         1016,   /* 83 */
  94         1018,   /* 84 */
  95         1020,   /* 85 */
  96         1022,   /* 86 */
  97         1023,   /* 87 */
  98         1023,   /* 88 */
  99         1024,   /* 89 */
 100         1024,   /* 90 */
 101 };
 102
 103 void
 104 icossin(int deg, int *cosp, int *sinp)
 105 {
 106         int sinsign, cossign;
 107         short *stp, *ctp;
 108
 109         deg %= 360;
 110         if(deg < 0)
 111                 deg += 360;
 112         sinsign = 1;
 113         cossign = 1;
 114         stp = 0;
 115         ctp = 0;
 116         switch(deg/90){
 117         case 2:
 118                 sinsign = -1;
 119                 cossign = -1;
 120                 deg -= 180;
 121                 /* fall through */
 122         case 0:
 123                 stp = &sinus[deg];
 124                 ctp = &sinus[90-deg];
 125                 break;
 126         case 3:
 127                 sinsign = -1;
 128                 cossign = -1;
 129                 deg -= 180;
 130                 /* fall through */
 131         case 1:
 132                 deg = 180-deg;
 133                 cossign = -cossign;
 134                 stp = &sinus[deg];
 135                 ctp = &sinus[90-deg];
 136                 break;
 137         }
 138         *sinp = sinsign*stp[0];
 139         *cosp = cossign*ctp[0];
 140 }