started tracking t_est[] and t_up[];

refactored code, moving math from ck_uptime() to BEGIN block;
switched to using t_up(new) = t_up(old) + dt to predict uptime;
added conditional use of "date +%s.%N" to get t as float on FreeBSD;
merged some changes between ckt0.awk and ckcptime.awk;

1 files changed, 7 insertions, 3 deletions

diff --git a/ckt0.sh b/ckt0.sh
index eeb0479..88fbb78 100644
--- a/ckt0.sh
+++ b/ckt0.sh
@@ -11,20 +11,24 @@ set -e
 #sysctl_t0_str=$( sysctl kern.boottime | mawk -F= '{print($NF)}' )
 #sysctl_t0=$( date -jf "%a %b %d %H:%M:%S %G" "${sysctl_t0_str}" +%s )
 
+
+# check for saved state
 if [ -f t0_est.log ] ; then {
     t0_est="$( cat t0_est.log )"
 } else {
     t0_est=""
 } fi
 
-printf "N t0_est t0_unc\n" | column -t
+printf "n t_est t_unc t_up_est t_up_unc t0_est t0_unc\n" | column -t
 for N in $(seq 0 10080) ; do
 
     printf "${N} ${t0_est}\n" | column -t
     t0_est=$( mawk -f ckt0.awk ${t0_est} | tail -n2 )
     #mawk -f ckt0.awk
-    sleep $( mawk 'BEGIN {print( sqrt(rand()^2.0 + 1.0) )}' )
-    #sleep $( echo $t0_est | mawk '{print(sqrt(1.0 + ((60.0*$NF)/(60.0 + $NF))^2.0 + rand()^2.0 ))}' )
+    sleep $( mawk 'BEGIN {print(sqrt( rand()^2.0 + 1.0 ))}' )
+    #sleep $( echo $t0_est | mawk '{print(sqrt( ((60.0*$4)/(60.0 + $4))^2.0 + rand()^2.0 ))}' )
+    #sleep $( echo $t0_est | mawk '{print(sqrt( ($2)^2.0 + ($4)^2.0 + ($6)^2.0 + rand()^2.0 ))}' )
+    #sleep $( echo $t0_est | mawk '{print( 10.0 - (systime() - $5)%10.0 + 0.5*rand()^2.0 )}' )
 
 done