[ViewVC] Contents of: cvs/aor/aorscan

#!/opt/bin/perl

# http://www.beepworld.de/members81/frequenzen/flughaefen.htm
# http://www.dia.unisa.it/professori/ads/corso-security/www/CORSO-9900/a5/Netsec/netsec.html
# http://www.scnt01426.pwp.blueyonder.co.uk/Articles/RS-232/AR-8000.htm

use List::Util;
use Array::Heap2;

use Fcntl;
use POSIX;
use Time::HiRes;

use Event;

use Coro;
use Coro::Event;
use Coro::Timer qw(sleep);
use Coro::Signal;

use PDL ();
#use PDL::FFTW;
#use PDL::Audio;
#use PDL::IO::Pic;

use PApp::SQL;

use strict;

sub BAUD()       { 9600 }
sub RESOLUTION() { 50 } # in Hz
sub round($)     { int ($_[0] / RESOLUTION + 0.5) * RESOLUTION }

sub SEARCH_INTERVAL_MIN() {  555 }
sub AFT_INTERVAL_MIN()    { 3333 }

$PApp::SQL::DBH = PApp::SQL::connect_cached __FILE__, "DBI:mysql:database=aor;mysql_read_default_file=/root/.my.cnf"
   or die "unable to open mysql database aor";

my $tty = $ARGV[0] || "/dev/ttyS0";

sub AGC_MAX() { 26000 }

my $record_rate  = 48000;
my $record_init  = "amixer sset Mic cap";
my $record_set   = "amixer sset Capture {}%";
my $record_cmd   = "arecord -traw -c1 -fS16_LE -r$record_rate 2>/dev/null";
my $playback_on  = "amixer sset Mic unmute";
my $playback_off = "amixer sset Mic mute";
my $playback_set = "amixer sset Mic {}%";

sub imag($) {
   my ($pdl) = @_;
   wpic $pdl, "/tmp/aor.pgm";
   system "cv", "/tmp/aor.pgm";
}

sub acmd($;$) {
   my ($cmd, $arg) = @_;

   $cmd =~ s/\{\}/$arg/g;
   system "$cmd >/dev/null";
}

my $record_buf;
my $record_vol = 100;

sub record($$) {
   my ($data_cb, $done_cb) = @_;

   $record_buf = "";

   open my $fh, "-|", $record_cmd
      or die "$record_cmd: $!";

   Event->io (fd => $fh, poll => 'r', cb => sub {
      if (sysread $fh, $record_buf, 65536, length $record_buf) {
         if (4096 < length $record_buf) {
            # do the agc adjustment
            my $pdl = PDL::short (unpack "v*", substr $record_buf, -4096);
            if (AGC_MAX < $pdl->abs->max) {
               acmd $record_set, --$record_vol;
#               printf "AGC %d $record_vol\n", max abs $pdl;
            }
         }

         unless ($data_cb->()) {
            $_[0]->w->cancel;
            $done_cb->();
         }
      } else {
         $_[0]->w->cancel;
         $done_cb->();
      }
   });
}

sub record_nsamples($) {
   my $bytes = $_[0] * 2;

   my $done = new Coro::Signal;

   record sub {
      $bytes > length $record_buf
   }, sub {
      $done->send;
   };

   $done->wait;

   substr $record_buf, -$bytes
}

if ($ARGV[0] eq "--addrange") {
   my (undef, $mode, $lo, $hi, $step, $level) = @ARGV;

   my $cnt;
   while () {
      my $f = $lo + $cnt++ * $step;

      last if $f > $hi;

      eval {
         sql_exec "insert into freq (mode, freq, width, activity_level) values (?, ?, ?, ?)",
                  $mode, round $f, $step, $level;
      };
   }

   exit;
} elsif ($ARGV[0] eq "--daemon") {
   # nop
} else {
   # exit
}

my $NOW = time;

Event->timer (after => 0, interval => 1, hard => 1, prio => 3, cb => sub { $NOW = time });

open my $aor, "+<:raw", $tty
   or die "$tty: $!";

