Commit | Line | Data |
---|---|---|
7af5ad59 PH |
1 | / |
2 | / | |
3 | / D S I N | |
4 | / - - - | |
5 | / | |
6 | /SUBROUTINE DSIN(X) | |
7 | / | |
8 | / VERSION 5A 4-27-77 PT | |
9 | / | |
10 | SECT DSIN | |
11 | JA #DSIN | |
12 | DPCHK | |
13 | TEXT +DSIN + | |
14 | SINXR, SETX XRSIN | |
15 | SETB BPSIN | |
16 | FPI2SN, 1 /PI DIVIDED BY 2 | |
17 | 3110 | |
18 | 3755 | |
19 | 2421 | |
20 | 0264 | |
21 | 3016 | |
22 | FPISIN, 2 /PI | |
23 | 3110 | |
24 | 3755 | |
25 | 2421 | |
26 | 0264 | |
27 | 3016 | |
28 | F2PISN, 3 /TWO PI | |
29 | 3110 | |
30 | 3755 | |
31 | 2421 | |
32 | 0264 | |
33 | 3016 | |
34 | BPSIN, F 0.0 | |
35 | XRSIN, F 0.0 | |
36 | X, F 0.0 | |
37 | F 0.0 | |
38 | ORG 10*3+BPSIN | |
39 | FNOP | |
40 | JA SINXR | |
41 | 0 | |
42 | SINRTN, JA . | |
43 | SIN1, F 0.0 | |
44 | F 0.0 | |
45 | F3PIB2, 0003 /4.71238898038468986 | |
46 | 2266 | |
47 | 2761 | |
48 | 7714 | |
49 | 6207 | |
50 | 2212 | |
51 | F1SIN, F 1.0 | |
52 | F 0.0 | |
53 | / | |
54 | SINC17, 7720 /1/17! | |
55 | 3124 | |
56 | 5435 | |
57 | 6014 | |
58 | 1265 | |
59 | 1236 | |
60 | SINC15, 7730 /-1/15! | |
61 | 4506 | |
62 | 0060 | |
63 | 3063 | |
64 | 0437 | |
65 | 4133 | |
66 | SINC13, 7740 /1/13! | |
67 | 2604 | |
68 | 4430 | |
69 | 2352 | |
70 | 0664 | |
71 | 1151 | |
72 | SINC11, 7747 /-1/11! | |
73 | 4506 | |
74 | 3352 | |
75 | 3002 | |
76 | 5354 | |
77 | 3710 | |
78 | SINC9, 7756 /1/9! | |
79 | 2707 | |
80 | 3616 | |
81 | 4525 | |
82 | 5434 | |
83 | 6350 | |
84 | / | |
85 | SINC7, 7764 /-1/7! | |
86 | 4577 | |
87 | 4577 | |
88 | 4577 | |
89 | 4577 | |
90 | 4636 | |
91 | / | |
92 | SINC5, 7772 /1/5! | |
93 | 2104 | |
94 | 2104 | |
95 | 2104 | |
96 | 2104 | |
97 | 2104 | |
98 | / | |
99 | SINC3, 7776 /-1/3! | |
100 | 5252 | |
101 | 5252 | |
102 | 5252 | |
103 | 5252 | |
104 | 5244 | |
105 | / | |
106 | SINTST, 7770 | |
107 | 2000 | |
108 | 0000 | |
109 | 0000 | |
110 | 0000 | |
111 | 0000 | |
112 | SEVTHR, 0007 | |
113 | 2217 | |
114 | 7777 | |
115 | 7777 | |
116 | 7777 | |
117 | 7777 | |
118 | BASE 0 | |
119 | #DSIN, STARTD | |
120 | FLDA 10*3 | |
121 | FSTA SINRTN | |
122 | FLDA 0 | |
123 | SETX XRSIN | |
124 | SETB BPSIN | |
125 | BASE BPSIN | |
126 | LDX 1,1 | |
127 | LDX 73,2 | |
128 | FSTA BPSIN | |
129 | FLDA% BPSIN,1 /ADDR OF X | |
130 | FSTA BPSIN | |
131 | STARTE | |
132 | FLDA% BPSIN /GET X | |
133 | LDX -1,0 /SET SIGN TO POSITIVE. | |
134 | JGT SINMOD /IF POSITIVE BYPASS FUDGE. | |
135 | JEQ SINRTN /IF ZERO EXIT. | |
136 | FNEG /NEGATIVE. NEGATE AC. SIN(-X)=-SIN(X) | |
137 | LDX 0,0 /SET SIGN TO MINUS. | |
138 | SINMOD, FSTA X | |
139 | FDIV F2PISN /X/2PI | |
140 | ALN 2 | |
141 | FNORM | |
142 | FMUL F2PISN /*2PI | |
143 | FSTA SIN1 | |
144 | FLDA X | |
145 | FSUB SIN1 | |
146 | FSTA SIN1 | |
147 | FSTA X | |
148 | /CHECK FOR QUADRANT | |
149 | /1 0-PI/2 SIN(PI/2)=1 | |
150 | /2 PI/2-PI SIN(PI)=0 | |
151 | /3 PI-3PI/2 SIN(3PI/2)=-1 | |
152 | /4 3PI/2-2PI SIN(2PI)=0 | |
153 | / | |
154 | FLDA FPI2SN /PI/2 | |
155 | FSUB SIN1 | |
156 | JGT TAYLOR | |
157 | JEQ SPIB2 /=PI/2 | |
158 | FLDA FPISIN /PI | |
159 | FSUB SIN1 | |
160 | JLT S1 | |
161 | JEQ SPI /=PI | |
162 | FLDA FPISIN | |
163 | FSUB SIN1 /SIN(X)=SIN(PI-X) | |
164 | FSTA X | |
165 | JA TAYLOR | |
166 | S1, FLDA F3PIB2 /3PI/2 | |
167 | FSUB SIN1 | |
168 | JLT S2 | |
169 | JEQ S3PIB2 /=3PI/2 | |
170 | FLDA SIN1 | |
171 | FSUB FPISIN | |
172 | FNEG | |
173 | FSTA X /SIN(X)=-SIN(X-PI) | |
174 | JA TAYLOR | |
175 | S2, FLDA F2PISN /2PI | |
176 | FSUB SIN1 | |
177 | JLT DSNER /ERROR | |
178 | JEQ SPI | |
179 | FNEG | |
180 | FSTA X | |
181 | /CALCULATE SIN VIA TAYLOR SERIES | |
182 | TAYLOR, FLDA X /RECALL NUMBER TO BE WORKED ON. | |
183 | FMUL X /MULTIPLY OUT. | |
184 | FSTA SIN1 | |
185 | FMUL SINC17 /NOW DO THE STANDARD ITERATION. | |
186 | FADD SINC15 | |
187 | FMUL SIN1 | |
188 | FADD SINC13 | |
189 | FMUL SIN1 | |
190 | FADD SINC11 | |
191 | FMUL SIN1 | |
192 | FADD SINC9 | |
193 | FMUL SIN1 | |
194 | FADD SINC7 | |
195 | FMUL SIN1 | |
196 | FADD SINC5 | |
197 | FMUL SIN1 | |
198 | FADD SINC3 | |
199 | FMUL SIN1 | |
200 | FADD F1SIN /ADD IN 1 FOR SERIES | |
201 | FMUL X /DO THE FINAL MULTIPLY. | |
202 | SING, JXN SINRTN,0 /SHALL WE NEGATE | |
203 | FNEG /YEP | |
204 | JA SINRTN /AND RETURN. | |
205 | SPIB2, FLDA F1SIN | |
206 | JA SING | |
207 | SPI, FCLA | |
208 | JA SINRTN | |
209 | S3PIB2, FLDA F1SIN | |
210 | FNEG | |
211 | JA SING | |
212 | EXTERN #ARGER | |
213 | DSNER, TRAP4 #ARGER | |
214 | \f |