/contrib/ntp/scripts/stats/etf.S
Assembly | 15 lines | 15 code | 0 blank | 0 comment | 0 complexity | a4689f09592a62cff31e6d780ba9ec3b MD5 | raw file
Possible License(s): MPL-2.0-no-copyleft-exception, BSD-3-Clause, LGPL-2.0, LGPL-2.1, BSD-2-Clause, 0BSD, JSON, AGPL-1.0, GPL-2.0
1options(digits=4) 2file2 <- "etf_summary" 3etf <- scan(file1, list(day=0, sec=0, offset=0, stab=0)) 4r <- lsfit(etf$sec, etf$offset) 5count<-length(etf$sec) 6mean<-r$coef[[1]] 7std<-sqrt(var(r$residuals)) 8slope<-r$coef[[2]] * 1000 9cat("\n", file=file2 , append=TRUE, fill=FALSE, sep="") 10cat(file1, "\n", file=file2, append=TRUE, fill=FALSE, sep="") 11cat("etf1 ", count, ", T ", mean, " ns, R ", slope, " ps/s, std ", std, " us\n", file=file2, append=TRUE, fill=FALSE, sep="") 12str <- paste("eps/", file1, ".eps", sep="") 13postscript(str, , , , 5, pointsize=18) 14par(mgp=c(1, 0, 0), tck=0.03, mar=c(2, 2, 1, 1)) 15plot(etf$sec, etf$offset, type="l", xlab=paste("MJD", etf$day, "Time (s)"), ylab="External Offset (ns)", ylim=c(-400, 400))