{
   my $tio = new POSIX::Termios;

   $tio->getattr (fileno $aor);

   $tio->setiflag (POSIX::IXOFF | POSIX::IXON);
   $tio->setoflag (0);
   $tio->setcflag (POSIX::CLOCAL | POSIX::CREAD | POSIX::CSTOPB | POSIX::CS8);
   $tio->setlflag (0);
   $tio->setispeed (&POSIX::B9600);
   $tio->setospeed (&POSIX::B9600);

   $tio->setattr (fileno $aor);
}

my %mode = (
   0 => "WFM",
   1 => "NFM",
   2 => "AM",
   3 => "USB",
   4 => "LSB",
   5 => "CW",
);

sub cmd($) {
#   print ">$_[0]<\n";#d#
#   fcntl $aor, F_SETFL, 0;
   syswrite $aor, "$_[0]\015\012";
}

sub delay($) {
   my $bytes = int ($_[0] * (BAUD / 11)) - 2;

   " " x ($bytes - 2)
}

cmd "MC1";   # mic off (0 = on, 1 = off, 2 = squelch)
cmd "ST000050"; # minimum stepsize
acmd $record_init;
acmd $record_set, $record_vol;
acmd $playback_off;

my @send_jobs;
my $send_requeue;
my ($curfreq, $curmode);
my $pipeline = 4;

my $rbuf;
my @resp_jobs;

sub send_scheduler {
   while (@send_jobs
          && @resp_jobs < $pipeline
          && (!@resp_jobs || !$resp_jobs[-1]{exclusive})
   ) {
      if ($send_requeue) {
         # re-sort queue
         $send_requeue = 0;

         @send_jobs = sort {
            $b->{nice} <=> $a->{nice}
            || $b->{mode} <=> $a->{mode}
            || $b->{freq} <=> $a->{freq};
         } @send_jobs;
      }

      my $job = pop @send_jobs;

      my ($cmd, @exp);

      if ($curmode != $job->{mode}) {
         $cmd .= "MD$job->{mode}\015\012";
         $curmode = $job->{mode};
      }

      if ($curfreq != $job->{freq}) {
         $cmd .= "RF$job->{freq}\015\012";
         $curfreq = $job->{freq};
      }

      printf "\rMD%d RF%9d %4d %4d         ", $curmode, $curfreq, $job->{nice}, scalar @send_jobs;

      $job->{exp} = [];
      $job->{res} = [];

      for (@{ $job->{cmd} }) {
         push @exp, $1 if s/^(..)=//;
         $cmd .= "$_\015\012";
      }

      $job->{exp} = \@exp;
      $job->{res} = [];

      syswrite $aor, $cmd;

      if (@exp) {
         push @resp_jobs, $job;
      } else {
         $job->{done}->send;
      }
   }
}

Event->io (fd => $aor, prio => 1, poll => 'r', cb => sub {
   sysread $aor, $rbuf, 4096, length $rbuf;

   while ($rbuf =~ s/^([^\015\012]*)[\012\015]+//s) {
      my $line = $1;
      $line =~ s/[\012\015]//g;

      next unless $line =~ /^\S/;

      @resp_jobs or die "out of sync: expected nothing, but got '$line'\n";

      my $exp = shift @{ $resp_jobs[0]{exp} };

      $exp eq substr $line, 0, 2
         or die "sync error: expected '$exp', got '$line'\n";

      push @{ $resp_jobs[0]{res} }, $line;

      unless (@{ $resp_jobs[0]{exp} }) {
         my $job = shift @resp_jobs;

         $job->{exclusive}->() if $job->{exclusive};
         $job->{done}->send;

         send_scheduler;
      }

      next;
   }
});

sub job::result {
   my ($self) = @_;

   $self->{done}->wait;

   wantarray ? @{ $self->{res} } : $self->{res}[0]
}

