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
icossin(int deg,int * cosp,int * sinp)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 }
141