summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat
-rw-r--r--conv.awk2
-rw-r--r--diff.awk6
-rw-r--r--diff1.awk4
-rw-r--r--hamming.awk2
-rw-r--r--kalman.awk4
-rw-r--r--kalman_init.awk10
-rw-r--r--lin_reg.awk4
-rwxr-xr-xlin_reg1.awk4
-rw-r--r--lin_reg2.awk4
-rw-r--r--mean.awk1
-rw-r--r--mean_avg.awk1
-rw-r--r--sum1.awk4
-rw-r--r--sum2.awk2
-rw-r--r--sum3.awk4
-rw-r--r--sum4.awk8
15 files changed, 31 insertions, 29 deletions
diff --git a/conv.awk b/conv.awk
index b1e9e62..f77cd67 100644
--- a/conv.awk
+++ b/conv.awk
@@ -4,7 +4,7 @@
# [PoC] linear convolution (with hardcoded IR window).
BEGIN {
- OFMT = "%.9g"
+ OFMT = "%.18g"
X = ARGV[1]
input_size = split(X, X_arr)
#H = "1.0 1.0 1.0" # rect
diff --git a/diff.awk b/diff.awk
index 03bc398..c79a5c3 100644
--- a/diff.awk
+++ b/diff.awk
@@ -4,7 +4,7 @@
# numerical diff along columns
BEGIN {
- OFMT = "%.9g"
+ OFMT = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -44,7 +44,7 @@ NF {
if ($y == header[y] || $y == dheader[y])
continue
if ($y ~ number)
- printf("%.12g", $y)
+ printf("%.18g", $y)
else
printf("")
printf(OFS)
@@ -58,7 +58,7 @@ NF {
data[y] = $y
(data_prev[y] ~ number) ? diff[y] = data[y] - data_prev[y] : diff[y] = ""
data_prev[y] = data[y]
- printf("%.12g", diff[y])
+ printf("%.18g", diff[y])
}
else
diff[y] = ""
diff --git a/diff1.awk b/diff1.awk
index e2c4a99..4132bf7 100644
--- a/diff1.awk
+++ b/diff1.awk
@@ -4,7 +4,7 @@
# numerical diff along columns
BEGIN {
- OFMT = "%.9g"
+ OFMT = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -41,7 +41,7 @@ NF {
data[y] = $y
(data_prev[y] ~ number) ? diff[y] = data[y] - data_prev[y] : diff[y] = ""
data_prev[y] = data[y]
- printf("%.12g", diff[y])
+ printf("%.18g", diff[y])
}
else
diff[y] = ""
diff --git a/hamming.awk b/hamming.awk
index 84dacd3..b048299 100644
--- a/hamming.awk
+++ b/hamming.awk
@@ -6,7 +6,7 @@
# 'a0' and 'a1' parameters of the raised cosine.
BEGIN {
- OMFT = "%.9g"
+ OMFT = "%.18g"
N = ARGV[1]
pi = 4*(4*atan2(1,5) - atan2(1,239))
diff --git a/kalman.awk b/kalman.awk
index 41936e2..b951d11 100644
--- a/kalman.awk
+++ b/kalman.awk
@@ -4,7 +4,7 @@
# simple linear regression between columns
BEGIN {
- OMFT = "%.9g"
+ OMFT = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -98,7 +98,7 @@ END {
for (y=1; y<=max_nf; y++) {
for (x=1; x<=max_nf; x++) {
if (x != y && r[x,y]) {
- printf("\n %g %g \t (%s) \t = (%g +/- %g)(%s) \t + (%g +/- %g)",
+ printf("\n %.18g %.18g \t (%s) \t = (%.18g +/- %.18g)(%s) \t + (%.18g +/- %.18g)",
meas_cov[x,y], est_err_cov[x,y], header[y], b[x,y], b_err[x,y], header[x],
a[x,y], a_err[x,y])
}
diff --git a/kalman_init.awk b/kalman_init.awk
index 9e0cd14..0ed72a6 100644
--- a/kalman_init.awk
+++ b/kalman_init.awk
@@ -4,7 +4,7 @@
# kalman csv experiemnt
BEGIN {
- OMFT = "%.9"
+ OMFT = "%.18"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -68,10 +68,10 @@ NF != 0 {
est_err[n] = (1.0 - KG[n])*est_err_last[n]
### visual check
- printf "meas: \t%g %g %g", $n, meas_err[n], count[n]
- printf "\nd1: \t%g", diff[n]
- printf "\nd2: \t%g", diff2[n]
- printf "\nest: \t%g %g %g", est[n], est_err[n], KG[n]
+ printf "meas: \t%.18g %.18g %.18g", $n, meas_err[n], count[n]
+ printf "\nd1: \t%.18g", diff[n]
+ printf "\nd2: \t%.18g", diff2[n]
+ printf "\nest: \t%.18g %.18g %.18g", est[n], est_err[n], KG[n]
last2[n] = last[n]
last[n] = $n
diff --git a/lin_reg.awk b/lin_reg.awk
index f82e2e8..8622289 100644
--- a/lin_reg.awk
+++ b/lin_reg.awk
@@ -4,7 +4,7 @@
# simple linear regression between columns
BEGIN {
- OFS = "%.9g"
+ OFS = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -84,7 +84,7 @@ END {
for (y=1; y<=nf_max; y++) {
for (x=1; x<=nf_max; x++) {
if (x != y && r[x,y]) {
- printf("\n %.9g \t (%s) \t = (%.9g +/- %.9g)(%s) \t + (%.9g +/- %.9g)",
+ printf("\n %.18g \t (%s) \t = (%.18g +/- %.18g)(%s) \t + (%.18g +/- %.18g)",
10.0*log(r[x,y]*r[x,y])/log(10), header[y], b[x,y], b_err[x,y], header[x],
a[x,y], a_err[x,y])
}
diff --git a/lin_reg1.awk b/lin_reg1.awk
index 2d21c29..cf3e69a 100755
--- a/lin_reg1.awk
+++ b/lin_reg1.awk
@@ -4,7 +4,7 @@
# simple linear regression between columns
BEGIN {
- OFS = "%.9g"
+ OFS = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -68,7 +68,7 @@ END {
for (x=1; x<=nf_max; x++) {
if (x != y && r[x,y]) {
r2[x,y] = r[x,y]*r[x,y]
- printf("\n %.9g \t (%s) \t = %.9g(%s) \t + %.9g",
+ printf("\n %.18g \t (%s) \t = %.18g(%s) \t + %.18g",
10.0*log(r2[x,y])/log(10), header[y], b[x,y], header[x], a[x,y])
}
}
diff --git a/lin_reg2.awk b/lin_reg2.awk
index 916b6cb..dea52b5 100644
--- a/lin_reg2.awk
+++ b/lin_reg2.awk
@@ -4,7 +4,7 @@
# simple linear regression between columns
BEGIN {
- OFS = "%.9g"
+ OFS = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
@@ -84,7 +84,7 @@ END {
for (y=1; y<=nf_max; y++) {
for (x=1; x<=nf_max; x++) {
if (x != y && r[x,y]) {
- printf("\n %.9g \t (%s) \t = (%.9g +/- %.9g)(%s) \t + (%.9g +/- %.9g)",
+ printf("\n %.18g \t (%s) \t = (%.18g +/- %.18g)(%s) \t + (%.18g +/- %.18g)",
10.0*log(r[x,y]*r[x,y])/log(10), header[y], b[x,y], b_err[x,y], header[x],
a[x,y], a_err[x,y])
}
diff --git a/mean.awk b/mean.awk
index c1d31a6..54b2e14 100644
--- a/mean.awk
+++ b/mean.awk
@@ -4,6 +4,7 @@
# calculate mean average
BEGIN {
+ OFS = "%.18g"
sign = "[+-±]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
diff --git a/mean_avg.awk b/mean_avg.awk
index 945ed3c..6f5a270 100644
--- a/mean_avg.awk
+++ b/mean_avg.awk
@@ -4,6 +4,7 @@
# average columns of numerical data
BEGIN {
+ OFS = "%.18g"
sign = "[+-±]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]*"
diff --git a/sum1.awk b/sum1.awk
index 2f5a041..8617e50 100644
--- a/sum1.awk
+++ b/sum1.awk
@@ -14,7 +14,7 @@
END {
for (i=1; i<=nf_max; i++) {
- printf("%g", sum[i])
- (i < nf_max) ? printf(" ") : printf("\n")
+ printf("%.18g", sum[i])
+ printf((i < nf_max) ? OFS : ORS)
}
}
diff --git a/sum2.awk b/sum2.awk
index 1df83ae..413f633 100644
--- a/sum2.awk
+++ b/sum2.awk
@@ -14,5 +14,5 @@ NR==1 { nf_max = NF }
END {
for (i=1; i<=NF; i++)
- printf("%g%s", sum[i], i < nf_max ? " " : "\n")
+ printf("%.18g%s", sum[i], i < nf_max ? OFS : ORS)
}
diff --git a/sum3.awk b/sum3.awk
index 8ceda8b..bca92e3 100644
--- a/sum3.awk
+++ b/sum3.awk
@@ -25,8 +25,8 @@ NR==1 {
END {
for (i=1; i<=nfld; i++) {
- (numcol[i]) ? printf("%g", sum[i]) : printf("--")
- printf(i < nfld ? " " : "\n")
+ printf(numcol[i] ? sum[i] : "--")
+ printf(i < nfld ? OFS : ORS)
}
}
diff --git a/sum4.awk b/sum4.awk
index 518f6e5..f9037da 100644
--- a/sum4.awk
+++ b/sum4.awk
@@ -5,6 +5,7 @@
# output: sums of numeric columns
function isnum(n) {
+ OFMT = "%.18g"
sign = "[+-]?"
decimal = "[0-9]+[.]?[0-9]*"
fraction = "[.][0-9]+"
@@ -16,8 +17,7 @@ function isnum(n) {
NR==1 {
nf_max = NF
for (i=1; i<=NF; i++) {
- if (!isnum($i))
- header[i] = $i
+ (!isnum($i)) ? header[i] = $i : header[i] = "col" i
}
}
@@ -32,7 +32,7 @@ END {
for (i=1; i<=nf_max; i++) {
if (header[i])
printf("%s: \t", header[i])
- printf("%g\n", sum[i])
- printf(i < nf_max ? "" : "\n")
+ printf("%.18g\n", sum[i])
+ printf(i < nf_max ? "" : OFS )
}
}
generated by cgit v1.2.3 (git 2.25.1) at 2026-03-18 05:51:21 +0000