blob: 7ce746064692d3b20cde0e3217af260f0b46380a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
%%
%% FILE calend1.tex
%%
\def\loadadvanced{\relax}
% Convert from Julian date and time in
% \date to long date (in local time)
\def\JDTtoL{\advance\date by500
\advance\date by\timezone\divide\date by1000}
%% Trigonometric functions
\def\sintable#1{\ifcase #1 0\or100\or199
\or296\or389\or479\or565\or644\or717
\or783\or841\or891\or932\or964\or985
\or997\or1000\or992\or974\or946\or909
\or863\or808\or746\or675\or598\or516
\or427\or335\or239\or141\or42\or-58
\or-158\fi}
% Reduces modulo 2\pi (requires positive
% argument theta):
% theta := theta MOD 2\pi, where
% theta = count1*10^(-3)
\def\twopimod{\count2 =\count1
\divide\count2 by6284 \count3 =1853
\count4 =6283\multiply\count3 by\count2
\multiply\count4 by\count2
\divide\count3 by10000
\advance\count3 by\count4
\advance\count1 by-\count3}
\newif\ifsign
% v := sin(theta), where
% v = count4*10^(-3);
% theta = count1*10^(-3)
% theta is reduced MOD 2\pi to be
% 0<=theta<2\pi by \TWOPIMOD,
% then linear interpolation is performed
% using \SINTABLE.
\def\Sin{
\ifnum\count1<0 \signtrue
\count1=-\count1\else \signfalse\fi
\loop\ifnum\count1>6284\twopimod\repeat
\ifnum\count1>3142
\advance\count1 by-3142
\ifsign\signfalse\else\signtrue\fi\fi
\multiply\count1 by10 \count3 =\count1
\divide\count3 by1000 \count2 =\count3
\multiply\count3 by1000
\advance\count3 by-\count1
\count5 =\sintable{\count2 }
\count4 =\count5\advance\count2 by1
\advance\count4 by -\sintable{\count2 }
\multiply\count4 by\count3
\divide\count4 by1000
\advance\count4 by\count5
\ifsign\count4 =-\count4 \fi}
% v := cos(theta), where
% v = count4*10^(-3);
% theta = count1*10^(-3).
% \SIN is evaluated on pi/2-theta.
\def\Cos{\advance\count1 by-1571
\multiply\count1 by-1\Sin}
% Linear transformation of T giving
% theta := a*T+b, where
% T = count0*10^(-6);
% theta = count1*10^(-3);
% a = #1*10^3+#2+#3*10^(-3);
% b = #4*10^(-3)
\def\lin#1.#2.#3+#4.{\count1 =#3
\count2 =#2 \count3 =#1
\multiply\count1 by\count0
\multiply\count2 by\count0
\multiply\count3 by\count0
\divide\count1 by1000
\advance\count1 by\count2
\divide\count1 by1000
\advance\count1 by\count3
\advance\count1 by #4}
% Accumulate value returned by a
% trigonometric function, scaled by
% factor f, into count6:
% ac := ac+f*v, where
% ac = count6*10^(-7)
% v = value of SIN or COS =count4*10^(-3)
% f = #1*10^(-4)
\def\fac#1{\multiply\count4 by #1
\advance\count6 by\count4}
\def\id{\count4=\count1} % Identity
|
