108 tgt.
x = (tmp *
vrt());
109 tgt.
y = (tmp *
vup());
110 tgt.
z = (tmp *
vpn());
119 double az = atan(tgt.
x/tgt.
z);
120 double el = atan(tgt.
y/tgt.
z);
130 double deflection =
vrt().
y;
131 while (fabs(deflection) > 0.001) {
132 double theta = asin(deflection/
vrt().length());
135 deflection =
vrt().
y;
152 tgt.
x = (tmp *
vrt());
153 tgt.
y = (tmp *
vup());
154 tgt.
z = (tmp *
vpn());
159 tgt.
x = (tmp *
vrt());
160 tgt.
y = (tmp *
vup());
161 tgt.
z = (tmp *
vpn());
168 double az = atan(tgt.
x/tgt.
z);
175 while (az >
PI) az -= 2*
PI;
176 while (az < -
PI) az += 2*
PI;
183 else if (az > alimit) {
192 tgt.
x = (tmp *
vrt());
193 tgt.
y = (tmp *
vup());
194 tgt.
z = (tmp *
vpn());
196 double el = atan(tgt.
y/tgt.
z);
198 if (e_lo > 0 && el < -e_lo) {
203 else if (e_hi > 0 && el > e_hi) {