blob: 9513e9919466cbfad44ae52ed1636aac87f797ac (
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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
%!ps
% PostScript prologue for pst-coil.tex.
% Version 1.07, 2015/05/13
% For distribution, see pstricks.tex.
%
/tx@CoilDict 40 dict def tx@CoilDict begin
/CoilLoop {
/t ED
t sin AspectSin mul t 180 div AspectCos mul add t cos lineto } def
/Coil {
/Inc ED
dup sin /AspectSin ED
cos /AspectCos ED
/ArmB ED
/ArmA ED
/h ED /w ED /y1 ED /x1 ED /y0 ED /x0 ED
x0 y0 translate y1 y0 sub x1 x0 sub 2 copy Pyth
/TotalLength ED
Atan rotate
/BeginAngle ArmA AspectCos Div w h mul Div 360 mul def
/EndAngle TotalLength ArmB sub AspectCos Div
w h mul Div 360 mul def
1 0 0 0 ArrowA ArmA 0 lineto
/mtrx CM def
w h mul 2 Div w 2 Div scale BeginAngle Inc
EndAngle { CoilLoop } for
EndAngle CoilLoop mtrx setmatrix TotalLength ArmB sub 0 lineto CP
% DG/SR modification begin - Jun. 2, 1998 - Patch 1 (from Michael Vulis)
% TotalLength 0 ArrowB lineto } def
TotalLength 0 ArrowB lineto
pop pop pop pop } def
% DG/SR modification end
%
/Sin {
25 dict begin
/Func ED
/PointsPerPeriod ED
/ArmB ED
/ArmA ED
/Amplitude ED
/roundValue ED
/relativePeriod ED
/Periods ED
/y1 ED /x1 ED /y0 ED /x0 ED
y1 y0 sub x1 x0 sub 2 copy Pyth
dup /TotalLength ED
ArmA sub ArmB sub /SinLength ED
Atan /RotAngle ED
relativePeriod
{}
{ SinLength Periods div
roundValue dup 0 eq { pop cvi } { 1 eq { round } if } ifelse
/Periods ED } ifelse
/dxStep SinLength 360 div def % the default step for x pos
/xStep 360 PointsPerPeriod div def % the step for the for loop
/mtrx CM def
x0 y0 translate
RotAngle rotate
1 0 0 0 ArrowA ArmA 0 lineto
1 xStep 360 { % PointsPerPeriod times
dup dxStep mul ArmA add exch % x pos and value for the func
Periods mul Func Amplitude mul lineto
} for
TotalLength ArmB sub 0 lineto CP
TotalLength 0 ArrowB lineto
mtrx setmatrix
pop pop pop pop
end
} def
%
/AltCoil {
/Inc ED
dup sin /AspectSin ED
cos /AspectCos ED /h ED /w ED /EndAngle ED /BeginAngle ED
/mtrx CM def
w h mul 2 Div w 2 Div scale BeginAngle sin AspectSin mul
BeginAngle 180 div AspectCos mul add BeginAngle cos
/lineto load stopped { moveto } if
BeginAngle Inc EndAngle { CoilLoop } for
EndAngle CoilLoop mtrx setmatrix } def
/ZigZag { 15 dict begin
/ArmB ED /ArmA ED
2 div /w ED
w mul /h ED /y1 ED /x1 ED /y0 ED /x0 ED
x1 y1 translate y0 y1 sub x0 x1 sub 2 copy Pyth
/TotalLength ED
Atan rotate TotalLength ArmA sub ArmB sub dup h div cvi /n ED n h mul
sub 2 div dup ArmA add
/ArmA ED ArmB add
/ArmB ED
/x ArmB h 2 div add def
mark 0 0 ArmB 0
n { x w /w w neg def /x x h add def } repeat
TotalLength ArmA sub 0 TotalLength 0
end } def
%
/ZigZagCirc { % x0 y0 x1 y1 bow h w ArmA ArmB on stack
/bow ED
/ArmB ED /ArmA ED
2 div /w ED
w mul /h ED /y1 ED /x1 ED /y0 ED /x0 ED
x1 y1 translate %
y0 y1 sub x0 x1 sub 2 copy Pyth 2 div
/HalfLength ED
Atan /ang ED %angle of A relative to B
%rotate % so end B is origin and BA is horizontal, A to right
/theta bow abs HalfLength Atan 2 mul def %halfangular sector for arc
/theta2 theta 2 mul def % full arc
%/psi ang 90 sub theta sub def % angle of B from center if bow>0
/rho HalfLength theta sin div def % radius of circle
bow 0 gt { /direc 1 def /thetaB ang 90 sub theta sub def } % travel B to A in pos drn
{ /direc -1 def /thetaB ang 90 add theta add def } ifelse % travel B to A in neg drn
%thetaB=angle from center to B
rho thetaB 180 add PtoC translate % origin now at center
/h h rho div RadtoDeg def
/ArmA ArmA rho div RadtoDeg def /ArmB ArmB rho div RadtoDeg def
theta2 ArmA sub ArmB sub dup h div cvi /n ED n h mul
%HalfLength 2 mul ArmA sub ArmB sub dup h div cvi /n ED n h mul
sub 2 div dup ArmA add
/ArmA ED ArmB add /ArmB ED
/h h direc mul def % h is now angle increment for half-period
/rhoo rho w add def /rhoi rho w sub def % outer and inner radii
/arg thetaB ArmB direc mul add def % argument for 2nd point
% thetaB = direc = bow = HalfLength = theta = rho = arg = h = n =
mark rho thetaB PtoC rho arg PtoC /arg arg h 2 div add def
n { w 0 gt { rhoo }{ rhoi } ifelse arg PtoC /w w neg def /arg arg h add def } repeat
rho thetaB theta2 ArmA sub direc mul add PtoC rho thetaB theta2 direc mul add PtoC
} def
%
end
% END pst-coil.pro
|
