bumped up number of printed digits

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 )
     }
 }