Warning: Attempt to read property "date" on null in /usr/local/www/websvn.planix.org/blame.php on line 247

Warning: Attempt to read property "msg" on null in /usr/local/www/websvn.planix.org/blame.php on line 247
WebSVN – planix.SVN – Blame – /os/branches/feature_posix/sys/src/cmd/unix/drawterm/libdraw/icossin2.c – Rev 2

Subversion Repositories planix.SVN

Rev

Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 - 1 
#include	<u.h>
2 
#include	<libc.h>
3 
#include	<draw.h>
4 
 
5 
/*
6 
 * Sine and Cosine of arctangents, calculated by
7 
 *   (sin(atan(index/100.0))*1024.+0.5)
8 
 *   (cos(atan(index/100.0))*1024.+0.5)
9 
 * To use, get rational tangent between 0<=tan<=1, scale by 100,
10 
 * and look up sin and cos, and use linear interpolation.  divide by 1024.
11 
 * Maximum error is 0.0020.  Without linear interpolation, it's 0.010.
12 
 */
13 
static
14 
short sinus[] = {
15 
	0,	/* 0.00 */
16 
	10,	/* 0.01 */
17 
	20,	/* 0.02 */
18 
	31,	/* 0.03 */
19 
	41,	/* 0.04 */
20 
	51,	/* 0.05 */
21 
	61,	/* 0.06 */
22 
	72,	/* 0.07 */
23 
	82,	/* 0.08 */
24 
	92,	/* 0.09 */
25 
	102,	/* 0.10 */
26 
	112,	/* 0.11 */
27 
	122,	/* 0.12 */
28 
	132,	/* 0.13 */
29 
	142,	/* 0.14 */
30 
	152,	/* 0.15 */
31 
	162,	/* 0.16 */
32 
	172,	/* 0.17 */
33 
	181,	/* 0.18 */
34 
	191,	/* 0.19 */
35 
	201,	/* 0.20 */
36 
	210,	/* 0.21 */
37 
	220,	/* 0.22 */
38 
	230,	/* 0.23 */
39 
	239,	/* 0.24 */
40 
	248,	/* 0.25 */
41 
	258,	/* 0.26 */
42 
	267,	/* 0.27 */
43 
	276,	/* 0.28 */
44 
	285,	/* 0.29 */
45 
	294,	/* 0.30 */
46 
	303,	/* 0.31 */
47 
	312,	/* 0.32 */
48 
	321,	/* 0.33 */
49 
	330,	/* 0.34 */
50 
	338,	/* 0.35 */
51 
	347,	/* 0.36 */
52 
	355,	/* 0.37 */
53 
	364,	/* 0.38 */
54 
	372,	/* 0.39 */
55 
	380,	/* 0.40 */
56 
	388,	/* 0.41 */
57 
	397,	/* 0.42 */
58 
	405,	/* 0.43 */
59 
	412,	/* 0.44 */
60 
	420,	/* 0.45 */
61 
	428,	/* 0.46 */
62 
	436,	/* 0.47 */
63 
	443,	/* 0.48 */
64 
	451,	/* 0.49 */
65 
	458,	/* 0.50 */
66 
	465,	/* 0.51 */
67 
	472,	/* 0.52 */
68 
	480,	/* 0.53 */
69 
	487,	/* 0.54 */
70 
	493,	/* 0.55 */
71 
	500,	/* 0.56 */
72 
	507,	/* 0.57 */
73 
	514,	/* 0.58 */
74 
	520,	/* 0.59 */
75 
	527,	/* 0.60 */
76 
	533,	/* 0.61 */
77 
	540,	/* 0.62 */
78 
	546,	/* 0.63 */
79 
	552,	/* 0.64 */
80 
	558,	/* 0.65 */
81 
	564,	/* 0.66 */
82 
	570,	/* 0.67 */
83 
	576,	/* 0.68 */
84 
	582,	/* 0.69 */
85 
	587,	/* 0.70 */
86 
	593,	/* 0.71 */
87 
	598,	/* 0.72 */
88 
	604,	/* 0.73 */
89 
	609,	/* 0.74 */
90 
	614,	/* 0.75 */
91 
	620,	/* 0.76 */
92 
	625,	/* 0.77 */
93 
	630,	/* 0.78 */
94 
	635,	/* 0.79 */
95 
	640,	/* 0.80 */
96 
	645,	/* 0.81 */
97 
	649,	/* 0.82 */
98 
	654,	/* 0.83 */
99 
	659,	/* 0.84 */
100 
	663,	/* 0.85 */
101 
	668,	/* 0.86 */
102 
	672,	/* 0.87 */
103 
	676,	/* 0.88 */
104 
	681,	/* 0.89 */
105 
	685,	/* 0.90 */
106 
	689,	/* 0.91 */
107 
	693,	/* 0.92 */
108 
	697,	/* 0.93 */
109 
	701,	/* 0.94 */
110 
	705,	/* 0.95 */
111 
	709,	/* 0.96 */
112 
	713,	/* 0.97 */
113 
	717,	/* 0.98 */
114 
	720,	/* 0.99 */
115 
	724,	/* 1.00 */
116 
	728,	/* 1.01 */
117 
};
118 
 
