From 4916e9b13765de970deff094abb3eb50c663834a Mon Sep 17 00:00:00 2001 From: wukong Date: Sat, 16 Jun 2018 23:38:02 -0700 Subject: bumped up number of printed digits --- conv.awk | 2 +- diff.awk | 6 +++--- diff1.awk | 4 ++-- hamming.awk | 2 +- kalman.awk | 4 ++-- kalman_init.awk | 10 +++++----- lin_reg.awk | 4 ++-- lin_reg1.awk | 4 ++-- lin_reg2.awk | 4 ++-- mean.awk | 1 + mean_avg.awk | 1 + sum1.awk | 4 ++-- sum2.awk | 2 +- sum3.awk | 4 ++-- sum4.awk | 8 ++++---- 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 ) } } -- cgit v1.2.3