sub job {
   my ($nice, $mode, $freq, @cmd) = @_;

   my $job = bless {
      freq => round $freq,
      mode => $mode,
      cmd  => \@cmd,
      done => new Coro::Signal,
      nice => $nice,
   }, job::;

   if (@cmd && ref $cmd[-1]) {
      $job->{exclusive} = pop @cmd;
   }

   push @send_jobs, $job;
   $send_requeue = 1;

   send_scheduler;

   $job
}

sub lm {
   map +(hex substr $_, 2) & 0x7f, @_
}

sub aft($$$$) {
   my ($nice, $mode, $center, $radius) = @_;

   $radius = List::Util::max RESOLUTION, $radius * 0.1;

   my @level = map [$_->{freq}, (List::Util::sum lm $_->result) / 5],
                  map +(job $nice, $mode, $center + $radius * $_, (delay 0.020, "LM=LM") x 5),
                     -10 .. 10;

   my ($tune, $weight);
   for (@level) {
      my ($f, $w) = @$_;

      $w **= 8;

      $tune   += $f * $w;
      $weight += $w;
   }

   my $lm = ($weight / @level) ** (1/8);

   (int ($tune / $weight / RESOLUTION) * RESOLUTION, $lm * 100)
}

sub sweep {
   my ($nice, $freqs, $cb) = @_;

   my @jobs = map {
      job $nice, $_->[0], $_->[1], "LM=LM";
   } @$freqs;

   $cb->($freqs->[$_], lm $jobs[$_]->result)
      for 0 .. $#jobs;
}

#############################################################################
# raw search

async {
   my $nice = 1000;

   while () {
      my @freq = sql_fetchall "select mode, freq, width from freq
                               where search_time < ?
                               order by search_time, mode, freq
                               limit 1000",
                              $NOW - SEARCH_INTERVAL_MIN;

      if (@freq) {
         sweep $nice, \@freq, sub {
            my ($info, $lm) = @_;

            sql_exec "update freq set search_time = ?, search_level = ?
                      where mode = ? and freq = ?",
                     $NOW, $lm,
                     $info->[0], $info->[1];
         };
      } else {
         Coro::Timer::sleep 60;
      }
   }
};

#############################################################################
# aft of active freqs

async {
   my $nice = 500;

   while () {
      my @info = sql_fetchall
                        "select mode, freq, width, activity_level
                         from freq
                         where (search_level >= activity_level
                                or aft_time > 0)
                           and aft_time < ?
                         order by mode, freq
                         limit 100",
                        $NOW;

      if (@info) {
         sweep $nice, \@info, sub {
            my ($info, $lm) = @_;
            my ($mode, $freq, $width, $activity_level) = @$info;

            if ($lm >= $activity_level) {
               my ($tune, $weight) = aft $nice, $mode, $freq, $width * 0.9;

               print "$mode $tune($freq) $lm>=$activity_level ";

               if ($freq - 0.5 * $width <= $tune && $tune <= $freq + 0.5 * $width) {
                  print "tuned\n";
                  sql_exec "update freq set aft_time = ?, aft_freq = ?, aft_level = ?
                            where mode = ? and freq = ?",
                           $NOW + 86400*5 - (rand 86400), $tune, $weight,
                           $mode, $freq;
               } else {
                  # outlier
                  print "outlier\n";
                  sql_exec "update freq set aft_time = ?, aft_level = 0
                            where mode = ? and freq = ?",
                           $NOW + 86400*7 - (rand 86400),
                           $mode, $freq;
               }
            } else {
               # currently inactive
               sql_exec "update freq set aft_time = ?
                         where mode = ? and freq = ?",
                        $NOW + 57,
                        $mode, $freq;
            }
         };
      } else {
         Coro::Timer::sleep 10;
      }
   }
};

#############################################################################
# scan active freqs