119 
static
120 
short cosinus[] = {
121 
	1024,	/* 0.00 */
122 
	1024,	/* 0.01 */
123 
	1024,	/* 0.02 */
124 
	1024,	/* 0.03 */
125 
	1023,	/* 0.04 */
126 
	1023,	/* 0.05 */
127 
	1022,	/* 0.06 */
128 
	1022,	/* 0.07 */
129 
	1021,	/* 0.08 */
130 
	1020,	/* 0.09 */
131 
	1019,	/* 0.10 */
132 
	1018,	/* 0.11 */
133 
	1017,	/* 0.12 */
134 
	1015,	/* 0.13 */
135 
	1014,	/* 0.14 */
136 
	1013,	/* 0.15 */
137 
	1011,	/* 0.16 */
138 
	1010,	/* 0.17 */
139 
	1008,	/* 0.18 */
140 
	1006,	/* 0.19 */
141 
	1004,	/* 0.20 */
142 
	1002,	/* 0.21 */
143 
	1000,	/* 0.22 */
144 
	998,	/* 0.23 */
145 
	996,	/* 0.24 */
146 
	993,	/* 0.25 */
147 
	991,	/* 0.26 */
148 
	989,	/* 0.27 */
149 
	986,	/* 0.28 */
150 
	983,	/* 0.29 */
151 
	981,	/* 0.30 */
152 
	978,	/* 0.31 */
153 
	975,	/* 0.32 */
154 
	972,	/* 0.33 */
155 
	969,	/* 0.34 */
156 
	967,	/* 0.35 */
157 
	963,	/* 0.36 */
158 
	960,	/* 0.37 */
159 
	957,	/* 0.38 */
160 
	954,	/* 0.39 */
161 
	951,	/* 0.40 */
162 
	947,	/* 0.41 */
163 
	944,	/* 0.42 */
164 
	941,	/* 0.43 */
165 
	937,	/* 0.44 */
166 
	934,	/* 0.45 */
167 
	930,	/* 0.46 */
168 
	927,	/* 0.47 */
169 
	923,	/* 0.48 */
170 
	920,	/* 0.49 */
171 
	916,	/* 0.50 */
172 
	912,	/* 0.51 */
173 
	909,	/* 0.52 */
174 
	905,	/* 0.53 */
175 
	901,	/* 0.54 */
176 
	897,	/* 0.55 */
177 
	893,	/* 0.56 */
178 
	890,	/* 0.57 */
179 
	886,	/* 0.58 */
180 
	882,	/* 0.59 */
181 
	878,	/* 0.60 */
182 
	874,	/* 0.61 */
183 
	870,	/* 0.62 */
184 
	866,	/* 0.63 */
185 
	862,	/* 0.64 */
186 
	859,	/* 0.65 */
187 
	855,	/* 0.66 */
188 
	851,	/* 0.67 */
189 
	847,	/* 0.68 */
190 
	843,	/* 0.69 */
191 
	839,	/* 0.70 */
192 
	835,	/* 0.71 */
193 
	831,	/* 0.72 */
194 
	827,	/* 0.73 */
195 
	823,	/* 0.74 */
196 
	819,	/* 0.75 */
197 
	815,	/* 0.76 */
198 
	811,	/* 0.77 */
199 
	807,	/* 0.78 */
200 
	804,	/* 0.79 */
201 
	800,	/* 0.80 */
202 
	796,	/* 0.81 */
203 
	792,	/* 0.82 */
204 
	788,	/* 0.83 */
205 
	784,	/* 0.84 */
206 
	780,	/* 0.85 */
207 
	776,	/* 0.86 */
208 
	773,	/* 0.87 */
209 
	769,	/* 0.88 */
210 
	765,	/* 0.89 */
211 
	761,	/* 0.90 */
212 
	757,	/* 0.91 */
213 
	754,	/* 0.92 */
214 
	750,	/* 0.93 */
215 
	746,	/* 0.94 */
216 
	742,	/* 0.95 */
217 
	739,	/* 0.96 */
218 
	735,	/* 0.97 */
219 
	731,	/* 0.98 */
220 
	728,	/* 0.99 */
221 
	724,	/* 1.00 */
222 
	720,	/* 1.01 */
223 
};
224 
 
225 
void
226 
icossin2(int x, int y, int *cosp, int *sinp)
227 
{
228 
	int sinsign, cossign, tan, tan10, rem;
229 
	short *stp, *ctp;
230 
 
231 
	if(x == 0){
232 
		if(y >= 0)
233 
			*sinp = ICOSSCALE, *cosp = 0;
234 
		else
235 
			*sinp = -ICOSSCALE, *cosp = 0;
236 
		return;
237 
	}
238 
	sinsign = cossign = 1;
239 
	if(x < 0){
240 
		cossign = -1;
241 
		x = -x;
242 
	}
243 
	if(y < 0){
244 
		sinsign = -1;
245 
		y = -y;
246 
	}
247 
	if(y > x){
248 
		tan = 1000*x/y;
249 
		tan10 = tan/10;
250 
		stp = &cosinus[tan10];
251 
		ctp = &sinus[tan10];
252 
	}else{
253 
		tan = 1000*y/x;
254 
		tan10 = tan/10;
255 
		stp = &sinus[tan10];
256 
		ctp = &cosinus[tan10];
257 
	}
258 
	rem = tan-(tan10*10);
259 
	*sinp = sinsign*(stp[0]+(stp[1]-stp[0])*rem/10);
260 
	*cosp = cossign*(ctp[0]+(ctp[1]-ctp[0])*rem/10);
261 
}