diff1.awk


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

#!/usr/bin/awk -f

### diff1.awk
# print numerical diff along columns

BEGIN {
    OFS = FS
    # OFMT = "%.9g"
    sign = "[+-]?"
    decimal = "[0-9]+[.]?[0-9]*"
    fraction = "[.][0-9]*"
    exponent = "([Ee]" sign "[0-9]+)?"
    number = "^" sign "(" decimal "|" fraction ")" exponent "$"
}

NR == 1 {
    # orig data columns
    for (n=1; n<=NF; n++) {
        if ($n !~ number) {
            header[n] = $n
            dheader[n] = "d" $n
        }
        else {
            header[n] = "col" n
            dheader[n] = "dcol" n
        }
        printf(dheader[n])
        printf(n < NF ? OFS : ORS)
    }
}

NF {
    if (NF > nf_max)
        nf_max = NF

    # data columns
    for (n=1; n<=nf_max; n++) {
        if ($n == header[n] || $n == dheader[n])
            continue
        if ($n ~ number) {
            data[n] = $n
            if (data_prev[n] ~ number)
                diff[n] = data[n] - data_prev[n]
            else
                diff[n] = ""
            data_prev[n] = data[n]
        }
        printf(diff[n])
        printf(n < nf_max ? OFS : ORS)
    }
}