det.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
53
54
55
56
57

#!/usr/bin/awk -f

### det.awk
# determinant via LU
# input: square array as delimited text
# output: (scalar) determinant

BEGIN {
    OFS = "\t"
    sign = "[+-]?"
    decimal = "[0-9]+[.]?[0-9]*"
    fraction = "[.][0-9]*"
    exponent = "([Ee]" sign "[0-9]+)?"
    number = "^" sign "(" decimal "|" fraction ")" exponent "$"
}

NR == 1 {
    for (y=1; y<=NF; y++)
        ($y ~ number) ? header[y] = "col" y : header[y] = $y
}

NF > 0 {
    (NF > nf_max) ? nf_max = NF : nf_max = nf_max

    ### columns
    for (y=1; y<=nf_max; y++) {
        ### rows
        if ($y !~ number) {
            continue
        }
        else {
            read_data[NR,y] = $y
            col_sum[y] += $y
            row_sum[NR] += $y
        }
    }
}

END {
    print length(row_sum), length(col_sum)
    ### columns
    for (y=1; y<=nf_max; y++) {
        if (y in col_sum) {
            ### rows
            for (x=1; x<=NR; x++) {
                if (x in row_sum) {
                    printf("[" OFMT "," OFMT "]" OFS OFMT OFS OFMT,
                        x, y, read_data[x,y], col_sum[y])
                    if (x < nf_max)
                        printf(OFS)
                }
            }
        }
        printf(ORS)
    }
}