blob: 470b27c8bef595adbe191a7dc761b24021e97ce0 (
plain) (
blame)
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
|
# awk program to scan loopstats files and report errors/statistics
#
# usage: awk -f loop.awk loopstats
#
# format of loopstats record
# MJD sec time (s) freq (ppm) tc
# 49235 3.943 0.000016 22.4716 0
#
BEGIN {
loop_tmax = loop_fmax = -1e9
loop_tmin = loop_fmin = 1e9
}
#
# scan all records in file
#
{
if (NF >= 5) {
loop_count++
if ($3 > loop_tmax)
loop_tmax = $3
if ($3 < loop_tmin)
loop_tmin = $3
if ($4 > loop_fmax)
loop_fmax = $4
if ($4 < loop_fmin)
loop_fmin = $4
loop_time += $3
loop_time_rms += $3 * $3
loop_freq += $4
loop_freq_rms += $4 * $4
}
} END {
if (loop_count > 0) {
loop_time /= loop_count
loop_time_rms = sqrt(loop_time_rms / loop_count - loop_time * loop_time)
loop_freq /= loop_count
loop_freq_rms = sqrt(loop_freq_rms / loop_count - loop_freq * loop_freq)
printf "loop %d, %.0f+/-%.1f, rms %.1f, freq %.2f+/-%0.3f, var %.3f\n", loop_count, (loop_tmax + loop_tmin) / 2 * 1e6, (loop_tmax - loop_tmin) / 2 * 1e6, loop_time_rms * 1e6, (loop_fmax + loop_fmin) / 2, (loop_fmax - loop_fmin) / 2, loop_freq_rms
}
}
|