async {
return;
   my $nice = 100;

   while () {
      my $st = sql_exec \my ($mode, $freq, $width, $activity_level),
                        "select mode, aft_freq, width, activity_level
                         from freq
                         where aft_level > 0 and aft_freq > 0
                         order by freq, mode";

      if ($st->rows) {
         my @jobs;

         while ($st->fetch) {
            push @jobs, [$activity_level, job $nice, $mode, $freq, delay 0.040, "LM=LM"];
         }

         for (@jobs) {
            my ($activity_level, $job) = @$_;

            my @lm = lm $job->result;
            
            if ($lm[0] >= $activity_level) {
               warn "scan $job->{mode} $job->{freq} @lm\n";
               (job $nice-1, $mode, $job->{freq}, "LM=LM", sub {
#               (job $nice-1, 0, 98400000, "LM=LM", sub {
                  acmd $playback_on;
                  Coro::Timer::sleep 1;
#                  my $data = 0.5 + (1 / 65536) * float unpack "v*", record_nsamples $record_rate * 0.1;
                  acmd $playback_off;
                  #$data = sin +(1/480) * xvals zeroes 48000;
#                  $data->reshape (100, $record_rate / 100);
                  #
#                  my $spectrum = 120 + cat map spectrum ($data->slice ("($_),"), "dB", "KAISER"), 0..99;
#                  printf "%s %s\n", (min $spectrum), (max $spectrum);
#                  imag short $spectrum;
               })->result;
            }
         }

         Coro::Timer::sleep 60;
      } else {
         Coro::Timer::sleep 60;
      }
   }
};

#############################################################################

$| = 1;

Coro::Event::loop;

__END__
mysql -e "delete from freq" aor
./aorscan --addrange 1   26565000   27405000 10000               25
./aorscan --addrange 1   34360000   35800000 20000               25
./aorscan --addrange 1   38460000   39840000 20000               25
./aorscan --addrange 1   74215000   77455000 20000               25
./aorscan --addrange 1   84015000   87255000 20000               24
./aorscan --addrange 2  108000000  144000000 25000               20
./aorscan --addrange 1  118000000  136000000  8333.3333333333333 24
./aorscan --addrange 1  108000000  144000000 25000               24
./aorscan --addrange 4  165210000  169380000 20000               24
./aorscan --addrange 3  169810000  173980000 20000               24
./aorscan --addrange 1  230000000  328000000  8333.3333333333333 10
./aorscan --addrange 2  230000000  328000000 50000               20
./aorscan --addrange 1  438650000  439425000 25000               10
./aorscan --addrange 4  443600000  444962500 12500               10
./aorscan --addrange 3  448600000  449962500 12500               10
./aorscan --addrange 1 1270200000 1270700000 25000               10
./aorscan --addrange 1 1298200000 1298700000 25000               10

./aorscan --addrange 3    2690000    2690000 25000               10
./aorscan --addrange 3    3413000    3413000 25000               10
./aorscan --addrange 4    3413000    3413000 25000               10
./aorscan --addrange 3    5640000    5640000 25000               10
./aorscan --addrange 4    5640000    5640000 25000               10
./aorscan --addrange 3    8957000    8957000 25000               10
./aorscan --addrange 4    8957000    8957000 25000               10
./aorscan --addrange 3   13264000   13264000 25000               10
./aorscan --addrange 4   13264000   13264000 25000               10

Revision:	1.4
Committed:	Sun Oct 2 09:36:37 2005 UTC (18 years, 8 months ago) by root
Branch:	MAIN
Changes since 1.3:	+53 -34 lines
Log Message:	* empty log message *
#	Content
1	#!/opt/bin/perl
2
3	# http://www.beepworld.de/members81/frequenzen/flughaefen.htm
4	# http://www.dia.unisa.it/professori/ads/corso-security/www/CORSO-9900/a5/Netsec/netsec.html
5	# http://www.scnt01426.pwp.blueyonder.co.uk/Articles/RS-232/AR-8000.htm
6
7	use List::Util;
8	use Array::Heap2;
9
10	use Fcntl;
11	use POSIX;
12	use Time::HiRes;
13
14	use Event;
15
16	use Coro;
17	use Coro::Event;
18	use Coro::Timer qw(sleep);
19	use Coro::Signal;
20
21	use PDL ();
22	#use PDL::FFTW;
23	#use PDL::Audio;
24	#use PDL::IO::Pic;
25
26	use PApp::SQL;
27
28	use strict;
29
30	sub BAUD() { 9600 }
31	sub RESOLUTION() { 50 } # in Hz
32	sub round($) { int ($_[0] / RESOLUTION + 0.5) * RESOLUTION }
33
34	sub SEARCH_INTERVAL_MIN() { 555 }
35	sub AFT_INTERVAL_MIN() { 3333 }
36
37	$PApp::SQL::DBH = PApp::SQL::connect_cached __FILE__, "DBI:mysql:database=aor;mysql_read_default_file=/root/.my.cnf"
38	or die "unable to open mysql database aor";
39
40	my $tty = $ARGV[0] \|\| "/dev/ttyS0";
41
42	sub AGC_MAX() { 26000 }
43
44	my $record_rate = 48000;
45	my $record_init = "amixer sset Mic cap";
46	my $record_set = "amixer sset Capture {}%";
47	my $record_cmd = "arecord -traw -c1 -fS16_LE -r$record_rate 2>/dev/null";
48	my $playback_on = "amixer sset Mic unmute";
49	my $playback_off = "amixer sset Mic mute";
50	my $playback_set = "amixer sset Mic {}%";
51
52	sub imag($) {
53	my ($pdl) = @_;
54	wpic $pdl, "/tmp/aor.pgm";
55	system "cv", "/tmp/aor.pgm";
56	}
57
58	sub acmd($;$) {
59	my ($cmd, $arg) = @_;
60
61	$cmd =~ s/\{\}/$arg/g;
62	system "$cmd >/dev/null";
63	}
64
65	my $record_buf;
66	my $record_vol = 100;
67
68	sub record($$) {
69	my ($data_cb, $done_cb) = @_;
70
71	$record_buf = "";
72
73	open my $fh, "-\|", $record_cmd
74	or die "$record_cmd: $!";
75
76	Event->io (fd => $fh, poll => 'r', cb => sub {
77	if (sysread $fh, $record_buf, 65536, length $record_buf) {
78	if (4096 < length $record_buf) {
79	# do the agc adjustment
80	my $pdl = PDL::short (unpack "v*", substr $record_buf, -4096);
81	if (AGC_MAX < $pdl->abs->max) {
82	acmd $record_set, --$record_vol;
83	# printf "AGC %d $record_vol\n", max abs $pdl;
84	}
85	}
86
87	unless ($data_cb->()) {
88	$_[0]->w->cancel;
89	$done_cb->();
90	}
91	} else {
92	$_[0]->w->cancel;
93	$done_cb->();
94	}
95	});
96	}
97
98	sub record_nsamples($) {
99	my $bytes = $_[0] * 2;
100
101	my $done = new Coro::Signal;
102
103	record sub {
104	$bytes > length $record_buf
105	}, sub {
106	$done->send;
107	};
108
109	$done->wait;
110
111	substr $record_buf, -$bytes
112	}
113
114	if ($ARGV[0] eq "--addrange") {
115	my (undef, $mode, $lo, $hi, $step, $level) = @ARGV;
116
117	my $cnt;
118	while () {
119	my $f = $lo + $cnt++ * $step;
120
121	last if $f > $hi;
122
123	eval {
124	sql_exec "insert into freq (mode, freq, width, activity_level) values (?, ?, ?, ?)",
125	$mode, round $f, $step, $level;
126	};
127	}
128
129	exit;
130	} elsif ($ARGV[0] eq "--daemon") {
131	# nop
132	} else {
133	# exit
134	}
135
136	my $NOW = time;
137
138	Event->timer (after => 0, interval => 1, hard => 1, prio => 3, cb => sub { $NOW = time });
139
140	open my $aor, "+<:raw", $tty
141	or die "$tty: $!";
142
143	{
144	my $tio = new POSIX::Termios;
145
146	$tio->getattr (fileno $aor);
147
148	$tio->setiflag (POSIX::IXOFF \| POSIX::IXON);
149	$tio->setoflag (0);
150	$tio->setcflag (POSIX::CLOCAL \| POSIX::CREAD \| POSIX::CSTOPB \| POSIX::CS8);
151	$tio->setlflag (0);
152	$tio->setispeed (&POSIX::B9600);
153	$tio->setospeed (&POSIX::B9600);
154
155	$tio->setattr (fileno $aor);
156	}
157
158	my %mode = (
159	0 => "WFM",
160	1 => "NFM",
161	2 => "AM",
162	3 => "USB",
163	4 => "LSB",
164	5 => "CW",
165	);
166
167	sub cmd($) {
168	# print ">$_[0]<\n";#d#
169	# fcntl $aor, F_SETFL, 0;
170	syswrite $aor, "$_[0]\015\012";
171	}
172
173	sub delay($) {
174	my $bytes = int ($_[0] * (BAUD / 11)) - 2;
175
176	" " x ($bytes - 2)
177	}
178
179	cmd "MC1"; # mic off (0 = on, 1 = off, 2 = squelch)
180	cmd "ST000050"; # minimum stepsize
181	acmd $record_init;
182	acmd $record_set, $record_vol;
183	acmd $playback_off;
184
185	my @send_jobs;
186	my $send_requeue;
187	my ($curfreq, $curmode);
188	my $pipeline = 4;
189
190	my $rbuf;
191	my @resp_jobs;
192
193	sub send_scheduler {
194	while (@send_jobs
195	&& @resp_jobs < $pipeline
196	&& (!@resp_jobs \|\| !$resp_jobs[-1]{exclusive})
197	) {
198	if ($send_requeue) {
199	# re-sort queue
200	$send_requeue = 0;
201
202	@send_jobs = sort {
203	$b->{nice} <=> $a->{nice}
204	\|\| $b->{mode} <=> $a->{mode}
205	\|\| $b->{freq} <=> $a->{freq};
206	} @send_jobs;
207	}
208
209	my $job = pop @send_jobs;
210
211	my ($cmd, @exp);
212
213	if ($curmode != $job->{mode}) {
214	$cmd .= "MD$job->{mode}\015\012";
215	$curmode = $job->{mode};
216	}
217
218	if ($curfreq != $job->{freq}) {
219	$cmd .= "RF$job->{freq}\015\012";
220	$curfreq = $job->{freq};
221	}
222
223	printf "\rMD%d RF%9d %4d %4d ", $curmode, $curfreq, $job->{nice}, scalar @send_jobs;
224
225	$job->{exp} = [];
226	$job->{res} = [];
227
228	for (@{ $job->{cmd} }) {
229	push @exp, $1 if s/^(..)=//;
230	$cmd .= "$_\015\012";
231	}
232
233	$job->{exp} = \@exp;
234	$job->{res} = [];
235
236	syswrite $aor, $cmd;
237
238	if (@exp) {
239	push @resp_jobs, $job;
240	} else {
241	$job->{done}->send;
242	}
243	}
244	}
245
246	Event->io (fd => $aor, prio => 1, poll => 'r', cb => sub {
247	sysread $aor, $rbuf, 4096, length $rbuf;
248
249	while ($rbuf =~ s/^([^\015\012]*)[\012\015]+//s) {
250	my $line = $1;
251	$line =~ s/[\012\015]//g;
252
253	next unless $line =~ /^\S/;
254
255	@resp_jobs or die "out of sync: expected nothing, but got '$line'\n";
256
257	my $exp = shift @{ $resp_jobs[0]{exp} };
258
259	$exp eq substr $line, 0, 2
260	or die "sync error: expected '$exp', got '$line'\n";
261
262	push @{ $resp_jobs[0]{res} }, $line;
263
264	unless (@{ $resp_jobs[0]{exp} }) {
265	my $job = shift @resp_jobs;
266
267	$job->{exclusive}->() if $job->{exclusive};
268	$job->{done}->send;
269
270	send_scheduler;
271	}
272
273	next;
274	}
275	});
276
277	sub job::result {
278	my ($self) = @_;
279
280	$self->{done}->wait;
281
282	wantarray ? @{ $self->{res} } : $self->{res}[0]
283	}
284
285	sub job {
286	my ($nice, $mode, $freq, @cmd) = @_;
287
288	my $job = bless {
289	freq => round $freq,
290	mode => $mode,
291	cmd => \@cmd,
292	done => new Coro::Signal,
293	nice => $nice,
294	}, job::;
295
296	if (@cmd && ref $cmd[-1]) {
297	$job->{exclusive} = pop @cmd;
298	}
299
300	push @send_jobs, $job;
301	$send_requeue = 1;
302
303	send_scheduler;
304
305	$job
306	}
307
308	sub lm {
309	map +(hex substr $_, 2) & 0x7f, @_
310	}
311
312	sub aft($$$$) {
313	my ($nice, $mode, $center, $radius) = @_;
314
315	$radius = List::Util::max RESOLUTION, $radius * 0.1;
316
317	my @level = map [$_->{freq}, (List::Util::sum lm $_->result) / 5],
318	map +(job $nice, $mode, $center + $radius * $_, (delay 0.020, "LM=LM") x 5),
319	-10 .. 10;
320
321	my ($tune, $weight);
322	for (@level) {
323	my ($f, $w) = @$_;
324
325	$w **= 8;
326
327	$tune += $f * $w;
328	$weight += $w;
329	}
330
331	my $lm = ($weight / @level) ** (1/8);
332
333	(int ($tune / $weight / RESOLUTION) * RESOLUTION, $lm * 100)
334	}
335
336	sub sweep {
337	my ($nice, $freqs, $cb) = @_;
338
339	my @jobs = map {
340	job $nice, $_->[0], $_->[1], "LM=LM";
341	} @$freqs;
342
343	$cb->($freqs->[$_], lm $jobs[$_]->result)
344	for 0 .. $#jobs;
345	}
346
347	#############################################################################
348	# raw search
349
350	async {
351	my $nice = 1000;
352
353	while () {
354	my @freq = sql_fetchall "select mode, freq, width from freq
355	where search_time < ?
356	order by search_time, mode, freq
357	limit 1000",
358	$NOW - SEARCH_INTERVAL_MIN;
359
360	if (@freq) {
361	sweep $nice, \@freq, sub {
362	my ($info, $lm) = @_;
363
364	sql_exec "update freq set search_time = ?, search_level = ?
365	where mode = ? and freq = ?",
366	$NOW, $lm,
367	$info->[0], $info->[1];
368	};
369	} else {
370	Coro::Timer::sleep 60;
371	}
372	}
373	};
374
375	#############################################################################
376	# aft of active freqs
377
378	async {
379	my $nice = 500;
380
381	while () {
382	my @info = sql_fetchall
383	"select mode, freq, width, activity_level
384	from freq
385	where (search_level >= activity_level
386	or aft_time > 0)
387	and aft_time < ?
388	order by mode, freq
389	limit 100",
390	$NOW;
391
392	if (@info) {
393	sweep $nice, \@info, sub {
394	my ($info, $lm) = @_;
395	my ($mode, $freq, $width, $activity_level) = @$info;
396
397	if ($lm >= $activity_level) {
398	my ($tune, $weight) = aft $nice, $mode, $freq, $width * 0.9;
399
400	print "$mode $tune($freq) $lm>=$activity_level ";
401
402	if ($freq - 0.5 * $width <= $tune && $tune <= $freq + 0.5 * $width) {
403	print "tuned\n";
404	sql_exec "update freq set aft_time = ?, aft_freq = ?, aft_level = ?
405	where mode = ? and freq = ?",
406	$NOW + 86400*5 - (rand 86400), $tune, $weight,
407	$mode, $freq;
408	} else {
409	# outlier
410	print "outlier\n";
411	sql_exec "update freq set aft_time = ?, aft_level = 0
412	where mode = ? and freq = ?",
413	$NOW + 86400*7 - (rand 86400),
414	$mode, $freq;
415	}
416	} else {
417	# currently inactive
418	sql_exec "update freq set aft_time = ?
419	where mode = ? and freq = ?",
420	$NOW + 57,
421	$mode, $freq;
422	}
423	};
424	} else {
425	Coro::Timer::sleep 10;
426	}
427	}
428	};
429
430	#############################################################################
431	# scan active freqs
432
433	async {
434	return;
435	my $nice = 100;
436
437	while () {
438	my $st = sql_exec \my ($mode, $freq, $width, $activity_level),
439	"select mode, aft_freq, width, activity_level
440	from freq
441	where aft_level > 0 and aft_freq > 0
442	order by freq, mode";
443
444	if ($st->rows) {
445	my @jobs;
446
447	while ($st->fetch) {
448	push @jobs, [$activity_level, job $nice, $mode, $freq, delay 0.040, "LM=LM"];
449	}
450
451	for (@jobs) {
452	my ($activity_level, $job) = @$_;
453
454	my @lm = lm $job->result;
455
456	if ($lm[0] >= $activity_level) {
457	warn "scan $job->{mode} $job->{freq} @lm\n";
458	(job $nice-1, $mode, $job->{freq}, "LM=LM", sub {
459	# (job $nice-1, 0, 98400000, "LM=LM", sub {
460	acmd $playback_on;
461	Coro::Timer::sleep 1;
462	# my $data = 0.5 + (1 / 65536) * float unpack "v", record_nsamples $record_rate 0.1;
463	acmd $playback_off;
464	#$data = sin +(1/480) * xvals zeroes 48000;
465	# $data->reshape (100, $record_rate / 100);
466	#
467	# my $spectrum = 120 + cat map spectrum ($data->slice ("($_),"), "dB", "KAISER"), 0..99;
468	# printf "%s %s\n", (min $spectrum), (max $spectrum);
469	# imag short $spectrum;
470	})->result;
471	}
472	}
473
474	Coro::Timer::sleep 60;
475	} else {
476	Coro::Timer::sleep 60;
477	}
478	}
479	};
480
481	#############################################################################
482
483	$\| = 1;
484
485	Coro::Event::loop;
486
487	__END__
488	mysql -e "delete from freq" aor
489	./aorscan --addrange 1 26565000 27405000 10000 25
490	./aorscan --addrange 1 34360000 35800000 20000 25
491	./aorscan --addrange 1 38460000 39840000 20000 25
492	./aorscan --addrange 1 74215000 77455000 20000 25
493	./aorscan --addrange 1 84015000 87255000 20000 24
494	./aorscan --addrange 2 108000000 144000000 25000 20
495	./aorscan --addrange 1 118000000 136000000 8333.3333333333333 24
496	./aorscan --addrange 1 108000000 144000000 25000 24
497	./aorscan --addrange 4 165210000 169380000 20000 24
498	./aorscan --addrange 3 169810000 173980000 20000 24
499	./aorscan --addrange 1 230000000 328000000 8333.3333333333333 10
500	./aorscan --addrange 2 230000000 328000000 50000 20
501	./aorscan --addrange 1 438650000 439425000 25000 10
502	./aorscan --addrange 4 443600000 444962500 12500 10
503	./aorscan --addrange 3 448600000 449962500 12500 10
504	./aorscan --addrange 1 1270200000 1270700000 25000 10
505	./aorscan --addrange 1 1298200000 1298700000 25000 10
506
507	./aorscan --addrange 3 2690000 2690000 25000 10
508	./aorscan --addrange 3 3413000 3413000 25000 10
509	./aorscan --addrange 4 3413000 3413000 25000 10
510	./aorscan --addrange 3 5640000 5640000 25000 10
511	./aorscan --addrange 4 5640000 5640000 25000 10
512	./aorscan --addrange 3 8957000 8957000 25000 10
513	./aorscan --addrange 4 8957000 8957000 25000 10
514	./aorscan --addrange 3 13264000 13264000 25000 10
515	./aorscan --addrange 4 13264000 13264000 25000 10
516